1
|
Gäding J, Della Balda V, Lan J, Konrad J, Iannuzzi M, Meißner RH, Tocci G. The role of the water contact layer on hydration and transport at solid/liquid interfaces. Proc Natl Acad Sci U S A 2024; 121:e2407877121. [PMID: 39259594 DOI: 10.1073/pnas.2407877121] [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: 04/22/2024] [Accepted: 08/04/2024] [Indexed: 09/13/2024] Open
Abstract
Understanding the structure in the nanoscopic region of water that is in direct contact with solid surfaces, so-called contact layer, is key to quantifying macroscopic properties that are of interest to e.g. catalysis, ice nucleation, nanofluidics, gas adsorption, and sensing. We explore the structure of the water contact layer on various technologically relevant solid surfaces, namely graphene, MoS[Formula: see text], Au(111), Au(100), Pt(111), and Pt(100), which have been previously hampered by time and length scale limitations of ab initio approaches or force field inaccuracies, by means of molecular dynamics simulations based on ab initio machine learning potentials built using an active learning scheme. Our results reveal that the in-plane intermolecular correlations of the water contact layer vary greatly among different systems: Whereas the contact layer on graphene and on Au(111) is predominantly homogeneous and isotropic, it is inhomogeneous and anisotropic on MoS[Formula: see text], on Au(100), and on the Pt surfaces, where it additionally forms two distinct sublayers. We apply hydrodynamics and the theory of the hydrophobic effect, to relate the energy corrugation and the characteristic length-scales of the contact layer with wetting, slippage, the hydration of small hydrophobic solutes and diffusio-osmotic transport. Thus, this work provides a microscopic picture of the water contact layer and links it to macroscopic properties of liquid/solid interfaces that are measured experimentally and that are relevant to wetting, hydrophobic solvation, nanofluidics, and osmotic transport.
Collapse
Affiliation(s)
- J Gäding
- Institute of Soft Matter Modeling, Hamburg University of Technology, Hamburg 21073, Germany
- Institute of Surface Science, Department of Atomistic Corrosion Informatics, Helmholtz-Zentrum Hereon, Geesthacht 21502, Germany
| | - V Della Balda
- Department of Chemistry, University of Zurich, Zurich 8057, Switzerland
| | - J Lan
- Department of Chemistry, New York University, New York, NY 10003
- Department of Chemistry, Simons Center for Computational Physical Chemistry at New York University, New York, NY 10003
| | - J Konrad
- Institute of Soft Matter Modeling, Hamburg University of Technology, Hamburg 21073, Germany
| | - M Iannuzzi
- Department of Chemistry, University of Zurich, Zurich 8057, Switzerland
| | - R H Meißner
- Institute of Soft Matter Modeling, Hamburg University of Technology, Hamburg 21073, Germany
- Institute of Surface Science, Department of Atomistic Corrosion Informatics, Helmholtz-Zentrum Hereon, Geesthacht 21502, Germany
| | - G Tocci
- Department of Chemistry, University of Zurich, Zurich 8057, Switzerland
| |
Collapse
|
2
|
Agosta L, Briels W, Hermansson K, Dzugutov M. The entropic origin of the enhancement of liquid diffusion close to a neutral confining surface. J Chem Phys 2024; 161:091102. [PMID: 39225520 DOI: 10.1063/5.0224016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2024] [Accepted: 08/13/2024] [Indexed: 09/04/2024] Open
Abstract
It is known that, in the proximity of a neutral wall, liquids experience diffusion enhancement relative to their bulk diffusion, but the origin of this phenomenon is still unknown. We report a molecular dynamics simulation investigating the dynamics of a simple liquid in the proximity to a non-interacting smooth confining wall, which exhibits a strong diffusion enhancement within the liquid layers adjacent to the wall. We present an analysis of these results, demonstrating that the observed diffusion enhancement can be accounted for, with numerical accuracy, using the universal scaling law that relates the liquid diffusion rate to the excess entropy. These results show that the scaling law, which has so far only been used for the description of the bulk liquid diffusion, can be successfully used to describe the diffusion in liquids under nano-scale confinement.
Collapse
Affiliation(s)
- Lorenzo Agosta
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
- Uppsala University, Ångström Laboratory, Department of Chemistry, 75121 Uppsala, Sweden
| | - Wim Briels
- University of Twente, Computational Chemical Physics, Postbus 217, Enschede 7500AE, Netherlands
- IBI-4:Biomacromolecular Systems and Processes, Forschungszentrum Jülich GmbH, Juelich D-52428, Germany
| | - Kersti Hermansson
- Uppsala University, Ångström Laboratory, Department of Chemistry, 75121 Uppsala, Sweden
| | - Mikhail Dzugutov
- Uppsala University, Ångström Laboratory, Department of Chemistry, 75121 Uppsala, Sweden
| |
Collapse
|
3
|
Bilichenko M, Iannuzzi M, Tocci G. Slip Opacity and Fast Osmotic Transport of Hydrophobes at Aqueous Interfaces with Two-Dimensional Materials. ACS NANO 2024; 18:24118-24127. [PMID: 39172927 DOI: 10.1021/acsnano.4c05118] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/24/2024]
Abstract
We investigate the interfacial transport of water and hydrophobic solutes on van der Waals bilayers and heterostructures formed by stacking graphene, hBN, and MoS2 using extensive ab initio molecular dynamics simulations. We compute water slippage and the diffusio-osmotic transport coefficient of hydrophobic particles at the interface by combining hydrodynamics and the theory of the hydrophobic effect. We find that slippage is dominated by the layer that is in direct contact with water and only marginally altered by the second layer, leading to a so-called "slip opacity". The screening of the lateral forces, where the liquid does not feel the forces coming from the second nearest layer, is one of the factors leading to the "slip opacity" in our systems. The diffusio-osmotic transport of small hydrophobes (with a radius below 2.5 Å) is also affected by the slip opacity, being dramatically enhanced by slippage. Furthermore, the direction of diffusio-osmotic flow is controlled by the solute size, with the flow in the opposite direction of the concentration gradient for smaller hydrophobes, and vice versa for larger ones. We connect our findings to the wetting properties of two-dimensional materials, and we propose that slippage and wetting can be controlled separately: whereas the slippage is mostly determined by the layer in closer proximity to water, wetting can be finely tuned by stacking different two-dimensional materials. Our study advances the computational design of two-dimensional materials and van der Waals heterostructures, enabling precise control over wetting and slippage properties for applications in coatings and water purification membranes.
Collapse
Affiliation(s)
- Maria Bilichenko
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Marcella Iannuzzi
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Gabriele Tocci
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| |
Collapse
|
4
|
Wang C, Tian W, Zhou K. Ab Initio Simulation of Liquid Water without Artificial High Temperature. J Chem Theory Comput 2024. [PMID: 39219067 DOI: 10.1021/acs.jctc.4c00650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
Comprehending the structure and dynamics of water is crucial in various fields, such as water desalination, ion separation, electrocatalysis, and biochemical processes. While reported works show that the ab initio molecular dynamics (AIMD) can accurately portray water's structure, the artificial high temperature (AHT) from 120 to 30 K is needed to mimic the quantum nature of hydrogen-bond network from GGA, metaGGA to hybrid functionals. The AHT proves to be an inadequate approach for systems involving aqueous multiphase mixtures, such as water-solid interfaces and aqueous solutions. This is due to the activation of additional phonons in other phases, which can lead to an overestimation of the dynamics of nearby water molecules. In this work, we find that the regularized SCAN (rSCAN) functional effectively captures both the structure and dynamics of liquid water at ambient conditions without AHT. Moreover, rSCAN closely matches experimental results for the hydration structures of alkali, alkali earth, and halide ions. We anticipate that the versatile and accurate rSCAN functional will emerge as a key tool based on ab initio simulation for investigating chemical processes in aqueous environments.
Collapse
Affiliation(s)
- Chenyu Wang
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou 215006, China
| | - Wei Tian
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou 215006, China
| | - Ke Zhou
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou 215006, China
| |
Collapse
|
5
|
Chen X, Cifuentes-Lopez A, Shao X, Lin L, Prokopchuk D, Pavanello M. Unraveling the Hydration Shell Structure and Dynamics of Group 10 Aqua Ions. J Phys Chem Lett 2024; 15:5517-5528. [PMID: 38749061 DOI: 10.1021/acs.jpclett.4c00464] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2024]
Abstract
We present ab initio simulations based on subsystem DFT of group 10 aqua ions accurately compared against experimental data on hydration structure. Our simulations provide insights into the molecular structures and dynamics of hydration shells, offering recalibrated interpretations of experimental results. We observe a soft, but distinct second hydration shell in Palladium (Pd) due to a balance between thermal fluctuations, metal-water interactions, and hydrogen bonding. Nickel (Ni) and platinum (Pt) exhibit more rigid hydration shells. Notably, our simulations align with experimental findings for Pd, showing axial hydration marked by a broad peak at about 3 Å in the Pd-O radial distribution function, revising the previously sharp "mesoshell" prediction. We introduce the "hydrogen bond dome" concept to describe a resilient network of hydrogen-bonded water molecules around the metal, which plays a critical role in the axial hydration dynamics.
Collapse
Affiliation(s)
- Xin Chen
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| | | | - Xuecheng Shao
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Lirong Lin
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Demyan Prokopchuk
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
6
|
Folkestad SD, Paul AC, Paul Née Matveeva R, Coriani S, Odelius M, Iannuzzi M, Koch H. Understanding X-ray absorption in liquid water using triple excitations in multilevel coupled cluster theory. Nat Commun 2024; 15:3551. [PMID: 38670938 PMCID: PMC11053016 DOI: 10.1038/s41467-024-47690-x] [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: 08/04/2023] [Accepted: 04/03/2024] [Indexed: 04/28/2024] Open
Abstract
X-ray absorption (XA) spectroscopy is an essential experimental tool to investigate the local structure of liquid water. Interpretation of the experiment poses a significant challenge and requires a quantitative theoretical description. High-quality theoretical XA spectra require reliable molecular dynamics simulations and accurate electronic structure calculations. Here, we present the first successful application of coupled cluster theory to model the XA spectrum of liquid water. We overcome the computational limitations on system size by employing a multilevel coupled cluster framework for large molecular systems. Excellent agreement with the experimental spectrum is achieved by including triple excitations in the wave function and using molecular structures from state-of-the-art path-integral molecular dynamics. We demonstrate that an accurate description of the electronic structure within the first solvation shell is sufficient to successfully model the XA spectrum of liquid water within the multilevel framework. Furthermore, we present a rigorous charge transfer analysis of the XA spectrum, which is reliable due to the accuracy and robustness of the electronic structure methodology. This analysis aligns with previous studies regarding the character of the prominent features of the XA spectrum of liquid water.
Collapse
Affiliation(s)
- Sarai Dery Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway
| | - Alexander C Paul
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway
| | - Regina Paul Née Matveeva
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway
| | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, DTU, 2800, Kongens Lyngby, Denmark
| | - Michael Odelius
- Department of Physics, Stockholm University, 10691, Stockholm, Sweden
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, 8057, Zürich, Switzerland
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway.
| |
Collapse
|
7
|
Villard J, Bircher MP, Rothlisberger U. Structure and dynamics of liquid water from ab initio simulations: adding Minnesota density functionals to Jacob's ladder. Chem Sci 2024; 15:4434-4451. [PMID: 38516095 PMCID: PMC10952088 DOI: 10.1039/d3sc05828j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Accepted: 02/12/2024] [Indexed: 03/23/2024] Open
Abstract
The accurate representation of the structural and dynamical properties of water is essential for simulating the unique behavior of this ubiquitous solvent. Here we assess the current status of describing liquid water using ab initio molecular dynamics, with a special focus on the performance of all the later generation Minnesota functionals. Findings are contextualized within the current knowledge on DFT for describing bulk water under ambient conditions and compared to experimental data. We find that, contrary to the prevalent idea that local and semilocal functionals overstructure water and underestimate dynamical properties, M06-L, revM06-L, and M11-L understructure water, while MN12-L and MN15-L overdistance water molecules due to weak cohesive effects. This can be attributed to a weakening of the hydrogen bond network, which leads to dynamical fingerprints that are over fast. While most of the hybrid Minnesota functionals (M06, M08-HX, M08-SO, M11, MN12-SX, and MN15) also yield understructured water, their dynamical properties generally improve over their semilocal counterparts. It emerges that exact exchange is a crucial component for accurately describing hydrogen bonds, which ultimately leads to corrections in both the dynamical and structural properties. However, an excessive amount of exact exchange strengthens hydrogen bonds and causes overstructuring and slow dynamics (M06-HF). As a compromise, M06-2X is the best performing Minnesota functional for water, and its D3 corrected variant shows very good structural agreement. From previous studies considering nuclear quantum effects (NQEs), the hybrid revPBE0-D3, and the rung-5 RPA (RPA@PBE) have been identified as the only two approximations that closely agree with experiments. Our results suggest that the M06-2X(-D3) functionals have the potential to further improve the reproduction of experimental properties when incorporating NQEs through path integral approaches. This work provides further proof that accurate modeling of water interactions requires the inclusion of both exact exchange and balanced (non-local) correlation, highlighting the need for higher rungs on Jacob's ladder to achieve predictive simulations of complex biological systems in aqueous environments.
Collapse
Affiliation(s)
- Justin Villard
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL) Lausanne CH-1015 Switzerland
| | - Martin P Bircher
- Computational and Soft Matter Physics, Universität Wien Wien A-1090 Austria
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL) Lausanne CH-1015 Switzerland
| |
Collapse
|
8
|
Gao Y, Wu J, Feng Y, Han J, Fang H. Structural effects of water clusters on viscosity at high shear rates. J Chem Phys 2024; 160:104502. [PMID: 38456533 DOI: 10.1063/5.0187906] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Accepted: 02/19/2024] [Indexed: 03/09/2024] Open
Abstract
In this study, we use molecular dynamics simulations of liquid water to investigate how shear thinning affects the viscosity of liquid water by structural changes of the hydrogen bond network. The effect of shear on viscosity can be divided into two parts: shear-induced destruction of the hydrogen bond network and the influence of the water structure on shear viscosity. First, strong shear destroys tetrahedral structures and thus reduces the connectivity of the hydrogen bond network. It is mainly because shear deformation, characterized by compression and expansion axes, respectively, triggers the destruction and formation of hydrogen bonds, resulting in anisotropic effects on water structures. At the same time, shear destroys large clusters and enhances the formation of small ones, resulting in a decrease in average cluster sizes. Second, the change of viscosity obeys a power law relationship with the change of hydrogen bond structures, highlighting a one-to-one correspondence between structure and property. Meanwhile, in order to explain why the structure affects viscosity, we define hydrogen-bond viscosity and find that the cooperative motion of the water structures can promote momentum transfer in the form of aggregations. Hydrogen-bond viscosity accounts for 5%-50% of the total viscosity. Our results elucidate that water structures are the important structural units to explain the change of water properties.
Collapse
Affiliation(s)
- Yitian Gao
- State Key Laboratory of Hydro-science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
| | - Jian Wu
- Department of Physics, Tsinghua University, Beijing 100084, China
| | - Yixuan Feng
- State Key Laboratory of Hydro-science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
| | - Jiale Han
- State Key Laboratory of Hydro-science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
| | - Hongwei Fang
- State Key Laboratory of Hydro-science and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
| |
Collapse
|
9
|
Quoika PK, Zacharias M. Liquid-Vapor Coexistence and Spontaneous Evaporation at Atmospheric Pressure of Common Rigid Three-Point Water Models in Molecular Simulations. J Phys Chem B 2024; 128:2457-2468. [PMID: 38427971 PMCID: PMC10945489 DOI: 10.1021/acs.jpcb.3c08183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2023] [Revised: 02/21/2024] [Accepted: 02/22/2024] [Indexed: 03/03/2024]
Abstract
Molecular dynamics (MD) simulations are widely used to investigate molecular systems at atomic resolution including biomolecular structures, drug-receptor interactions, and novel materials. Frequently, MD simulations are performed in an aqueous solution with explicit models of water molecules. Commonly, such models are parameterized to reproduce the liquid phase of water under ambient conditions. However, often, simulations at significantly higher temperatures are also of interest. Hence, it is important to investigate the equilibrium of the liquid and vapor phases of molecular models of water at elevated temperatures. Here, we evaluate the behavior of 11 common rigid three-point water models over a wide range of temperatures. From liquid-vapor coexistence simulations, we estimated the critical points and studied the spontaneous evaporation of these water models. Moreover, we investigated the influence of the system size, choice of the pressure-coupling algorithm, and rate of heating on the process and compared them with the experimental data. We found that modern rigid three-point water models reproduce the critical point surprisingly well. Furthermore, we discovered that the critical temperature correlates with the quadrupole moment of the respective water model. This indicates that the spatial arrangement of the partial charges is important for reproducing the liquid-vapor phase transition. Our findings may guide the selection of water models for simulations conducted at high temperatures.
Collapse
Affiliation(s)
- Patrick K. Quoika
- Center for Functional Protein
Assemblies, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, Garching 85748, Germany
| | - Martin Zacharias
- Center for Functional Protein
Assemblies, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, Garching 85748, Germany
| |
Collapse
|
10
|
Silvestrelli PL. Transport properties in liquids from first-principles: The case of liquid water and liquid argon. J Chem Phys 2023; 158:134503. [PMID: 37031126 DOI: 10.1063/5.0144353] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/05/2023] Open
Abstract
Shear and bulk viscosities of liquid water and argon are evaluated from first-principles in the density functional theory (DFT) framework, by performing molecular dynamics simulations in the NVE ensemble and using the Kubo–Greenwood equilibrium approach. The standard DFT functional is corrected in such a way to allow for a reasonable description of van der Waals effects. For liquid argon, the thermal conductivity has been also calculated. Concerning liquid water, to our knowledge, this is the first estimate of the bulk viscosity and of the shear-viscosity/bulk-viscosity ratio from first-principles. By analyzing our results, we can conclude that our first-principles simulations, performed at a nominal average temperature of 366 to guarantee that the systems are liquid-like, actually describe the basic dynamical properties of liquid water at about 330 K. In comparison with liquid water, the normal, monatomic liquid Ar is characterized by a much smaller bulk-viscosity/shear-viscosity ratio (close to unity) and this feature is well reproduced by our first-principles approach, which predicts a value of the ratio in better agreement with experimental reference data than that obtained using the empirical Lennard-Jones potential. The computed thermal conductivity of liquid argon is also in good agreement with the experimental value.
Collapse
Affiliation(s)
- Pier Luigi Silvestrelli
- Dipartimento di Fisica e Astronomia “G. Galilei,” Università di Padova, via Marzolo 8, I-35131 Padova, Italy
| |
Collapse
|