1
|
Moonkaen P, Finney JM, McCoy AB. Isotope Effects on Ground and Excited States of Ethyl Cation, H +(C 2H 4). J Phys Chem A 2023; 127:1196-1205. [PMID: 36705480 DOI: 10.1021/acs.jpca.2c07334] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The structure and spectra of ethyl cation, H+(C2H4), and its deuterated analogues are investigated using diffusion Monte Carlo (DMC). These calculations all show that the ground state wave function for H+(C2H4) is localized near the minimum energy configuration in which the excess proton is in a bridging configuration, although the amplitude of the vibrational motions of the bridging proton is large. Deuteration of the bridging proton reduces the amplitude of this motion, while deuteration of only the ethylenic hydrogen atoms in H+(C2D4) has little effect on the amplitude of the motion of the bridging proton. Excited states that are accessed by spectroscopically observed transitions in H+(C2H4) are calculated using fixed-node DMC. The calculated and measured frequencies for the states with one quantum of excitation in the ethylenic CH stretching vibrations show good agreement. We also explore the excited state with one quantum of excitation in the proton transfer vibration of the bridging proton and obtain a frequency of 616 cm-1 for H+(C2H4). This frequency increases to 629 cm-1 in H+(C2D4). Deuteration decreases this frequency to 491 and 495 cm-1 in D+(C2H4) and D+(C2D4), respectively. The effects of partial deuteration on the frequencies of the CH stretching vibrations, and the corresponding probability amplitudes are also explored. Finally, we report the vibrationally averaged rotational constants for the four isotopologues of ethyl cation considered in this study.
Collapse
Affiliation(s)
- Pattarapon Moonkaen
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Jacob M Finney
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
2
|
DiRisio RJ, Finney JM, McCoy AB. Diffusion Monte Carlo approaches for studying nuclear quantum effects in fluxional molecules. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1615] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Ryan J. DiRisio
- Department of Chemistry University of Washington Seattle Washington USA
| | - Jacob M. Finney
- Department of Chemistry University of Washington Seattle Washington USA
| | - Anne B. McCoy
- Department of Chemistry University of Washington Seattle Washington USA
| |
Collapse
|
3
|
Puzzarini C, Bloino J, Tasinato N, Barone V. Accuracy and Interpretability: The Devil and the Holy Grail. New Routes across Old Boundaries in Computational Spectroscopy. Chem Rev 2019; 119:8131-8191. [DOI: 10.1021/acs.chemrev.9b00007] [Citation(s) in RCA: 114] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
4
|
Fore ME, McCoy AB. Statistical Analysis of the Effect of Deuteration on Quantum Delocalization in CH5+. J Phys Chem A 2019; 123:4623-4631. [DOI: 10.1021/acs.jpca.9b02685] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Meredith E. Fore
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B. McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
5
|
Curotto E, Mella M. Diffusion Monte Carlo simulations of gas phase and adsorbed D 2-(H 2) n clusters. J Chem Phys 2018; 148:102315. [PMID: 29544319 DOI: 10.1063/1.5000372] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We have computed ground state energies and analyzed radial distributions for several gas phase and adsorbed D2(H2)n and HD(H2)n clusters. An external model potential designed to mimic ionic adsorption sites inside porous materials is used [M. Mella and E. Curotto, J. Phys. Chem. A 121, 5005 (2017)]. The isotopic substitution lowers the ground state energies by the expected amount based on the mass differences when these are compared with the energies of the pure clusters in the gas phase. A similar impact is found for adsorbed aggregates. The dissociation energy of D2 from the adsorbed clusters is always much higher than that of H2 from both pure and doped aggregates. Radial distributions of D2 and H2 are compared for both the gas phase and adsorbed species. For the gas phase clusters, two types of hydrogen-hydrogen interactions are considered: one based on the assumption that rotations and translations are adiabatically decoupled and the other based on nonisotropic four-dimensional potential. In the gas phase clusters of sufficiently large size, we find the heavier isotopomer more likely to be near the center of mass. However, there is a considerable overlap among the radial distributions of the two species. For the adsorbed clusters, we invariably find the heavy isotope located closer to the attractive interaction source than H2, and at the periphery of the aggregate, H2 molecules being substantially excluded from the interaction with the source. This finding rationalizes the dissociation energy results. For D2-(H2)n clusters with n≥12, such preference leads to the desorption of D2 from the aggregate, a phenomenon driven by the minimization of the total energy that can be obtained by reducing the confinement of (H2)12. The same happens for (H2)13, indicating that such an effect may be quite general and impact on the absorption of quantum species inside porous materials.
Collapse
Affiliation(s)
- E Curotto
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | - M Mella
- Dipartimento di Scienza ed Alta Tecnologia, Università degli Studi dell'Insubria, Via Valleggio 11, 22100 Como, Italy
| |
Collapse
|
6
|
Bulik IW, Frisch MJ, Vaccaro PH. Fixed-Node, Importance-Sampling Diffusion Monte Carlo for Vibrational Structure with Accurate and Compact Trial States. J Chem Theory Comput 2018; 14:1554-1563. [DOI: 10.1021/acs.jctc.8b00016] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ireneusz W. Bulik
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Michael J. Frisch
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, United States
| | - Patrick H. Vaccaro
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| |
Collapse
|
7
|
Bulik IW, Frisch MJ, Vaccaro PH. Vibrational self-consistent field theory using optimized curvilinear coordinates. J Chem Phys 2017; 147:044110. [DOI: 10.1063/1.4995440] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ireneusz W. Bulik
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, USA
| | - Michael J. Frisch
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| | - Patrick H. Vaccaro
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, USA
| |
Collapse
|
8
|
Hochlaf M. Advances in spectroscopy and dynamics of small and medium sized molecules and clusters. Phys Chem Chem Phys 2017; 19:21236-21261. [DOI: 10.1039/c7cp01980g] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Investigations of the spectroscopy and dynamics of small- and medium-sized molecules and clusters represent a hot topic in atmospheric chemistry, biology, physics, atto- and femto-chemistry and astrophysics.
Collapse
Affiliation(s)
- Majdi Hochlaf
- Université Paris-Est
- Laboratoire Modélisation et Simulation Multi Echelle
- MSME UMR 8208 CNRS
- 77454 Marne-la-Vallée
- France
| |
Collapse
|
9
|
Christensen HM, Jake LC, Curotto E. Smart darting diffusion Monte Carlo: Applications to lithium ion-Stockmayer clusters. J Chem Phys 2016; 144:174115. [PMID: 27155633 DOI: 10.1063/1.4948562] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In a recent investigation [K. Roberts et al., J. Chem. Phys. 136, 074104 (2012)], we have shown that, for a sufficiently complex potential, the Diffusion Monte Carlo (DMC) random walk can become quasiergodic, and we have introduced smart darting-like moves to improve the sampling. In this article, we systematically characterize the bias that smart darting moves introduce in the estimate of the ground state energy of a bosonic system. We then test a simple approach to eliminate completely such bias from the results. The approach is applied for the determination of the ground state of lithium ion-n-dipoles clusters in the n = 8-20 range. For these, the smart darting diffusion Monte Carlo simulations find the same ground state energy and mixed-distribution as the traditional approach for n < 14. In larger systems we find that while the ground state energies agree quantitatively with or without smart darting moves, the mixed-distributions can be significantly different. Some evidence is offered to conclude that introducing smart darting-like moves in traditional DMC simulations may produce a more reliable ground state mixed-distribution.
Collapse
Affiliation(s)
- H M Christensen
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | - L C Jake
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | - E Curotto
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| |
Collapse
|
10
|
Calculating rovibrationally excited states of H2D+ and HD2+ by combination of fixed node and multi-state rotational diffusion Monte Carlo. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2015.12.030] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
11
|
Lin Z, McCoy AB. Probing the Relationship Between Large-Amplitude Motions in H5(+) and Proton Exchange Between H3(+) and H2. J Phys Chem A 2015; 119:12109-18. [PMID: 26244451 DOI: 10.1021/acs.jpca.5b05774] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding the spectroscopy and dynamics of H5(+) is central in gaining insights into the H3(+) + H2 → H5(+) → H2 + H3(+) proton transfer reaction. This molecular ion exhibits large-amplitude vibrations, which allow for the transfer of a proton between H3(+) and H2 even in its ground vibrational state. With vibrational excitation, the number of open channels for permutations of protons increase. In this work, the minimized energy path variant of diffusion Monte Carlo is used to investigate how the energetically accessible proton permutations evolve as H5(+) is dissociated into H3(+) + H2. Two mechanisms for proton permutation are investigated. The first is the proton hop, which correlates to large-amplitude vibrations of the central proton in H5(+). The second is the exchange of a pair of hydrogen atoms between H3(+) and H2. This mechanism requires several proton hops along with a 120° rotation of H3(+) within the H5(+) molecular ion. This analysis shows that while there is a narrow region of configuration space over which both isomerization processes are energetically accessible, full permutation of the five protons in H5(+) more likely occurs through a stepwise mechanism. Such full permutation of the protons becomes accessible when the shared proton stretch is excited to the vpt = 2 or 3 excited state. The effects of deuteration and rotational excitation of the H2 and H3(+) products are also investigated. Deuteration inhibits permutation of protons, while rotational excitation has only a small impact on these processes.
Collapse
Affiliation(s)
- Zhou Lin
- Department of Chemistry and Biochemistry, The Ohio State University , Columbus, Ohio 43210, United States
| | - Anne B McCoy
- Department of Chemistry and Biochemistry, The Ohio State University , Columbus, Ohio 43210, United States
| |
Collapse
|
12
|
Mallory JD, Brown SE, Mandelshtam VA. Assessing the Performance of the Diffusion Monte Carlo Method As Applied to the Water Monomer, Dimer, and Hexamer. J Phys Chem A 2015; 119:6504-15. [PMID: 26001418 DOI: 10.1021/acs.jpca.5b02511] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The diffusion Monte Carlo (DMC) method is applied to the water monomer, dimer, and hexamer using q-TIP4P/F, one of the most simple empirical water models with flexible monomers. The bias in the time step (Δτ) and population size (Nw) is investigated. For the binding energies, the bias in Δτ cancels nearly completely, whereas a noticeable bias in Nw remains. However, for the isotope shift (e.g, in the dimer binding energies between (H2O)2 and (D2O)2), the systematic errors in Nw do cancel. Consequently, very accurate results for the latter (within ∼0.01 kcal/mol) are obtained with moderate numerical effort (Nw ∼ 10(3)). For the water hexamer and its (D2O)6 isotopomer, the DMC results as a function of Nw are examined for the cage and prism isomers. For a given isomer, the issue of the walker population leaking out of the corresponding basin of attraction is addressed by using appropriate geometric constraints. The population size bias for the hexamer is more severe, and to maintain accuracy similar to that of the dimer, Nw must be increased by ∼2 orders of magnitude. Fortunately, when the energy difference between the cage and prism is taken, the biases cancel, thereby reducing the systematic errors to within ∼0.01 kcal/mol when using a population of Nw = 4.8 × 10(5) walkers. Consequently, a very accurate result for the isotope shift is also obtained. Notably, both the quantum and isotope effects for the prism-cage energy difference are small.
Collapse
Affiliation(s)
- Joel D Mallory
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697, United States
| | - Sandra E Brown
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697, United States
| | - Vladimir A Mandelshtam
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697, United States
| |
Collapse
|
13
|
Hermes MR, Hirata S. Stochastic many-body perturbation theory for anharmonic molecular vibrations. J Chem Phys 2015; 141:084105. [PMID: 25173003 DOI: 10.1063/1.4892614] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new quantum Monte Carlo (QMC) method for anharmonic vibrational zero-point energies and transition frequencies is developed, which combines the diagrammatic vibrational many-body perturbation theory based on the Dyson equation with Monte Carlo integration. The infinite sums of the diagrammatic and thus size-consistent first- and second-order anharmonic corrections to the energy and self-energy are expressed as sums of a few m- or 2m-dimensional integrals of wave functions and a potential energy surface (PES) (m is the vibrational degrees of freedom). Each of these integrals is computed as the integrand (including the value of the PES) divided by the value of a judiciously chosen weight function evaluated on demand at geometries distributed randomly but according to the weight function via the Metropolis algorithm. In this way, the method completely avoids cumbersome evaluation and storage of high-order force constants necessary in the original formulation of the vibrational perturbation theory; it furthermore allows even higher-order force constants essentially up to an infinite order to be taken into account in a scalable, memory-efficient algorithm. The diagrammatic contributions to the frequency-dependent self-energies that are stochastically evaluated at discrete frequencies can be reliably interpolated, allowing the self-consistent solutions to the Dyson equation to be obtained. This method, therefore, can compute directly and stochastically the transition frequencies of fundamentals and overtones as well as their relative intensities as pole strengths, without fixed-node errors that plague some QMC. It is shown that, for an identical PES, the new method reproduces the correct deterministic values of the energies and frequencies within a few cm(-1) and pole strengths within a few thousandths. With the values of a PES evaluated on the fly at random geometries, the new method captures a noticeably greater proportion of anharmonic effects.
Collapse
Affiliation(s)
- Matthew R Hermes
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| |
Collapse
|
14
|
On the fly nodal searches in importance sampled fixed-node diffusion Monte Carlo using a parallel, fine-grained, genetic algorithm. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2014.11.054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
15
|
Samanta K, Beames JM, Lester MI, Subotnik JE. Quantum dynamical investigation of the simplest Criegee intermediate CH2OO and its O–O photodissociation channels. J Chem Phys 2014; 141:134303. [DOI: 10.1063/1.4894746] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Affiliation(s)
- Kousik Samanta
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA
| | - Joseph M. Beames
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA
| | - Marsha I. Lester
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA
| |
Collapse
|
16
|
Lin Z, McCoy AB. The role of large-amplitude motions in the spectroscopy and dynamics of H5+. J Chem Phys 2014; 140:114305. [DOI: 10.1063/1.4868098] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
|