1
|
Creazzo F, Luber S. Water-air interface revisited by means of path-integral ab initio molecular dynamics. Phys Chem Chem Phys 2024. [PMID: 39078670 DOI: 10.1039/d4cp02500h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/31/2024]
Abstract
Although nuclear quantum effects (NQEs) have been considered on bulk liquid water, the impact of these latter on the air-water interface has not yet been reported. Herein, by performing and comparing ab initio molecular dynamics (AIMD) and path integral AIMD (PI-AIMD) simulations, we reveal the impact of NQEs on structural, dynamical and electronic properties as well as IR spectra of the air-water interface at room temperature. NQEs, being able to describe a more accurate proton delocalization in H-bonded system than AIMD, reveal a different structural arrangement and dynamical behaviour of both bulk and interfacial water molecules in comparison to AIMD results. A more de-structured and de-bound water arrangement and coordination are identified when the quantum nature of nuclei are considered for both bulk and interfacial water molecules. Structural properties, such as inter-/intra-molecular bond lengths, coordination numbers and H-bonding angles of bulk and interfacial water molecules here calculated, are affected by NQEs mitigating the overstructured description given by AIMD. Further evidences of an AIMD overstructured description of bulk water are in the computed IR spectra, where an increased absorption peak intensity and an increased strength of the hydrogen-bond network are alleviated by NQEs. In addition, NQEs show a valuable impact on the electronic structure of the air-water interface, reducing the total valence bandwidth and the electronic energy band-gap when passing from bulk to interfacial water. This work proves how NQEs significantly affect properties and features of the air-water interface, that are essential to accurately describe H-bonded interfacial systems.
Collapse
Affiliation(s)
- Fabrizio Creazzo
- Department of Chemistry, University of Zurich, Zurich, Switzerland.
| | - Sandra Luber
- Department of Chemistry, University of Zurich, Zurich, Switzerland.
| |
Collapse
|
2
|
Warias JE, Petersdorf L, Hövelmann SC, Giri RP, Lemke C, Festersen S, Greve M, Mandin P, LeBideau D, Bertram F, Magnussen OM, Murphy BM. The laser pump X-ray probe system at LISA P08 PETRA III. JOURNAL OF SYNCHROTRON RADIATION 2024; 31:779-790. [PMID: 38843001 PMCID: PMC11226150 DOI: 10.1107/s1600577524003400] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Accepted: 04/17/2024] [Indexed: 07/06/2024]
Abstract
Understanding and controlling the structure and function of liquid interfaces is a constant challenge in biology, nanoscience and nanotechnology, with applications ranging from molecular electronics to controlled drug release. X-ray reflectivity and grazing incidence diffraction provide invaluable probes for studying the atomic scale structure at liquid-air interfaces. The new time-resolved laser system at the LISA liquid diffractometer situated at beamline P08 at the PETRA III synchrotron radiation source in Hamburg provides a laser pump with X-ray probe. The femtosecond laser combined with the LISA diffractometer allows unique opportunities to investigate photo-induced structural changes at liquid interfaces on the pico- and nanosecond time scales with pump-probe techniques. A time resolution of 38 ps has been achieved and verified with Bi. First experiments include laser-induced effects on salt solutions and liquid mercury surfaces with static and varied time scales measurements showing the proof of concept for investigations at liquid surfaces.
Collapse
Affiliation(s)
- Jonas Erik Warias
- Institute of Experimental and Applied PhysicsKiel UniversityLeibnizstrasse 1924118KielGermany
| | - Lukas Petersdorf
- Institute of Experimental and Applied PhysicsKiel UniversityLeibnizstrasse 1924118KielGermany
- Ruprecht-Haensel Laboratory, Olshausenstrasse 40, 24098Kiel, Germany
| | - Svenja Carolin Hövelmann
- Institute of Experimental and Applied PhysicsKiel UniversityLeibnizstrasse 1924118KielGermany
- Ruprecht-Haensel Laboratory, Olshausenstrasse 40, 24098Kiel, Germany
- Deutsches Elektronen-Synchrotron DESYNotkestrasse 8522607HamburgGermany
| | - Rajendra Prasad Giri
- Institute of Experimental and Applied PhysicsKiel UniversityLeibnizstrasse 1924118KielGermany
| | - Christoph Lemke
- Institute of Experimental and Applied PhysicsKiel UniversityLeibnizstrasse 1924118KielGermany
| | - Sven Festersen
- Institute of Experimental and Applied PhysicsKiel UniversityLeibnizstrasse 1924118KielGermany
| | - Matthias Greve
- Institute of Experimental and Applied PhysicsKiel UniversityLeibnizstrasse 1924118KielGermany
| | | | | | - Florian Bertram
- Deutsches Elektronen-Synchrotron DESYNotkestrasse 8522607HamburgGermany
| | - Olaf Magnus Magnussen
- Institute of Experimental and Applied PhysicsKiel UniversityLeibnizstrasse 1924118KielGermany
- Ruprecht-Haensel Laboratory, Olshausenstrasse 40, 24098Kiel, Germany
| | - Bridget Mary Murphy
- Institute of Experimental and Applied PhysicsKiel UniversityLeibnizstrasse 1924118KielGermany
- Ruprecht-Haensel Laboratory, Olshausenstrasse 40, 24098Kiel, Germany
| |
Collapse
|
3
|
Ojha D, Henao A, Zysk F, Kühne TD. Nuclear quantum effects on the vibrational dynamics of the water-air interface. J Chem Phys 2024; 160:204114. [PMID: 38804494 DOI: 10.1063/5.0204071] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Accepted: 04/11/2024] [Indexed: 05/29/2024] Open
Abstract
We have applied path-integral molecular dynamics simulations to investigate the impact of nuclear quantum effects on the vibrational dynamics of water molecules at the water-air interface. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time series analysis, while the time scales of vibrational spectral diffusion are determined from frequency-time correlation functions and joint probability distributions. We find that the inclusion of nuclear quantum effects leads not only to a redshift in the vibrational frequency distribution by about 120 cm-1 for both the bulk and interfacial water molecules but also to an acceleration of the vibrational dynamics at the water-air interface by as much as 35%. In addition, a blueshift of about 45 cm-1 is seen in the vibrational frequency distribution of interfacial water molecules compared to that of the bulk. Furthermore, the dynamics of water molecules beyond the topmost molecular layer was found to be rather similar to that of bulk water.
Collapse
Affiliation(s)
- Deepak Ojha
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Andrés Henao
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Frederik Zysk
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Thomas D Kühne
- Center for Advanced Systems Understanding (CASUS), Untermarkt 20, D-02826 Görlitz, Germany, Helmholtz Zentrum Dresden-Rossendorf, Bautzner Landstraße 400, D-01328 Dresden, Germany, and TU Dresden, Institute of Artificial Intelligence, Chair of Computational System Sciences, Nöthnitzer Straße 46, D-01187 Dresden, Germany
| |
Collapse
|
4
|
Ming X, Si W, Yu Q, Sun Z, Qiu G, Cao M, Li Y, Li Z. Molecular insight into the initial hydration of tricalcium aluminate. Nat Commun 2024; 15:2929. [PMID: 38575602 PMCID: PMC10995194 DOI: 10.1038/s41467-024-47164-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Accepted: 03/21/2024] [Indexed: 04/06/2024] Open
Abstract
Portland cement (PC) is ubiquitously used in construction for centuries, yet the elucidation of its early-age hydration remains a challenge. Understanding the initial hydration progress of tricalcium aluminate (C3A) at molecular scale is thus crucial for tackling this challenge as it exhibits a proclivity for early-stage hydration and plays a pivotal role in structural build-up of cement colloids. Herein, we implement a series of ab-initio calculations to probe the intricate molecular interactions of C3A during its initial hydration process. The C3A surface exhibits remarkable chemical activity in promoting water dissociation, which in turn facilitates the gradual desorption of Ca ions through a metal-proton exchange reaction. The dissolution pathways and free energies of these Ca ions follow the ligand-exchange mechanism with multiple sequential reactions to form the ultimate products where Ca ions adopt fivefold or sixfold coordination. Finally, these Ca complexes reprecipitate on the remaining Al-rich layer through the interface-coupled dissolution-reprecipitation mechanism, demonstrating dynamically stable inner-sphere adsorption states. The above results are helpful in unmasking the early-age hydration of PC and advancing the rational design of cement-based materials through the bottom-up approach.
Collapse
Affiliation(s)
- Xing Ming
- Faculty of Innovation Engineering, Macau University of Science and Technology, Avenida Wai Long, Taipa, Macao SAR, China
| | - Wen Si
- School of Civil Engineering, Dalian University of Technology, Dalian, China
| | - Qinglu Yu
- Institute of Applied Physics and Materials Engineering, University of Macau, Avenida da Universidade, Taipa, Macao SAR, China
| | - Zhaoyang Sun
- Institute of Applied Physics and Materials Engineering, University of Macau, Avenida da Universidade, Taipa, Macao SAR, China
| | - Guotao Qiu
- Institute of Applied Physics and Materials Engineering, University of Macau, Avenida da Universidade, Taipa, Macao SAR, China
| | - Mingli Cao
- School of Civil Engineering, Dalian University of Technology, Dalian, China
| | - Yunjian Li
- Faculty of Innovation Engineering, Macau University of Science and Technology, Avenida Wai Long, Taipa, Macao SAR, China.
| | - Zongjin Li
- Faculty of Innovation Engineering, Macau University of Science and Technology, Avenida Wai Long, Taipa, Macao SAR, China.
| |
Collapse
|
5
|
Devlin SW, Bernal F, Riffe EJ, Wilson KR, Saykally RJ. Spiers Memorial Lecture: Water at interfaces. Faraday Discuss 2024; 249:9-37. [PMID: 37795954 DOI: 10.1039/d3fd00147d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/06/2023]
Abstract
In this article we discuss current issues in the context of the four chosen subtopics for the meeting: dynamics and nano-rheology of interfacial water, electrified/charged aqueous interfaces, ice interfaces, and soft matter/water interfaces. We emphasize current advances in both theory and experiment, as well as important practical manifestations and areas of unresolved controversy.
Collapse
Affiliation(s)
- Shane W Devlin
- Department of Chemistry, University of California, Berkeley, CA 94720, USA.
- Chemical Sciences Division, Lawrence Berkeley National Lab, Berkeley, CA 94720, USA
| | - Franky Bernal
- Department of Chemistry, University of California, Berkeley, CA 94720, USA.
- Chemical Sciences Division, Lawrence Berkeley National Lab, Berkeley, CA 94720, USA
| | - Erika J Riffe
- Department of Chemistry, University of California, Berkeley, CA 94720, USA.
- Chemical Sciences Division, Lawrence Berkeley National Lab, Berkeley, CA 94720, USA
| | - Kevin R Wilson
- Chemical Sciences Division, Lawrence Berkeley National Lab, Berkeley, CA 94720, USA
| | - Richard J Saykally
- Department of Chemistry, University of California, Berkeley, CA 94720, USA.
- Chemical Sciences Division, Lawrence Berkeley National Lab, Berkeley, CA 94720, USA
| |
Collapse
|
6
|
Becker M, Loche P, Rezaei M, Wolde-Kidan A, Uematsu Y, Netz RR, Bonthuis DJ. Multiscale Modeling of Aqueous Electric Double Layers. Chem Rev 2024; 124:1-26. [PMID: 38118062 PMCID: PMC10785765 DOI: 10.1021/acs.chemrev.3c00307] [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/07/2023] [Revised: 11/17/2023] [Accepted: 11/30/2023] [Indexed: 12/22/2023]
Abstract
From the stability of colloidal suspensions to the charging of electrodes, electric double layers play a pivotal role in aqueous systems. The interactions between interfaces, water molecules, ions and other solutes making up the electrical double layer span length scales from Ångströms to micrometers and are notoriously complex. Therefore, explaining experimental observations in terms of the double layer's molecular structure has been a long-standing challenge in physical chemistry, yet recent advances in simulations techniques and computational power have led to tremendous progress. In particular, the past decades have seen the development of a multiscale theoretical framework based on the combination of quantum density functional theory, force-field based simulations and continuum theory. In this Review, we discuss these theoretical developments and make quantitative comparisons to experimental results from, among other techniques, sum-frequency generation, atomic-force microscopy, and electrokinetics. Starting from the vapor/water interface, we treat a range of qualitatively different types of surfaces, varying from soft to solid, from hydrophilic to hydrophobic, and from charged to uncharged.
Collapse
Affiliation(s)
| | - Philip Loche
- Fachbereich
Physik, Freie Universität Berlin, 14195 Berlin, Germany
- Laboratory
of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Majid Rezaei
- Fachbereich
Physik, Freie Universität Berlin, 14195 Berlin, Germany
- Institute
of Theoretical Chemistry, Ulm University, 89081 Ulm, Germany
| | | | - Yuki Uematsu
- Department
of Physics and Information Technology, Kyushu
Institute of Technology, 820-8502 Iizuka, Japan
- PRESTO,
Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Roland R. Netz
- Fachbereich
Physik, Freie Universität Berlin, 14195 Berlin, Germany
| | - Douwe Jan Bonthuis
- Institute
of Theoretical and Computational Physics, Graz University of Technology, 8010 Graz, Austria
| |
Collapse
|
7
|
Ricard TC, Zhu X, Iyengar SS. Capturing Weak Interactions in Surface Adsorbate Systems at Coupled Cluster Accuracy: A Graph-Theoretic Molecular Fragmentation Approach Improved through Machine Learning. J Chem Theory Comput 2023. [PMID: 38019639 DOI: 10.1021/acs.jctc.3c00955] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2023]
Abstract
The accurate and efficient study of the interactions of organic matter with the surface of water is critical to a wide range of applications. For example, environmental studies have found that acidic polyfluorinated alkyl substances, especially perfluorooctanoic acid (PFOA), have spread throughout the environment and bioaccumulate into human populations residing near contaminated watersheds, leading to many systemic maladies. Thus, the study of the interactions of PFOA with water surfaces became important for the mitigation of their activity as pollutants and threats to public health. However, theoretical study of the interactions of such organic adsorbates on the surface of water, and their bulk concerted properties, often necessitates the use of ab initio methods to properly incorporate the long-range electronic properties that govern these extended systems. Notable theoretical treatments of "on-water" reactions thus far have employed hybrid DFT and semilocal DFT, but the interactions involved are weak interactions that may be best described using post-Hartree-Fock theory. Here, we aim to demonstrate the utility of a graph-theoretic approach to molecular fragmentation that accurately captures the critical "weak" interactions while maintaining an efficient ab initio treatment of the long-range periodic interactions that underpin the physics of extended systems. We apply this graph-theoretical treatment to study PFOA on the surface of water as a model system for the study of weak interactions seen in the wide range of surface interactions and reactions. The approach divides a system into a set of vertices, that are then connected through edges, faces, and higher order graph theoretic objects known as simplexes, to represent a collection of locally interacting subsystems. These subsystems are then used to construct ab initio molecular dynamics simulations and for computing multidimensional potential energy surfaces. To further improve the computational efficiency of our graph theoretic fragmentation method, we use a recently developed transfer learning protocol to construct the full system potential energy from a family of neural networks each designed to accurately model the behavior of individual simplexes. We use a unique multidimensional clustering algorithm, based on the k-means clustering methodology, to define our training space for each separate simplex. These models are used to extrapolate the energies for molecular dynamics trajectories at PFOA water interfaces, at less than one-tenth the cost as compared to a regular molecular fragmentation-based dynamics calculation with excellent agreement with couple cluster level of full system potential energies.
Collapse
Affiliation(s)
- Timothy C Ricard
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Xiao Zhu
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
8
|
Shen H, Shen X, Wu Z. Simulating the isotropic Raman spectra of O-H stretching mode in liquid H 2O based on a machine learning potential: the influence of vibrational couplings. Phys Chem Chem Phys 2023; 25:28180-28188. [PMID: 37819214 DOI: 10.1039/d3cp03035k] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/13/2023]
Abstract
In this study, we trained a deep potential (DP) for H2O, an accurate machine learning (ML) potential. We performed molecular dynamics (MD) simulations of liquid water using the DP model (or DeePMD simulations). Our results showed that the DP model exhibits DFT-level accuracy, and the DeePMD simulation is a promising approach for modeling the structural properties of liquid water. Based on the DeePMD simulation trajectories, we calculated the isotropic Raman spectra of the O-H stretching mode using the surface-specific velocity-velocity correlation function (ssVVCF), showing that the DeePMD/ssVVCF approach can correctly capture the bimodal characteristics of the experimental Raman spectra, with one peak located near 3400 cm-1 and the other near 3250 cm-1. The success of the DeePMD/ssVVCF approach should be credited to (1) the DFT-level accuracy of the DP model for H2O, (2) the ssVVCF formulation considering the coupling between vibrational modes, and (3) non-Condon effects. Furthermore, the DeePMD simulations revealed that the anharmonic interactions between the coupled water molecules in the first and second hydration shells should play an essential role in the strong mixing of the H-O-H bending mode and the O-H stretching mode, leading to the delocalization of the O-H stretching band. In particular, increasing the strength of hydrogen bonds would enhance the bend-stretch coupling, leading to the red-shifting of the O-H vibrational spectra and the increase in the intensity of the shoulder around 3250 cm-1.
Collapse
Affiliation(s)
- Hujun Shen
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Education University, Guiyang 550018, China.
| | - Xu Shen
- National Center of Technology Innovation for Intelligent Design and Numerical Control, Huazhong University of Science and Technology, Wuhan, 430074, China
| | - Zhenhua Wu
- Department of Big Data and Artificial Intelligence, Guizhou Vocational Technology College of Electronics & Information, Kaili, 556000, China
| |
Collapse
|
9
|
Sanchez-Burgos I, Muniz MC, Espinosa JR, Panagiotopoulos AZ. A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water. J Chem Phys 2023; 158:2889532. [PMID: 37158636 DOI: 10.1063/5.0144500] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Accepted: 03/30/2023] [Indexed: 05/10/2023] Open
Abstract
Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology-a machine learning approach-to study this ubiquitous phase transition, starting from the phase diagram in the liquid-vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. We find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid-vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid-vapor coexistence and water cavitation.
Collapse
Affiliation(s)
- Ignacio Sanchez-Burgos
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue,Cambridge CB3 0HE, United Kingdom
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Maria Carolina Muniz
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Jorge R Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue,Cambridge CB3 0HE, United Kingdom
- Departamento de Química Fisica, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | | |
Collapse
|
10
|
Donkor ED, Laio A, Hassanali A. Do Machine-Learning Atomic Descriptors and Order Parameters Tell the Same Story? The Case of Liquid Water. J Chem Theory Comput 2023. [PMID: 36920997 DOI: 10.1021/acs.jctc.2c01205] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
Machine-learning (ML) has become a key workhorse in molecular simulations. Building an ML model in this context involves encoding the information on chemical environments using local atomic descriptors. In this work, we focus on the Smooth Overlap of Atomic Positions (SOAP) and their application in studying the properties of liquid water both in the bulk and at the hydrophobic air-water interface. By using a statistical test aimed at assessing the relative information content of different distance measures defined on the same data space, we investigate if these descriptors provide the same information as some of the common order parameters that are used to characterize local water structure such as hydrogen bonding, density, or tetrahedrality to name a few. Our analysis suggests that the ML description and the standard order parameters of the local water structure are not equivalent. In particular, a combination of these order parameters probing local water environments can predict SOAP similarity only approximately, and vice versa, the environments that are similar according to SOAP are not necessarily similar according to the standard order parameters. We also elucidate the role of some of the metaparameters in the SOAP definition in encoding chemical information.
Collapse
Affiliation(s)
- Edward Danquah Donkor
- The Abdus Salam International Center for Theoretical Physics (ICTP), Strada Costiera 11, 34151 Trieste, Italy.,Scuola Internazionale Superiore di Studi Avanzati (SISSA), via Bonomea 265, 34136 Trieste, Italy
| | - Alessandro Laio
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), via Bonomea 265, 34136 Trieste, Italy
| | - Ali Hassanali
- The Abdus Salam International Center for Theoretical Physics (ICTP), Strada Costiera 11, 34151 Trieste, Italy
| |
Collapse
|
11
|
Han Y, Zhang T, Guo X, Jiao T. Insights into the mechanism of electrostatic field promoting ozone mass transfer in water: A molecular dynamics perspective. THE SCIENCE OF THE TOTAL ENVIRONMENT 2022; 848:157710. [PMID: 35908697 DOI: 10.1016/j.scitotenv.2022.157710] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/16/2022] [Revised: 07/17/2022] [Accepted: 07/26/2022] [Indexed: 06/15/2023]
Abstract
Ozone is the main role of ozone-based advanced oxidation process for organic wastewater treatment, which is usually added in water by aeration. However, the low solubility of ozone in water seriously affects the degradation efficiency. In this article, the external electrostatic field (EF) was proposed to improve the ozone solubility in water. The mass transfer characteristics of ozone in a gas-liquid two-phase system under EF were studied by molecular dynamics (MD) simulation. The microscopic mechanism of ozone mass transfer in water promoted by external EF was revealed by analyzing Gibbs dividing surface (GDS) interface structure, interfacial water molecular orientation, surface tension, liquid phase viscosity, hydrogen bond network and ozone self-diffusion coefficient. Our findings reveal that EF can enhance the thickness of GDS region (from 0.4648 nm to 0.4996 nm when EF is 0.2 V/nm) as well as its ozone content. The dipole moment orientation of water molecules also tends to point in the EF direction due to the influence of EF, making the difference in dipole moment orientation of water molecules in the first and second layers of GDS region gradually disappear. In addition, compared with the absence of EF, the existence of external EF can decrease the surface tension (from 77.6162 mN/m to 73.3480 mN/m when EF is 0.2 V/nm) at the gas-liquid interface and the viscosity of liquid phase (from 0.293 mPa·s to 0.162 mPa·s when EF is 0.2 V/nm), break the network of hydrogen bond in liquid phase, and increase the mobility of ozone (self-diffusion coefficient of ozone changes from 1.3232 × 10-6 cm2·s-1 to 1.8812 × 10-6 cm2·s-1 when EF is 0.2 V/nm). All these properties changes indicate that the presence of external EF enhances the ability of ozone to penetrate the interface of the two-phase system, and then improves the ozone mass transfer efficiency in water.
Collapse
Affiliation(s)
- Yong Han
- State Key Laboratory of Metastable Materials Science and Technology, Yanshan University, Qinhuangdao, Hebei 066004, PR China; School of Electrical Engineering, Yanshan University, Qinhuangdao, Hebei 066004, PR China.
| | - Teng Zhang
- School of Electrical Engineering, Yanshan University, Qinhuangdao, Hebei 066004, PR China
| | - Xiaoqiang Guo
- School of Electrical Engineering, Yanshan University, Qinhuangdao, Hebei 066004, PR China
| | - Tifeng Jiao
- State Key Laboratory of Metastable Materials Science and Technology, Yanshan University, Qinhuangdao, Hebei 066004, PR China.
| |
Collapse
|
12
|
Odendahl NL, Geissler PL. Local Ice-like Structure at the Liquid Water Surface. J Am Chem Soc 2022; 144:11178-11188. [PMID: 35696525 DOI: 10.1021/jacs.2c01827] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Experiments and computer simulations have established that liquid water's surfaces can deviate in important ways from familiar bulk behavior. Even in the simplest case of an air-water interface, distinctive layering, orientational biases, and hydrogen bond arrangements have been reported, but an overarching picture of their origins and relationships has been incomplete. Here we show that a broad set of such observations can be understood through an analogy with the basal face of crystalline ice. Using simulations, we demonstrate a number of structural similarities between water and ice surfaces, suggesting the presence of domains at the air-water interface with ice-like features that persist over 2-3 molecular diameters. Most prominent is a shared characteristic layering of molecular density and orientation perpendicular to the interface. Lateral correlations of hydrogen bond network geometry point to structural similarities in the parallel direction as well. Our results bolster and significantly extend previous conceptions of ice-like structure at the liquid's boundary and suggest that the much-discussed quasi-liquid layer on ice evolves subtly above the melting point into a quasi-ice layer at the surface of liquid water.
Collapse
Affiliation(s)
- Nathan L Odendahl
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Phillip L Geissler
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
13
|
Capelli R, Muniz-Miranda F, Pavan GM. Ephemeral ice-like local environments in classical rigid models of liquid water. J Chem Phys 2022; 156:214503. [DOI: 10.1063/5.0088599] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Despite great efforts over the past 50 years, the simulation of water still presents significant challenges and open questions. At room temperature and pressure, the collective molecular interactions and dynamics of water molecules may form local structural arrangements that are non-trivial to classify. Here, we employ a data-driven approach built on Smooth Overlap of Atomic Position (SOAP) that allows us to compare and classify how widely used classical models represent liquid water. Macroscopically, the obtained results are rationalized based on water thermodynamic observables. Microscopically, we directly observe how transient ice-like ordered environments may dynamically/statistically form in liquid water, even above freezing temperature, by comparing the SOAP spectra for different ice structures with those of the simulated liquid systems. This confirms recent ab initio-based calculations but also reveals how the emergence of ephemeral local ice-like environments in liquid water at room conditions can be captured by classical water models.
Collapse
Affiliation(s)
- Riccardo Capelli
- Department of Applied Science and Technology, Politecnico di Torino, Corso Duca degli Abruzzi 24, I-10129 Torino, Italy
| | - Francesco Muniz-Miranda
- Department of Chemical and Geological Sciences, University of Modena and Reggio-Emilia, Via Campi 103, I-41125 Modena, Italy
| | - Giovanni M. Pavan
- Department of Applied Science and Technology, Politecnico di Torino, Corso Duca degli Abruzzi 24, I-10129 Torino, Italy
- Department of Innovative Technologies, University of Applied Sciences and Arts of Southern Switzerland, Polo Universitario Lugano, Campus Est, Via la Santa 1, CH-6962 Lugano-Viganello, Switzerland
| |
Collapse
|
14
|
Thomsen B, Shiga M. Structures of liquid and aqueous water isotopologues at ambient temperature from ab initio path integral simulations. Phys Chem Chem Phys 2022; 24:10851-10859. [PMID: 35504275 DOI: 10.1039/d2cp00499b] [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
The heavy hydrogen isotopes D and T are found in trace amounts in water, and when their concentration increases they can play an intricate role in modulating the physical properties of the liquid. We present an analysis of the microscopic structures of ambient light water (H2O(l)), heavy water (D2O(l)), T2O(l), HDO(aq) and HTO(aq) studied by ab initio path integral molecular dynamics (PIMD). Unlike previous ab initio PIMD investigations of H2O(l) and D2O(l) [Chen et al., Phys. Rev. Lett., 2003, 91, 215503] [Machida et al., J. Chem. Phys., 2017, 148, 102324] we find that D2O(l) is more structured than H2O(l), as is predicted by the experiment. The agreement between the experiment and our simulation for H2O(l) and D2O(l) allows us to accurately predict the intra- and intermolecular structures of T2O(l) HDO(aq) and HTO(aq). T2O(l) is found to have a similar intermolecular structure to that of D2O(l), while the intramolecular structure is more compact, giving rise to a smaller dipole moment than those of H2O(l) and D2O(l). For the mixed isotope species, HDO(aq) and HTO(aq), we find smaller dipole moments and fewer hydrogen bonds when compared with the pure species H2O and D2O. We can attribute this effect to the relative compactness of the mixed isotope species, which results in a lower dipole moment than that of the pure species.
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
|
15
|
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
|
16
|
Goujon F, Ghoufi A, Malfreyt P. Associated molecular liquids at the graphene monolayer interface. J Chem Phys 2021; 154:104504. [PMID: 33722040 DOI: 10.1063/5.0042438] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
We report molecular simulations of the interaction between a graphene sheet and different liquids such as water, ethanol, and ethylene glycol. We describe the structural arrangements at the graphene interface in terms of density profiles, number of hydrogen bonds (HBs), and local structuration in neighboring layers close to the surface. We establish the formation of a two-dimensional HB network in the layer closest to the graphene. We also calculate the interfacial tension of liquids with a graphene monolayer and its profile along the direction normal to the graphene to rationalize and quantify the strengthening of the intermolecular interactions in the liquid due to the presence of the surface.
Collapse
Affiliation(s)
- Florent Goujon
- Université Clermont Auvergne, CNRS, SIGMA Clermont, Institut de Chimie de Clermont-Ferrand (ICCF), F-63000 Clermont-Ferrand, France
| | - Aziz Ghoufi
- Université de Rennes, CNRS, IPR (Institut de Physique de Rennes), UMR 6251, F-35000 Rennes, France
| | - Patrice Malfreyt
- Université Clermont Auvergne, CNRS, SIGMA Clermont, Institut de Chimie de Clermont-Ferrand (ICCF), F-63000 Clermont-Ferrand, France
| |
Collapse
|
17
|
Ojha D, Kühne TD. Hydrogen bond dynamics of interfacial water molecules revealed from two-dimensional vibrational sum-frequency generation spectroscopy. Sci Rep 2021; 11:2456. [PMID: 33510246 PMCID: PMC7844302 DOI: 10.1038/s41598-021-81635-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2020] [Accepted: 01/06/2021] [Indexed: 11/17/2022] Open
Abstract
Vibrational sum-frequency generation (vSFG) spectroscopy allows the study of the structure and dynamics of interfacial systems. In the present work, we provide a simple recipe, based on a narrowband IR pump and broadband vSFG probe technique, to computationally obtain the two-dimensional vSFG spectrum of water molecules at the air-water interface. Using this technique, to study the time-dependent spectral evolution of hydrogen-bonded and free water molecules, we demonstrate that at the interface, the vibrational spectral dynamics of the free OH bond is faster than that of the bonded OH mode.
Collapse
Affiliation(s)
- Deepak Ojha
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, Paderborn University, Warburger Str. 100, 33098, Paderborn, Germany
| | - Thomas D Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, Paderborn University, Warburger Str. 100, 33098, Paderborn, Germany.
- Paderborn Center for Parallel Computing and Institute for Lightweight Design, Paderborn University, Warburger Str. 100, 33098, Paderborn, Germany.
| |
Collapse
|
18
|
Chen M, Li L, Zhu R, Zhu J, He H. Intrinsic water layering next to soft, solid, hydrophobic, and hydrophilic substrates. J Chem Phys 2020; 153:224702. [DOI: 10.1063/5.0030021] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Meng Chen
- CAS Key Laboratory of Mineralogy and Metallogeny/Guangdong Provincial Key Laboratory of Mineral Physics and Materials, Guangzhou Institute of Geochemistry, Institutions of Earth Science, Chinese Academy of Sciences (CAS), Guangzhou 510640, China
| | - Lin Li
- CAS Key Laboratory of Mineralogy and Metallogeny/Guangdong Provincial Key Laboratory of Mineral Physics and Materials, Guangzhou Institute of Geochemistry, Institutions of Earth Science, Chinese Academy of Sciences (CAS), Guangzhou 510640, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Runliang Zhu
- CAS Key Laboratory of Mineralogy and Metallogeny/Guangdong Provincial Key Laboratory of Mineral Physics and Materials, Guangzhou Institute of Geochemistry, Institutions of Earth Science, Chinese Academy of Sciences (CAS), Guangzhou 510640, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Jianxi Zhu
- CAS Key Laboratory of Mineralogy and Metallogeny/Guangdong Provincial Key Laboratory of Mineral Physics and Materials, Guangzhou Institute of Geochemistry, Institutions of Earth Science, Chinese Academy of Sciences (CAS), Guangzhou 510640, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hongping He
- CAS Key Laboratory of Mineralogy and Metallogeny/Guangdong Provincial Key Laboratory of Mineral Physics and Materials, Guangzhou Institute of Geochemistry, Institutions of Earth Science, Chinese Academy of Sciences (CAS), Guangzhou 510640, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
19
|
Ansari N, Karmakar T, Parrinello M. Molecular Mechanism of Gas Solubility in Liquid: Constant Chemical Potential Molecular Dynamics Simulations. J Chem Theory Comput 2020; 16:5279-5286. [PMID: 32551636 DOI: 10.1021/acs.jctc.0c00450] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Accurate prediction of gas solubility in a liquid is crucial in many areas of chemistry, and a detailed understanding of the molecular mechanism of the gas solvation continues to be an active area of research. Here, we extend the idea of the constant chemical potential molecular dynamics (CμMD) approach to the calculation of the gas solubility in the liquid under constant gas chemical potential conditions. As a representative example, we utilize this method to calculate the isothermal solubility of carbon dioxide in water. Additionally, we provide microscopic insight into the mechanism of solvation that preferentially occurs in areas of the surface where the hydrogen network is broken.
Collapse
Affiliation(s)
- Narjes Ansari
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland.,Facoltà di informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, CH-6900 Lugano, Switzerland
| | - Tarak Karmakar
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland.,Facoltà di informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, CH-6900 Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland.,Facoltà di informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, CH-6900 Lugano, Switzerland.,Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy
| |
Collapse
|
20
|
Kühne TD, Iannuzzi M, Del Ben M, Rybkin VV, Seewald P, Stein F, Laino T, Khaliullin RZ, Schütt O, Schiffmann F, Golze D, Wilhelm J, Chulkov S, Bani-Hashemian MH, Weber V, Borštnik U, Taillefumier M, Jakobovits AS, Lazzaro A, Pabst H, Müller T, Schade R, Guidon M, Andermatt S, Holmberg N, Schenter GK, Hehn A, Bussy A, Belleflamme F, Tabacchi G, Glöß A, Lass M, Bethune I, Mundy CJ, Plessl C, Watkins M, VandeVondele J, Krack M, Hutter J. CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations. J Chem Phys 2020; 152:194103. [PMID: 33687235 DOI: 10.1063/5.0007045] [Citation(s) in RCA: 896] [Impact Index Per Article: 224.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.
Collapse
Affiliation(s)
- Thomas D Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Mauro Del Ben
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Vladimir V Rybkin
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Patrick Seewald
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Frederick Stein
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Teodoro Laino
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | - Rustam Z Khaliullin
- Department of Chemistry, McGill University, CH-801 Sherbrooke St. West, Montreal, Quebec H3A 0B8, Canada
| | - Ole Schütt
- Department of Materials, ETH Zürich, CH-8092 Zürich, Switzerland
| | | | - Dorothea Golze
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| | - Jan Wilhelm
- Institute of Theoretical Physics, University of Regensburg, Universitätsstraße 31, D-93053 Regensburg, Germany
| | - Sergey Chulkov
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | | | - Valéry Weber
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | | | | | | | | | - Hans Pabst
- Intel Extreme Computing, Software and Systems, Zürich, Switzerland
| | - Tiziano Müller
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Robert Schade
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Manuel Guidon
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Samuel Andermatt
- Integrated Systems Laboratory, ETH Zürich, CH-8092 Zürich, Switzerland
| | - Nico Holmberg
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076 Aalto, Finland
| | - Gregory K Schenter
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Anna Hehn
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Augustin Bussy
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Fabian Belleflamme
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Gloria Tabacchi
- Department of Science and High Technology, University of Insubria and INSTM, via Valleggio 9, I-22100 Como, Italy
| | - Andreas Glöß
- BASF SE, Carl-Bosch-Straße 38, D-67056 Ludwigshafen am Rhein, Germany
| | - Michael Lass
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Iain Bethune
- Hartree Centre, Science and Technology Facilities Council, Sci-Tech Daresbury, Warrington WA4 4AD, United Kingdom
| | - Christopher J Mundy
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Christian Plessl
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Matt Watkins
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | - Joost VandeVondele
- Swiss National Supercomputing Centre (CSCS), ETH Zürich, Zürich, Switzerland
| | - Matthias Krack
- Laboratory for Scientific Computing and Modelling, Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
21
|
Tang F, Ohto T, Sun S, Rouxel JR, Imoto S, Backus EHG, Mukamel S, Bonn M, Nagata Y. Molecular Structure and Modeling of Water-Air and Ice-Air Interfaces Monitored by Sum-Frequency Generation. Chem Rev 2020; 120:3633-3667. [PMID: 32141737 PMCID: PMC7181271 DOI: 10.1021/acs.chemrev.9b00512] [Citation(s) in RCA: 65] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2019] [Indexed: 12/26/2022]
Abstract
From a glass of water to glaciers in Antarctica, water-air and ice-air interfaces are abundant on Earth. Molecular-level structure and dynamics at these interfaces are key for understanding many chemical/physical/atmospheric processes including the slipperiness of ice surfaces, the surface tension of water, and evaporation/sublimation of water. Sum-frequency generation (SFG) spectroscopy is a powerful tool to probe the molecular-level structure of these interfaces because SFG can specifically probe the topmost interfacial water molecules separately from the bulk and is sensitive to molecular conformation. Nevertheless, experimental SFG has several limitations. For example, SFG cannot provide information on the depth of the interface and how the orientation of the molecules varies with distance from the surface. By combining the SFG spectroscopy with simulation techniques, one can directly compare the experimental data with the simulated SFG spectra, allowing us to unveil the molecular-level structure of water-air and ice-air interfaces. Here, we present an overview of the different simulation protocols available for SFG spectra calculations. We systematically compare the SFG spectra computed with different approaches, revealing the advantages and disadvantages of the different methods. Furthermore, we account for the findings through combined SFG experiments and simulations and provide future challenges for SFG experiments and simulations at different aqueous interfaces.
Collapse
Affiliation(s)
- Fujie Tang
- Max
Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
- Department
of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Tatsuhiko Ohto
- Graduate
School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
| | - Shumei Sun
- Max
Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
- Department
of Physical Chemistry, University of Vienna, Währinger Strasse 42, 1090 Vienna, Austria
| | - Jérémy R. Rouxel
- Department
of Chemistry and Department of Physics and Astronomy, University of California, Irvine, Irvine, California 92697-2025, United States
| | - Sho Imoto
- Max
Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Ellen H. G. Backus
- Max
Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
- Department
of Physical Chemistry, University of Vienna, Währinger Strasse 42, 1090 Vienna, Austria
| | - Shaul Mukamel
- Department
of Chemistry and Department of Physics and Astronomy, University of California, Irvine, Irvine, California 92697-2025, United States
| | - Mischa Bonn
- Max
Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Yuki Nagata
- Max
Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
- Department
of Physics, State Key Laboratory of Surface Physics and Key Laboratory
of Micro- and Nano-Photonic Structures (MOE), Fudan University, Shanghai 200433, China
| |
Collapse
|
22
|
Elgabarty H, Kühne TD. Tumbling with a limp: local asymmetry in water's hydrogen bond network and its consequences. Phys Chem Chem Phys 2020; 22:10397-10411. [DOI: 10.1039/c9cp06960g] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Ab initio molecular dynamics simulations of ambient liquid water and energy decomposition analysis have recently shown that water molecules exhibit significant asymmetry between the strengths of the two donor and/or the two acceptor interactions.
Collapse
Affiliation(s)
- Hossam Elgabarty
- Dynamics of Condensed Matter and Center for Sustainable Systems Design
- Chair of Theoretical Chemistry
- University of Paderborn
- Paderborn
- Germany
| | - Thomas D. Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design
- Chair of Theoretical Chemistry
- University of Paderborn
- Paderborn
- Germany
| |
Collapse
|
23
|
Altan I, Khan AR, James S, Quinn MK, McManus JJ, Charbonneau P. Using Schematic Models to Understand the Microscopic Basis for Inverted Solubility in γD-Crystallin. J Phys Chem B 2019; 123:10061-10072. [PMID: 31557434 DOI: 10.1021/acs.jpcb.9b07774] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Inverted solubility-melting a crystal by cooling-is observed in a handful of proteins, such as carbomonoxy hemoglobin C and γD-crystallin. In human γD-crystallin, the phenomenon is associated with the mutation of the 23rd residue, a proline, to a threonine, serine, or valine. One proposed microscopic mechanism entails an increase in surface hydrophobicity upon mutagenesis. Recent crystal structures of a double mutant that includes the P23T mutation allow for a more careful investigation of this proposal. Here, we first measure the surface hydrophobicity of various mutant structures of γD-crystallin and discern no notable increase in hydrophobicity upon mutating the 23rd residue. We then investigate the solubility inversion regime with a schematic patchy particle model that includes one of three variants of temperature-dependent patch energies: two of the hydrophobic effect, and one of a more generic nature. We conclude that, while solubility inversion due to the hydrophobic effect may be possible, microscopic evidence to support it in γD-crystallin is weak. More generally, we find that solubility inversion requires a fine balance between patch strengths and their temperature-dependent component, which may explain why inverted solubility is not commonly observed in proteins. We also find that the temperature-dependent interaction has only a negligible impact on liquid-liquid phase boundaries of γD-crystallin, in line with previous experimental observations.
Collapse
Affiliation(s)
| | - Amir R Khan
- School of Biochemistry and Immunology , Trinity College Dublin , Dublin , Ireland
| | - Susan James
- Department of Chemistry , Maynooth University , Maynooth , Ireland
| | - Michelle K Quinn
- Department of Chemistry , Maynooth University , Maynooth , Ireland
| | | | | |
Collapse
|
24
|
Das S, Bonn M, Backus EHG. The Surface Activity of the Hydrated Proton Is Substantially Higher than That of the Hydroxide Ion. Angew Chem Int Ed Engl 2019; 58:15636-15639. [PMID: 31418999 PMCID: PMC6856863 DOI: 10.1002/anie.201908420] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2019] [Indexed: 11/30/2022]
Abstract
The behavior of hydroxide and hydrated protons, the auto-ionization products of water, at surfaces is important for a wide range of applications and disciplines. However, it is unknown at which bulk concentration these ions start to become surface active at the water-air interface. Here, we report changes in the D2 O-air interface in the presence of excess D+ hyd /OD- hyd determined using surface-sensitive vibrational sum-frequency generation (SFG) spectroscopy. The onset of the perturbation of the D2 O surface occurs at a bulk concentration as low as 2.7±0.2 mm D+ hyd . In contrast, a concentration of several hundred mm OD- hyd is required to change the D2 O surface. The hydrated proton is thus orders of magnitude more surface-active than hydroxide at the water-air interface.
Collapse
Affiliation(s)
- Sudipta Das
- Department of Molecular SpectroscopyMax Planck Institute for Polymer ResearchAckermannweg 1055128MainzGermany
| | - Mischa Bonn
- Department of Molecular SpectroscopyMax Planck Institute for Polymer ResearchAckermannweg 1055128MainzGermany
| | - Ellen H. G. Backus
- Department of Physical ChemistryUniversity of ViennaWähringer Strasse 421090ViennaAustria
- Department of Molecular SpectroscopyMax Planck Institute for Polymer ResearchAckermannweg 1055128MainzGermany
| |
Collapse
|
25
|
Ojha D, Kaliannan NK, Kühne TD. Time-dependent vibrational sum-frequency generation spectroscopy of the air-water interface. Commun Chem 2019. [DOI: 10.1038/s42004-019-0220-6] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Abstract
Vibrational sum-frequency generation spectroscopy is a powerful method to study the microscopic structure and dynamics of interfacial systems. Here we demonstrate a simple computational approach to calculate the time-dependent, frequency-resolved vibrational sum-frequency generation spectrum (TD-vSFG) of the air-water interface. Using this approach, we show that at the air-water interface, the transition of water molecules with bonded OH modes to free OH modes occurs at a time scale of $$\sim$$
~
3 ps, whereas water molecules with free OH modes rapidly make a transition to a hydrogen-bonded state within $$\sim$$
~
2 ps. Furthermore, we also elucidate the origin of the observed differential dynamics based on the time-dependent evolution of water molecules in the different local solvent environments.
Collapse
|
26
|
Das S, Bonn M, Backus EHG. Das hydratisierte Proton besitzt eine deutlich höhere Oberflächenaktivität als das Hydroxidion. Angew Chem Int Ed Engl 2019. [DOI: 10.1002/ange.201908420] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Sudipta Das
- Department of Molecular Spectroscopy Max Planck Institute for Polymer Research Ackermannweg 10 55128 Mainz Deutschland
| | - Mischa Bonn
- Department of Molecular Spectroscopy Max Planck Institute for Polymer Research Ackermannweg 10 55128 Mainz Deutschland
| | - Ellen H. G. Backus
- Department of Physical Chemistry University of Vienna Währinger Strasse 42 1090 Vienna Österreich
- Department of Molecular Spectroscopy Max Planck Institute for Polymer Research Ackermannweg 10 55128 Mainz Deutschland
| |
Collapse
|
27
|
Robalo JR, Streacker LM, Mendes de Oliveira D, Imhof P, Ben-Amotz D, Verde AV. Hydrophobic but Water-Friendly: Favorable Water–Perfluoromethyl Interactions Promote Hydration Shell Defects. J Am Chem Soc 2019; 141:15856-15868. [DOI: 10.1021/jacs.9b06862] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- João R. Robalo
- Department of Theory & Bio-systems, Max Planck Institute for Colloids and Interfaces, Science Park, Potsdam 14476, Germany
| | - Louis M. Streacker
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | | | - Petra Imhof
- Institute for Theoretical Physics, Free University of Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Dor Ben-Amotz
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Ana Vila Verde
- Department of Theory & Bio-systems, Max Planck Institute for Colloids and Interfaces, Science Park, Potsdam 14476, Germany
| |
Collapse
|
28
|
Ohto T, Dodia M, Xu J, Imoto S, Tang F, Zysk F, Kühne TD, Shigeta Y, Bonn M, Wu X, Nagata Y. Accessing the Accuracy of Density Functional Theory through Structure and Dynamics of the Water-Air Interface. J Phys Chem Lett 2019; 10:4914-4919. [PMID: 31393136 PMCID: PMC6748669 DOI: 10.1021/acs.jpclett.9b01983] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Accepted: 08/08/2019] [Indexed: 05/31/2023]
Abstract
Density functional theory-based molecular dynamics simulations are increasingly being used for simulating aqueous interfaces. Nonetheless, the choice of the appropriate density functional, critically affecting the outcome of the simulation, has remained arbitrary. Here, we assess the performance of various exchange-correlation (XC) functionals, based on the metrics relevant to sum-frequency generation spectroscopy. The structure and dynamics of water at the water-air interface are governed by heterogeneous intermolecular interactions, thereby providing a critical benchmark for XC functionals. We find that the XC functionals constrained by exact functional conditions (revPBE and revPBE0) with the dispersion correction show excellent performance. The poor performance of the empirically optimized density functional (M06-L) indicates the importance of satisfying the exact functional condition. Understanding the performance of different XC functionals can aid in resolving the controversial interpretation of the interfacial water structure and direct the design of novel, improved XC functionals better suited to describing the heterogeneous interactions in condensed phases.
Collapse
Affiliation(s)
- Tatsuhiko Ohto
- Graduate
School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
| | - Mayank Dodia
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Jianhang Xu
- Department
of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Sho Imoto
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Fujie Tang
- Department
of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Frederik Zysk
- Dynamics
of Condensed Matter and Center for Sustainable Systems Design, Chair
of Theoretical Chemistry, University of
Paderborn, Warburger Strasse 100, 33098 Paderborn, Germany
| | - Thomas D. Kühne
- Dynamics
of Condensed Matter and Center for Sustainable Systems Design, Chair
of Theoretical Chemistry, University of
Paderborn, Warburger Strasse 100, 33098 Paderborn, Germany
| | - Yasuteru Shigeta
- Graduate
School of Pure and Applied Sciences, University
of Tsukuba, Tennodai 1-1-1, Tsukuba, Ibaraki 305-8571, Japan
- Center
for Computational Sciences, University of
Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| | - Mischa Bonn
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Xifan Wu
- Department
of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Yuki Nagata
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
29
|
Kaliannan NK, Henao Aristizabal A, Wiebeler H, Zysk F, Ohto T, Nagata Y, Kühne TD. Impact of intermolecular vibrational coupling effects on the sum-frequency generation spectra of the water/air interface. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1620358] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Naveen Kumar Kaliannan
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Paderborn, Germany
| | - Andres Henao Aristizabal
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Paderborn, Germany
| | - Hendrik Wiebeler
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Paderborn, Germany
| | - Frederik Zysk
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Paderborn, Germany
| | - Tatsuhiko Ohto
- Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka, Japan
| | - Yuki Nagata
- Max-Planck Institute for Polymer Research, Mainz, Germany
| | - Thomas D. Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Paderborn, Germany
- Paderborn Center for Parallel Computing and Institute for Lightweight Design, University of Paderborn, Paderborn, Germany
| |
Collapse
|
30
|
Spectroscopic BIL-SFG Invariance Hides the Chaotropic Effect of Protons at the Air-Water Interface. ATMOSPHERE 2018. [DOI: 10.3390/atmos9100396] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
The knowledge of the water structure at the interface with the air in acidic pH conditions is of utmost importance for chemistry in the atmosphere. We shed light on the acidic air-water (AW) interfacial structure by DFT-MD simulations of the interface containing one hydronium ion coupled with theoretical SFG (Sum Frequency Generation) spectroscopy. The interpretation of SFG spectra at charged interfaces requires a deconvolution of the signal into BIL (Binding Interfacial Layer) and DL (Diffuse Layer) SFG contributions, which is achieved here, and hence reveals that even though H 3 O + has a chaotropic effect on the BIL water structure (by weakening the 2D-HBond-Network observed at the neat air-water interface) it has no direct probing in SFG spectroscopy. The changes observed experimentally in the SFG of the acidic AW interface from the SFG at the neat AW are shown here to be solely due to the DL-SFG contribution to the spectroscopy. Such BIL-SFG and DL-SFG deconvolution rationalizes the experimental SFG data in the literature, while the hydronium chaotropic effect on the water 2D-HBond-Network in the BIL can be put in perspective of the decrease in surface tension at acidic AW interfaces.
Collapse
|
31
|
Richters D, Kühne TD. Linear-scaling self-consistent field theory based molecular dynamics: application to C60buckyballs colliding with graphite. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1511899] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Dorothee Richters
- Institute of Physical Chemistry and Center of Computational Sciences, Johannes Gutenberg University Mainz, Mainz, Germany
| | - Thomas D. Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Paderborn, Germany
- Paderborn Center for Parallel Computing and Institute for Lightweight Design with Hybrid Systems, Paderborn, Germany
| |
Collapse
|
32
|
Serva A, Pezzotti S, Bougueroua S, Galimberti DR, Gaigeot MP. Combining ab-initio and classical molecular dynamics simulations to unravel the structure of the 2D-HB-network at the air-water interface. J Mol Struct 2018. [DOI: 10.1016/j.molstruc.2018.03.074] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
33
|
Pezzotti S, Serva A, Gaigeot MP. 2D-HB-Network at the air-water interface: A structural and dynamical characterization by means of ab initio and classical molecular dynamics simulations. J Chem Phys 2018; 148:174701. [DOI: 10.1063/1.5018096] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Simone Pezzotti
- Laboratoire Analyse et Modélisation pour la Biologie et l’Environnement, LAMBE CNRS UMR8587, Université d’Evry val d’Essonne, Blvd. F. Mitterrand, Bat Maupertuis, 91025 Evry, France and Université Paris-Saclay, Orsay, France
| | - Alessandra Serva
- Laboratoire Analyse et Modélisation pour la Biologie et l’Environnement, LAMBE CNRS UMR8587, Université d’Evry val d’Essonne, Blvd. F. Mitterrand, Bat Maupertuis, 91025 Evry, France and Université Paris-Saclay, Orsay, France
| | - Marie-Pierre Gaigeot
- Laboratoire Analyse et Modélisation pour la Biologie et l’Environnement, LAMBE CNRS UMR8587, Université d’Evry val d’Essonne, Blvd. F. Mitterrand, Bat Maupertuis, 91025 Evry, France and Université Paris-Saclay, Orsay, France
| |
Collapse
|
34
|
Machida M, Kato K, Shiga M. Nuclear quantum effects of light and heavy water studied by all-electron first principles path integral simulations. J Chem Phys 2018; 148:102324. [DOI: 10.1063/1.5000091] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Affiliation(s)
- Masahiko Machida
- CCSE, Japan Atomic Energy Agency (JAEA), 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| | - Koichiro Kato
- Mizuho Information and Research Institute, Inc., 2-3, Kandanishiki-cho, Chiyoda-ku, Tokyo 101-8443, Japan
| | - Motoyuki Shiga
- CCSE, Japan Atomic Energy Agency (JAEA), 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| |
Collapse
|
35
|
Affiliation(s)
- Xinzijian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| | - Jian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| |
Collapse
|
36
|
Roy Choudhuri J, Chandra A. Effects of dispersion interactions on the structure, polarity, and dynamics of liquid-vapor interface of an aqueous NaCl solution: Results of first principles simulations at room temperature. J Chem Phys 2018; 148:024702. [DOI: 10.1063/1.5005951] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jyoti Roy Choudhuri
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Amalendu Chandra
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| |
Collapse
|
37
|
Mamatkulov SI, Allolio C, Netz RR, Bonthuis DJ. Adsorption von hydrierten Protonen an der Luft-Wasser-Grenzfläche durch Orientierung. Angew Chem Int Ed Engl 2017. [DOI: 10.1002/ange.201707391] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Shavkat I. Mamatkulov
- Fachbereich Physik; Freie Universität Berlin; 14195 Berlin Deutschland
- Center of Higher Technologies; 100174 Tashkent Uzbekistan
| | - Christoph Allolio
- Fachbereich Physik; Freie Universität Berlin; 14195 Berlin Deutschland
- Institute of Chemistry and the Fritz Haber Research Center; The Hebrew University of Jerusalem; Jerusalem 91904 Israel
| | - Roland R. Netz
- Fachbereich Physik; Freie Universität Berlin; 14195 Berlin Deutschland
| | | |
Collapse
|
38
|
Tang F, Ohto T, Hasegawa T, Xie WJ, Xu L, Bonn M, Nagata Y. Definition of Free O–H Groups of Water at the Air–Water Interface. J Chem Theory Comput 2017; 14:357-364. [DOI: 10.1021/acs.jctc.7b00566] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Fujie Tang
- International
Center for Quantum Materials, School of Physics, Peking University, 5
Yiheyuan Road, Haidian, Beijing 100871, China
- Max Planck Institute
for Polymer Research, Ackermannweg
10, D-55128 Mainz, Germany
| | - Tatsuhiko Ohto
- Graduate
School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
| | - Taisuke Hasegawa
- Max Planck Institute
for Polymer Research, Ackermannweg
10, D-55128 Mainz, Germany
- Department
of Chemistry, Graduate School of Science, Kyoto University, Sakyoku, Kyoto 606-8502, Japan
| | - Wen Jun Xie
- Max Planck Institute
for Polymer Research, Ackermannweg
10, D-55128 Mainz, Germany
- College
of Chemistry and Molecular Engineering, Peking University, 5
Yiheyuan Road, Haidian, Beijing 100871, China
| | - Limei Xu
- International
Center for Quantum Materials, School of Physics, Peking University, 5
Yiheyuan Road, Haidian, Beijing 100871, China
- Collaborative Innovation
Center of Quantum Matter, Beijing 100871, China
| | - Mischa Bonn
- Max Planck Institute
for Polymer Research, Ackermannweg
10, D-55128 Mainz, Germany
| | - Yuki Nagata
- Max Planck Institute
for Polymer Research, Ackermannweg
10, D-55128 Mainz, Germany
- Institute for
Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
39
|
Mamatkulov SI, Allolio C, Netz RR, Bonthuis DJ. Orientation-Induced Adsorption of Hydrated Protons at the Air-Water Interface. Angew Chem Int Ed Engl 2017; 56:15846-15851. [DOI: 10.1002/anie.201707391] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2017] [Indexed: 11/06/2022]
Affiliation(s)
- Shavkat I. Mamatkulov
- Fachbereich Physik; Freie Universität Berlin; 14195 Berlin Germany
- Center of Higher Technologies; 100174 Tashkent Uzbekistan
| | - Christoph Allolio
- Fachbereich Physik; Freie Universität Berlin; 14195 Berlin Germany
- Institute of Chemistry and the Fritz Haber Research Center; The Hebrew University of Jerusalem; Jerusalem 91904 Israel
| | - Roland R. Netz
- Fachbereich Physik; Freie Universität Berlin; 14195 Berlin Germany
| | | |
Collapse
|
40
|
Kayal A, Chandra A. Orientational order and dynamics of interfacial water near a hexagonal boron-nitride sheet: An ab initio molecular dynamics study. J Chem Phys 2017; 147:164704. [DOI: 10.1063/1.4991594] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Abhijit Kayal
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Amalendu Chandra
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| |
Collapse
|
41
|
Pezzotti S, Galimberti DR, Gaigeot MP. 2D H-Bond Network as the Topmost Skin to the Air-Water Interface. J Phys Chem Lett 2017; 8:3133-3141. [PMID: 28644626 DOI: 10.1021/acs.jpclett.7b01257] [Citation(s) in RCA: 78] [Impact Index Per Article: 11.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
We provide a detailed description of the structure of water at the interface with the air (liquid-vapor LV interface) from state-of-the-art DFT-based molecular dynamics simulations. For the first time, a two-dimensional (2D) H-bond extended network has been identified and fully characterized, demonstrating that interfacial water is organized into a 2D sheet with H-bonds oriented parallel to the instantaneous surface and following its spatial and temporal oscillations. By analyzing the nonlinear vSFG (vibrational sum frequency generation) spectrum of the LV interface in terms of layer-by-layer signal, we demonstrate that the 2D water sheet is solely responsible for the spectral signatures, hence providing the interfacial 3.5 Å thickness effectively probed in nonlinear interfacial spectroscopy. The 2D H-bond network unraveled here is the essential key to rationalize macroscopic properties of water-air interfaces, as demonstrated here for spectroscopy and the surface potential.
Collapse
Affiliation(s)
- Simone Pezzotti
- LAMBE CNRS UMR8587, Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, Université d'Evry val d'Essonne , Boulevard F. Mitterrand, Bat Maupertuis, 91025 Evry, France
- Université Paris-Saclay , 91190 Saint-Aubin, France
| | - Daria Ruth Galimberti
- LAMBE CNRS UMR8587, Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, Université d'Evry val d'Essonne , Boulevard F. Mitterrand, Bat Maupertuis, 91025 Evry, France
- Université Paris-Saclay , 91190 Saint-Aubin, France
| | - Marie-Pierre Gaigeot
- LAMBE CNRS UMR8587, Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, Université d'Evry val d'Essonne , Boulevard F. Mitterrand, Bat Maupertuis, 91025 Evry, France
- Université Paris-Saclay , 91190 Saint-Aubin, France
| |
Collapse
|
42
|
Galib M, Duignan TT, Misteli Y, Baer MD, Schenter GK, Hutter J, Mundy CJ. Mass density fluctuations in quantum and classical descriptions of liquid water. J Chem Phys 2017; 146:244501. [DOI: 10.1063/1.4986284] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Mirza Galib
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Timothy T. Duignan
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Yannick Misteli
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Marcel D. Baer
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Gregory K. Schenter
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Jürg Hutter
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Christopher J. Mundy
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| |
Collapse
|
43
|
Giberti F, Hassanali AA. The excess proton at the air-water interface: The role of instantaneous liquid interfaces. J Chem Phys 2017; 146:244703. [DOI: 10.1063/1.4986082] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Federico Giberti
- Institute for Molecular Engineering, The University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, USA
| | - Ali A. Hassanali
- Condensed Matter and Statistical Physics, The Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
| |
Collapse
|
44
|
del Rio BG, Souto J, Alemany MMG, González LE. Ab initio study of intrinsic profiles of liquid metals and their reflectivity. EPJ WEB OF CONFERENCES 2017. [DOI: 10.1051/epjconf/201715105002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
45
|
Muñoz-Santiburcio D, Marx D. On the complex structural diffusion of proton holes in nanoconfined alkaline solutions within slit pores. Nat Commun 2016; 7:12625. [PMID: 27550616 PMCID: PMC4996981 DOI: 10.1038/ncomms12625] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2016] [Accepted: 07/19/2016] [Indexed: 11/24/2022] Open
Abstract
The hydroxide anion OH−(aq) in homogeneous bulk water, that is, the solvated proton hole, is known to feature peculiar properties compared with excess protons solvated therein. In this work, it is disclosed that nanoconfinement of such alkaline aqueous solutions strongly affects the key structural and dynamical properties of OH−(aq) compared with the bulk limit. The combined effect of the preferred hypercoordinated solvation pattern of OH−(aq), its preferred perpendicular orientation relative to the confining surfaces, the pronounced layering of nanoconfined water and the topology of the hydrogen bond network required for proton hole transfer lead to major changes of the charge transport mechanism, in such a way that the proton hole migration mechanism depends exquisitely on the width of the confined space that hosts the water film. Moreover, the anionic Zundel complex, which is of transient nature in homogeneous bulk solutions, can be dynamically trapped as a shallow intermediate species by suitable nanoconfinement conditions. Confined liquids can exhibit very different properties compared with the bulk. Here, the authors show that the migration mechanism of the hydroxide anion in water is greatly affected by bi-dimensional nano-confinement and strongly depends on the width of the confined space.
Collapse
Affiliation(s)
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
46
|
Reif MM, Hünenberger PH. Origin of Asymmetric Solvation Effects for Ions in Water and Organic Solvents Investigated Using Molecular Dynamics Simulations: The Swain Acity-Basity Scale Revisited. J Phys Chem B 2016; 120:8485-517. [PMID: 27173101 DOI: 10.1021/acs.jpcb.6b02156] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The asymmetric solvation of ions can be defined as the tendency of a solvent to preferentially solvate anions over cations or cations over anions, at identical ionic charge magnitudes and effective sizes. Taking water as a reference, these effects are quantified experimentally for many solvents by the relative acity (A) and basity (B) parameters of the Swain scale. The goal of the present study is to investigate the asymmetric solvation of ions using molecular dynamics simulations, and to connect the results to this empirical scale. To this purpose, the charging free energies of alkali and halide ions, and of their hypothetical oppositely charged counterparts, are calculated in a variety of solvents. In a first set of calculations, artificial solvent models are considered that present either a charge or a shape asymmetry at the molecular level. The solvation asymmetry, probed by the difference in charging free energy between the two oppositely charged ions, is found to encompass a term quadratic in the ion charge, related to the different solvation structures around the anion and cation, and a term linear in the ion charge, related to the solvation structure around the uncharged ion-sized cavity. For these simple solvent models, the two terms are systematically counteracting each other, and it is argued that only the quadratic term should be retained when comparing the results of simulations involving physical solvents to experimental data. In a second set of calculations, 16 physical solvents are considered. The theoretical estimates for the acity A are found to correlate very well with the Swain parameters, whereas the correlation for B is very poor. Based on this observation, the Swain scale is reformulated into a new scale involving an asymmetry parameter Σ, positive for acitic solvents and negative for basitic ones, and a polarity parameter Π. This revised scale has the same predictive power as the original scale, but it characterizes asymmetry in an absolute sense, the atomistic simulations playing the role of an extra-thermodynamic assumption, and is optimally compatible with the simulation results. Considering the 55 solvents in the Swain set, it is observed that a moderate basity (Σ between -0.9 and -0.3, related to electronic polarization) represents the baseline for most solvents, while a highly variable acity (Σ between 0.0 and 3.0, related to hydrogen-bond donor capacity modulated by inductive effects) represents a landmark of protic solvents.
Collapse
Affiliation(s)
- Maria M Reif
- Physics Department (T38), Technische Universität München , D-85748 Garching, Germany
| | | |
Collapse
|
47
|
Nagata Y, Ohto T, Bonn M, Kühne TD. Surface tension of ab initio liquid water at the water-air interface. J Chem Phys 2016; 144:204705. [DOI: 10.1063/1.4951710] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
|
48
|
Perakis F, De Marco L, Shalit A, Tang F, Kann ZR, Kühne TD, Torre R, Bonn M, Nagata Y. Vibrational Spectroscopy and Dynamics of Water. Chem Rev 2016; 116:7590-607. [DOI: 10.1021/acs.chemrev.5b00640] [Citation(s) in RCA: 225] [Impact Index Per Article: 28.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Fivos Perakis
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Luigi De Marco
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
- Department of Chemistry and James Franck Institute, University of Chicago, 929 East 57th Street, Chicago, Illinois 60637, United States
| | - Andrey Shalit
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Fujie Tang
- International Center for Quantum Materials, Peking University, 5 Yiheyuan Road, Haidian, Beijing 100871, China
| | - Zachary R. Kann
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States,
| | - Thomas D. Kühne
- Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| | - Renato Torre
- European Lab for Nonlinear Spectroscopy and Dipartimento di Fisica e Astronomia, Università di Firenze, Via Nello Carrara 1, Sesto Fiorentino (Firenze) I-50019, Italy
| | - 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
|
49
|
Nagata Y, Ohto T, Backus EHG, Bonn M. Molecular Modeling of Water Interfaces: From Molecular Spectroscopy to Thermodynamics. J Phys Chem B 2016; 120:3785-96. [DOI: 10.1021/acs.jpcb.6b01012] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Tatsuhiko Ohto
- Graduate
School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
| | - Ellen H. G. Backus
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Mischa Bonn
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
50
|
John C, Spura T, Habershon S, Kühne TD. Quantum ring-polymer contraction method: Including nuclear quantum effects at no additional computational cost in comparison to ab initio molecular dynamics. Phys Rev E 2016; 93:043305. [PMID: 27176426 DOI: 10.1103/physreve.93.043305] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2016] [Indexed: 06/05/2023]
Abstract
We present a simple and accurate computational method which facilitates ab initio path-integral molecular dynamics simulations, where the quantum-mechanical nature of the nuclei is explicitly taken into account, at essentially no additional computational cost in comparison to the corresponding calculation using classical nuclei. The predictive power of the proposed quantum ring-polymer contraction method is demonstrated by computing various static and dynamic properties of liquid water at ambient conditions using density functional theory. This development will enable routine inclusion of nuclear quantum effects in ab initio molecular dynamics simulations of condensed-phase systems.
Collapse
Affiliation(s)
- Christopher John
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| | - Thomas Spura
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Thomas D Kühne
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany and Paderborn Center for Parallel Computing and Institute for Lightweight Design, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| |
Collapse
|