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.

Analytic evaluation of energy gradients for the single and double excitation coupled cluster (CCSD) wave function: Theory and application

Analytic evaluation of energy gradients for the single, double and linearized triple excitation coupled-cluster CCSDT-1 wavefunction: Theory and applications

Analytic evaluation of energy gradients for the singles and doubles coupled cluster method including perturbative triple excitations: Theory and applications to FOOF and Cr2

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

Coupled cluster energy derivatives. Analytic Hessian for the closed‐shell coupled cluster singles and doubles wave function: Theory and applications

Analytical evaluation of ground state gradients in quantum electrodynamics coupled cluster theory

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

Analytical Energy Gradient Evaluation in Relativistic and Nonrelativistic Density Functional Calculations

Cavity Quantum-Electrodynamical Time-Dependent Density Functional Theory Within Gaussian Atomic Basis. II. Analytic Energy Gradient.

Analytical Energy Gradients for the Cluster-in-Molecule MP2 Method and Its Application to Geometry Optimizations of Large Systems

Analytic Energy Gradient of Excited Electronic State Within TDDFT/MMpol Framework: Benchmark Tests and Parallel Implementation.

Analytical Energy Gradients for Ionized States Using Equation-of-motion Coupled-Cluster Theory with Spin Orbit Coupling

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

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

The optimization of molecular orbitals for coupled cluster wavefunctions

Converging High-Level Coupled-Cluster Energetics by Monte Carlo Sampling and Moment Expansions

Coupled cluster theory for nonadiabatic dynamics: nuclear gradients and nonadiabatic couplings in similarity constrained coupled cluster theory

Development and Application of the Analytical Energy Gradient for the Normalized Elimination of the Small Component Method.

Analytic Energy Gradients and Hessians of Exact Two-Component Relativistic Methods: Efficient Implementation and Extensive Applications.

A systematic theoretical study of harmonic vibrational frequencies: The single and double excitation coupled cluster (CCSD) method