1
|
Eraković M, Cvitaš MT. Tunneling splittings in the vibrationally excited states of water trimer. Phys Chem Chem Phys 2024; 26:12965-12981. [PMID: 38634688 DOI: 10.1039/d4cp00013g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/19/2024]
Abstract
Tunneling splitting (TS) patterns in vibrationally excited states of the water trimer are calculated, taking into account six tunneling pathways that describe the flips of free OH bonds and five bifurcation mechanisms that break and reform hydrogen bonds in the trimer ring. A tunneling matrix (TM) model is used to derive the energy shifts due to tunneling in terms of the six distinct TM elements in symbolic form. TM elements are calculated using the recently-developed modified WKB (Wentzel-Kramers-Brillouin) method in full dimensionality. Convergence was achieved for the lowest six excited vibrational modes. Bifurcation widths of the pseudorotational quartets are shown to be of comparable size to the ground-state widths, obtained using instanton theory, or increased for some particular modes of vibration. The largest increase is obtained for the excited out-of-phase flip of two adjacent water monomers with free OH bonds pointing in opposite directions relative to the ring plane. Bifurcation widths in (D2O)3 are found to be two orders of magnitude smaller than in (H2O)3. Geometrical arguments were used to explain the order of states in some TS multiplets in vibrationally excited water trimers.
Collapse
Affiliation(s)
- Mihael Eraković
- Department of Physical Chemistry, Ruđer Bošković Institute, Zagreb, Croatia.
| | - Marko T Cvitaš
- Department of Physics, University of Zagreb Faculty of Science, Zagreb, Croatia.
| |
Collapse
|
2
|
Eraković M, Cvitaš MT. Tunneling splittings using modified WKB method in Cartesian coordinates: The test case of vinyl radical. J Chem Phys 2024; 160:154112. [PMID: 38639313 DOI: 10.1063/5.0204986] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2024] [Accepted: 03/31/2024] [Indexed: 04/20/2024] Open
Abstract
Modified WKB theory for calculating tunneling splittings in symmetric multi-well systems in full dimensionality is re-derived using Cartesian coordinates. It is explicitly shown that the theory rests on the wavefunction that is exact for harmonic potentials. The theory was applied to calculate tunneling splittings in vinyl radical and some of its deuterated isotopologues in their vibrational ground states and the low-lying vibrationally excited states and compared to exact variational results. The exact results are reproduced within a factor of 2 in most states. Remarkably, all large enhancements of tunneling splittings relative to the ground state, up to three orders in magnitude in some excited mode combinations, are well reproduced. It is also shown that in the asymmetrically deuterated vinyl radical, the theory correctly predicts the states that are localized in a single well and the delocalized tunneling states. Modified WKB theory on the minimum action path is computationally inexpensive and can also be applied without modification to much larger systems in full dimensionality; the results of this test case serve to give insight into the expected accuracy of the method.
Collapse
Affiliation(s)
- Mihael Eraković
- Department of Physical Chemistry, Ruđer Bošković Institute, Bijenička Cesta 54, 10000 Zagreb, Croatia
| | - Marko T Cvitaš
- Department of Physics, University of Zagreb Faculty of Science, Bijenička Cesta 32, 10000 Zagreb, Croatia
| |
Collapse
|
3
|
Trenins G, Meuser L, Bertschi H, Vavourakis O, Flütsch R, Richardson JO. Exact tunneling splittings from symmetrized path integrals. J Chem Phys 2023; 159:034108. [PMID: 37466233 DOI: 10.1063/5.0158879] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2023] [Accepted: 06/13/2023] [Indexed: 07/20/2023] Open
Abstract
We develop a new simulation technique based on path-integral molecular dynamics for calculating ground-state tunneling splitting patterns from ratios of symmetrized partition functions. In particular, molecular systems are rigorously projected onto their J = 0 rotational state by an "Eckart spring" that connects two adjacent beads in a ring polymer. Using this procedure, the tunneling splitting can be obtained from thermodynamic integration at just one (sufficiently low) temperature. Converged results are formally identical to the values that would have been obtained by solving the full rovibrational Schrödinger equation on a given Born-Oppenheimer potential energy surface. The new approach is showcased with simulations of hydronium and methanol, which are in good agreement with wavefunction-based calculations and experimental measurements. The method will be of particular use for the study of low-barrier methyl rotations and other floppy modes, where instanton theory is not valid.
Collapse
Affiliation(s)
- George Trenins
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Lars Meuser
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Hannah Bertschi
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Odysseas Vavourakis
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Reto Flütsch
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Jeremy O Richardson
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| |
Collapse
|
4
|
Eraković M, Cvitaš MT. Vibrational Tunneling Spectra of Molecules with Asymmetric Wells: A Combined Vibrational Configuration Interaction and Instanton Approach. J Chem Theory Comput 2022; 18:2785-2802. [PMID: 35439012 PMCID: PMC9097297 DOI: 10.1021/acs.jctc.2c00124] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
A combined approach
that uses the vibrational configuration interaction
(VCI) and semiclassical instanton theory was developed to study vibrational
tunneling spectra of molecules with multiple wells in full dimensionality.
The method can be applied to calculate low-lying vibrational states in the systems with an
arbitrary number of minima, which are not necessarily equal in energy
or shape. It was tested on a two-dimensional double-well model system
and on malonaldehyde, and the calculations reproduced the exact quantum
mechanical (QM) results with high accuracy. The method was subsequently
applied to calculate the vibrational spectrum of the asymmetrically
deuterated malonaldehyde with nondegenerate vibrational frequencies
in the two wells. The spectrum is obtained at a cost of single-well
VCI calculations used to calculate the local energies. The interactions
between states of different wells are computed semiclassically using
the instanton theory at a comparatively negligible computational cost.
The method is particularly suited to systems in which the wells are
separated by large potential barriers and tunneling splittings are
small, for example, in some water clusters, when the exact QM methods
come at a prohibitive computational cost.
Collapse
Affiliation(s)
- Mihael Eraković
- Department of Physical Chemistry, Rud̵er Bošković Institute, Bijenička Cesta 54, 10000 Zagreb, Croatia
| | - Marko T Cvitaš
- Department of Physics, Faculty of Science, University of Zagreb, Bijenička Cesta 32, 10000 Zagreb, Croatia
| |
Collapse
|
5
|
Eraković M, Cvitaš MT. Tunnelling splitting patterns in some partially deuterated water trimers. Phys Chem Chem Phys 2021; 23:4240-4254. [PMID: 33586727 DOI: 10.1039/d0cp06135b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We apply our recently developed semiclassical method for calculating tunnelling splittings (TS) in asymmetric systems to make the first characterization of the ground-state TS pattern of some partially deuterated water trimers. Similarly to homoisotopic water trimers, the ground-state TS patterns are explained in terms of six distinct rearrangement mechanisms. TS patterns in (D2O)(H2O)2 and (H2O)(D2O)2 are composed of sextets induced by the dynamics of flips, and each of its levels is further finely split into a quartet of doublets and a doublet of quartets, respectively, due to various bifurcation dynamics. The TS pattern is obtained using 18 distinct tunnelling matrix elements. TS patterns of (HOD)(H2O)2 and (HOD)(D2O)2 each consists of two sextets, belonging to in-bond and out-of-bond substituted isomers. These sextet levels are further split into quartets by bifurcations. The TS pattern is computed in terms of 13 matrix elements. We also derive analytic expressions for bifurcation tunnelling splittings in terms of tunnelling matrix elements using symmetry. The present approach can be applied to other water clusters and also to the low-lying vibrationally excited states and should help in the interpretation and assignment of experimental spectra in the future.
Collapse
Affiliation(s)
- Mihael Eraković
- Department of Physical Chemistry, Ruđer Bošković Institute, 10000 Zagreb, Croatia.
| | - Marko T Cvitaš
- Department of Physical Chemistry, Ruđer Bošković Institute, 10000 Zagreb, Croatia.
| |
Collapse
|
6
|
Eraković M, Cvitaš MT. Tunneling splittings of vibrationally excited states using general instanton paths. J Chem Phys 2020; 153:134106. [PMID: 33032414 DOI: 10.1063/5.0024210] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
A multidimensional semiclassical method for calculating tunneling splittings in vibrationally excited states of molecules using Cartesian coordinates is developed. It is an extension of the theory by Mil'nikov and Nakamura [J. Chem. Phys. 122, 124311 (2005)] to asymmetric paths that are necessary for calculating tunneling splitting patterns in multi-well systems, such as water clusters. Additionally, new terms are introduced in the description of the semiclassical wavefunction that drastically improves the splitting estimates for certain systems. The method is based on the instanton theory and builds the semiclassical wavefunction of the vibrationally excited states from the ground-state instanton wavefunction along the minimum action path and its harmonic neighborhood. The splittings of excited states are thus obtained at a negligible added numerical effort. The cost is concentrated, as for the ground-state splittings, in the instanton path optimization and the hessian evaluation along the path. The method can thus be applied without modification to many mid-sized molecules in full dimensionality and in combination with on-the-fly evaluation of electronic potentials. The tests were performed on several model potentials and on the water dimer.
Collapse
Affiliation(s)
- Mihael Eraković
- Department of Physical Chemistry, Ruder Bošković Institute, Bijenička Cesta 54, 10000 Zagreb, Croatia
| | - Marko T Cvitaš
- Department of Physical Chemistry, Ruder Bošković Institute, Bijenička Cesta 54, 10000 Zagreb, Croatia
| |
Collapse
|
7
|
Jahr E, Laude G, Richardson JO. Instanton theory of tunneling in molecules with asymmetric isotopic substitutions. J Chem Phys 2020; 153:094101. [PMID: 32891112 DOI: 10.1063/5.0021831] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We consider quantum tunneling in asymmetric double-well systems for which the local minima in the two wells have the same energy, but the frequencies differ slightly. In a molecular context, this situation can arise if the symmetry is broken by isotopic substitutions. We derive a generalization of instanton theory for these asymmetric systems, leading to a semiclassical expression for the tunneling matrix element and hence the energy-level splitting. We benchmark the method using a set of one- and two-dimensional models, for which the results compare favorably with numerically exact quantum calculations. Using the ring-polymer instanton approach, we apply the method to compute the level splittings in various isotopomers of malonaldehyde in full dimensionality and analyze the relative contributions from the zero-point energy difference and tunneling effects.
Collapse
Affiliation(s)
- Elena Jahr
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Gabriel Laude
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
8
|
Eraković M, Vaillant CL, Cvitaš MT. Instanton theory of ground-state tunneling splittings with general paths. J Chem Phys 2020; 152:084111. [PMID: 32113369 DOI: 10.1063/1.5145278] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We derive a multidimensional instanton theory for calculating ground-state tunneling splittings in Cartesian coordinates for general paths. It is an extension of the method by Mil'nikov and Nakamura [J. Chem. Phys. 115, 6881 (2001)] to include asymmetric paths that are necessary for calculating tunneling splitting patterns in multi-well systems, such as water clusters. The approach avoids multiple expensive matrix diagonalizations to converge the fluctuation prefactor in the ring-polymer instanton (RPI) method, and instead replaces them by an integration of a Riccati differential equation. When combined with the string method for locating instantons, we avoid the need to converge the calculation with respect to the imaginary time period of the semiclassical orbit, thereby reducing the number of convergence parameters of the optimized object to just one: the number of equally spaced system replicas used to represent the instanton path. The entirety of the numerical effort is thus concentrated in optimizing the shape of the path and evaluating hessians along the path, which is a dramatic improvement over RPI. In addition to the standard instanton approximations, we neglect the coupling of vibrational modes to external rotations. The method is tested on the model potential of malonaldehyde and on the water dimer and trimer, giving close agreement with RPI at a much-reduced cost.
Collapse
Affiliation(s)
- Mihael Eraković
- Department of Physical Chemistry, Rudđđer Bošković Institute, Bijenička Cesta 54, 10000 Zagreb, Croatia
| | - Christophe L Vaillant
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Marko T Cvitaš
- Department of Physical Chemistry, Rudđđer Bošković Institute, Bijenička Cesta 54, 10000 Zagreb, Croatia
| |
Collapse
|
9
|
Cvitaš MT, Richardson JO. Quantum tunnelling pathways of the water pentamer. Phys Chem Chem Phys 2020; 22:1035-1044. [DOI: 10.1039/c9cp05561d] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Five tunnelling rearrangement pathways in water pentamer are responsible for the ground-state tunnelling splitting pattern of 320 states.
Collapse
Affiliation(s)
- Marko T. Cvitaš
- Department of Physical Chemistry
- Ruđer Bošković Institute
- Croatia
| | | |
Collapse
|
10
|
Vaillant CL, Wales DJ, Althorpe SC. Tunneling Splittings in Water Clusters from Path Integral Molecular Dynamics. J Phys Chem Lett 2019; 10:7300-7304. [PMID: 31682130 DOI: 10.1021/acs.jpclett.9b02951] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
We present calculations of tunneling splittings in selected small water clusters, based on a recently developed path integral molecular dynamics (PIMD) method. The ground-rotational-state tunneling motions associated with the largest splittings in the water dimer, trimer, and hexamer are considered, and we show that the PIMD predictions are in very good agreement with benchmark quantum and experimental results. As the tunneling spectra are highly sensitive to both the details of the quantum dynamics and the potential energy surface, our calculations are a validation of the MB-Pol surface as well as the accuracy of PIMD. The favorable scaling of PIMD with system size paves the way for calculations of tunneling splittings in large, nonrigid molecular systems with motions that cannot be treated accurately by other methods, such as the semiclassical instanton.
Collapse
Affiliation(s)
- C L Vaillant
- Laboratory of Theoretical Physical Chemistry , Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne , Switzerland
| | - D J Wales
- Department of Chemistry , University of Cambridge , Lensfield Road , Cambridge CB2 1EW , United Kingdom
| | - S C Althorpe
- Department of Chemistry , University of Cambridge , Lensfield Road , Cambridge CB2 1EW , United Kingdom
| |
Collapse
|
11
|
Thapa MJ, Fang W, Richardson JO. Nonadiabatic quantum transition-state theory in the golden-rule limit. I. Theory and application to model systems. J Chem Phys 2019; 150:104107. [DOI: 10.1063/1.5081108] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Manish J. Thapa
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Wei Fang
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
12
|
Vaillant CL, Althorpe SC, Wales DJ. Path Integral Energy Landscapes for Water Clusters. J Chem Theory Comput 2019; 15:33-42. [PMID: 30550261 DOI: 10.1021/acs.jctc.8b00675] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
The energy landscapes for a discretized path integral representation of the water dimer, trimer and pentamer are characterized in terms of the localized (classical) and delocalized minima and transition states. The transition states are finite-temperature approximations to the exact instanton path, and they are typically used to calculate the tunneling splittings or reaction rates. The features of the path integral landscape are explored, thus elucidating procedures that could usefully be automated when searching for instantons in larger systems. Our work not only clarifies the role of minima and transition states in path integral calculations but also enables us to analyze the quantum-to-classical transition.
Collapse
Affiliation(s)
- Christophe L Vaillant
- University Chemical Laboratories , Lensfield Road , Cambridge CB2 1EW , United Kingdom
| | - Stuart C Althorpe
- University Chemical Laboratories , Lensfield Road , Cambridge CB2 1EW , United Kingdom
| | - David J Wales
- University Chemical Laboratories , Lensfield Road , Cambridge CB2 1EW , United Kingdom
| |
Collapse
|
13
|
Vaillant CL, Wales DJ, Althorpe SC. Tunneling splittings from path-integral molecular dynamics using a Langevin thermostat. J Chem Phys 2018; 148:234102. [DOI: 10.1063/1.5029258] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- C. L. Vaillant
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - D. J. Wales
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - S. C. Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
14
|
|
15
|
|
16
|
Vaillant CL, Cvitaš MT. Rotation-tunneling spectrum of the water dimer from instanton theory. Phys Chem Chem Phys 2018; 20:26809-26813. [DOI: 10.1039/c8cp04991b] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A linearly-scaling path-integral method is developed to calculate rotation-tunneling spectra, and is applied to the water dimer.
Collapse
Affiliation(s)
| | - Marko T. Cvitaš
- Department of Physical Chemistry, Ruđer Bošković Institute
- 10000 Zagreb
- Croatia
| |
Collapse
|