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
|
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
|
3
|
Higashi M, Truhlar DG. Combined Electrostatically Embedded Multiconfiguration Molecular Mechanics and Molecular Mechanical Method: Application to Molecular Dynamics Simulation of a Chemical Reaction in Aqueous Solution with Hybrid Density Functional Theory. J Chem Theory Comput 2015; 4:1032-9. [PMID: 26636356 DOI: 10.1021/ct8000816] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We here combine the electrostatically embedded multiconfiguration molecular mechanics (EE-MCMM) method for generating global potential energy surfaces in the presence of an electrostatic potential with molecular mechanics (MM). The resulting EE-MCMM/MM method is illustrated by applying it to carry out a molecular dynamics simulation for the symmetric bimolecular reaction Cl(-) + CH3Cl' → ClCH3 + Cl'(-) in aqueous solution with hybrid density functional theory as the quantum mechanical level. The potential of mean force is calculated, and the free energy barrier is found to be 25.3 kcal/mol, which is in good agreement with previous work. The advantage of the combined EE-MCMM and MM method is that the number of quantum mechanical calculations required for the active subsystem is very small compared to straight direct dynamics.
Collapse
Affiliation(s)
- Masahiro Higashi
- Department of Chemistry and Supercomputing Institute, 207 Pleasant Street SE, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, 207 Pleasant Street SE, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
4
|
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]
|
5
|
Cybulski H, Sadlej J. Calculated Nuclear Magnetic Resonance Parameters for Multiproton-Exchange and Nonbonded-Hydrogen Rotation Processes in Cyclic Water Clusters. J Phys Chem A 2010; 115:5774-84. [DOI: 10.1021/jp107595n] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hubert Cybulski
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Joanna Sadlej
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
6
|
Tishchenko O, Truhlar DG. Gradient-based multiconfiguration Shepard interpolation for generating potential energy surfaces for polyatomic reactions. J Chem Phys 2010; 132:084109. [PMID: 20192292 DOI: 10.1063/1.3310296] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper describes and illustrates a way to construct multidimensional representations of reactive potential energy surfaces (PESs) by a multiconfiguration Shepard interpolation (MCSI) method based only on gradient information, that is, without using any Hessian information from electronic structure calculations. MCSI, which is called multiconfiguration molecular mechanics (MCMM) in previous articles, is a semiautomated method designed for constructing full-dimensional PESs for subsequent dynamics calculations (classical trajectories, full quantum dynamics, or variational transition state theory with multidimensional tunneling). The MCSI method is based on Shepard interpolation of Taylor series expansions of the coupling term of a 2 x 2 electronically diabatic Hamiltonian matrix with the diagonal elements representing nonreactive analytical PESs for reactants and products. In contrast to the previously developed method, these expansions are truncated in the present version at the first order, and, therefore, no input of electronic structure Hessians is required. The accuracy of the interpolated energies is evaluated for two test reactions, namely, the reaction OH+H(2)-->H(2)O+H and the hydrogen atom abstraction from a model of alpha-tocopherol by methyl radical. The latter reaction involves 38 atoms and a 108-dimensional PES. The mean unsigned errors averaged over a wide range of representative nuclear configurations (corresponding to an energy range of 19.5 kcal/mol in the former case and 32 kcal/mol in the latter) are found to be within 1 kcal/mol for both reactions, based on 13 gradients in one case and 11 in the other. The gradient-based MCMM method can be applied for efficient representations of multidimensional PESs in cases where analytical electronic structure Hessians are too expensive or unavailable, and it provides new opportunities to employ high-level electronic structure calculations for dynamics at an affordable cost.
Collapse
Affiliation(s)
- Oksana Tishchenko
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA.
| | | |
Collapse
|
7
|
Han JA, Kim YH. Multiconfiguration Molecular Mechanics Studies for the Potential Energy Surfaces of the Excited State Double Proton Transfer in the 1:1 7-Azaindole:H2OComplex. B KOREAN CHEM SOC 2010. [DOI: 10.5012/bkcs.2010.31.02.365] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
8
|
Duong MPT, Kim Y. Theoretical Studies for the Rates and Kinetic Isotope Effects of the Excited-State Double Proton Transfer in the 1:1 7-Azaindole:H2O Complex Using Variational Transition State Theory Including Multidimensional Tunneling. J Phys Chem A 2010; 114:3403-10. [DOI: 10.1021/jp910533m] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- My Phu Thi Duong
- Department of Chemistry, Kyung Hee University, 1 Seochun-dong, Giheung-Gu, Yongin-si, Gyeonggi-do, 449-701, Korea
| | - Yongho Kim
- Department of Chemistry, Kyung Hee University, 1 Seochun-dong, Giheung-Gu, Yongin-si, Gyeonggi-do, 449-701, Korea
| |
Collapse
|
9
|
Higashi M, Truhlar DG. Electrostatically Embedded Multiconfiguration Molecular Mechanics Based on the Combined Density Functional and Molecular Mechanical Method. J Chem Theory Comput 2008; 4:790-803. [DOI: 10.1021/ct800004y] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Masahiro Higashi
- Department of Chemistry and Supercomputing Institute, 207 Pleasant Street SE, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, 207 Pleasant Street SE, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
10
|
Tishchenko O, Truhlar DG. Global Potential Energy Surfaces with Correct Permutation Symmetry by Multiconfiguration Molecular Mechanics. J Chem Theory Comput 2007; 3:938-48. [DOI: 10.1021/ct600315h] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Oksana Tishchenko
- Chemistry Department and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Chemistry Department and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
11
|
Fernandez-Ramos A, Miller JA, Klippenstein SJ, Truhlar DG. Modeling the kinetics of bimolecular reactions. Chem Rev 2007; 106:4518-84. [PMID: 17091928 DOI: 10.1021/cr050205w] [Citation(s) in RCA: 474] [Impact Index Per Article: 27.9] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Antonio Fernandez-Ramos
- Departamento de Quimica Fisica, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| | | | | | | |
Collapse
|
12
|
Mrázek J, Burda JV. Can the pH value of water solutions be estimated by quantum chemical calculations of small water clusters? J Chem Phys 2006; 125:194518. [PMID: 17129134 DOI: 10.1063/1.2363383] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In the study, various water clusters were explored from the point of view of the proton transfer between H-bonded neighbors. A relatively modest approach--the MP26-31++G(d,p) level--was chosen as acceptable considering the fact that also larger systems have to be included. The tight ion-pair model (with usually three fixed O-O distances) was adopted for the autodissociation process. First, cluster-estimated pH values rapidly decrease as cluster size increases from 2 to 6. For larger clusters in gas phase, the topology of H bonds plays an important role, varying pH from 7 to 13 in hexamers and from 5 to 15 in octamer clusters. The relationship energy/distance was quantified, too. Enhancing our model with the conductorlike screening model (COSMO) approach brought significant improvement in description of the autodissociation reaction with a stable zwitterionic structure. However, when the O-O restrictions were released, the small barrier for backward reaction disappeared, reforming neutral cluster spontaneously. Also Boltzmann weighting procedure was applied, and for the explored clusters in vacuo, the series of pH 25-18-14-13-10 was obtained for cluster sizes n=2, 3, 4, 6, 8. Using the COSMO approach, the analogous series is 15-14-12-10-9. The limit of the series is still about two to three units above the experimentally known pH. In order to reach the size-independent (bulk) value, larger clusters are needed. However, the situation is far from hopeless since (as it was proven in the study) four-coordinated molecules are not involved in the proton transfer process directly; they can only be a part of the surrounding environment.
Collapse
Affiliation(s)
- Jirí Mrázek
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| | | |
Collapse
|
13
|
Antoniotti S, Golebiowski J, Cabrol-Bass D, Duñach E. Density functional theory investigations on acid-catalysed epoxide oxidative ring-opening by DMSO. Competition between oxidation processes. ACTA ACUST UNITED AC 2006. [DOI: 10.1016/j.theochem.2006.01.029] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|