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.

Molecular-Orbital-Free Algorithm For The Excited-State Force In Time-Dependent Density Functional Theory

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

Excited state orbital optimization via minimizing the square of the gradient: General approach and application to singly and doubly excited states via density functional theory

Excitation Energies from Thermally-Assisted-Occupation Density Functional Theory: Theory and Computational Implementation

A New Fragment-Based Approach For Calculating Electronic Excitation Energies Of Large Systems

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

Time-Dependent Density Functional Theory Excited-State Forces with the Gaussian and Augmented Plane Wave Method

Direct Unconstrained Optimization of Molecular Orbital Coefficients in Density Functional Theory

Orbital-optimized versus time-dependent density functional calculations of intramolecular charge transfer excited states

Efficient time-dependent orbital-free density functional theory: Semilocal adiabatic response

Orbital Optimized Density Functional Theory for Electronic Excited States

Time-Dependent Quantum Transport: An Efficient Method Based On Liouville-Von-Neumann Equation For Single-Electron Density Matrix

A Constraint-Based Orbital-Optimized Excited State Method (COOX)

Optimization of highly excited matrix product states with an application to vibrational spectroscopy

A Simplified Scheme for Relativistic Density Functional Computation in the Zeroth-Order Regular Approximation

Nonadiabatic excited-state dynamics and energy gradients in the framework of FMO-LC-TDDFTB

Foundation for the ΔSCF Approach in Density Functional Theory

Orbital-optimized Density Functional Calculations of Molecular Rydberg Excited States with Real Space Grid Representation and Self-Interaction Correction

Orbital Energies Are Chemical Potentials in Ground-State Density Functional Theory and Excited-State $Δ$SCF Theory

Developing Orbital-Dependent Corrections for the Non-Additive Kinetic Energy in Subsystem Density Functional Theory