1
|
Méndez E, Videla PE, Laria D. Equilibrium and Dynamical Characteristics of Hydrogen Bond Bifurcations in Water-Water and Water-Ammonia Dimers: A Path Integral Molecular Dynamics Study. J Phys Chem A 2022; 126:4721-4733. [PMID: 35834556 DOI: 10.1021/acs.jpca.2c02525] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present path integral molecular dynamics results that describe the effects of nuclear quantum fluctuations on equilibrium and dynamical characteristics pertaining to bifurcation pathways in hydrogen bonded dimers combining water and ammonia, at cryogenic temperatures of the order of 20 K. Along these isomerizations, the hydrogen atoms in the molecules acting as hydrogen-bond donors interchange their original dangling/connective characters. Our results reveal that the resulting quantum transition paths comprise three stages: the initial and final ones involve overall rotations during which the two protons retain their classical-like characteristics. Effects from quantum fluctuation are clearly manifested in the changes operated at the intermediate passages over transition states, as the spatial extents of the protons stretch over typical lengths comparable to the distances between connective and dangling basins of attractions. Consequently, the classical over-the-hill path is replaced by a tunneling controlled mechanism which, within the path integral perspective, can be cast in terms of concerted inter-basin migrations of polymer beads from dangling-to-connective and from connective-to-dangling, at practically no energy costs. We also estimated the characteristic timescales describing such interconversions within the approximate ring polymer rate theory. Effects derived from full and partial deuteration are also discussed.
Collapse
Affiliation(s)
- Emilio Méndez
- Departamento de Química Inorgánica, Analítica y Química-Física and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Pablo E Videla
- Department of Chemistry and Energy Sciences Institute, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Daniel Laria
- Departamento de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|
2
|
Rodríguez-Segundo R, Gijón A, Prosmiti R. Quantum molecular simulations of micro-hydrated halogen anions. Phys Chem Chem Phys 2022; 24:14964-14974. [PMID: 35686995 DOI: 10.1039/d2cp01396g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report the results of a detailed and accurate investigation focused on structures and energetics of poly-hydrated halides employing first-principles polarizable halide-water potentials to describe the underlying forces. Following a bottom-up data-driven potential approach, we initially looked into the classical behavior of higher-order X-(H2O)N clusters. We have located several low-lying energies, such as global and local minima, structures for each cluster, with various water molecules (up to N = 8) surrounding the halide anion (X- = F-, Cl-, Br-, I-), employing an evolutionary programming method. It is found that the F--water clusters exhibit different structural configurations than the heavier halides, however independently of the halide anion, all clusters show in general a selective growth with the anion preferring to be connected to the outer shell of the water molecule arrangements. In turn, path-integral molecular dynamics simulations are performed to incorporate explicitly nuclear quantum and thermal effects in describing the nature of halide ion microsolvation in such prototypical model systems. Our data reveal that at low finite temperatures, nuclear quantum effects affect certain structural properties, such as weakening hydrogen bonding between the halide anion and water molecules, with minor distortions in the water network beyond the first hydration shell, indicating local structure rearrangements. Such structural characteristics and the promising cluster size trends observed in the single-ion solvation energies motivated us to draw connections of small size cluster data to the limits of continuum bulk values, toward the investigation of the challenging computational modeling of bulk single ion hydration.
Collapse
Affiliation(s)
- Raúl Rodríguez-Segundo
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain. .,Atelgraphics S.L., Mota de Cuervo 42, 28043, Madrid, Spain
| | - Alfonso Gijón
- Materials Science Institute of Madrid (ICMM-CSIC), CSIC, Sor Juana Inés de la Cruz 3, 28049 Madrid, Spain
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| |
Collapse
|
3
|
Opoku RA, Toubin C, Gomes ASP. Simulating core electron binding energies of halogenated species adsorbed on ice surfaces and in solution via relativistic quantum embedding calculations. Phys Chem Chem Phys 2022; 24:14390-14407. [PMID: 35647703 DOI: 10.1039/d1cp05836c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work, we investigate the effects of the environment on the X-ray photoelectron spectra of hydrogen chloride and chloride ions adsorbed on ice surfaces, as well as of chloride ions in water droplets. In our approach, we combine a density functional theory (DFT) description of the ice surface with that of halogen species using the recently developed relativistic core-valence separation equation of motion coupled cluster (CVS-EOM-IP-CCSD) via the frozen density embedding formalism (FDE), to determine the K and L1,2,3 edges of chlorine. Our calculations, which incorporate temperature effects through snapshots from classical molecular dynamics simulations, are shown to reproduce experimental trends in the change of the core binding energies of Cl- upon moving from a liquid (water droplets) to an interfacial (ice quasi-liquid layer) environment. Our simulations yield water valence band binding energies in good agreement with experiment, which vary little between the droplets and the ice surface. For halide core binding energies there is an overall trend for overestimating experimental values, though good agreement between theory and experiment is found for Cl- in water droplets and on ice. For HCl on the other hand there are significant discrepancies between experimental and calculated core binding energies when we consider structural models that maintain the H-Cl bond more or less intact. An analysis of models that allow for pre-dissociated and dissociated structures suggests that experimentally observed chemical shifts in binding energies between Cl- and HCl would require that H+ (in the form of H3O+) and Cl- are separated by roughly 4-6 Å.
Collapse
Affiliation(s)
- Richard A Opoku
- Univ. Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, F-59000 Lille, France.
| | - Céline Toubin
- Univ. Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, F-59000 Lille, France.
| | | |
Collapse
|
4
|
Punyain W, Takahashi K. Evaluation of Ar tagging toward the vibrational spectra and zero point energy of X -HOH, X -DOH, and X -HOD, for X = F, Cl, Br. Phys Chem Chem Phys 2021; 23:9492-9499. [PMID: 33885081 DOI: 10.1039/d0cp06339h] [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
In this study, we theoretically evaluated the effect of argon tagging toward the binding energy and vibrational spectra of water halide anion complexes Ar.X-HOH, Ar.X-HOD, and Ar.X-DOH (X = F, Cl, Br). The ionic hydrogen bond (IHB) OH stretching mode was calculated to have a strong peak in the vibrational spectra, and coupling to intermolecular modes as well as bending modes was observed as combination bands and Fermi resonances. We found that the argon tagging affected the IHB OH stretching peak position in Ar.F-H2O, but not in Ar.Cl-H2O and Ar.Br-H2O. Furthermore, D-binding is favored for Cl and Br based on zero point energies, but for F our calculated zero point energies did not show a preference between H- and D-binding. We show that the competition of the energy lowering in the zero point energy of the anharmonic IHB OH (OD) stretching mode versus the intermolecular out-of-plane IHB OH (OD) wagging mode is important for determining the preference between H- and D-binding for these monohydrated halide clusters. We also found that for X-HOD the HOD bending fundamental peak is blue shifted compared to bare HOD, but is redshifted for F-DOH.
Collapse
Affiliation(s)
- Wikorn Punyain
- Institute of Atomic and Molecular Sciences, Academia Sinica, PO Box 23-166, Taipei 10617, Taiwan.
| | | |
Collapse
|
5
|
Méndez E, Laria D. Nuclear quantum effects on the hydrogen bond donor-acceptor exchange in water-water and water-methanol dimers. J Chem Phys 2020; 153:054302. [PMID: 32770908 DOI: 10.1063/5.0016122] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present results from path integral molecular dynamics simulations that describe effects from the explicit incorporation of nuclear quantum fluctuations on the topology of the free energy associated with the geared exchange of hydrogen bonds in the water-water dimer. Compared to the classical treatment, our results reveal important reductions in the free energy barriers and changes at a qualitative level in the overall profile. Most notable are those manifested by a plateau behavior, ascribed to nuclear tunneling, which bridges reactant and product states, contrasting with the usual symmetric double-well profile. The characteristics of the proton localizations along the pathway are examined. An imaginary time analysis of the rotational degrees of freedom of the partners in the dimer at the vicinities of transition states shows a clear "anticorrelation" between intermolecular interactions coupling beads localized in connective and dangling basins of attractions. As such, the transfer is operated by gradual concerted inter-basin migrations in opposite directions, at practically no energy costs. Modifications operated by partial deuteration and by the asymmetries in the hydrogen bonding characteristics prevailing in water-methanol heterodimers are also examined.
Collapse
Affiliation(s)
- Emilio Méndez
- Departamento de Química Inorgánica Analítica y Química-Física e INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Daniel Laria
- Departamento de Química Inorgánica Analítica y Química-Física e INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| |
Collapse
|
6
|
Sexton TM, Tschumper GS. 2-body:Many-body QM:QM study of structures, energetics, and vibrational frequencies for microhydrated halide ions. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1554827] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Thomas More Sexton
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS, USA
| | - Gregory S. Tschumper
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS, USA
| |
Collapse
|
7
|
Jung KA, Videla PE, Batista VS. Inclusion of nuclear quantum effects for simulations of nonlinear spectroscopy. J Chem Phys 2018; 148:244105. [PMID: 29960352 DOI: 10.1063/1.5036768] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The computation and interpretation of nonlinear vibrational spectroscopy is of vital importance for understanding a wide range of dynamical processes in molecular systems. Here, we introduce an approach to evaluate multi-time response functions in terms of multi-time double symmetrized Kubo transformed thermal correlation functions. Furthermore, we introduce a multi-time extension of ring polymer molecular dynamics to evaluate these Kubo transforms. Benchmark calculations show that the approximations are useful for short times even for nonlinear operators, providing a consistent improvement over classical simulations of multi-time correlation functions. The introduced methodology thus provides a practical way of including nuclear quantum effects in multi-time response functions of non-linear optical spectroscopy.
Collapse
Affiliation(s)
- Kenneth A Jung
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Pablo E Videla
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| |
Collapse
|
8
|
Mallory JD, Mandelshtam VA. Nuclear Quantum Effects and Thermodynamic Properties for Small (H2O)1–21X– Clusters (X– = F–, Cl–, Br–, I–). J Phys Chem A 2018; 122:4167-4180. [DOI: 10.1021/acs.jpca.8b00917] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Joel D. Mallory
- Department of Chemistry, University of California, Irvine, California 92697, United States
| | | |
Collapse
|
9
|
Videla PE, Rossky PJ, Laria D. Isotope effects in aqueous solvation of simple halides. J Chem Phys 2018; 148:102306. [DOI: 10.1063/1.4986231] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Affiliation(s)
- Pablo E. Videla
- Departamento de Química Inorgánica Analítica y Química-Física e INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Peter J. Rossky
- Department of Chemistry, Rice University, Houston, Texas 77251-1892, USA
| | - D. Laria
- Departamento de Química Inorgánica Analítica y Química-Física e INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
- Departamento de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|
10
|
|
11
|
Litman YE, Videla PE, Rodriguez J, Laria D. Positional Isotope Exchange in HX·(H2O)n (X = F, I) Clusters at Low Temperatures. J Phys Chem A 2016; 120:7213-24. [DOI: 10.1021/acs.jpca.6b06681] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yair E. Litman
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Pablo E. Videla
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Javier Rodriguez
- Departamento
de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
- ECyT, UNSAM, Martín
de Irigoyen 3100, 1650 San Martín, Pcia. de Buenos Aires, Argentina
| | - Daniel Laria
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
- Departamento
de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|