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
Assessment of electron–proton correlation functionals for vibrational spectra of shared-proton systems by constrained nuclear-electronic orbital density functional theory
Proton transfer plays a crucial role in various chemical and biological processes. A major theoretical challenge in simulating proton transfer arises from the quantum nature of the proton. The constrained nuclear-electronic orbital (CNEO) framework was recently developed to efficiently and accurately account for nuclear quantum effects, particularly quantum nuclear delocalization effects, in quantum chemistry calculations and molecular dynamics simulations. In this paper, we systematically investigate challenging proton transfer modes in a series of shared-proton systems using CNEO density functional theory (CNEO-DFT), focusing on evaluating existing electron–proton correlation functionals. Our results show that CNEO-DFT accurately describes proton transfer vibrational modes and significantly outperforms conventional DFT. The inclusion of the epc17-2 electron–proton correlation functional in CNEO-DFT produces similar performance to that without electron–proton correlations, while the epc17-1 functional yields less accurate results, comparable with conventional DFT. These findings hold true for both asymmetrical and symmetrical shared-proton systems. Therefore, until a more accurate electron–proton correlation functional is developed, we currently recommend performing vibrational spectrum calculations using CNEO-DFT without electron–proton correlation functionals.
more »
« less
- Award ID(s):
- 2238473
- PAR ID:
- 10650164
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 161
- Issue:
- 24
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Phenol–benzimidazole and phenol–pyridine proton-coupled electron transfer (PCET) dyad systems are computationally investigated to resolve the origins of the asymmetrically broadened H-bonded OH stretch transitions that have been previously reported using cryogenic ion vibrational spectroscopy in the ground electronic state. Two-dimensional (2D) potentials describing the strongly shared H atom are predicted to be very shallow along the H atom transfer coordinate, enabling dislocation of the H atom between the donor and acceptor groups upon excitation of the OH vibrational modes. These soft H atom potentials result in strong coupling between the OH modes, which exhibit significant bend-stretch mixing, and a large number of normal mode coordinates. Vibrational spectra are calculated using a Hamiltonian that linearly and quadratically couples the H atom potentials to over two dozen of the most strongly coupled normal modes treated at the harmonic level. The calculated vibrational spectra qualitatively reproduce the asymmetric shape and breadth of the experimentally observed bands in the 2300–3000 cm–1 range. Interestingly, these transitions fall well above the predicted OH stretch fundamentals, which are computed to be surprisingly red-shifted (<2000 cm–1). Time-dependent calculations predict rapid (<100 fs) relaxation of the excited OH modes and instant response from the lower-frequency normal modes, corroborating the strong coupling predicted by the model Hamiltonian. The results highlight a unique broadening mechanism and complicated anharmonic effects present within these biologically relevant PCET model systems.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
-
We present a graph-theory-based reformulation of all ONIOM-based molecular fragmentation methods. We discuss applications to (a) accurate post-Hartree–Fock AIMD that can be conducted at DFT cost for medium-sized systems, (b) hybrid DFT condensed-phase studies at the cost of pure density functionals, (c) reduced cost on-the-fly large basis gas-phase AIMD and condensed-phase studies, (d) post-Hartree–Fock-level potential surfaces at DFT cost to obtain quantum nuclear effects, and (e) novel transfer machine learning protocols derived from these measures. Additionally, in previous work, the unifying strategy discussed here has been used to construct new quantum computing algorithms. Thus, we conclude that this reformulation is robust and accurate.more » « less
An official website of the United States government

