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.

Analytical Hartree-Fock gradients for periodic systems

Analytical Hartree-Fock gradients with respect to the cell parameter for systems periodic in three dimensions

General technique for analytical derivatives of post-projected Hartree-Fock

Sparse tensor based nuclear gradients for periodic Hartree-Fock and low-scaling correlated wave function methods in the CP2K software package: A massively parallel and GPU accelerated implementation

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

Analytic Gradient for Time-Dependent Density Functional Theory Combined with the Fragment Molecular Orbital Method

Analytical Energy Gradient Evaluation in Relativistic and Nonrelativistic Density Functional Calculations

Analytic Gradients for the Long-Range-Corrected XYG3 Type of Doubly Hybrid Density Functionals: Theory, Implementation, and Assessment.

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

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Closed-shell Coupled-Cluster Theory with Spin-Orbit Coupling

Efficient Treatment of Relativistic Effects with Periodic Density Functional Methods: Energies, Gradients, and Stress Tensors

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

A Parallel Implementation of the Analytic Nuclear Gradient for Time-Dependent Density Functional Theory Within the Tamm–Dancoff Approximation

Analytic gradients for equation-of-motion coupled cluster with single, double, and perturbative triple excitations

Efficient periodic resolution-of-the-identity Hartree–Fock exchange method with k -point sampling and Gaussian basis sets

Analytic gradients for compressed multistate pair-density functional theory

Density fitting in periodic systems: application to TDHF in diamond and oxides

Current density functional framework for spin-orbit coupling: Extension to periodic systems

Analytical Gradients for Electron-Attached and Ionized States for the Algebraic-Diagrammatic Construction Scheme for the Electron Propagator up to Third Order

Hybrid functionals for periodic systems in the density functional tight-binding method