An official website of the United States government
Proton transfer is crucial in various chemical and biological processes. Because of significant nuclear quantum effects, accurate and efficient description of proton transfer remains a great challenge. In this Communication, we apply constrained nuclear–electronic orbital density functional theory (CNEO-DFT) and constrained nuclear–electronic orbital molecular dynamics (CNEO-MD) to three prototypical shared proton systems and investigate their proton transfer modes. We find that with a good description of nuclear quantum effects, CNEO-DFT and CNEO-MD can well describe the geometries and vibrational spectra of the shared proton systems. Such a good performance is in significant contrast to DFT and DFT-based ab initio molecular dynamics, which often fail for shared proton systems. As an efficient method based on classical simulations, CNEO-MD is promising for future investigations of larger and more complex proton transfer systems.
more »
« less
Triple electron–electron–proton excitations and second-order approximations in nuclear–electronic orbital coupled cluster methods
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
- Award ID(s):
- 1954348
- PAR ID:
- 10415455
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 157
- Issue:
- 7
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- 074104
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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.more » « less
-
Proton transfer and hydrogen tunneling play key roles in many processes of chemical and biological importance. The generalized nuclear-electronic orbital multistate density functional theory (NEO-MSDFT) method was developed in order to capture hydrogen tunneling effects in systems involving the transfer and tunneling of one or more protons. The generalized NEO-MSDFT method treats the transferring protons quantum mechanically on the same level as the electrons and obtains the delocalized vibronic states associated with hydrogen tunneling by mixing localized NEO-DFT states in a nonorthogonal configuration interaction scheme. Herein, we present the derivation and implementation of analytical gradients for the generalized NEO-MSDFT vibronic state energies and the nonadiabatic coupling vectors between these vibronic states. We use this methodology to perform adiabatic and nonadiabatic dynamics simulations of the double proton transfer reactions in the formic acid dimer and the heterodimer of formamidine and formic acid. The generalized NEO-MSDFT method is shown to capture the strongly coupled synchronous or asynchronous tunneling of the two protons in these processes. Inclusion of vibronically nonadiabatic effects is found to significantly impact the double proton transfer dynamics. This work lays the foundation for a variety of nonadiabatic dynamics simulations of multiple proton transfer systems, such as proton relays and hydrogen-bonding networks.more » « less
-
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
-
The interaction between free electrons and photons in electron microscopes offers unique opportunities for microscopy and quantum science. For example, modulating electron beams with multiple laser excitations, researchers have demonstrated a novel near-field electron microscope, capable of probing electromagnetic excitations on the nanometer spatial scale and in the attosecond (10 −18 s) temporal range [see D. Nabben et al., Nature, 619, 63 (2023)]. Additionally, it has recently been demonstrated that the interaction between free electrons and photons in an electron microscope can be quantum coherent, and furthermore, this quantum coherence could potentially be leveraged for heralded sources of single electrons and photons [see A. Feist et al., Science, 377, 777 (2022)]. Although promising, these innovations in free-electron-photon interactions have thus far suffered a significant limitation: they require high-energy (>100-ke V) electron beams. Accordingly, these demonstrations have taken place in energetic (and expensive) transmission electron microscopes (TEMs). TEMs are a logical setting for these experiments, as their high-energy electrons can be velocity-matched to co-propagating photons in dielectric waveguides. However, achieving such velocity-matching between photons in conventional dielectric waveguides and electrons is not feasible for the low electron energies (<30-keV) in more common scanning electron microscope (SEMs).more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0106173
