1
|
Wang Z, Ju S, Wang Y, Zhang R, Ma L, Song J, Lin K. The isosbestic point in the Raman spectra of the hydration shell. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 317:124413. [PMID: 38728849 DOI: 10.1016/j.saa.2024.124413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Revised: 04/16/2024] [Accepted: 05/05/2024] [Indexed: 05/12/2024]
Abstract
Isosbestic point is often observed in a series of spectra, but their interpretation is still controversial, such as whether the continuum model can produce an isosbestic point. In order to answer this question, the Raman spectra of hydration shell with continuous distribution structure in different ionic aqueous solutions were separated by Raman ratio spectra, and an isosbestic point was successfully observed. Our experimental results show that the continuum model can indeed produce the isosbestic point. In order to deepen the understanding of the isosbestic point, we calculate the first moment of the Raman spectra and conduct molecular dynamics (MD) simulations. Both experimental and theoretical findings indicate that elevated temperatures lead to increased disorder among water molecules within the hydration shell.
Collapse
Affiliation(s)
- Zhiqiang Wang
- School of Physics, Xidian University, Xi'an, 710071, P. R. China
| | - Siwen Ju
- School of Physics, Xidian University, Xi'an, 710071, P. R. China
| | - Yuxi Wang
- School of Flexible Electronics (SOFE) & State Key Laboratory of Optoelectronic Materials and Technologies (OEMT), Sun Yat-sen University, Shenzhen, 5181071, P. R. China
| | - Ruiting Zhang
- School of Physics, Xidian University, Xi'an, 710071, P. R. China.
| | - Lin Ma
- School of Physics, Xidian University, Xi'an, 710071, P. R. China
| | - Jiangluqi Song
- School of Physics, Xidian University, Xi'an, 710071, P. R. China
| | - Ke Lin
- School of Physics, Xidian University, Xi'an, 710071, P. R. China; Interdisciplinary Research Center of Smart Sensor, Xidian University, Xi'an, 710071, P. R. China.
| |
Collapse
|
2
|
Hunter KM, Paesani F. Monitoring water harvesting in metal-organic frameworks, one water molecule at a time. Chem Sci 2024; 15:5303-5310. [PMID: 38577368 PMCID: PMC10988614 DOI: 10.1039/d3sc06162k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Accepted: 03/05/2024] [Indexed: 04/06/2024] Open
Abstract
Metal-organic frameworks (MOFs) have gained prominence as potential materials for atmospheric water harvesting, a vital solution for arid regions and areas experiencing severe water shortages. However, the molecular factors influencing the performance of MOFs in capturing water from the air remain elusive. Among all MOFs, Ni2X2BTDD (X = F, Cl, Br) stands out as a promising water harvester due to its ability to adsorb substantial amounts of water at low relative humidity (RH). Here, we use advanced molecular dynamics simulations carried out with the state-of-the-art MB-pol data-driven many-body potential to monitor water adsorption in the three Ni2X2BTDD variants as a function of RH. Our simulations reveal that the type of halide atom in the three Ni2X2BTDD frameworks significantly influences the corresponding molecular mechanisms of water adsorption: while water molecules form strong hydrogen bonds with the fluoride atoms in Ni2F2BTDD, they tend to form hydrogen bonds with the nitrogen atoms of the triazolate linkers in Ni2Cl2BTDD and Ni2Br2BTDD. Importantly, the large size of the bromide atoms reduces the void volume in the Ni2Br2BTDD pores, which enable water molecules to initiate an extended hydrogen-bond network at lower RH. These findings not only underscore the prospect for precisely tuning structural and chemical modifications of the frameworks to optimize their interaction with water, but also highlight the predictive power of simulations with the MB-pol data-driven many-body potential. By providing a realistic description of water under different thermodynamic conditions and environments, these simulations yield unique, molecular-level insights that can guide the design and optimization of energy-efficient water harvesting materials.
Collapse
Affiliation(s)
- Kelly M Hunter
- Department of Chemistry and Biochemistry, University of California La Jolla San Diego California 92093 USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California La Jolla San Diego California 92093 USA
- Materials Science and Engineering, University of California La Jolla San Diego California 92093 USA
- Halicioğlu Data Science Institute, University of California La Jolla San Diego California 92093 USA
- San Diego Supercomputer Center, University of California La Jolla San Diego California 92093 USA
| |
Collapse
|
3
|
Savoj R, Agnew H, Zhou R, Paesani F. Molecular Insights into the Influence of Ions on the Water Structure. I. Alkali Metal Ions in Solution. J Phys Chem B 2024; 128:1953-1962. [PMID: 38373140 DOI: 10.1021/acs.jpcb.3c08150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/21/2024]
Abstract
In this study, we explore the impact of alkali metal ions (Li+, Na+, K+, Rb+, and Cs+) on the hydration structure of water using molecular dynamics simulations carried out with MB-nrg potential energy functions (PEFs). Our analyses include radial distribution functions, coordination numbers, dipole moments, and infrared spectra of water molecules, calculated as a function of solvation shells. The results collectively indicate a highly local influence of all of the alkali metal ions on the hydrogen-bond network established by the surrounding water molecules, with the smallest and most densely charged Li+ ion exerting the most pronounced effect. Remarkably, the MB-nrg PEFs demonstrate excellent agreement with available experimental data for the position and size of the first solvation shells, underscoring their potential as predictive models for realistic simulations of ionic aqueous solutions across various thermodynamic conditions and environments.
Collapse
Affiliation(s)
- Roya Savoj
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Henry Agnew
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Ruihan Zhou
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- Halicioğlu Data Science Institute, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
4
|
LaCour RA, Heindel JP, Head-Gordon T. Predicting the Raman Spectra of Liquid Water with a Monomer-Field Model. J Phys Chem Lett 2023; 14:11742-11749. [PMID: 38116782 DOI: 10.1021/acs.jpclett.3c02873] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2023]
Abstract
The Raman spectrum of liquid water is quite complex, reflecting its strong sensitivity to the local environment of the individual waters. The OH-stretch region of the spectrum, which captures the influence of hydrogen bonding, has only just begun to be unraveled. Here we develop a model for predicting the Raman spectra of the OH-stretch region by considering how local electric fields distort the energy surface of each water monomer. We find that our model is capable of reproducing the bimodal nature of the main peak, with the shoulder at 3250 cm-1 resulting almost entirely from Fermi resonance. Furthermore, we capture the temperature and polarization dependence of the shoulder, which has proven to be difficult to obtain with previous methods, and analyze the origin of this dependence. We expect our model to be generally useful for understanding and predicting how Raman spectra change under different conditions and with different probe reporters beyond water.
Collapse
Affiliation(s)
- R Allen LaCour
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Joseph P Heindel
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Departments of Bioengineering and Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, United States
| |
Collapse
|
5
|
Lwin E, Fischer TL, Suhm MA. Microhydration of Tertiary Amines: Robust Resonances in Red-Shifted Water. J Phys Chem Lett 2023; 14:10194-10199. [PMID: 37930195 PMCID: PMC10658632 DOI: 10.1021/acs.jpclett.3c02517] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 10/29/2023] [Accepted: 11/01/2023] [Indexed: 11/07/2023]
Abstract
Tertiary amines are strong hydrogen bond acceptors. When a water molecule donates one of the OH groups, its in-phase stretching vibration wavenumber is decreased to such an extent that it comes close to the water bending overtone. This gives rise to anharmonic phenomena such as classical Fermi resonance, resonance with multiple-quantum dark states, or combination transitions with low-frequency intermolecular modes. These effects, which contribute to the spectral breadth of room-temperature hydrogen-bonded amine complexes, are disentangled by Fourier transform infrared spectroscopy in pulsed supersonic slit jet expansions. Monohydrates of the amines quinuclidine, N-methylpyrrolidine, N-methylpiperidine, and dimethylcyclohexylamine exhibit systematic mode coupling signatures. These suggest relatively fast energy flow out of the excited OH stretching fundamental into intra- and intermolecular degrees of freedom of the hydrogen-bonded water molecule. Trimeric complexes are spectroscopically separated from the amine monohydrates.
Collapse
Affiliation(s)
- Eaindra Lwin
- Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Taija L. Fischer
- Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Martin A. Suhm
- Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| |
Collapse
|
6
|
Mendonça BHS, de Moraes EE, Kirch A, Batista RJC, de Oliveira AB, Barbosa MC, Chacham H. Flow through Deformed Carbon Nanotubes Predicted by Rigid and Flexible Water Models. J Phys Chem B 2023; 127:8634-8643. [PMID: 37754781 DOI: 10.1021/acs.jpcb.3c02889] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/28/2023]
Abstract
In this study, using nonequilibrium molecular dynamics simulation, the flow of water in deformed carbon nanotubes is studied for two water models TIP4P/2005 and simple point charge/FH (SPC/FH). The results demonstrated a nonuniform dependence of the flow on the tube deformation and the flexibility imposed on the water molecules, leading to an unexpected increase in the flow in some cases. The effects of the tube diameter and pressure gradient are investigated to explain the abnormal flow behavior with different degrees of structural deformation.
Collapse
Affiliation(s)
- Bruno H S Mendonça
- Departamento de Física, ICEX, Universidade Federal de Minas Gerais, CP 702, Belo Horizonte 30123-970, MG, Brazil
| | - Elizane E de Moraes
- Instituto de Física, Universidade Federal da Bahia, Campus Universitário de Ondina, Salvador 40210-340, BA, Brazil
| | - Alexsandro Kirch
- Instituto de Física, Universidade de São Paulo, CP 66318, São Paulo 05315-970, SP, Brazil
| | - Ronaldo J C Batista
- Departamento de Física, Universidade Federal de Ouro Preto, Campus Morro do Cruzeiro, Ouro Preto 35400-000, MG, Brazil
| | - Alan B de Oliveira
- Departamento de Física, Universidade Federal de Ouro Preto, Campus Morro do Cruzeiro, Ouro Preto 35400-000, MG, Brazil
| | - Marcia C Barbosa
- Instituto de Física, Universidade Federal do Rio Grande do Sul, Porto Alegre 91501-970, RS, Brazil
| | - Hélio Chacham
- Departamento de Física, ICEX, Universidade Federal de Minas Gerais, CP 702, Belo Horizonte 30123-970, MG, Brazil
| |
Collapse
|
7
|
Riera M, Knight C, Bull-Vulpe EF, Zhu X, Agnew H, Smith DGA, Simmonett AC, Paesani F. MBX: A many-body energy and force calculator for data-driven many-body simulations. J Chem Phys 2023; 159:054802. [PMID: 37526156 PMCID: PMC10550339 DOI: 10.1063/5.0156036] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Accepted: 07/11/2023] [Indexed: 08/02/2023] Open
Abstract
Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the "many-body energy" (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.
Collapse
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Christopher Knight
- Argonne National Laboratory, Computational Science Division, Lemont, Illinois 60439, USA
| | - Ethan F. Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Xuanyu Zhu
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Henry Agnew
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | | | - Andrew C. Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | | |
Collapse
|
8
|
Pastorczak M, Duk K, Shahab S, Kananenka AA. Combinational Vibration Modes in H 2O/HDO/D 2O Mixtures Detected Thanks to the Superior Sensitivity of Femtosecond Stimulated Raman Scattering. J Phys Chem B 2023. [PMID: 37201478 DOI: 10.1021/acs.jpcb.3c01334] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Overtones and combinational modes frequently play essential roles in ultrafast vibrational energy relaxation in liquid water. However, these modes are very weak and often overlap with fundamental modes, particularly in isotopologues mixtures. We measured VV and HV Raman spectra of H2O and D2O mixtures with femtosecond stimulated Raman scattering (FSRS) and compared the results with calculated spectra. Specifically, we observed the mode at around 1850 cm-1 and assigned it to H-O-D bend + rocking libration. Second, we found that the H-O-D bend overtone band and the OD stretch + rocking libration combination band contribute to the band located between 2850 and 3050 cm-1. Furthermore, we assigned the broad band located between 4000 and 4200 cm-1 to be composed of combinational modes of high-frequency OH stretching modes with predominantly twisting and rocking librations. These results should help in a proper interpretation of Raman spectra of aqueous systems as well as in the identification of vibrational relaxation pathways in isotopically diluted water.
Collapse
Affiliation(s)
- Marcin Pastorczak
- Institute of Physical Chemistry, Polish Academy of Sciences, Laser Centre, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Katsiaryna Duk
- Institute of Physical Chemistry, Polish Academy of Sciences, Laser Centre, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Samaneh Shahab
- Institute of Physical Chemistry, Polish Academy of Sciences, Laser Centre, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Alexei A Kananenka
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
9
|
Inoue K, Litman Y, Wilkins DM, Nagata Y, Okuno M. Is Unified Understanding of Vibrational Coupling of Water Possible? Hyper-Raman Measurement and Machine Learning Spectra. J Phys Chem Lett 2023; 14:3063-3068. [PMID: 36947156 DOI: 10.1021/acs.jpclett.3c00398] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
The impact of the vibrational coupling of the OH stretch mode on the spectra differs significantly between IR and Raman spectra of water. Unified understanding of the vibrational couplings is not yet achieved. By using a different class of vibrational spectroscopy, hyper-Raman (HR) spectroscopy, together with machine-learning-assisted HR spectra calculation, we examine the impact of the vibrational couplings of water through the comparison of isotopically diluted H2O and pure H2O. We found that the isotopic dilution reduces the HR bandwidths, but the impact of the vibrational coupling is smaller than in the IR and parallel-polarized Raman. Machine learning HR spectra indicate that the intermolecular coupling plays a major role in broadening the bandwidth, while the intramolecular coupling is negligibly small, which is consistent with the IR and Raman spectra. Our result clearly demonstrates a limited impact of the intramolecular vibration, independent of the selection rules of vibrational spectroscopies.
Collapse
Affiliation(s)
- Kazuki Inoue
- Department of Basic Science, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan
| | - Yair Litman
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - David M Wilkins
- Atomistic Simulation Centre, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, Northern Ireland, United Kingdom
| | - Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Masanari Okuno
- Department of Basic Science, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan
| |
Collapse
|
10
|
Mayerhöfer TG, Pahlow S, Ivanovski V, Popp J. Dispersion related coupling effects in IR spectra on the example of water and Amide I bands. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2023; 288:122115. [PMID: 36436263 DOI: 10.1016/j.saa.2022.122115] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Revised: 11/02/2022] [Accepted: 11/08/2022] [Indexed: 06/16/2023]
Abstract
We discuss coupling effects in infrared spectra which are caused by dispersion and local field effects. The first effect is instigated by changes of the refractive index due to absorption which have an impact on the strength of adjacent absorptions. The second effect is a consequence of the light-induced polarization of one molecule affecting neighboring ones. These coupling effects do not only effect band positions, but also influence relative intensities. They are particularly strong in case of overlapping bands and complicate their deconvolution by band fitting. We investigated the corresponding challenges for the HO-stretching vibrations in water and the Amide I band in proteins. Our findings show that the effects are significant and of high interest for protein and water structure determination. Especially, for the water stretching vibrations we conclude that it is of utmost importance to consider such coupling effects in quantum mechanical calculations of water spectra. Otherwise, progress in understanding band positions and profiles is likely to be hampered. Also, in case of the Amide I band we found a distinct impact of such coupling effects. Accordingly, we strongly recommend consideration of dispersion and local field effects to ensure the possibility of an accurate, quantitative determination of α-helix and β-sheet structures.
Collapse
Affiliation(s)
- Thomas G Mayerhöfer
- Leibniz Institute of Photonic Technology (IPHT), Jena, 07745, Albert-Einstein-Str. 9, Germany; Institute of Physical Chemistry and Abbe Center of Photonics, Friedrich Schiller University, Jena, 07743, Helmholtzweg 4, Germany.
| | - Susanne Pahlow
- Leibniz Institute of Photonic Technology (IPHT), Jena, 07745, Albert-Einstein-Str. 9, Germany; Institute of Physical Chemistry and Abbe Center of Photonics, Friedrich Schiller University, Jena, 07743, Helmholtzweg 4, Germany
| | - Vladimir Ivanovski
- Institute of Chemistry, Faculty of Natural Sciences and Mathematics, Ss. Cyril and Methodius University in Skopje, Arhimedova 5, 1000 Skopje, The Former Yugolav Republic of Macedonia
| | - Jürgen Popp
- Leibniz Institute of Photonic Technology (IPHT), Jena, 07745, Albert-Einstein-Str. 9, Germany; Institute of Physical Chemistry and Abbe Center of Photonics, Friedrich Schiller University, Jena, 07743, Helmholtzweg 4, Germany
| |
Collapse
|
11
|
Takahashi H, Tanimura Y. Discretized hierarchical equations of motion in mixed Liouville-Wigner space for two-dimensional vibrational spectroscopies of liquid water. J Chem Phys 2023; 158:044115. [PMID: 36725520 DOI: 10.1063/5.0135725] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
A model of a bulk water system describing the vibrational motion of intramolecular and intermolecular modes is constructed, enabling analysis of its linear and nonlinear vibrational spectra as well as the energy transfer processes between the vibrational modes. The model is described as a system of four interacting anharmonic oscillators nonlinearly coupled to their respective heat baths. To perform a rigorous numerical investigation of the non-Markovian and nonperturbative quantum dissipative dynamics of the model, we derive discretized hierarchical equations of motion in mixed Liouville-Wigner space, with Lagrange-Hermite mesh discretization being employed in the Liouville space of the intramolecular modes and Lagrange-Hermite mesh discretization and Hermite discretization in the Wigner space of the intermolecular modes. One-dimensional infrared and Raman spectra and two-dimensional terahertz-infrared-visible and infrared-infrared-Raman spectra are computed as demonstrations of the quantum dissipative description provided by our model.
Collapse
Affiliation(s)
- Hideaki Takahashi
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
12
|
Yu X, Chiang KY, Yu CC, Bonn M, Nagata Y. On the Fresnel factor correction of sum-frequency generation spectra of interfacial water. J Chem Phys 2023; 158:044701. [PMID: 36725499 DOI: 10.1063/5.0133428] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Insights into the microscopic structure of aqueous interfaces are essential for understanding the chemical and physical processes on the water surface, including chemical synthesis, atmospheric chemistry, and events in biomolecular systems. These aqueous interfaces have been probed by heterodyne-detected sum-frequency generation (HD-SFG) spectroscopy. To obtain the molecular response from the measured HD-SFG spectra, one needs to correct the measured ssp spectra for local electromagnetic field effects at the interface due to a spatially varying dielectric function. This so-called Fresnel factor correction can change the inferred response substantially, and different ways of performing this correction lead to different conclusions about the interfacial water response. Here, we compare the simulated and experimental spectra at the air/water interface. We use three previously developed models to compare the experiment with theory: an advanced approach taking into account the detailed inhomogeneous interfacial dielectric profile and the Lorentz and slab models to approximate the interfacial dielectric function. Using the advanced model, we obtain an excellent quantitative agreement between theory and experiment, in both spectral shape and amplitude. Remarkably, we find that for the Fresnel factor correction of the ssp spectra, the Lorentz model for the interfacial dielectric function is equally accurate in the hydrogen (H)-bonded region of the response, while the slab model underestimates this response significantly. The Lorentz model, thus, provides a straightforward method to obtain the molecular response from the measured spectra of aqueous interfaces in the H-bonded region.
Collapse
Affiliation(s)
- Xiaoqing Yu
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Kuo-Yang Chiang
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Chun-Chieh Yu
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Mischa Bonn
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
13
|
Zhuang D, Riera M, Zhou R, Deary A, Paesani F. Hydration Structure of Na + and K + Ions in Solution Predicted by Data-Driven Many-Body Potentials. J Phys Chem B 2022; 126:9349-9360. [PMID: 36326071 DOI: 10.1021/acs.jpcb.2c05674] [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/05/2022]
Abstract
The hydration structure of Na+ and K+ ions in solution is systematically investigated using a hierarchy of molecular models that progressively include more accurate representations of many-body interactions. We found that a conventional empirical pairwise additive force field that is commonly used in biomolecular simulations is unable to reproduce the extended X-ray absorption fine structure (EXAFS) spectra for both ions. In contrast, progressive inclusion of many-body effects rigorously derived from the many-body expansion of the energy allows the MB-nrg potential energy functions (PEFs) to achieve nearly quantitative agreement with the experimental EXAFS spectra, thus enabling the development of a molecular-level picture of the hydration structure of both Na+ and K+ in solution. Since the MB-nrg PEFs have already been shown to accurately describe isomeric equilibria and vibrational spectra of small ion-water clusters in the gas phase, the present study demonstrates that the MB-nrg PEFs effectively represent the long-sought-after models able to correctly predict the properties of ionic aqueous systems from the gas to the liquid phase, which has so far remained elusive.
Collapse
Affiliation(s)
- Debbie Zhuang
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Ruihan Zhou
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Alexander Deary
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California92093, United States.,Materials Science and Engineering, University of California San Diego, La Jolla, California92093, United States.,San Diego Supercomputer Center, University of California San Diego, La Jolla, California92093, United States
| |
Collapse
|
14
|
Stefaniuk A, Gawinkowski S, Golec B, Gorski A, Szutkowski K, Waluk J, Poznański J. Isotope effects observed in diluted D 2O/H 2O mixtures identify HOD-induced low-density structures in D 2O but not H 2O. Sci Rep 2022; 12:18732. [PMID: 36333587 PMCID: PMC9636167 DOI: 10.1038/s41598-022-23551-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2022] [Accepted: 11/02/2022] [Indexed: 11/06/2022] Open
Abstract
Normal and heavy water are solvents most commonly used to study the isotope effect. The isotope effect of a solvent significantly influences the behavior of a single molecule in a solution, especially when there are interactions between the solvent and the solute. The influence of the isotope effect becomes more significant in D2O/H2O since the hydrogen bond in H2O is slightly weaker than its counterpart (deuterium bond) in D2O. Herein, we characterize the isotope effect in a mixture of normal and heavy water on the solvation of a HOD molecule. We show that the HOD molecule affects the proximal solvent molecules, and these disturbances are much more significant in heavy water than in normal water. Moreover, in D2O, we observe the formation of low-density structures indicative of an ordering of the solvent around the HOD molecule. The qualitative differences between HOD interaction with D2O and H2O were consistently confirmed with Raman spectroscopy and NMR diffusometry.
Collapse
Affiliation(s)
- Anna Stefaniuk
- grid.418825.20000 0001 2216 0871Institute of Biochemistry and Biophysics Polish Academy of Sciences, Pawińskiego 5a, 02-106 Warsaw, Poland
| | - Sylwester Gawinkowski
- grid.425290.80000 0004 0369 6111Institute of Physical Chemistry Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Barbara Golec
- grid.425290.80000 0004 0369 6111Institute of Physical Chemistry Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Aleksander Gorski
- grid.425290.80000 0004 0369 6111Institute of Physical Chemistry Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Kosma Szutkowski
- grid.5633.30000 0001 2097 3545Adam Mickiewicz University, NanoBioMedical Centre, Wszechnicy Piastowskiej 3, 61-614 Poznan, Poland
| | - Jacek Waluk
- grid.425290.80000 0004 0369 6111Institute of Physical Chemistry Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland ,grid.440603.50000 0001 2301 5211Faculty of Mathematics and Science, Cardinal Stefan Wyszyński University, Dewajtis 5, 01-815 Warsaw, Poland
| | - Jarosław Poznański
- grid.418825.20000 0001 2216 0871Institute of Biochemistry and Biophysics Polish Academy of Sciences, Pawińskiego 5a, 02-106 Warsaw, Poland
| |
Collapse
|
15
|
Bull-Vulpe EF, Riera M, Bore SL, Paesani F. Data-Driven Many-Body Potential Energy Functions for Generic Molecules: Linear Alkanes as a Proof-of-Concept Application. J Chem Theory Comput 2022. [PMID: 36113028 DOI: 10.1021/acs.jctc.2c00645] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a generalization of the many-body energy (MB-nrg) theoretical/computational framework that enables the development of data-driven potential energy functions (PEFs) for generic covalently bonded molecules, with arbitrary quantum mechanical accuracy. The "nearsightedness of electronic matter" is exploited to define monomers as "natural building blocks" on the basis of their distinct chemical identity. The energy of generic molecules is then expressed as a sum of individual many-body energies of incrementally larger subsystems. The MB-nrg PEFs represent the low-order n-body energies, with n = 1-4, using permutationally invariant polynomials derived from electronic structure data carried out at an arbitrary quantum mechanical level of theory, while all higher-order n-body terms (n > 4) are represented by a classical many-body polarization term. As a proof-of-concept application of the general MB-nrg framework, we present MB-nrg PEFs for linear alkanes. The MB-nrg PEFs are shown to accurately reproduce reference energies, harmonic frequencies, and potential energy scans of alkanes, independently of their length. Since, by construction, the MB-nrg framework introduced here can be applied to generic covalently bonded molecules, we envision future computer simulations of complex molecular systems using data-driven MB-nrg PEFs, with arbitrary quantum mechanical accuracy.
Collapse
Affiliation(s)
- Ethan F. Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Sigbjørn L. Bore
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
16
|
Ryu WH, Voth GA. Coarse-Graining of Imaginary Time Feynman Path Integrals: Inclusion of Intramolecular Interactions and Bottom-up Force-Matching. J Phys Chem A 2022; 126:6004-6019. [PMID: 36007243 PMCID: PMC9466601 DOI: 10.1021/acs.jpca.2c04349] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Revised: 08/16/2022] [Indexed: 11/29/2022]
Abstract
Feynman's imaginary time path integral formalism of quantum statistical mechanics and the corresponding quantum-classical isomorphism provide a tangible way of incorporating nuclear quantum effect (NQE) in the simulation of condensed matter systems using well-developed classical simulation techniques. Our previous work has presented the many-body coarse-graining of path integral (CG-PI) theory that builds an isomorphism between the quantum partition function of N distinguishable particles and the classical partition function of 2N pseudoparticles. In this present work, we develop a generalized version of the many-body CG-PI theory that incorporates many-body interactions in the force field. Based on the new derivation, we provide a numerical CG-PI (n-CG-PI) modeling strategy parametrized from the underlying path integral molecular dynamics (PIMD) trajectories using force matching and Boltzmann inversion. The n-CG-PI models for two liquid systems are shown to capture well both the intramolecular and intermolecular structural correlations of the reference PIMD simulations. The generalized derivation of the many-body CG-PI theory and the n-CG-PI model presented in this work extend the scope of the CG-PI formalism by generalizing the previously limited theory to incorporate force fields of realistic molecular systems.
Collapse
Affiliation(s)
- Won Hee Ryu
- Department of Chemistry,
James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department of Chemistry,
James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
17
|
Liu J, Lan J, He X. Toward High-level Machine Learning Potential for Water Based on Quantum Fragmentation and Neural Networks. J Phys Chem A 2022; 126:3926-3936. [PMID: 35679610 DOI: 10.1021/acs.jpca.2c00601] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Accurate and efficient simulation of liquids, such as water and salt solutions, using high-level wave function theories is still a formidable task for computational chemists owing to the high computational costs. In this study, we develop a deep machine learning potential based on fragment-based second-order Møller-Plesset perturbation theory (DP-MP2) for water through neural networks. We show that the DP-MP2 potential predicts the structural, dynamical, and thermodynamic properties of liquid water in better agreement with the experimental data than previous studies based on density functional theory (DFT). The nuclear quantum effects (NQEs) on the properties of liquid water are also examined, which are noticeable in affecting the structural and dynamical properties of liquid water under ambient conditions. This work provides a general framework for quantitative predictions of the properties of condensed-phase systems with the accuracy of high-level wave function theory while achieving significant computational savings compared to ab initio simulations.
Collapse
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China.,Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Jinggang Lan
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,New York University-East China Normal University Center for Computational Chemistry, NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
18
|
Valtera S, Jašík J, Vaidulych M, Olszówka JE, Zlámalová M, Tarábková H, Kavan L, Vajda Š. Atom by atom built subnanometer copper cluster catalyst for the highly selective oxidative dehydrogenation of cyclohexene. J Chem Phys 2022; 156:114302. [DOI: 10.1063/5.0065350] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The effect of particle size and support on the catalytic performance of supported subnanometer copper clusters was investigated in the oxidative dehydrogenation of cyclohexene. From among the investigated seven size-selected subnanometer copper particles between a single atom and clusters containing 2–7 atoms, the highest activity was observed for the titania-supported copper tetramer with 100% selectivity toward benzene production and being about an order of magnitude more active than not only all the other investigated cluster sizes on the same support but also the same tetramer on the other supports, Al2O3, SiO2, and SnO2. In addition to the profound effect of cluster size on activity and with Cu4 outstanding from the studied series, Cu4 clusters supported on SiO2 provide an example of tuning selectivity through support effects when this particular catalyst also produces cyclohexadiene with about 30% selectivity. Titania-supported Cu5 and Cu7 clusters supported on TiO2 produce a high fraction of cyclohexadiene in contrast to their neighbors, while Cu4 and Cu6 solely produce benzene without any combustion, thus representing odd–even oscillation of selectivity with the number of atoms in the cluster.
Collapse
Affiliation(s)
- Stanislav Valtera
- Department of Nanocatalysis, J. Heyrovský Institute of Physical Chemistry, v.v.i., Czech Academy of Sciences, Dolejškova 2155/3, CZ-182 23 Prague 8, Czech Republic
| | - Juraj Jašík
- Department of Nanocatalysis, J. Heyrovský Institute of Physical Chemistry, v.v.i., Czech Academy of Sciences, Dolejškova 2155/3, CZ-182 23 Prague 8, Czech Republic
| | - Mykhailo Vaidulych
- Department of Nanocatalysis, J. Heyrovský Institute of Physical Chemistry, v.v.i., Czech Academy of Sciences, Dolejškova 2155/3, CZ-182 23 Prague 8, Czech Republic
| | - Joanna Elżbieta Olszówka
- Department of Nanocatalysis, J. Heyrovský Institute of Physical Chemistry, v.v.i., Czech Academy of Sciences, Dolejškova 2155/3, CZ-182 23 Prague 8, Czech Republic
| | - Magda Zlámalová
- Department of Electrochemical Materials, J. Heyrovský Institute of Physical Chemistry, v.v.i., Czech Academy of Sciences, Dolejškova 2155/3, CZ-182 23 Prague 8, Czech Republic
- Department of Inorganic Chemistry, Faculty of Science, Charles University, Hlavova 2030, CZ-128 40 Prague, Czech Republic
| | - Hana Tarábková
- Department of Electrochemical Materials, J. Heyrovský Institute of Physical Chemistry, v.v.i., Czech Academy of Sciences, Dolejškova 2155/3, CZ-182 23 Prague 8, Czech Republic
| | - Ladislav Kavan
- Department of Electrochemical Materials, J. Heyrovský Institute of Physical Chemistry, v.v.i., Czech Academy of Sciences, Dolejškova 2155/3, CZ-182 23 Prague 8, Czech Republic
| | - Štefan Vajda
- Department of Nanocatalysis, J. Heyrovský Institute of Physical Chemistry, v.v.i., Czech Academy of Sciences, Dolejškova 2155/3, CZ-182 23 Prague 8, Czech Republic
| |
Collapse
|
19
|
Uptake of N 2O 5 by aqueous aerosol unveiled using chemically accurate many-body potentials. Nat Commun 2022; 13:1266. [PMID: 35273144 PMCID: PMC8913772 DOI: 10.1038/s41467-022-28697-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Accepted: 01/28/2022] [Indexed: 11/24/2022] Open
Abstract
The reactive uptake of N2O5 to aqueous aerosol is a major loss channel for nitrogen oxides in the troposphere. Despite its importance, a quantitative picture of the uptake mechanism is missing. Here we use molecular dynamics simulations with a data-driven many-body model of coupled-cluster accuracy to quantify thermodynamics and kinetics of solvation and adsorption of N2O5 in water. The free energy profile highlights that N2O5 is selectively adsorbed to the liquid–vapor interface and weakly solvated. Accommodation into bulk water occurs slowly, competing with evaporation upon adsorption from gas phase. Leveraging the quantitative accuracy of the model, we parameterize and solve a reaction–diffusion equation to determine hydrolysis rates consistent with experimental observations. We find a short reaction–diffusion length, indicating that the uptake is dominated by interfacial features. The parameters deduced here, including solubility, accommodation coefficient, and hydrolysis rate, afford a foundation for which to consider the reactive loss of N2O5 in more complex solutions. The reactive uptake of N2O5 to aqueous aerosol is a major loss channel for nitrogen oxides in the troposphere. Here authors report a theoretical investigation on the N2O5 uptake into aqueous aerosol and determine the hydrolysis rates by numerically solving a molecularly detailed reaction–diffusion equation.
Collapse
|
20
|
Laity PR, Holland C. Seeking Solvation: Exploring the Role of Protein Hydration in Silk Gelation. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27020551. [PMID: 35056868 PMCID: PMC8781151 DOI: 10.3390/molecules27020551] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Revised: 12/31/2021] [Accepted: 01/11/2022] [Indexed: 02/05/2023]
Abstract
The mechanism by which arthropods (e.g., spiders and many insects) can produce silk fibres from an aqueous protein (fibroin) solution has remained elusive, despite much scientific investigation. In this work, we used several techniques to explore the role of a hydration shell bound to the fibroin in native silk feedstock (NSF) from Bombyx mori silkworms. Small angle X-ray and dynamic light scattering (SAXS and DLS) revealed a coil size (radius of gyration or hydrodynamic radius) around 12 nm, providing considerable scope for hydration. Aggregation in dilute aqueous solution was observed above 65 °C, matching the gelation temperature of more concentrated solutions and suggesting that the strength of interaction with the solvent (i.e., water) was the dominant factor. Infrared (IR) spectroscopy indicated decreasing hydration as the temperature was raised, with similar changes in hydration following gelation by freezing or heating. It was found that the solubility of fibroin in water or aqueous salt solutions could be described well by a relatively simple thermodynamic model for the stability of the protein hydration shell, which suggests that the affected water is enthalpically favoured but entropically penalised, due to its reduced (vibrational or translational) dynamics. Moreover, while the majority of this investigation used fibroin from B. mori, comparisons with published work on silk proteins from other silkworms and spiders, globular proteins and peptide model systems suggest that our findings may be of much wider significance.
Collapse
|
21
|
Vazquez de Vasquez MG, Carter-Fenk KA, McCaslin LM, Beasley EE, Clark JB, Allen HC. Hydration and Hydrogen Bond Order of Octadecanoic Acid and Octadecanol Films on Water at 21 and 1 °C. J Phys Chem A 2021; 125:10065-10078. [PMID: 34761931 DOI: 10.1021/acs.jpca.1c06101] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The temperature-dependent hydration structure of long-chain fatty acids and alcohols at air-water interfaces has great significance in the fundamental interactions underlying ice nucleation in the atmosphere. We present an integrated theoretical and experimental study of the temperature-dependent vibrational structure and electric field character of the immediate hydration shells of fatty alcohol and acid headgroups. We use a combination of surface-sensitive infrared reflection-absorption spectroscopy (IRRAS), surface potentiometry, and ab initio molecular dynamics simulations to elucidate detailed molecular structures of the octadecanoic acid and octadecanol (stearic acid and stearyl alcohol) headgroup hydration shells at room temperature and near freezing. In experiments, the alcohol at high surface concentration exhibits the largest surface potential; yet we observe a strengthening of the hydrogen-bonding for the solvating water molecules near freezing for both the alcohol and the fatty acid IRRAS experiments. Results reveal that the hydration shells for both compounds screen their polar headgroup dipole moments reducing the surface potential at low surface coverages; at higher surface coverage, the polar headgroups become dehydrated, which reduces the screening, correlating to higher observed surface potential values. Lowering the temperature promotes tighter chain packing and an increase in surface potential. IRRAS reveals that the intra- and intermolecular vibrational coupling mechanisms are highly sensitive to changes in temperature. We find that intramolecular coupling dominates the vibrational relaxation pathways for interfacial water determined by comparing the H2O and the HOD spectra. Using ab initio molecular dynamics (AIMD) calculations on cluster systems of propanol + 6H2O and propionic acid + 10H2O, a spectral decomposition scheme was used to correlate the OH stretching motion with the IRRAS spectral features, revealing the effects of intra- and intermolecular coupling on the spectra. Spectra calculated with AIMD reproduce the red shift and increase in intensity observed in experimental spectra corresponding to the OH stretching region of the first solvation shell. These findings suggest that intra- and intermolecular vibrational couplings strongly impact the OH stretching region at fatty acid and fatty alcohol water interfaces. Overall, results are consistent with ice templating behavior for both the fatty acid and the alcohol, yet the surface potential signature is strongest for the fatty alcohol. These findings develop a better understanding of the complex surface potential and spectral signatures involved in ice templating.
Collapse
Affiliation(s)
- Maria G Vazquez de Vasquez
- Department of Chemistry & Biochemistry, The Ohio State University, 100 West 18th Avenue, Columbus, Ohio 43210, United States
| | - Kimberly A Carter-Fenk
- Department of Chemistry & Biochemistry, The Ohio State University, 100 West 18th Avenue, Columbus, Ohio 43210, United States
| | - Laura M McCaslin
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, United States
| | - Emma E Beasley
- Department of Chemistry & Biochemistry, The Ohio State University, 100 West 18th Avenue, Columbus, Ohio 43210, United States
| | - Jessica B Clark
- Department of Chemistry & Biochemistry, The Ohio State University, 100 West 18th Avenue, Columbus, Ohio 43210, United States
| | - Heather C Allen
- Department of Chemistry & Biochemistry, The Ohio State University, 100 West 18th Avenue, Columbus, Ohio 43210, United States
| |
Collapse
|
22
|
Kwac K, Freedman H, Cho M. Machine Learning Approach for Describing Water OH Stretch Vibrations. J Chem Theory Comput 2021; 17:6353-6365. [PMID: 34498885 DOI: 10.1021/acs.jctc.1c00540] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A machine learning approach employing neural networks is developed to calculate the vibrational frequency shifts and transition dipole moments of the symmetric and antisymmetric OH stretch vibrations of a water molecule surrounded by water molecules. We employed the atom-centered symmetry functions (ACSFs), polynomial functions, and Gaussian-type orbital-based density vectors as descriptor functions and compared their performances in predicting vibrational frequency shifts using the trained neural networks. The ACSFs perform best in modeling the frequency shifts of the OH stretch vibration of water among the types of descriptor functions considered in this paper. However, the differences in performance among these three descriptors are not significant. We also tried a feature selection method called CUR matrix decomposition to assess the importance and leverage of the individual functions in the set of selected descriptor functions. We found that a significant number of those functions included in the set of descriptor functions give redundant information in describing the configuration of the water system. We here show that the predicted vibrational frequency shifts by trained neural networks successfully describe the solvent-solute interaction-induced fluctuations of OH stretch frequencies.
Collapse
Affiliation(s)
- Kijeong Kwac
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science (IBS), Seoul 02841, Republic of Korea
| | - Holly Freedman
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science (IBS), Seoul 02841, Republic of Korea
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science (IBS), Seoul 02841, Republic of Korea.,Department of Chemistry, Korea University, Seoul 02841, Republic of Korea
| |
Collapse
|
23
|
Bull-Vulpe EF, Riera M, Götz AW, Paesani F. MB-Fit: Software infrastructure for data-driven many-body potential energy functions. J Chem Phys 2021; 155:124801. [PMID: 34598567 DOI: 10.1063/5.0063198] [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/14/2022] Open
Abstract
Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.
Collapse
Affiliation(s)
- Ethan F Bull-Vulpe
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
24
|
Gittus OR, Bresme F. Thermophysical properties of water using reactive force fields. J Chem Phys 2021; 155:114501. [PMID: 34551553 DOI: 10.1063/5.0057868] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The widescale importance and rich phenomenology of water continue to motivate the development of computational models. ReaxFF force fields incorporate many characteristics desirable for modeling aqueous systems: molecular flexibility, polarization, and chemical reactivity (bond formation and breaking). However, their ability to model the general properties of water has not been evaluated in detail. We present comprehensive benchmarks of the thermophysical properties of water for two ReaxFF models, the water-2017 and CHON-2017_weak force fields. These include structural, electrostatic, vibrational, thermodynamic, coexistence, and transport properties at ambient conditions (300 K and 0.997 g cm-3) and along the standard pressure (1 bar) isobar. Overall, CHON-2017_weak predicts more accurate thermophysical properties than the water-2017 force field. Based on our results, we recommend potential avenues for improvement: the dipole moment to quadrupole moment ratio, the self-diffusion coefficient, especially for water-2017, and the gas phase vibrational frequencies with the aim to improve the vibrational properties of liquid water.
Collapse
Affiliation(s)
- Oliver R Gittus
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London W12 0BZ, United Kingdom
| | - Fernando Bresme
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London W12 0BZ, United Kingdom
| |
Collapse
|
25
|
Benson RL, Althorpe SC. On the "Matsubara heating" of overtone intensities and Fermi splittings. J Chem Phys 2021; 155:104107. [PMID: 34525826 DOI: 10.1063/5.0056829] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Classical molecular dynamics (MD) and imaginary-time path-integral dynamics methods underestimate the infrared absorption intensities of overtone and combination bands by typically an order of magnitude. Plé et al. [J. Chem. Phys. 155, 2863 (2021)] have shown that this is because such methods fail to describe the coupling of the centroid to the Matsubara dynamics of the fluctuation modes; classical first-order perturbation theory (PT) applied to the Matsubara dynamics is sufficient to recover most of the lost intensity in simple models and gives identical results to quantum (Rayleigh-Schrödinger) PT. Here, we show numerically that the results of this analysis can be used as post-processing correction factors, which can be applied to realistic (classical MD or path-integral dynamics) simulations of infrared spectra. We find that the correction factors recover most of the lost intensity in the overtone and combination bands of gas-phase water and ammonia and much of it for liquid water. We then re-derive and confirm the earlier PT analysis by applying canonical PT to Matsubara dynamics, which has the advantage of avoiding secular terms and gives a simple picture of the perturbed Matsubara dynamics in terms of action-angle variables. Collectively, these variables "Matsubara heat" the amplitudes of the overtone and combination vibrations of the centroid to what they would be in a classical system with the oscillators (of frequency Ωi) held at their quantum effective temperatures [of ℏΩi coth(βℏΩi/2)/2kB]. Numerical calculations show that a similar neglect of "Matsubara heating" causes path-integral methods to underestimate Fermi resonance splittings.
Collapse
Affiliation(s)
- Raz L Benson
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stuart C Althorpe
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
26
|
Caruso A, Paesani F. Data-driven many-body models enable a quantitative description of chloride hydration from clusters to bulk. J Chem Phys 2021; 155:064502. [PMID: 34391363 DOI: 10.1063/5.0059445] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
We present a new data-driven potential energy function (PEF) describing chloride-water interactions, which is developed within the many-body-energy (MB-nrg) theoretical framework. Besides quantitatively reproducing low-order many-body energy contributions, the new MB-nrg PEF is able to correctly predict the interaction energies of small chloride-water clusters calculated at the coupled cluster level of theory. Importantly, classical and quantum molecular dynamics simulations of a single chloride ion in water demonstrate that the new MB-nrg PEF predicts x-ray spectra in close agreement with the experimental results. Comparisons with an popular empirical model and a polarizable PEF emphasize the importance of an accurate representation of short-range many-body effect while demonstrating that pairwise additive representations of chloride-water and water-water interactions are inadequate for correctly representing the hydration structure of chloride in both gas-phase clusters and solution. We believe that the analyses presented in this study provide additional evidence for the accuracy and predictive ability of the MB-nrg PEFs, which can then enable more realistic simulations of ionic aqueous systems in different environments.
Collapse
Affiliation(s)
- Alessandro Caruso
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
27
|
Cruzeiro VWD, Wildman A, Li X, Paesani F. Relationship between Hydrogen-Bonding Motifs and the 1b 1 Splitting in the X-ray Emission Spectrum of Liquid Water. J Phys Chem Lett 2021; 12:3996-4002. [PMID: 33877847 DOI: 10.1021/acs.jpclett.1c00486] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The split of the 1b1 peak observed in the X-ray emission (XE) spectrum of liquid water has been the focus of intense research. Although several hypotheses have been proposed to explain the origin of this split, a consensus has not yet been reached. Here, we introduce a novel theoretical/computation approach which, combining path-integral molecular dynamics simulations with the MB-pol model and time-dependent density functional theory calculations, predicts the 1b1 splitting in liquid water and not in crystalline ice, in agreement with the experimental observations. A systematic analysis of the underlying local structure of liquid water at ambient conditions indicates that several different hydrogen-bonding motifs contribute to the overall XE line shape in the energy range corresponding to emissions from the 1b1 orbitals. This suggests that it is not possible to unambiguously attribute the split of the 1b1 peak to only two specific structural arrangements of the underlying hydrogen-bonding network.
Collapse
Affiliation(s)
- Vinícius Wilian D Cruzeiro
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, United States
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Andrew Wildman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
- Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Francesco Paesani
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, United States
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California, San Diego, La Jolla, California 92093, United States
| |
Collapse
|
28
|
Fischer TL, Wagner T, Gottschalk HC, Nejad A, Suhm MA. A Rather Universal Vibrational Resonance in 1:1 Hydrates of Carbonyl Compounds. J Phys Chem Lett 2021; 12:138-144. [PMID: 33315407 DOI: 10.1021/acs.jpclett.0c03197] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
When the lower frequency OH stretching fundamental of a water molecule is shifted to the 3500 cm-1 spectral range by the solvation of a carbonyl compound, in this case a ketone, its infrared intensity is shared with a dark state. It is shown by chemical and isotope substitution for more than a dozen systems that the location of this resonance is remarkably substitution-independent. Harmonic and anharmonic model calculations support its assignment to a combination of the water bending overtone and in-plane water libration. This previously unrecognized intramolecular-intermolecular coupling in single solvent water has a strength of 7-10 cm-1. It may have been sporadically observed before in a few other carbonyl compounds such as amides, without any previous exploration of its potential universality. The resulting generic picosecond energy redistribution channel for aqueous solutions may represent a slow counterpart and doorway model of what happens on a subpicosecond time scale when the hydrogen bonds become stronger, such as in carboxylic acid dimers or protonated water clusters.
Collapse
Affiliation(s)
- Taija L Fischer
- Institute of Physical Chemistry, University of Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| | - Till Wagner
- Institute of Physical Chemistry, University of Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| | - Hannes C Gottschalk
- Institute of Physical Chemistry, University of Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| | - Arman Nejad
- Institute of Physical Chemistry, University of Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| | - Martin A Suhm
- Institute of Physical Chemistry, University of Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
| |
Collapse
|
29
|
Lambros E, Paesani F. How good are polarizable and flexible models for water: Insights from a many-body perspective. J Chem Phys 2020; 153:060901. [DOI: 10.1063/5.0017590] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
30
|
Riera M, Talbot JJ, Steele RP, Paesani F. Infrared signatures of isomer selectivity and symmetry breaking in the Cs+(H2O)3 complex using many-body potential energy functions. J Chem Phys 2020; 153:044306. [DOI: 10.1063/5.0013101] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Justin J. Talbot
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Ryan P. Steele
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
31
|
Guo J, Cao D, Chen J, Bian K, Xu LM, Wang EG, Jiang Y. Probing the intermolecular coupled vibrations in a water cluster with inelastic electron tunneling spectroscopy. J Chem Phys 2020; 152:234301. [PMID: 32571057 DOI: 10.1063/5.0009385] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The hydrogen-bonding networks of water have strong intra- and intermolecular vibrational coupling which influences the energy dissipation and proton transfer in water. Disentangling and quantitative characterization of different coupling effects in water at a single-molecular level still remains a great challenge. Using tip-enhanced inelastic electron tunneling spectroscopy (IETS) based on low-temperature scanning tunneling microscopy, we report the direct quantitative assessment of the intermolecular coupling constants of the OH-stretch vibrational bands of an isolated water tetramer adsorbed on a Au(111)-supported NaCl(001) bilayer film. This is achieved by distinguishing various coupled modes of the H-bonded O-H stretching vibrations through tip-height dependent IET spectra. In contrast, such vibrational coupling is negligible in the half-deuterated water tetramer owing to the large energy mismatch between the OH and OD stretching modes. Not only do these findings advance our understanding on the effects of local environment on the intermolecular vibrational coupling in water, but also open up a new route for vibrational spectroscopic studies of extended H-bonded network at the single-molecular level.
Collapse
Affiliation(s)
- Jing Guo
- College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Duanyun Cao
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Ji Chen
- School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Ke Bian
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Li-Mei Xu
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - En-Ge Wang
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Ying Jiang
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, People's Republic of China
| |
Collapse
|
32
|
Zhai Y, Caruso A, Gao S, Paesani F. Active learning of many-body configuration space: Application to the Cs+–water MB-nrg potential energy function as a case study. J Chem Phys 2020; 152:144103. [DOI: 10.1063/5.0002162] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Yaoguang Zhai
- Department of Computer Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
| | - Alessandro Caruso
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Sicun Gao
- Department of Computer Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
33
|
Benson RL, Trenins G, Althorpe SC. Which quantum statistics–classical dynamics method is best for water? Faraday Discuss 2020; 221:350-366. [DOI: 10.1039/c9fd00077a] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We make a systematic comparison of TRPMD, CMD, QCMD and LSC-IVR by calculating the infrared spectrum of water in its gas, liquid and ice phases.
Collapse
Affiliation(s)
- Raz L. Benson
- Department of Chemistry
- University of Cambridge
- Cambridge
- UK
| | - George Trenins
- Department of Chemistry
- University of Cambridge
- Cambridge
- UK
| | | |
Collapse
|
34
|
Multistate Reactive Molecular Dynamics Simulations of Proton Diffusion in Water Clusters and in the Bulk. J Phys Chem B 2019; 123:9846-9861. [PMID: 31647873 DOI: 10.1021/acs.jpcb.9b03258] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The molecular mechanics with proton transfer (MMPT) force field is combined with multistate adiabatic reactive molecular dynamics (MS-ARMD) to describe proton transport in the condensed phase. Parametrization for small protonated water clusters based on electronic structure calculations at the MP2/6-311+G(2d,2p) level of theory and refinement by comparing with infrared spectra for a protonated water tetramer yields a force field which faithfully describes the minimum energy structures of small protonated water clusters. In protonated water clusters up to (H2O)100H+, the proton hopping rate is around 100 hops/ns. This rate converges for 21 ≤ n ≤ 31, and no further speedup in bulk water is found. This indicates that bulklike behavior requires the solvation of a Zundel motif by ∼25 water molecules, which corresponds to the second solvation sphere. For smaller cluster sizes, the number of available states (i.e., the number of proton acceptors) is too small and slows down proton-transfer rates. The cluster simulations confirm that the excess proton is typically located on the surface. The free-energy surface as a function of the weights of the two lowest states and a configurational parameter suggests that the "special pair" plays a central role in rapid proton transport. The barriers between this minimum-energy structure and the Zundel and Eigen minima are sufficiently low (∼1 kcal/mol, consistent with recent experiments and commensurate with a hopping rate of ∼100/ns or 1 every 10 ps), leading to a highly dynamic environment. These findings are also consistent with recent experiments which find that Zundel-type hydration geometries are prevalent in bulk water.
Collapse
|
35
|
Rieth AJ, Wright AM, Skorupskii G, Mancuso JL, Hendon CH, Dincă M. Record-Setting Sorbents for Reversible Water Uptake by Systematic Anion Exchanges in Metal-Organic Frameworks. J Am Chem Soc 2019; 141:13858-13866. [PMID: 31398286 PMCID: PMC6748661 DOI: 10.1021/jacs.9b06246] [Citation(s) in RCA: 62] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The reversible capture of water vapor at low humidity can enable transformative applications such as atmospheric water harvesting and heat transfer that uses water as a refrigerant, replacing environmentally detrimental hydro- and chloro-fluorocarbons. The driving force for these applications is governed by the relative humidity at which the pores of a porous material fill with water. Here, we demonstrate modulation of the onset of pore-filling in a family of metal-organic frameworks with record water sorption capacities by employing anion exchange. Unexpectedly, the replacement of the structural bridging Cl- with the more hydrophilic anions F- and OH- does not induce pore-filling at lower relative humidity, whereas the introduction of the larger Br- results in a substantial shift toward lower relative humidity. We rationalize these results in terms of pore size modifications as well as the water hydrogen bonding structure based on detailed infrared spectroscopic measurements. Fundamentally, our data suggest that, in the presence of strong nucleation sites, the thermodynamic favorability of water pore-filling depends more strongly on the pore diameter and the interface between water in the center of the pore and water bound to the pore walls than the hydrophilicity of the pore wall itself. On the basis of these results, we report two materials that exhibit record water uptake capacities in their respective humidity regions and extended stability over 400 water adsorption-desorption cycles.
Collapse
Affiliation(s)
- Adam J Rieth
- Department of Chemistry , Massachusetts Institute of Technology , 77 Massachusetts Avenue , Cambridge , Massachusetts 02139 , United States
| | - Ashley M Wright
- Department of Chemistry , Massachusetts Institute of Technology , 77 Massachusetts Avenue , Cambridge , Massachusetts 02139 , United States
| | - Grigorii Skorupskii
- Department of Chemistry , Massachusetts Institute of Technology , 77 Massachusetts Avenue , Cambridge , Massachusetts 02139 , United States
| | - Jenna L Mancuso
- Materials Science Institute, Department of Chemistry and Biochemistry , University of Oregon , Eugene , Oregon 97403 , United States
| | - Christopher H Hendon
- Materials Science Institute, Department of Chemistry and Biochemistry , University of Oregon , Eugene , Oregon 97403 , United States
| | - Mircea Dincă
- Department of Chemistry , Massachusetts Institute of Technology , 77 Massachusetts Avenue , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
36
|
Trenins G, Willatt MJ, Althorpe SC. Path-integral dynamics of water using curvilinear centroids. J Chem Phys 2019. [DOI: 10.1063/1.5100587] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- George Trenins
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Michael J. Willatt
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Stuart C. Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
37
|
Kananenka AA, Hestand NJ, Skinner JL. OH-Stretch Raman Multivariate Curve Resolution Spectroscopy of HOD/H2O Mixtures. J Phys Chem B 2019; 123:5139-5146. [DOI: 10.1021/acs.jpcb.9b02686] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Alexei A. Kananenka
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| | - Nicholas J. Hestand
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| | - J. L. Skinner
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
38
|
Sengupta S, Moberg DR, Paesani F, Tyrode E. Neat Water-Vapor Interface: Proton Continuum and the Nonresonant Background. J Phys Chem Lett 2018; 9:6744-6749. [PMID: 30407831 DOI: 10.1021/acs.jpclett.8b03069] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Whether the surface of neat water is "acidic" or "basic" remains an active and controversial field of research. Most of the experimental evidence supporting the preferential adsorption of H3O+ ions stems from nonlinear optical spectroscopy methods typically carried out at extreme pH conditions (pH < 1). Here, we use vibrational sum frequency spectroscopy (VSFS) to target the "proton continuum", an unexplored frequency range characteristic of hydrated protons and hydroxide ions. The VSFS spectra of neat water show a broad and nonzero signal intensity between 1700 and 3000 cm-1 in the three different polarization combinations examined. By comparing the SF response of water with that from dilute HCl and NaOH aqueous solutions, we conclude the intensity does not originate from either adsorbed H3O+ or OH- ions. Contributions from the nonresonant background are then critically considered by comparing the experimental results with many-body molecular dynamics (MB-MD) simulated spectra.
Collapse
Affiliation(s)
- Sanghamitra Sengupta
- Department of Chemistry , KTH Royal Institute of Technology , SE-10044 Stockholm , Sweden
| | | | | | - Eric Tyrode
- Department of Chemistry , KTH Royal Institute of Technology , SE-10044 Stockholm , Sweden
| |
Collapse
|