Abstract:The computational investigation of photochemical processes often entails the calculation of excited state geometries, energies, and energy gradients. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons, thereby including the associated nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations. The multicomponent density functional theory (NEO-DFT) and time-dependent DFT (NEO-TDDFT) methods allow efficient calculations of ground and excited states, respectively. Herein, the analytical gradients are derived and implemented for the NEO-TDDFT method and the associated Tamm-Dancoff approximation (NEO-TDA). The programmable equations for these analytical gradients, as well as the NEO-DFT analytical Hessian, are provided. The NEO approach includes the anharmonic zero-point energy and density delocalization associated with the quantum protons, as well as vibronic mixing, in geometry optimizations and energy calculations of ground and excited states. The harmonic zero-point energy associated with the other nuclei can be computed via the NEO Hessian. This approach is used to compute the 0-0 adiabatic excitation energies for a set of nine small molecules with all protons quantized, exhibiting slight improvement over the conventional electronic approach. Geometry optimizations of two excited state intramolecular proton transfer systems are performed with one and two quantized protons, respectively. The NEO calculations for these systems produce electronically excited state geometries with stronger intramolecular hydrogen bonds and similar relative stabilities compared to conventional electronic methods. This work provides the foundation for nonadiabatic dynamics simulations of fundamental processes such as photoinduced proton transfer and proton-coupled electron transfer.

Strategy for constructing compact numerical atomic orbital basis sets by incorporating the gradients of reference wavefunctions

Systematically Improvable Optimized Atomic Basis Sets Forab Initiocalculations

All-electron periodic G 0 W 0 implementation with numerical atomic orbital basis functions: Algorithm and benchmarks

Electronic Structure Interpolation Via Atomic Orbitals

Implementation of Screened Hybrid Density Functional for Periodic Systems with Numerical Atomic Orbitals: Basis Function Fitting and Integral Screening.

A General Second Order Complete Active Space Self-Consistent-field Solver for Large-Scale Systems

Efficient implementation of analytical gradients for periodic hybrid functional calculations within fitted numerical atomic orbitals from NAO2GTO

Resolution-of-identity approach to Hartree–Fock, hybrid density functionals, RPA, MP2 and GW with numeric atom-centered orbital basis functions

Accurate stress calculations based on numerical atomic orbital bases: Implementation and benchmarks

Optimization of atomic density-fitting basis functions for molecular two-electron integral approximations

An Efficient Implementation of Analytical Nuclear Gradients for Linear-Response Time-Dependent Density Functional Theory in the Plane Wave Basis

Optimized Local Basis Set for Kohn-Sham Density Functional Theory.

Accurate Stress Calculations Based On Numerical Atomic Orbital Bases: And Benchmarks

Analytical Gradients for Projection-Based Wavefunction-in-DFT Embedding

Introduction to First-Principles Simulation Package Abacus Based on Systematically Improvable Atomic Orbitals

Efficient Calculation of Derivatives of Integrals in a Basis of Non-Separable Gaussians Through Exploitation of Sparsity

Adaptive Algorithm for Electronic Structure Calculations Using Reduction of Gaussian Mixtures

Analytic Gradients for Density Fitting MP2 Using Natural Auxiliary Functions

Calculation of electric field gradients with the exact two-component (X2C) quasi-relativistic method and its local approximations

Robust linear-scaling optimization of compact localized orbitals in density functional theory

Analytical Gradients for Nuclear-Electronic Orbital Time-Dependent Density Functional Theory: Excited State Geometry Optimizations and Adiabatic Excitation Energies