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

Zhen Tao,Saswata Roy,Patrick E. Schneider,Fabijan Pavošević,Sharon Hammes-Schiffer
DOI: https://doi.org/10.48550/arXiv.2105.05191
2021-05-12
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.
Chemical Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to derive and implement the analytical gradient in the nuclear - electron orbital time - dependent density functional theory (NEO - TDDFT) for excited - state geometry optimization and the calculation of adiabatic excitation energy. Specifically, the researchers hope: 1. **Derive and implement the analytical gradient of NEO - TDDFT**: This includes the derivation and programming implementation of the analytical gradients of the NEO - TDDFT method and its related Tamm - Dancoff approximation (NEO - TDA). 2. **Calculate the excited - state geometric structure**: Through the analytical gradient, the excited - state geometry optimization can be carried out, thereby obtaining a more accurate excited - state geometric structure. 3. **Calculate the adiatic excitation energy**: Using the analytical gradient, the adiatic excitation energy (ΔEadia) from the ground state to the excited state can be calculated, as well as the 0 - 0 adiatic excitation energy (ΔE0 - 0) including the zero - point energy correction. 4. **Verify the accuracy of the method**: Through the calculations of a series of small molecules, verify the accuracy and effectiveness of the NEO method in dealing with the quantum proton effect and the non - Born - Oppenheimer effect. 5. **Study the proton transfer process in specific systems**: In particular, perform excited - state geometry optimization on two systems with internal intermolecular proton transfer (ESIPT) - [2,2′ - bipyridine] - 3 - ol and [2,2′ - bipyridine] - 3,3′ - diol - to study the influence of the quantum nuclear effect on their geometric structures and energies. Through these works, the researchers hope to provide a basis for the non - adiabatic dynamics simulation of basic processes such as photo - induced proton transfer and proton - coupled electron transfer.