1
|
Xue Y, Sexton TM, Yang J, Tschumper GS. Systematic analysis of electronic barrier heights and widths for concerted proton transfer in cyclic hydrogen bonded clusters: (HF) n, (HCl) n and (H 2O) n where n = 3, 4, 5. Phys Chem Chem Phys 2024; 26:12483-12494. [PMID: 38619858 DOI: 10.1039/d4cp00422a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/16/2024]
Abstract
The MP2 and CCSD(T) methods are paired with correlation consistent basis sets as large as aug-cc-pVQZ to optimize the structures of the cyclic minima for (HF)n, (HCl)n and (H2O)n where n = 3-5, as well as the corresponding transition states (TSs) for concerted proton transfer (CPT). MP2 and CCSD(T) harmonic vibrational frequencies confirm the nature of each minimum and TS. Both conventional and explicitly correlated CCSD(T) computations are employed to assess the electronic dissociation energies and barrier heights for CPT near the complete basis (CBS) limit for all 9 clusters. Results for (HF)n are consistent with prior studies identifying Cnh and Dnh point group symmetry for the minima and TSs, respectively. Our computations also confirm that CPT proceeds through Cs TS structures for the C1 minima of (H2O)3 and (H2O)5, whereas the process goes through a TS with D2d symmetry for the S4 global minimum of (H2O)4. This work corroborates earlier findings that the minima for (HCl)3, (HCl)4 and (HCl)5 have C3h, S4 and C1 point group symmetry, respectively, and that the Cnh structures are not minima for n = 4 and 5. Moreover, our computations show the TSs for CPT in (HCl)3, (HCl)4 and (HCl)5 have D3h, D2d, and C2 point group symmetry, respectively. At the CCSD(T) CBS limit, (HF)4 and (HF)5 have the smallest electronic barrier heights for CPT (≈15 kcal mol-1 for both), followed by the HF trimer (≈21 kcal mol-1). The barriers are appreciably higher for the other clusters (around 27 kcal mol-1 for (H2O)4 and (HCl)3; roughly 30 kcal mol-1 for (H2O)3, (H2O)5 and (HCl)4; up to 38 kcal mol-1 for (HCl)5). At the CBS limit, MP2 significantly underestimates the CCSD(T) barrier heights (e.g., by ca. 2, 4 and 7 kcal mol-1 for the pentamers of HF, H2O and HCl, respectively), whereas CCSD overestimates these barriers by roughly the same magnitude. Scaling the barrier heights and dissociation energies by the number of fragments in the cluster reveals strong linear relationships between the two quantities and with the magnitudes of the imaginary vibrational frequency for the TSs.
Collapse
Affiliation(s)
- Yuan Xue
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS 38677-1848, USA.
| | - Thomas More Sexton
- School of Arts and Sciences, Chemistry University of Mary, Bismark, ND 58504, USA.
| | - Johnny Yang
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS 38677-1848, USA.
| | - Gregory S Tschumper
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS 38677-1848, USA.
| |
Collapse
|
2
|
Zane AR, Curotto E. Electrolyte clusters as hydrogen sponges: diffusion Monte Carlo simulations. Phys Chem Chem Phys 2022; 24:26094-26101. [PMID: 36268906 DOI: 10.1039/d2cp03658d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
We carry out Diffusion Monte Carlo simulations of up to five hydrogen molecules aggregated with two Stockmayer clusters that solvate a single lithium ion. The first one contains six point dipole solvent particles with parameters tuned to emulate nitromethane. The second cluster is a relative large system investigated recently [G. DiEmma, S. Kalette, and E. Curotto, Chem. Phys. Lett. 2019, 725, 80-86]. In both cases we find that the aggregated hydrogen molecules perturb significantly the ground state of the host cluster and form a distorted tetrahedral cage around the Li+ ion. The fifth hydrogen molecule is absorbed by the larger Stockmayer cluster while remaining in the proximity of the solvated charge.
Collapse
Affiliation(s)
- A R Zane
- Department of Chemistry & Physics Arcadia University, 450 S. Easton Rd., Glenside, PA 19038, USA.
| | - E Curotto
- Department of Chemistry & Physics Arcadia University, 450 S. Easton Rd., Glenside, PA 19038, USA.
| |
Collapse
|
3
|
Patkar D, Ahirwar MB, Shrivastava SP, Deshmukh MM. Assessment of hydrogen bond strengths and cooperativity in self- and cross-associating cyclic (HF)m(H2O)n (m + n = 2 to 8) clusters. NEW J CHEM 2022. [DOI: 10.1039/d1nj05431g] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
In this work, we investigated the strengths of various self- and cross-associating hydrogen bonds (HBs) in mixed hydrogen fluoride–water cyclic (HF)m(H2O)n (m + n = 2 to 8) clusters, employing a molecular tailoring approach (MTA)-based method.
Collapse
Affiliation(s)
- Deepak Patkar
- Department of Chemistry, Dr Harisingh Gour Vishwavidyalaya (A Central University), Sagar, 470003, India
| | - Mini Bharati Ahirwar
- Department of Chemistry, Dr Harisingh Gour Vishwavidyalaya (A Central University), Sagar, 470003, India
| | - Satya Prakash Shrivastava
- Department of Chemistry, Dr Harisingh Gour Vishwavidyalaya (A Central University), Sagar, 470003, India
| | - Milind M. Deshmukh
- Department of Chemistry, Dr Harisingh Gour Vishwavidyalaya (A Central University), Sagar, 470003, India
| |
Collapse
|
4
|
Patkar D, Ahirwar MB, Gadre SR, Deshmukh MM. Unusually Large Hydrogen-Bond Cooperativity in Hydrogen Fluoride Clusters, (HF) n, n = 3 to 8, Revealed by the Molecular Tailoring Approach. J Phys Chem A 2021; 125:8836-8845. [PMID: 34612647 DOI: 10.1021/acs.jpca.1c06478] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this work, our recently proposed molecular tailoring approach (MTA)-based method is employed for the evaluation of individual hydrogen-bond (HB) energies in linear (L) and cyclic (C) hydrogen fluoride clusters, (HF)n (n = 3 to 8). The estimated individual HB energies calculated at the MP2(full)/aug-cc-pVTZ level for the L-(HF)n are between 6.2 to 9.5 kcal/mol and those in the C-(HF)n lie between 7.9 to 11.4 kcal/mol. The zero-point energy corrections and basis set superposition corrections are found to be very small (less than 0.6 and 1.2 kcal/mol, respectively). The cooperativity contribution toward individual HBs is seen to fall between 1.0 to 4.8 kcal/mol and 3.2 to 6.9 kcal/mol for linear and cyclic clusters, respectively. Interestingly, the HB energies in dimers, cleaved from these clusters, lie in a narrow range (4.4 to 5.2 kcal/mol) suggesting that the large HB strength in (HF)n clusters is mainly due to the large cooperativity contribution, especially for n ≥ 5 (50 to 62% of the HBs energy). Furthermore, the HB energies in these clusters show a good qualitative correlation with geometrical parameters (H···F distance and F-H···F angles), stretching frequencies of F-H bonds, and electron density values at the (3, -1) bond critical points.
Collapse
Affiliation(s)
- Deepak Patkar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), Sagar, 470003, India
| | - Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), Sagar, 470003, India
| | - Shridhar R Gadre
- Department of Scientific Computing, Modelling and Simulation, Savitribai Phule Pune University, Pune, 411 007, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), Sagar, 470003, India
| |
Collapse
|
5
|
von Domaros M, Perlt E. Anharmonic effects in the quantum cluster equilibrium method. J Chem Phys 2017; 146:124114. [PMID: 28388115 DOI: 10.1063/1.4978958] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The well-established quantum cluster equilibrium (QCE) model provides a statistical thermodynamic framework to apply high-level ab initio calculations of finite cluster structures to macroscopic liquid phases using the partition function. So far, the harmonic approximation has been applied throughout the calculations. In this article, we apply an important correction in the evaluation of the one-particle partition function and account for anharmonicity. Therefore, we implemented an analytical approximation to the Morse partition function and the derivatives of its logarithm with respect to temperature, which are required for the evaluation of thermodynamic quantities. This anharmonic QCE approach has been applied to liquid hydrogen chloride and cluster distributions, and the molar volume, the volumetric thermal expansion coefficient, and the isobaric heat capacity have been calculated. An improved description for all properties is observed if anharmonic effects are considered.
Collapse
Affiliation(s)
- Michael von Domaros
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 6, D-53115 Bonn, Germany
| | - Eva Perlt
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 6, D-53115 Bonn, Germany
| |
Collapse
|
6
|
Samanta AK, Wang Y, Mancini JS, Bowman JM, Reisler H. Energetics and Predissociation Dynamics of Small Water, HCl, and Mixed HCl–Water Clusters. Chem Rev 2016; 116:4913-36. [DOI: 10.1021/acs.chemrev.5b00506] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Amit K. Samanta
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| | - Yimin Wang
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - John S. Mancini
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Hanna Reisler
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| |
Collapse
|
7
|
Curotto E. Ion-Stockmayer clusters: Minima, classical thermodynamics, and variational ground state estimates of Li+(CH3NO2)n(n= 1–20). J Chem Phys 2015; 143:214301. [DOI: 10.1063/1.4936587] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- E. Curotto
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| |
Collapse
|
8
|
Curotto E, Mella M. On the convergence of diffusion Monte Carlo in non-Euclidean spaces. II. Diffusion with sources and sinks. J Chem Phys 2015; 142:114111. [PMID: 25796235 DOI: 10.1063/1.4914516] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We test the second order Milstein method adapted to simulate diffusion in general compact Riemann manifolds on a number of systems characterized by nonconfining potential energy surfaces of increasing complexity. For the 2-sphere and more complex spaces derived from it, we compare the Milstein method with a number of other first and second order approaches. In each case tested, we find evidence that demonstrate the versatility and relative ease of implementation of the Milstein method derived in Part I.
Collapse
Affiliation(s)
- E Curotto
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | - Massimo Mella
- Dipartimento di Scienza ed Alta Tecnologia, via Valleggio 9, Università degli studi dell'Insubria, 22100 Como, Italy
| |
Collapse
|
9
|
Wolf S, Curotto E. Ring polymer dynamics for rigid tops with an improved integrator. J Chem Phys 2014; 141:024116. [DOI: 10.1063/1.4887460] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Mancini JS, Bowman JM. A New Many-Body Potential Energy Surface for HCl Clusters and Its Application to Anharmonic Spectroscopy and Vibration–Vibration Energy Transfer in the HCl Trimer. J Phys Chem A 2014; 118:7367-74. [DOI: 10.1021/jp412264t] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- John S. Mancini
- Cherry
L. Emerson Center
for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Cherry
L. Emerson Center
for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
11
|
|
12
|
Mancini JS, Bowman JM. On-the-fly ab intito calculations of anharmonic vibrational frequencies: Local-monomer theory and application to HCl clusters. J Chem Phys 2013; 139:164115. [DOI: 10.1063/1.4826351] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
13
|
Samanta AK, Ch’ng LC, Reisler H. Imaging bond breaking and vibrational energy transfer in small water containing clusters. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2013.05.003] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
14
|
Curotto E, Mella M. Quantum Monte Carlo simulations of selected ammonia clusters (n = 2–5): Isotope effects on the ground state of typical hydrogen bonded systems. J Chem Phys 2010; 133:214301. [DOI: 10.1063/1.3506027] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
|
15
|
Asare E, Musah AR, Curotto E, Freeman DL, Doll JD. The thermodynamic and ground state properties of the TIP4P water octamer. J Chem Phys 2009; 131:184508. [DOI: 10.1063/1.3259047] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
Lubombo C, Curotto E, Janeiro Barral PE, Mella M. Thermodynamic properties of ammonia clusters (NH3)n n=2–11: Comparing classical and quantum simulation results for hydrogen bonded species. J Chem Phys 2009; 131:034312. [DOI: 10.1063/1.3159398] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|