1
|
Tachikawa H. Mechanism of ionic dissociation of HCl in the smallest water clusters. Phys Chem Chem Phys 2024; 26:3623-3631. [PMID: 38224187 DOI: 10.1039/d3cp05715a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2024]
Abstract
The dissociation of strong acids into water is a fundamental process in chemistry and biology. Determining the minimum number of water molecules that can result in an ionic dissociation of hydrochloric acid (HCl → H+ + Cl-) remains a challenging subject. In this study, the reactions of H2O with HCl(H2O)n-1 (HCl-H2O cluster), i.e., HCl(H2O)n-1 + H2O (n = 3-7), were investigated by using the direct ab initio molecular dynamics (AIMD) method. Direct AIMD calculations were performed to set the collision energy of H2O to zero for all trajectories. For n = 3, no reaction occurred. In contrast, HCl dissociated to H+ + Cl- at n = 4, forming a contact ion pair (cIP) and solvent-separated ion pair (ssIP) as products. The reactions were expressed as HCl(H2O)3 + H2O → H3O+(H2O)2Cl- (ssIP), and HCl(H2O)3 + H2O → H3O+(Cl-)(H2O)2 (cIP). The ion pair (IP) products were dependent on the collision site of H2O relative to HCl(H2O)3. For n = 5-7, both IPs were formed through the reaction between H2O and HCl(H2O)n-1 (n = 5-7). The reaction between HCl and (H2O)4 (HCl + (H2O)4 → HCl(H2O)4) was non-reactive in IP formation. The reaction mechanism was discussed based on the theoretical results.
Collapse
Affiliation(s)
- Hiroto Tachikawa
- Division of Applied Chemistry, Faculty of Engineering, Hokkaido University, Sapporo 060-8628, Japan.
| |
Collapse
|
2
|
Bresnahan CG, David R, Milet A, Kumar R. Ion Pairing in HCl-Water Clusters: From Electronic Structure Investigations to Multiconfigurational Force-Field Development. J Phys Chem A 2019; 123:9371-9381. [PMID: 31589444 DOI: 10.1021/acs.jpca.9b07775] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In the bulk, condensed-phase HCl exists as a dissociated Cl- ion and a proton that is delocalized over solvating water molecules. However, in the gas phase, HCl is covalent, and even on the introduction of hydrating water molecules, the HCl covalent state dominates small clusters and is relevant at larger clusters including 21 water molecules. Electronic structure calculations (at the MP2 level) and ab initio metadynamics simulations (at the DFT level) have been carried out on HCl-(H2O)n clusters with n = 2-22 to investigate distinct solvation environments in clusters from covalent HCl structure, to contact ion pairs and solvent-separated ion pairs. The data were further used to train and validate a multiconfigurational force-field for HCl-water clusters that incorporates covalent HCl states into the MS-EVB3.2 formalism. Additionally, the many-body interaction of the Cl- ion with water and the excess proton was modeled by the introduction of two geometric three-body terms that incorporates the dominant many-body interaction in an efficient noniterative manner.
Collapse
Affiliation(s)
- Caitlin G Bresnahan
- Department of Chemistry , 232 Choppin Hall , Louisiana State University , Baton Rouge , Louisiana 70803 , United States
| | - Rolf David
- Department of Chemistry , 232 Choppin Hall , Louisiana State University , Baton Rouge , Louisiana 70803 , United States.,Univ. Grenoble Alpes, CNRS, DCM , 38000 Grenoble , France
| | - Anne Milet
- Univ. Grenoble Alpes, CNRS, DCM , 38000 Grenoble , France
| | - Revati Kumar
- Department of Chemistry , 232 Choppin Hall , Louisiana State University , Baton Rouge , Louisiana 70803 , United States
| |
Collapse
|
3
|
Christensen EG, Steele RP. Probing the Partial Activation of Water by Open-Shell Interactions, Cl(H 2O) 1-4. J Phys Chem A 2019; 123:8657-8673. [PMID: 31513400 DOI: 10.1021/acs.jpca.9b07235] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The partial chemical activation of water by reactive radicals was examined computationally for small clusters of chlorine and water, Cl•(H2O)n=1-4. Using an automated isomer-search procedure, dozens of unique, stable structures were computed. Among the resulting structural classes were intact, hydrated-chlorine isomers, as well as hydrogen-abstracted (HCl)(OH)(H2O)n-1 configurations. The latter showed increased stability as the degree of hydration increased, until n = 4, where a new class of structures was discovered with a chloride ion bound to an oxidized water network. The electronic structure of these three structural classes was investigated, and spectral signatures of this hydration-based evolution were connected to these electronic properties. An ancillary outcome of this detailed computational analysis, including coupled-cluster benchmarks, was the calibration of cost-effective quantum chemistry methods for future studies of these radical-water complexes.
Collapse
Affiliation(s)
- Elizabeth G Christensen
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry , University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112 , United States
| | - Ryan P Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry , University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112 , United States
| |
Collapse
|
4
|
Horká-Zelenková V, Seyfang G, Dietiker P, Quack M. Nuclear Spin Symmetry Conservation Studied for Symmetric Top Molecules (CH3D, CHD3, CH3F, and CH3Cl) in Supersonic Jet Expansions. J Phys Chem A 2019; 123:6160-6174. [DOI: 10.1021/acs.jpca.9b02580] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Veronika Horká-Zelenková
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
- J. Heyrovský Institute of Physical Chemistry of the CAS, Dolejskova 2155/3, 182 23 Prague 8, Czech Republic
| | - Georg Seyfang
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Peter Dietiker
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Martin Quack
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
| |
Collapse
|
5
|
Kwasniewski D, Butler M, Reisler H. Vibrational predissociation of the phenol-water dimer: a view from the water. Phys Chem Chem Phys 2019; 21:13968-13976. [PMID: 30511053 DOI: 10.1039/c8cp06581k] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The vibrational predissociation (VP) dynamics of the phenol-water (PhOH-H2O) dimer were studied by detecting H2O fragments and using velocity map imaging (VMI) to infer the internal energy distributions of PhOH cofragments, pair-correlated with selected rotational levels of the H2O fragments. Following infrared (IR) laser excitation of the hydrogen-bonded OH stretch fundamental of PhOH (Pathway 1) or the asymmetric OH stretch localized on H2O (Pathway 2), dissociation to H2O + PhOH was observed. H2O fragments were monitored state-selectively by using 2+1 Resonance-Enhanced Multiphoton Ionization (REMPI) combined with time-of-flight mass spectrometry (TOF-MS). VMI of H2O in selected rotational levels was used to derive center-of-mass (c.m.) translational energy (ET) distributions. The pair-correlated internal energy distributions of the PhOH cofragments derived via Pathway 1 were well described by a statistical prior distribution. On the other hand, the corresponding distributions obtained via Pathway 2 show a propensity to populate higher-energy rovibrational levels of PhOH than expected from a statistical distribution and agree better with an energy-gap model. The REMPI spectra of the H2O fragments from both pathways could be fit by Boltzmann plots truncated at the maximum allowed energy, with a higher temperature for Pathway 2 than that for Pathway 1. We conclude that the VP dynamics depends on the OH stretch level initially excited.
Collapse
Affiliation(s)
- Daniel Kwasniewski
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA.
| | - Mitchell Butler
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA.
| | - Hanna Reisler
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA.
| |
Collapse
|
6
|
Zuraski K, Wang QK, Kwasniewski D, Bowman JM, Reisler H. Predissociation dynamics of the HCl-(H 2O) 3 tetramer: An experimental and theoretical investigation. J Chem Phys 2018; 148:204303. [PMID: 29865837 DOI: 10.1063/1.5026585] [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
The cyclic HCl-(H2O)3 tetramer is the largest observed neutral HCl-(H2O)n cluster. The vibrational predissociation of HCl-(H2O)3 is investigated by theory, quasiclassical trajectory (QCT) calculations, and experiment, following the infrared excitation of the hydrogen-bonded OH-stretch fundamental. The energetically possible dissociation pathways are HCl + (H2O)3 (Pathway 1) and H2O + HCl-(H2O)2 (Pathway 2). The HCl and H2O monomer fragments are observed by 2 + 1 resonance enhanced multiphoton ionization combined with time-of-flight mass spectrometry, and their rotational energy distributions are inferred and compared to the theoretical results. Velocity map images of the monomer fragments in selected rotational levels are used for each pathway to obtain pair-correlated speed distributions. The fragment speed distributions obtained by experiment and QCT calculations are broad and structureless, encompassing the entire range of allowed speeds for each pathway. Bond dissociation energies, D0, are estimated experimentally from the endpoints of the speed distributions: 2100 ± 300 cm-1 and 2400 ± 100 cm-1 for Pathway 1 and Pathway 2, respectively. These values are lower but in the same order as the corresponding calculated D0: 2426 ± 23 cm-1 and 2826 ± 19 cm-1. The differences are attributed to contributions from vibrational hot bands of the clusters that appear in the high-speed tails of the experimental pair-correlated distributions. Satisfactory agreement between theory and experiment is achieved when comparing the monomer fragments' rotational energies, the shapes of the speed distributions, and the average fragment speeds and center-of-mass translational energies. Insights into the dissociation mechanism and lifetime are gained from QCT calculations performed on a previously reported many-body potential energy surface. It is concluded that the dissociation lifetime is on the order of 10 ps and that the final trimer products are in their lowest energy cyclic forms.
Collapse
Affiliation(s)
- Kristen Zuraski
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| | - Qingfeng Kee Wang
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Daniel Kwasniewski
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Hanna Reisler
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| |
Collapse
|
7
|
Pérez de Tudela R, Marx D. Acid Dissociation in HCl-Water Clusters is Temperature Dependent and Cannot be Detected Based on Dipole Moments. PHYSICAL REVIEW LETTERS 2017; 119:223001. [PMID: 29286767 DOI: 10.1103/physrevlett.119.223001] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/24/2017] [Indexed: 05/09/2023]
Abstract
The dissociation of acids in aqueous environments at low temperatures in the presence of a limited amount of water is underlying a wealth of processes from atmospheric to interstellar science. For the paradigmatic case of HCl(H_{2}O)_{n} clusters, our extensive ab initio path integral simulations quantify in terms of free energy differences and barriers that n=4 water molecules are indeed required to dissociate HCl at low temperatures. Increasing the temperature, however, reverses the process and thus counteracts dissociation by fluctuation-driven recombination. The size of the electric dipole moment is shown to not correlate with the acid being in its dissociated or molecular state, thus rendering its measurement as a function of n unable to detect the dissociation transition.
Collapse
Affiliation(s)
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
8
|
Tachikawa H. Proton Transfer Rates in Ionized Hydrogen Chloride–Water Clusters: A Direct Ab Initio Molecular Dynamics Study. J Phys Chem A 2017; 121:5237-5244. [DOI: 10.1021/acs.jpca.7b05112] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hiroto Tachikawa
- Division of Applied Chemistry, Graduate
School of Engineering, Hokkaido University, Sapporo 060-8628, Japan
| |
Collapse
|