1
|
Dickinson JA, Hammes-Schiffer S. Nonadiabatic Hydrogen Tunneling Dynamics for Multiple Proton Transfer Processes with Generalized Nuclear-Electronic Orbital Multistate Density Functional Theory. J Chem Theory Comput 2024. [PMID: 39259939 DOI: 10.1021/acs.jctc.4c00737] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/13/2024]
Abstract
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.
Collapse
Affiliation(s)
- Joseph A Dickinson
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
2
|
Dickinson JA, Yu Q, Hammes-Schiffer S. Generalized Nuclear-Electronic Orbital Multistate Density Functional Theory for Multiple Proton Transfer Processes. J Phys Chem Lett 2023:6170-6178. [PMID: 37379485 DOI: 10.1021/acs.jpclett.3c01422] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/30/2023]
Abstract
Proton transfer and hydrogen tunneling play pivotal roles in many chemical and biological processes. The nuclear-electronic orbital multistate density functional theory (NEO-MSDFT) approach was developed to describe hydrogen tunneling systems within the multicomponent NEO framework, where the transferring proton is quantized and treated with molecular orbital techniques on the same level as the electrons. Herein, the NEO-MSDFT framework is generalized to an arbitrary number of quantum protons to allow applications to systems involving the transfer and tunneling of multiple protons. The generalized NEO-MSDFT approach is shown to produce delocalized, bilobal proton densities and accurate tunneling splittings for fixed geometries of the formic acid dimer and asymmetric substituted variants, as well as the porphycene molecule. Investigation of a protonated water chain highlights the applicability of this approach to proton relay systems. This work provides the foundation for nuclear-electronic quantum dynamics simulations of a wide range of multiple proton transfer processes.
Collapse
Affiliation(s)
- Joseph A Dickinson
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
3
|
Li F, Yang X, Liu X, Cao J, Bian W. An Ab Initio Neural Network Potential Energy Surface for the Dimer of Formic Acid and Further Quantum Tunneling Dynamics. ACS OMEGA 2023; 8:17296-17303. [PMID: 37214673 PMCID: PMC10193396 DOI: 10.1021/acsomega.3c02169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Accepted: 04/19/2023] [Indexed: 05/24/2023]
Abstract
We construct a full-dimensional ab initio neural network potential energy surface (PES) for the isomerization system of the formic acid dimer (FAD). This is based upon ab initio calculations using the DLPNO-CCSD(T) approach with the aug-cc-pVTZ basis set, performed at over 14000 symmetry-unique geometries. An accurate fit to the obtained energies is generated using a general neural network fitting procedure combined with the fundamental invariant method, and the overall energy-weighted root-mean-square fitting error is about 6.4 cm-1. Using this PES, we present a multidimensional quantum dynamics study on tunneling splittings with an efficient theoretical scheme developed by our group. The ground-state tunneling splitting of FAD calculated with a four-mode coupled method is in good agreement with the most recent experimental measurements. The PES can be applied for further dynamics studies. The effectiveness of the present scheme for constructing a high-dimensional PES is demonstrated, and this scheme is expected to be feasible for larger molecular systems.
Collapse
Affiliation(s)
- Fengyi Li
- Beijing
National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, People’s Republic of China
- School
of Chemical Sciences, University of Chinese
Academy of Sciences, Beijing 100049, People’s
Republic of China
| | - Xingyu Yang
- Beijing
National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, People’s Republic of China
- School
of Chemical Sciences, University of Chinese
Academy of Sciences, Beijing 100049, People’s
Republic of China
| | - Xiaoxi Liu
- Beijing
National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, People’s Republic of China
- School
of Chemical Sciences, University of Chinese
Academy of Sciences, Beijing 100049, People’s
Republic of China
| | - Jianwei Cao
- Beijing
National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, People’s Republic of China
| | - Wensheng Bian
- Beijing
National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, People’s Republic of China
- School
of Chemical Sciences, University of Chinese
Academy of Sciences, Beijing 100049, People’s
Republic of China
| |
Collapse
|
4
|
Liu H, Cao J, Bian W. Efficient Quantum Mechanical Calculations of Mode-Specific Tunneling Splittings upon Fundamental Excitation in the Dimer of Formic Acid. J Phys Chem A 2020; 124:6536-6543. [PMID: 32662997 DOI: 10.1021/acs.jpca.0c05471] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The formic acid dimer (FAD) is an important benchmark system for understanding the double hydrogen transfer process. Most recently, Zhang et al. measured a few tunneling splittings upon fundamental excitation of FAD precisely (Zhang, Y. et al. J. Chem. Phys. 2017, 146, 244306); however, relevant theoretical studies are very limited. Here, we present a multidimensional quantum dynamics study on mode-specific tunneling splittings upon fundamental excitation in FAD with an efficient theoretical scheme developed by our group in which the process-oriented basis function customization strategy is combined with the preconditioned inexact spectral transform method. Various mode-specific tunneling splittings upon fundamental excitation are systematically calculated, and interesting mode-specific excitation effects on tunneling rate are identified. In particular, the calculated tunneling splittings for the ν22 and ν21 states are in good agreement with experiment, and the remarkable mode-specific suppression effects upon excitation should result from that the antisymmetric vibrational modes hinder the concerted double H-transfer. The present work is helpful to acquire a better understanding of the mode-specific excitation effects on tunneling processes.
Collapse
Affiliation(s)
- Hao Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, 100190 Beijing, China.,School of Chemical Sciences, University of Chinese Academy of Sciences, 100049 Beijing, China
| | - Jianwei Cao
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, 100190 Beijing, China
| | - Wensheng Bian
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, 100190 Beijing, China.,School of Chemical Sciences, University of Chinese Academy of Sciences, 100049 Beijing, China
| |
Collapse
|
5
|
Abstract
AbstractFormic acid dimer as the prototypical doubly hydrogen-bonded gas-phase species is discussed from the perspective of the three translational and the three rotational degrees of freedom which are lost when two formic acid molecules form a stable complex. The experimental characterisation of these strongly hindered translations and rotations is reviewed, as are attempts to describe the associated fundamental vibrations, their combinations, and their thermal shifts by different electronic structure calculations and vibrational models. A remarkable match is confirmed for the combination of a CCSD(T)-level harmonic treatment and an MP2-level anharmonic VPT2 correction. Qualitatively correct thermal shifts of the vibrational spectra can be obtained from classical molecular dynamics in CCSD(T)-quality force fields. A detailed analysis suggests that this agreement between experiment and composite theoretical treatment is not strongly affected by fortuitous error cancellation but fully converged variational treatments of the six pair or intermolecular modes and their overtones and combinations in this model system would be welcome.
Collapse
|
6
|
Liu H, Cao J, Bian W. Double Proton Transfer in the Dimer of Formic Acid: An Efficient Quantum Mechanical Scheme. Front Chem 2019; 7:676. [PMID: 31750286 PMCID: PMC6842929 DOI: 10.3389/fchem.2019.00676] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2019] [Accepted: 09/30/2019] [Indexed: 11/13/2022] Open
Abstract
Double proton transfer plays an important role in biology and chemistry, such as with DNA base pairs, proteins and molecular clusters, and direct information about these processes can be obtained from tunneling splittings. Carboxylic acid dimers are prototypes for multiple proton transfer, of which the formic acid dimer is the simplest one. Here, we present efficient quantum dynamics calculations of ground-state and fundamental excitation tunneling splittings in the formic acid dimer and its deuterium isotopologues. These are achieved with a multidimensional scheme developed by us, in which the saddle-point normal coordinates are chosen, the basis functions are customized for the proton transfer process, and the preconditioned inexact spectral transform method is used to solve the resultant eigenvalue problem. Our computational results are in excellent agreement with the most recent experiments (Zhang et al., 2017; Li et al., 2019).
Collapse
Affiliation(s)
- Hao Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing, China.,School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing, China
| | - Jianwei Cao
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing, China
| | - Wensheng Bian
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing, China.,School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
7
|
Qu C, Bowman JM. Quantum approaches to vibrational dynamics and spectroscopy: is ease of interpretation sacrificed as rigor increases? Phys Chem Chem Phys 2019; 21:3397-3413. [DOI: 10.1039/c8cp04990d] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The subject of this Perspective is quantum approaches, beyond the harmonic approximation, to vibrational dynamics and IR spectroscopy.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University
- Atlanta
- USA
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University
- Atlanta
- USA
| |
Collapse
|
8
|
Qu C, Bowman JM. High-dimensional fitting of sparse datasets of CCSD(T) electronic energies and MP2 dipole moments, illustrated for the formic acid dimer and its complex IR spectrum. J Chem Phys 2018; 148:241713. [DOI: 10.1063/1.5017495] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chen Qu
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| | - Joel M. Bowman
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
9
|
Smedarchina Z, Siebrand W, Fernández-Ramos A. Entanglement and co-tunneling of two equivalent protons in hydrogen bond pairs. J Chem Phys 2018; 148:102307. [DOI: 10.1063/1.5000681] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Zorka Smedarchina
- Centro Singular de Investigación en Química Biolóxica e Materiais Moleculares (CIQUS), e Departamento de Química Física, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
- National Research Council of Canada, Ottawa, Ontario K1A 0R6 Canada
| | - Willem Siebrand
- National Research Council of Canada, Ottawa, Ontario K1A 0R6 Canada
| | - Antonio Fernández-Ramos
- Centro Singular de Investigación en Química Biolóxica e Materiais Moleculares (CIQUS), e Departamento de Química Física, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| |
Collapse
|
10
|
Oswald S, Meyer E, Suhm MA. Dinitrogen as a Sensor for Metastable Carboxylic Acid Dimers and a Weak Hydrogen Bond Benchmarking Tool. J Phys Chem A 2018. [DOI: 10.1021/acs.jpca.8b00334] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Sönke Oswald
- Institut für Physikalische Chemie, Georg-August-Universität Göttingen, Tammannstraße 6, 37077 Göttingen, Germany
| | - Enno Meyer
- Institut für Physikalische Chemie, Georg-August-Universität Göttingen, Tammannstraße 6, 37077 Göttingen, Germany
| | - Martin A. Suhm
- Institut für Physikalische Chemie, Georg-August-Universität Göttingen, Tammannstraße 6, 37077 Göttingen, Germany
| |
Collapse
|
11
|
Qu C, Bowman J. Quantum and classical IR spectra of (HCOOH)2, (DCOOH)2 and (DCOOD)2 using ab initio potential energy and dipole moment surfaces. Faraday Discuss 2018; 212:33-49. [DOI: 10.1039/c8fd00077h] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Full-dimensional (24 modes) quantum calculation of the IR spectrum of (DCOOD)2, and comparison with classical MD one.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry
- Cherry L. Emerson Center for Scientific Computations
- Emory University
- Atlanta
- USA
| | - Joel M. Bowman
- Department of Chemistry
- Cherry L. Emerson Center for Scientific Computations
- Emory University
- Atlanta
- USA
| |
Collapse
|
12
|
Richardson JO. Full- and reduced-dimensionality instanton calculations of the tunnelling splitting in the formic acid dimer. Phys Chem Chem Phys 2017; 19:966-970. [DOI: 10.1039/c6cp07808g] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Nearly all degrees of freedom need to be included for accurate theoretical predictions of quantum dynamics.
Collapse
|
13
|
Fillaux F, Cousson A. A neutron diffraction study of the crystal of benzoic acid from 6 to 293 K and a macroscopic-scale quantum theory of the lattice of hydrogen-bonded dimers. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.09.016] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
14
|
Ivanov SD, Grant IM, Marx D. Quantum free energy landscapes from ab initio path integral metadynamics: Double proton transfer in the formic acid dimer is concerted but not correlated. J Chem Phys 2016; 143:124304. [PMID: 26429008 DOI: 10.1063/1.4931052] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
With the goal of computing quantum free energy landscapes of reactive (bio)chemical systems in multi-dimensional space, we combine the metadynamics technique for sampling potential energy surfaces with the ab initio path integral approach to treating nuclear quantum motion. This unified method is applied to the double proton transfer process in the formic acid dimer (FAD), in order to study the nuclear quantum effects at finite temperatures without imposing a one-dimensional reaction coordinate or reducing the dimensionality. Importantly, the ab initio path integral metadynamics technique allows one to treat the hydrogen bonds and concomitant proton transfers in FAD strictly independently and thus provides direct access to the much discussed issue of whether the double proton transfer proceeds via a stepwise or concerted mechanism. The quantum free energy landscape we compute for this H-bonded molecular complex reveals that the two protons move in a concerted fashion from initial to product state, yet world-line analysis of the quantum correlations demonstrates that the protons are as quantum-uncorrelated at the transition state as they are when close to the equilibrium structure.
Collapse
Affiliation(s)
- Sergei D Ivanov
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Ian M Grant
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
15
|
Qu C, Bowman JM. An ab initio potential energy surface for the formic acid dimer: zero-point energy, selected anharmonic fundamental energies, and ground-state tunneling splitting calculated in relaxed 1–4-mode subspaces. Phys Chem Chem Phys 2016; 18:24835-24840. [DOI: 10.1039/c6cp03073d] [Citation(s) in RCA: 59] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report a full-dimensional, permutationally invariant potential energy surface (PES) for the cyclic formic acid dimer.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation
- Emory University
- Atlanta
- USA
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation
- Emory University
- Atlanta
- USA
| |
Collapse
|
16
|
Ren Y, Bian W. Mode-Specific Tunneling Splittings for a Sequential Double-Hydrogen Transfer Case: An Accurate Quantum Mechanical Scheme. J Phys Chem Lett 2015; 6:1824-1829. [PMID: 26263255 DOI: 10.1021/acs.jpclett.5b00672] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We present the first accurate quantum dynamics calculations of mode-specific tunneling splittings in a sequential double-hydrogen transfer process. This is achieved in the vinylidene-acetylene system, the simplest molecular system of this kind, and by large-scale parallel computations with an efficient theoretical scheme developed by us. In our scheme, basis functions are customized for the hydrogen transfer process; a 4-dimensional basis contraction strategy is combined with the preconditioned inexact spectral transform method; efficient parallel implementation is achieved. Mode-specific permutation tunneling splittings of vinylidene states are reported and tremendous mode-specific promotion effects are revealed; in particular, the CH2 rock mode enhances the ground-state splitting by a factor of 10(3). We find that the ground-state vinylidene has a reversible-isomerization time of 622 ps, much longer than all previous estimates. Our calculations also shed light on the importance of the deep intermediate well and vibrational excitation in the double-hydrogen transfer processes.
Collapse
Affiliation(s)
- Yinghui Ren
- †Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- ‡University of Chinese Academy of Sciences, Beijing 100049, China
| | - Wensheng Bian
- †Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|