1
|
Hou R, Li C, Pan D. Raman and IR spectra of water under graphene nanoconfinement at ambient and extreme pressure-temperature conditions: a first-principles study. Faraday Discuss 2024; 249:181-194. [PMID: 37791622 DOI: 10.1039/d3fd00111c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The nanoconfinement of water can result in dramatic differences in its physical and chemical properties compared to bulk water. However, a detailed molecular-level understanding of these properties is still lacking. Vibrational spectroscopy, such as Raman and infrared, is a popular experimental tool for studying the structure and dynamics of water, and is often complemented by atomistic simulations to interpret experimental spectra, but there have been few theoretical spectroscopy studies of nanoconfined water using first-principles methods at ambient conditions, let alone under extreme pressure-temperature conditions. Here, we compute the Raman and IR spectra of water nanoconfined by graphene at ambient and extreme pressure-temperature conditions using ab initio simulations. Our results revealed alterations in the Raman stretching and low-frequency bands due to the graphene confinement. We also found spectroscopic evidence indicating that nanoconfinement considerably changes the tetrahedral hydrogen bond network, which is typically found in bulk water. Furthermore, we observed an unusual bending band in the Raman spectrum at ∼10 GPa and 1000 K, which is attributed to the unique molecular structure of confined ionic water. Additionally, we found that at ∼20 GPa and 1000 K, confined water transformed into a superionic fluid, making it challenging to identify the IR stretching band. Finally, we computed the ionic conductivity of confined water in the ionic and superionic phases. Our results highlight the efficacy of Raman and IR spectroscopy in studying the structure and dynamics of nanoconfined water in a large pressure-temperature range. Our predicted Raman and IR spectra can serve as a valuable guide for future experiments.
Collapse
Affiliation(s)
- Rui Hou
- Department of Physics, Hong Kong University of Science and Technology, Hong Kong, China.
- HKUST Shenzhen-Hong Kong Collaborative Innovation Research Institute, Shenzhen, China
| | - Chu Li
- Department of Physics, Hong Kong University of Science and Technology, Hong Kong, China.
- HKUST Shenzhen-Hong Kong Collaborative Innovation Research Institute, Shenzhen, China
| | - Ding Pan
- Department of Physics, Hong Kong University of Science and Technology, Hong Kong, China.
- Department of Chemistry, Hong Kong University of Science and Technology, Hong Kong, China
- HKUST Shenzhen-Hong Kong Collaborative Innovation Research Institute, Shenzhen, China
| |
Collapse
|
2
|
Zhang C, Puligheddu M, Zhang L, Car R, Galli G. Thermal Conductivity of Water at Extreme Conditions. J Phys Chem B 2023; 127:7011-7017. [PMID: 37524047 PMCID: PMC10424233 DOI: 10.1021/acs.jpcb.3c02972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2023] [Revised: 07/06/2023] [Indexed: 08/02/2023]
Abstract
Measuring the thermal conductivity (κ) of water at extreme conditions is a challenging task, and few experimental data are available. We predict κ for temperatures and pressures relevant to the conditions of the Earth mantle, between 1,000 and 2,000 K and up to 22 GPa. We employ close to equilibrium molecular dynamics simulations and a deep neural network potential fitted to density functional theory data. We then interpret our results by computing the equation of state of water on a fine grid of points and using a simple model for κ. We find that the thermal conductivity is weakly dependent on temperature and monotonically increases with pressure with an approximate square-root behavior. In addition, we show how the increase of κ at high pressure, relative to ambient conditions, is related to the corresponding increase in the sound velocity. Although the relationships between the thermal conductivity, pressure and sound velocity established here are not rigorous, they are sufficiently accurate to allow for a robust estimate of the thermal conductivity of water in a broad range of temperatures and pressures, where experiments are still difficult to perform.
Collapse
Affiliation(s)
- Cunzhi Zhang
- Pritzker
School of Molecular Engineering, University
of Chicago, Chicago, Illinois 60637, United States
| | - Marcello Puligheddu
- Pritzker
School of Molecular Engineering, University
of Chicago, Chicago, Illinois 60637, United States
- Materials
Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Linfeng Zhang
- Program
in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, United States
| | - Roberto Car
- Program
in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, United States
- Department
of Chemistry, Department of Physics, and Princeton Institute for the
Science and Technology of Materials, Princeton
University, Princeton, New Jersey 08544, United States
| | - Giulia Galli
- Pritzker
School of Molecular Engineering, University
of Chicago, Chicago, Illinois 60637, United States
- Materials
Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, United States
- Department
of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
3
|
Ye Z, Zhang C, Galli G. Photoelectron spectra of water and simple aqueous solutions at extreme conditions. Faraday Discuss 2022; 236:352-363. [DOI: 10.1039/d2fd00003b] [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
Determining the electronic structure of aqueous solutions at extreme conditions is an important step towards understanding chemical bonding and reactions in water under pressure (P) and at high temperature (T)....
Collapse
|
4
|
Resta R. Faraday law, oxidation numbers, and ionic conductivity: The role of topology. J Chem Phys 2021; 155:244503. [PMID: 34972381 DOI: 10.1063/5.0077718] [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/15/2022] Open
Abstract
Faraday's experiment measures-within a modern view-the charge adiabatically transported over a macroscopic distance by a given nuclear species in insulating liquids: the reason why it is an integer is deeply rooted in topology. Whole numbers enter chemistry in a different form: atomic oxidation states. They are not directly measurable and are determined instead from an agreed set of rules. Insulating liquids are a remarkable exception; Faraday's experiment indeed measures the oxidation numbers of each dissociated component in the liquid phase, whose topological values are unambiguous. Ionic conductivity in insulating liquids is expressed in terms of the autocorrelation function of the fluctuating charge current at a given temperature in a zero electric field; topology plays a major role in this important observable as well. The existing literature deals with the above issues by adopting the independent-electron framework; here, I provide the many-body generalization of all the above findings, which, furthermore, allows for compact and very transparent notations and formulas.
Collapse
Affiliation(s)
- Raffaele Resta
- Istituto Officina dei Materiali IOM-CNR, Strada Costiera 11, 34151 Trieste, Italy and Donostia International Physics Center, 20018 San Sebastián, Spain
| |
Collapse
|
5
|
Cassone G, Sponer J, Saija F. Molecular dissociation and proton transfer in aqueous methane solution under an electric field. Phys Chem Chem Phys 2021; 23:25649-25657. [PMID: 34782902 DOI: 10.1039/d1cp04202e] [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
Methane-water mixtures are ubiquitous in our solar system and they have been the subject of a wide variety of experimental, theoretical, and computational studies aimed at understanding their behaviour under disparate thermodynamic scenarios, up to extreme planetary ice conditions of pressures and temperatures [Lee and Scandolo, Nat. Commun., 2011, 2, 185]. Although it is well known that electric fields, by interacting with condensed matter, can produce a range of catalytic effects which can be similar to those observed when material systems are pressurised, to the best of our knowledge, no quantum-based computational investigations of methane-water mixtures under an electric field have been reported so far. Here we present a study relying upon state-of-the-art ab initio molecular dynamics simulations where a liquid aqueous methane solution is exposed to strong oriented static and homogeneous electric fields. It turns out that a series of field-induced effects on the dipoles, polarisation, and the electronic structure of both methane and water molecules are recorded. Moreover, upon increasing the field strength, increasing fractions of water molecules are not only re-oriented towards the field direction, but are also dissociated by the field, leading to the release of oxonium and hydroxyde ions in the mixture. However, in contrast to what is observed upon pressurisation (∼50 GPa), where the presence of the water counterions triggers methane ionisation and other reactions, methane molecules preserve their integrity up to the strongest field explored (i.e., 0.50 V Å-1). Interestingly, neither the field-induced molecular dissociation of neat water (i.e., 0.30 V Å-1) nor the proton conductivity typical of pure aqueous samples at these field regimes (i.e., 1.3 S cm-1) are affected by the presence of hydrophobic interactions, at least in a methane-water mixture containing a molar fraction of 40% methane.
Collapse
Affiliation(s)
- Giuseppe Cassone
- Institute for Chemical-Physical Processes, National Research Council of Italy (IPCF-CNR), Viale F. Stagno d'Alcontres 37, 98158 Messina, Italy.
| | - Jiri Sponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolska 135, 61265 Brno, Czech Republic
| | - Franz Saija
- Institute for Chemical-Physical Processes, National Research Council of Italy (IPCF-CNR), Viale F. Stagno d'Alcontres 37, 98158 Messina, Italy.
| |
Collapse
|
6
|
Thomsen B, Shiga M. Ab initio study of nuclear quantum effects on sub- and supercritical water. J Chem Phys 2021; 155:194107. [PMID: 34800944 DOI: 10.1063/5.0071857] [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/15/2022] Open
Abstract
The structures of water in the ambient, subcritical, and supercritical conditions at various densities were studied systematically by ab initio path integral molecular dynamics simulations. It was found that the nuclear quantum effects (NQEs) have a significant impact on the structure of hydrogen bonds in close contact, not only in the ambient condition but also in the sub- and supercritical conditions. The NQEs on the structure beyond the hydrogen bond contact are important in ambient water, but not much for water in the sub- and supercritical conditions. The NQEs are furthermore important for determining the number of hydrogen bonds in the ambient conditions, and this role is, however, diminished in the sub- and supercritical conditions. The NQEs do, nevertheless, show their importance in determining the intramolecular structure of water and the close contact structures of the hydrogen bonds, even at sub- and supercritical conditions. Using the RPBE-D3 functional, the computed radial distribution functions for ambient water are in excellent agreement with experimental data, upgrading our previous results using the BLYP-D2 functional [Machida et al., J. Chem. Phys. 148, 102324 (2018)]. The computed radial distribution functions for water in the sub- and supercritical conditions were carefully compared with experiment. In particular, we found that the first peak in hydrogen pair distribution functions matches only when the NQEs are taken into account.
Collapse
Affiliation(s)
- Bo Thomsen
- CCSE, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| | - Motoyuki Shiga
- CCSE, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| |
Collapse
|
7
|
Cassone G, Sponer J, Saija F. Ab Initio Molecular Dynamics Studies of the Electric-Field-Induced Catalytic Effects on Liquids. Top Catal 2021. [DOI: 10.1007/s11244-021-01487-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
8
|
Grasselli F, Baroni S. Invariance principles in the theory and computation of transport coefficients. THE EUROPEAN PHYSICAL JOURNAL. B 2021; 94:160. [PMID: 34776779 PMCID: PMC8550620 DOI: 10.1140/epjb/s10051-021-00152-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/04/2021] [Accepted: 07/05/2021] [Indexed: 06/13/2023]
Abstract
ABSTRACT In this work, we elaborate on two recently discovered invariance principles, according to which transport coefficients are, to a large extent, independent of the microscopic definition of the densities and currents of the conserved quantities being transported (energy, momentum, mass, charge). The first such principle, gauge invariance, allows one to define a quantum adiabatic energy current from density-functional theory, from which the heat conductivity can be uniquely defined and computed using equilibrium ab initio molecular dynamics. When combined with a novel topological definition of atomic oxidation states, gauge invariance also sheds new light onto the mechanisms of charge transport in ionic conductors. The second principle, convective invariance, allows one to extend the analysis to multi-component systems. These invariance principles can be combined with new spectral analysis methods for the current time series to be fed into the Green-Kubo formula to obtain accurate estimates of transport coefficients from relatively short molecular dynamics simulations.
Collapse
Affiliation(s)
- Federico Grasselli
- COSMO–Laboratory of Computational Science and Modelling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Stefano Baroni
- SISSA–Scuola Internazionale Superiore di Studi Avanzati, 34136 EU Trieste, Italy
- CNR-IOM DEMOCRITOS Simulation Center, SISSA, 34136 Trieste, EU Italy
| |
Collapse
|
9
|
de Almeida JM, Nguyen NL, Colonna N, Chen W, Rodrigues Miranda C, Pasquarello A, Marzari N. Electronic Structure of Water from Koopmans-Compliant Functionals. J Chem Theory Comput 2021; 17:3923-3930. [PMID: 34137253 DOI: 10.1021/acs.jctc.1c00063] [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/28/2022]
Abstract
Obtaining a precise theoretical description of the spectral properties of liquid water poses challenges for both molecular dynamics (MD) and electronic structure methods. The lower computational cost of the Koopmans-compliant functionals with respect to Green's function methods allows the simulations of many MD trajectories, with a description close to the state-of-art quasi-particle self-consistent GW plus vertex corrections method (QSGW + fxc). Thus, we explore water spectral properties when different MD approaches are used, ranging from classical MD to first-principles MD, and including nuclear quantum effects. We have observed that different MD approaches lead to up to 1 eV change in the average band gap; thus, we focused on the band gap dependence with the geometrical properties of a system to explain such spread. We have evaluated the changes in the band gap due to variations in the intramolecular O-H bond distance and HOH angle, as well as the intermolecular hydrogen bond O···O distance and the OHO angles. We have observed that the dominant contribution comes from the O-H bond length; the O···O distance plays a secondary role, and the other geometrical properties do not significantly influence the gap. Furthermore, we analyze the electronic density of states (DOS), where the KIPZ functional shows good agreement with the DOS obtained with state-of-art approaches employing quasi-particle self-consistent GW plus vertex corrections. The O-H bond length also significantly influences the DOS. When nuclear quantum effects are considered, broadening of the peaks driven by the broader distribution of the O-H bond lengths is observed, leading to a closer agreement with the experimental photoemission spectra.
Collapse
Affiliation(s)
- James Moraes de Almeida
- Universidade Federal do ABC, Centro de Ciências Naturais e Humanas, Santo André, 09210-580 SP, Brazil.,Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Ngoc Linh Nguyen
- Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Nicola Colonna
- Laboratory for Neutron Scattering and Imaging, Paul Scherrer Institute, CH-5232 Villigen-PSI, Switzerland.,National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Wei Chen
- Institute of Condensed Matter and Nanoscience, Université Catholique de Louvain, B-1348 Louvain-la-Neuve, Belgium
| | | | - Alfredo Pasquarello
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Nicola Marzari
- Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
10
|
Zhang L, Wang H, Car R, E W. Phase Diagram of a Deep Potential Water Model. PHYSICAL REVIEW LETTERS 2021; 126:236001. [PMID: 34170175 DOI: 10.1103/physrevlett.126.236001] [Citation(s) in RCA: 89] [Impact Index Per Article: 29.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Accepted: 04/28/2021] [Indexed: 06/13/2023]
Abstract
Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. The fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII^{''} and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.
Collapse
Affiliation(s)
- Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA
| | - Han Wang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Fenghao East Road 2, Beijing 100094, People's Republic of China
| | - Roberto Car
- Department of Chemistry, Department of Physics, Program in Applied and Computational Mathematics, Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08544, USA
| | - Weinan E
- Department of Mathematics and Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA and Beijing Institute of Big Data Research, Beijing 100871, People's Republic of China
| |
Collapse
|
11
|
Tripathi R, Durán Caballero L, Pérez de Tudela R, Hölzl C, Marx D. Unveiling Zwitterionization of Glycine in the Microhydration Limit. ACS OMEGA 2021; 6:12676-12683. [PMID: 34056419 PMCID: PMC8154221 DOI: 10.1021/acsomega.1c00869] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/17/2021] [Accepted: 03/24/2021] [Indexed: 06/12/2023]
Abstract
Charge separation under solvation stress conditions is a fundamental process that comes in many forms in doped water clusters. Yet, the mechanism of intramolecular charge separation, where constraints due to the molecular structure might be intricately tied to restricted solvation structures, remains largely unexplored. Microhydrated amino acids are such paradigmatic molecules. Ab initio simulations are carried out at 300 K in the frameworks of metadynamics sampling and thermodynamic integration to map the thermal mechanisms of zwitterionization using Gly(H2O) n with n = 4 and 10. In both cases, a similar water-mediated proton transfer chain mechanism is observed; yet, detailed analyses of thermodynamics and kinetics demonstrate that the charge-separated zwitterion is the preferred species only for n = 10 mainly due to kinetic stabilization. Structural analyses disclose that bifurcated H-bonded water bridges, connecting the cationic and anionic sites in the fluctuating microhydration network at room temperature, are enhanced in the transition-state ensemble exclusively for n = 10 and become overwhelmingly abundant in the stable zwitterion. The findings offer potential insights into charge separation under solvation stress conditions beyond the present example.
Collapse
|
12
|
Rozsa V, Galli G. Solvation of simple ions in water at extreme conditions. J Chem Phys 2021; 154:144501. [PMID: 33858154 DOI: 10.1063/5.0046193] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The interaction of ions and water at high pressure and temperature plays a critical role in Earth and planetary science yet remains poorly understood. Aqueous fluids affect geochemical properties ranging from water phase stability to mineral solubility and reactivity. Here, we report first-principles molecular dynamics simulations of mono-valent ions (Li+, K+, Cl-) as well as NaCl in liquid water at temperatures and pressures relevant to the Earth's upper mantle (11 GPa, 1000 K) and concentrations in the dilute limit (0.44-0.88 m), in the regime of ocean salinity. We find that, at extreme conditions, the average structural and vibrational properties of water are weakly affected by the presence of ions, beyond the first solvation shell, similar to what was observed at ambient conditions. We also find that the ionic conductivity of the liquid increases in the presence of ions by less than an order of magnitude and that the dielectric constant is moderately reduced by at most ∼10% at these conditions. Our findings may aid in the parameterization of deep earth water models developed to describe water-rock reactions.
Collapse
Affiliation(s)
- Viktor Rozsa
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| | - Giulia Galli
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
13
|
Henao A, Salazar-Rios JM, Guardia E, Pardo LC. Structure and dynamics of water plastic crystals from computer simulations. J Chem Phys 2021; 154:104501. [PMID: 33722053 DOI: 10.1063/5.0038762] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Water has a rich phase diagram with several crystals, as confirmed by experiments. High-pressure and high-temperature water is of interest for Earth's mantle and exoplanetary investigations. It is in this region of the phase diagram of water that new plastic crystal phases of water have been revealed via computer simulations by both classical forcefields and ab initio calculations. However, these plastic phases still remain elusive in experiments. Here, we present a complete characterization of the structure, dynamics, and thermodynamics of the computational plastic crystal phases of water using molecular dynamics and the two-phase thermodynamic method and uncover the interplay between them. The relaxation times of different reorientational correlation functions are obtained for the hypothetical body-centered-cubic and face-centered-cubic plastic crystal phases of water at T = 440 K and P = 8 GPa. Results are compared to a high pressure liquid and ice VII phases to improve the understanding of the plastic crystal phases. Entropy results indicate that the fcc crystal is more stable compared to the bcc structure under the studied conditions.
Collapse
Affiliation(s)
- Andrés Henao
- Grup de Simulació per Ordinador en Matèria Condensada, Departament de Física, B4-B5 Campus Nord, Universitat Politècnica de Catalunya, E-08034 Barcelona, Catalonia, Spain
| | | | - Elvira Guardia
- Grup de Simulació per Ordinador en Matèria Condensada, Departament de Física, B4-B5 Campus Nord, Universitat Politècnica de Catalunya, E-08034 Barcelona, Catalonia, Spain
| | - Luis C Pardo
- Grup de Caracterització de Materials, Departament de Física, EEBE, Universitat Politècnica de Catalunya, Avgda. Eduard Maristany 16, E-080197 Barcelona, Catalonia, Spain
| |
Collapse
|
14
|
Rozsa VF, Galli G. Molecular Polarizabilities in Aqueous Systems from First-Principles. J Phys Chem B 2021; 125:2183-2192. [DOI: 10.1021/acs.jpcb.0c10732] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Viktor F. Rozsa
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Giulia Galli
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Materials Science Division, Argonne National Laboratory, Chicago, Illinois 60439, United States
| |
Collapse
|
15
|
Cassone G. Nuclear Quantum Effects Largely Influence Molecular Dissociation and Proton Transfer in Liquid Water under an Electric Field. J Phys Chem Lett 2020; 11:8983-8988. [PMID: 33035059 DOI: 10.1021/acs.jpclett.0c02581] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Proton transfer in liquid water controls acid-base chemistry, crucial enzyme reactions, and the functioning of fuel cells. Externally applied static electric fields in water are capable of dissociating molecules and transferring protons across the H-bond network. However, the impact of nuclear quantum effects (NQEs) on these fundamental field-induced phenomena has not yet been reported. By comparing state-of-the-art ab initio molecular dynamics (AIMD) and path integral AIMD simulations of water under electric fields, I show that quantum delocalization of the proton lowers the molecular ionization threshold to approximately one-third. Moreover, also the water behavior as a protonic semiconductor is considerably modified by the inclusion of NQEs. In fact, when the quantum nature of the nuclei is taken into account, the proton conductivity is ∼50% larger. This work proves that NQEs sizably affect the protolysis phenomenon and proton transfer in room-temperature liquid water.
Collapse
Affiliation(s)
- Giuseppe Cassone
- Institute for Chemical-Physical Processes, National Research Council, Viale F. Stagno d'Alcontres 37, 98158 Messina, Italy
| |
Collapse
|
16
|
Heat and charge transport in H 2O at ice-giant conditions from ab initio molecular dynamics simulations. Nat Commun 2020; 11:3605. [PMID: 32681002 PMCID: PMC7367872 DOI: 10.1038/s41467-020-17275-5] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2020] [Accepted: 06/11/2020] [Indexed: 11/09/2022] Open
Abstract
The impact of the inner structure and thermal history of planets on their observable features, such as luminosity or magnetic field, crucially depends on the poorly known heat and charge transport properties of their internal layers. The thermal and electric conductivities of different phases of water (liquid, solid, and super-ionic) occurring in the interior of ice giant planets, such as Uranus or Neptune, are evaluated from equilibrium ab initio molecular dynamics, leveraging recent progresses in the theory and data analysis of transport in extended systems. The implications of our findings on the evolution models of the ice giants are briefly discussed.
Collapse
|
17
|
Zhang C, Giberti F, Sevgen E, de Pablo JJ, Gygi F, Galli G. Dissociation of salts in water under pressure. Nat Commun 2020; 11:3037. [PMID: 32546791 PMCID: PMC7298052 DOI: 10.1038/s41467-020-16704-9] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2019] [Accepted: 05/15/2020] [Indexed: 11/09/2022] Open
Abstract
The investigation of salts in water at extreme conditions is crucial to understanding the properties of aqueous fluids in the Earth. We report first principles (FP) and classical molecular dynamics simulations of NaCl in the dilute limit, at temperatures and pressures relevant to the Earth’s upper mantle. Similar to ambient conditions, we observe two metastable states of the salt: the contact (CIP) and the solvent-shared ion-pair (SIP), which are entropically and enthalpically favored, respectively. We find that the free energy barrier between the CIP and SIP minima increases at extreme conditions, and that the stability of the CIP is enhanced in FP simulations, consistent with the decrease of the dielectric constant of water. The minimum free energy path between the CIP and SIP becomes smoother at high pressure, and the relative stability of the two configurations is affected by water self-dissociation, which can only be described properly by FP simulations. Salts in water at extreme conditions play a fundamental role in determining the properties of the Earthʼs mantle constituents. Here the authors shed light on ion-water and ion-ion interactions for NaCl dissolved in water at conditions relevant to the Earthʼs upper mantle by molecular dynamics simulations.
Collapse
Affiliation(s)
- Cunzhi Zhang
- Department of Materials Science and Engineering, COE, Peking University, 100871, Beijing, China
| | - Federico Giberti
- University of Chicago, 5640 S. Ellis Ave., Chicago, IL, 60637, USA
| | - Emre Sevgen
- University of Chicago, 5640 S. Ellis Ave., Chicago, IL, 60637, USA
| | - Juan J de Pablo
- University of Chicago, 5640 S. Ellis Ave., Chicago, IL, 60637, USA.,Materials Science Division, Argonne National Laboratory, Argonne, IL, 60439, USA
| | - Francois Gygi
- University of California Davis, Davis, CA, 95616, USA
| | - Giulia Galli
- University of Chicago, 5640 S. Ellis Ave., Chicago, IL, 60637, USA. .,Materials Science Division, Argonne National Laboratory, Argonne, IL, 60439, USA.
| |
Collapse
|
18
|
Dettori R, Donadio D. Carbon dioxide, bicarbonate and carbonate ions in aqueous solutions under deep Earth conditions. Phys Chem Chem Phys 2020; 22:10717-10725. [PMID: 32103223 DOI: 10.1039/c9cp06904f] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
We investigate the effect of pressure, temperature and acidity on the composition of water-rich carbon-bearing fluids under thermodynamic conditions that correspond to the Earth's deep crust and upper mantle. Our first-principles molecular dynamics simulations provide mechanistic insight into the hydration shell of carbon dioxide, bicarbonate and carbonate ions, and into the pathways of the acid/base reactions that convert these carbon species into one another in aqueous solutions. At temperatures of 1000 K and higher, our simulations can sample the chemical equilibrium of these acid/base reactions, thus allowing us to estimate the chemical composition of diluted carbon dioxide and (bi)carbonate ions as a function of acidity and thermodynamic conditions. We find that, especially at the highest temperature, the acidity of the solution is essential to determine the stability domain of CO2vs. HCO3-vs. CO32-.
Collapse
Affiliation(s)
- Riccardo Dettori
- Department of Chemistry, University of California Davis, One Shields Avenue, Davis, California 95616, USA.
| | | |
Collapse
|
19
|
Shao Y, Hellström M, Yllö A, Mindemark J, Hermansson K, Behler J, Zhang C. Temperature effects on the ionic conductivity in concentrated alkaline electrolyte solutions. Phys Chem Chem Phys 2020; 22:10426-10430. [PMID: 31895378 DOI: 10.1039/c9cp06479f] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Alkaline electrolyte solutions are important components in rechargeable batteries and alkaline fuel cells. As the ionic conductivity is thought to be a limiting factor in the performance of these devices, which are often operated at elevated temperatures, its temperature dependence is of significant interest. Here we use NaOH as a prototypical example of alkaline electrolytes, and for this system we have carried out reactive molecular dynamics simulations with an experimentally verified high-dimensional neural network potential derived from density-functional theory calculations. It is found that in concentrated NaOH solutions elevated temperatures enhance both the contributions of proton transfer to the ionic conductivity and deviations from the Nernst-Einstein relation. These findings are expected to be of practical relevance for electrochemical devices based on alkaline electrolyte solutions.
Collapse
Affiliation(s)
- Yunqi Shao
- Department of Chemistry -Ångström Laboratory, Uppsala University, Box 538, 751 21 Uppsala, Sweden.
| | | | | | | | | | | | | |
Collapse
|
20
|
Cassone G, Sponer J, Trusso S, Saija F. Ab initio spectroscopy of water under electric fields. Phys Chem Chem Phys 2019; 21:21205-21212. [PMID: 31368466 DOI: 10.1039/c9cp03101d] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Whereas a broad range of literature exists on the spectroscopy of water in disparate conditions, infrared (IR) and Raman spectra of water subjected to electric fields have never extensively been investigated so far. Based on ab initio molecular dynamics simulations, here we present IR and Raman spectra of bulk liquid water under the effect of static electric fields. A contraction of the entire frequency range is recorded upon increasing the field intensity both in the IR and in the Raman spectra. Whilst the OH stretching band is progressively shifted toward lower frequencies - indicating a field-induced strengthening of the H-bond network - all the other bands are up-shifted by the field. Furthermore, an evident modification of the librational mode band appears in all the spectra. Finally, the order-maker action of the field emerges also from the increase of the water orientational tetrahedral order. Upon field exposure, the water structure becomes more "ice like".
Collapse
Affiliation(s)
- Giuseppe Cassone
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 61265 Brno, Czech Republic.
| | | | | | | |
Collapse
|
21
|
Stolte N, Pan D. Large Presence of Carbonic Acid in CO 2-Rich Aqueous Fluids under Earth's Mantle Conditions. J Phys Chem Lett 2019; 10:5135-5141. [PMID: 31411889 DOI: 10.1021/acs.jpclett.9b01919] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The chemistry of carbon in aqueous fluids at extreme pressure and temperature conditions is of great importance to Earth's deep carbon cycle, which substantially affects the carbon budget at Earth's surface and global climate change. At ambient conditions, the concentration of carbonic acid in water is negligible; therefore, aqueous carbonic acid was simply ignored in previous geochemical models. However, by applying extensive ab initio molecular dynamics simulations at pressure and temperature conditions similar to those in Earth's upper mantle, we found that carbonic acid can be the most abundant carbon species in aqueous CO2 solutions at ∼10 GPa and 1000 K. The mole percent of carbonic acid in total dissolved carbon species increases with increasing pressure along an isotherm, while its mole percent decreases with increasing temperature along an isobar. In CO2-rich solutions, we found significant proton transfer between carbonic acid molecules and bicarbonate ions, which may enhance the conductivity of the solutions. The effects of pH buffering by carbonic acid may play an important role in water-rock interactions in Earth's interior. Our findings suggest that carbonic acid is an important carbon carrier in the deep carbon cycle.
Collapse
Affiliation(s)
- Nore Stolte
- Department of Physics , Hong Kong University of Science and Technology , Hong Kong , China
| | - Ding Pan
- Department of Physics , Hong Kong University of Science and Technology , Hong Kong , China
- Department of Chemistry , Hong Kong University of Science and Technology , Hong Kong , China
- HKUST Fok Ying Tung Research Institute , Guangzhou 511458 , China
| |
Collapse
|
22
|
Hestand NJ, Strong SE, Shi L, Skinner JL. Mid-IR spectroscopy of supercritical water: From dilute gas to dense fluid. J Chem Phys 2019; 150:054505. [DOI: 10.1063/1.5079232] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nicholas J. Hestand
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| | - Steven E. Strong
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| | - Liang Shi
- School of Natural Sciences, University of California, Merced, California 95344, USA
| | - J. L. Skinner
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|