Hydrogen tunneling plays a critical role in many biologically and chemically important processes. The nuclear–electronic orbital multistate density functional theory (NEO-MSDFT) method was developed to describe hydrogen transfer systems. In this approach, the transferring proton is treated quantum mechanically on the same level as the electrons within multicomponent DFT, and a nonorthogonal configuration interaction scheme is used to produce delocalized vibronic states from localized vibronic states. The NEO-MSDFT method has been shown to provide accurate hydrogen tunneling splittings for fixed molecular systems. Herein, the NEO-MSDFT analytical gradients for both ground and excited vibronic states are derived and implemented. The analytical gradients and semi-numerical Hessians are used to optimize and characterize equilibrium and transition state geometries and to generate minimum energy paths (MEPs), for proton transfer in the deprotonated acetylene dimer and malonaldehyde. The barriers along the resulting MEPs are lower when the transferring proton is quantized because the NEO-MSDFT method inherently includes the zero-point energy of the transferring proton. Analysis of the proton densities along the MEPs illustrates that the proton density can exhibit symmetric or asymmetric bilobal character associated with symmetric or slightly asymmetric double-well potential energy surfaces and hydrogen tunneling. Analysis of the contributions to the intrinsic reaction coordinate reveals that changes in the C–O bond lengths drive proton transfer in malonaldehyde. This work provides the foundation for future reaction path studies and direct nonadiabatic dynamics simulations of a wide range of hydrogen transfer reactions.
Within the nuclear-electronic orbital (NEO) framework, the real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) approach enables the simulation of coupled electronic-nuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of long-time nuclear quantum dynamics. Herein, the electronic Born–Oppenheimer (BO) approximation within the NEO framework is presented. In this approach, the electronic density is quenched to the ground state at each time step, and the real-time nuclear quantum dynamics is propagated on an instantaneous electronic ground state defined by both the classical nuclear geometry and the nonequilibrium quantum nuclear density. Because the electronic dynamics is no longer propagated, this approximation enables the use of an order-of-magnitude larger time step, thus greatly reducing the computational cost. Moreover, invoking the electronic BO approximation also fixes the unphysical asymmetric Rabi splitting observed in previous semiclassical RT-NEO-TDDFT simulations of vibrational polaritons even for small Rabi splitting, instead yielding a stable, symmetric Rabi splitting. For the intramolecular proton transfer in malonaldehyde, both RT-NEO-Ehrenfest dynamics and its BO counterpart can describe proton delocalization during the real-time nuclear quantum dynamics. Thus, the BO RT-NEO approach provides the foundation for a wide range of chemical and biological applications.
more » « less- Award ID(s):
- 1954348
- NSF-PAR ID:
- 10440610
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 11
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
Laser-induced molecular alignment is well understood within the framework of the Born–Oppenheimer (BO) approximation. Without the BO approximation, however, the concept of molecular structure is lost, making it hard to precisely define alignment. In this work, we demonstrate the emergence of alignment from the first-ever non-BO quantum dynamics simulations, using the HD molecule exposed to ultrashort laser pulses as a few-body test case. We extract the degree of alignment from the non-BO wave function by means of an operator expressed in terms of pseudo-proton coordinates that mimics the BO-based definition of alignment. The only essential approximation, in addition to the semiclassical electric-dipole approximation for the matter–field interaction, is the choice of time-independent explicitly correlated Gaussian basis functions. We use a variational, electric-field-dependent basis-set construction procedure, which allows us to keep the basis-set dimension low while capturing the main effects of electric polarization on the nuclear and electronic degrees of freedom. The basis-set construction procedure is validated by comparing with virtually exact grid-based simulations for two one-dimensional model systems: laser-driven electron dynamics in a soft attractive Coulomb potential and nuclear rovibrational dynamics in a Morse potential.more » « less
-
The accurate description of nuclear quantum effects, such as zero-point energy, is important for modeling a wide range of chemical and biological processes. Within the nuclear–electronic orbital (NEO) approach, such effects are incorporated in a computationally efficient way by treating electrons and select nuclei, typically protons, quantum mechanically with molecular orbital techniques. Herein, we implement and test a NEO coupled cluster method that explicitly includes the triple electron–electron–proton excitations, where two electrons and one proton are excited simultaneously, using automatic differentiation. Our calculations show that this NEO-CCSDT eep method provides highly accurate proton densities and proton affinities, outperforming any previously studied NEO method. These examples highlight the importance of the triple electron–electron–proton excitations for an accurate description of nuclear quantum effects. Additionally, we also implement and test the second-order approximate coupled cluster with singles and doubles (NEO-CC2) method as well as its scaled-opposite-spin (SOS) versions. The NEO-SOS′-CC2 method, which scales the electron–proton correlation energy as well as the opposite-spin and same-spin components of the electron–electron correlation energy, achieves nearly the same accuracy as the NEO-CCSDT eep method for the properties studied. Because of its low computational cost, this method will enable a wide range of chemical and photochemical applications for large molecular systems. This work sets the stage for a variety of developments and applications within the NEO framework.more » « less
-
more » « less
Abstract The Hohenberg-Kohn theorem of density-functional theory establishes the existence of a bijection between the ground-state electron density and the external potential of a many-body system. This guarantees a one-to-one map from the electron density to all observables of interest including electronic excited-state energies. Time-Dependent Density-Functional Theory (TDDFT) provides one framework to resolve this map; however, the approximations inherent in practical TDDFT calculations, together with their computational expense, motivate finding a cheaper, more direct map for electronic excitations. Here, we show that determining density and energy functionals via machine learning allows the equations of TDDFT to be bypassed. The framework we introduce is used to perform the first excited-state molecular dynamics simulations with a machine-learned functional on malonaldehyde and correctly capture the kinetics of its excited-state intramolecular proton transfer, allowing insight into how mechanical constraints can be used to control the proton transfer reaction in this molecule. This development opens the door to using machine-learned functionals for highly efficient excited-state dynamics simulations.
-
A computational strategy to simulate two-dimensional electronic spectra (2DES) is introduced, which allows us to analyse ground state dynamics and to sample and measure different conformations attained by flexible molecular systems in solution. An explicit mixed quantum mechanics/molecular mechanics (QM/MM) approach is employed for the evaluation of the necessary electronic excited state energies and transition dipole moments. The method is applied towards a study of the highly flexible water-solvated adenine–adenine monophosphate (ApA), a system featuring two interacting adenine moieties that display various intermolecular arrangements, known to deeply affect their photochemical outcome. Molecular dynamics simulations and cluster analysis have been used to select the molecular conformations, reducing the complexity of the flexible ApA conformational space. By using our sum-over-states (SOS) approach to obtain the 2DES spectra for each of these selected conformations, we can discern spectral changes and relate them to specific nuclear arrangements: close lying π-stacked bases exhibit a splitting of their respective 1 L a signal traces; T-stacked bases exhibit the appearance of charge transfer states in the low-energy Vis probing window while displaying no 1 L a splitting, being particularly favoured when promoting amino to 5-ring interactions; unstacked and distant adenine moieties exhibit signals similar to those of the adenine monomer, as is expected for non-interacting nucleobases. 2DES maps reveal the spectral fingerprints associated with specific molecular conformations, and are thus a promising option to enable their quantitative spectroscopic detection beyond standard 1D pump-probe techniques. This is expected to aid the understanding of how nucleobase aggregation controls and modulates the photostability and photo-damage of extended DNA/RNA systems.more » « less
