The surface tension of TIP4P/2005 water model using the Ewald sums for the dispersion interactions

Second inflection point of supercooled water surface tension induced by hydrogen bonds: A molecular-dynamics study

Surface tension of supercooled water is a fundamental property in various scientific processes. In this study, we perform molecular dynamics simulations with the TIP4P-2005 model to investigate the surface tension of supercooled water down to 220 K. Our results show a second inflection point (SIP) in the surface tension at temperature TSIP ≈ 267.5 ± 2.3 K. Using an extended IAPWS-E functional fit for the water surface tension, we calculate the surface excess internal-energy and entropy terms of the excess Helmholtz free energy. Similar to prior studies [Wang et al., Phys. Chem. Chem. Phys. 21, 3360 (2019); Gorfer et al., J. Chem. Phys. 158, 054503 (2023)], our results show that the surface tension is governed by two driving forces: a surface excess entropy change above the SIP and a surface excess internal-energy change below it. We study hydrogen-bonding near the SIP because it is the main cause of water's anomalous properties. With decreasing temperature, our results show that the entropy contribution to the surface tension reaches a maximum slightly below the SIP and then decreases. This is because the number of hydrogen bonds increases more slowly below the SIP. Moreover, the strengths and lifetimes of the hydrogen bonds also rise dramatically below the SIP, causing the internal-energy term to dominate the excess surface free energy. Thus, the SIP in the surface tension of supercooled TIP4P-2005 water is associated with an increase in the strengths and lifetimes of hydrogen bonds, along with a decrease in the formation rate (#/K) of new hydrogen bonds.

Phase transitions of fluorotelomer alcohols at the water|alkane interface studied via molecular dynamics simulation

Fluorosurfactants are long-lasting environmental pollutants that accumulate at interfaces ranging from aerosol droplet surfaces to cell membranes. Modeling of adsorption-based removal technologies for fluorosurfactants requires accurate simulation methods which can predict their adsorption isotherm and monolayer structure. Fluorotelomer alcohols with one or two methylene groups adjacent to the alcohol (7 : 1 FTOH and 6 : 2 FTOH, respectively) are investigated using the OPLS-AA force field at the water|hexane interface, varying the interfacial area per surfactant. The acquired interfacial pressure isotherms and monolayer phase behavior are compared with previous experimental results. The results are consistent with the experimental data inasmuch as, at realistic adsorption densities, only 7 : 1 FTOH shows a phase transition between liquid-expanded (LE) and 2D crystalline phases. Structures of the LE and crystalline phases are in good agreement with the sticky disc and Langmuir defective crystal models, respectively, used previously to interpret experimental data. Interfacial pressure of the LE phase agrees well with experiment, and sticky disc interaction parameters indicate no 2D LE-gas transition is present for either molecule. Conformation analysis reveals 7 : 1 FTOH favors conformers where the OH dipole is perpendicular to the molecular backbone, such that the crystalline phase is stabilized when these dipoles align.

Relation between Double Layer Structure, Capacitance, and Surface Tension in Electrowetting of Graphene and Aqueous Electrolytes

Deciphering the mechanisms of charge storage on carbon-based materials is pivotal for the development of next-generation electrochemical energy storage systems. Graphene, the building block of graphitic electrodes, is an ideal model for probing such processes on a fundamental level. Herein, we investigate the thermodynamics of the graphene/aqueous electrolyte interface by utilizing a multiscale quantum mechanics-classical molecular dynamics (QM/MD) approach to provide insights into the effect of alkali metal ion (Li+) concentration on the interfacial tension (γSL) of the charged graphene/electrolyte interface. We demonstrate that the dependence of γSL on the applied surface charge exhibits an asymmetric behavior relative to the neutral surface. At the positively charged graphene sheet, the electrowetting response is amplified by electrolyte concentration, resulting in a strongly hydrophilic surface. On the contrary, at negative potential bias, γSL shows a weaker response to the charging of the electrode. Changes in γSL greatly affect the total areal capacitance predicted by the Young-Lippmann equation but have a negligible impact on the simulated total areal capacitance, indicating that the EDL structure is not directly correlated with the wettability of the surface and different interfacial mechanisms drive the two phenomena. The proposed model is validated experimentally by studying the electrowetting response of highly oriented pyrolytic graphite over a wide range of electrolyte concentrations. Our work presents the first combined theoretical and experimental study on electrowetting using carbon surfaces, introducing new conceptual routes for the investigation of wetting phenomena under potential bias.

Hydrophobicity of molecular-scale textured surfaces: The case of zeolitic imidazolate frameworks, an atomistic perspective

Hydrophobicity has proven fundamental in an inexhaustible amount of everyday applications. Material hydrophobicity is determined by chemical composition and geometrical characteristics of its macroscopic surface. Surface roughness or texturing enhances intrinsic hydrophilic or hydrophobic characteristics of a material. Here we consider crystalline surfaces presenting molecular-scale texturing typical of crystalline porous materials, e.g., metal-organic frameworks. In particular, we investigate one such material with remarkable hydrophobic qualities, ZIF-8. We show that ZIF-8 hydrophobicity is driven not only by its chemical composition but also its sub-nanoscale surface corrugations, a physical enhancement rare amongst hydrophobes. Studying ZIF-8's hydrophobic properties is challenging as experimentally it is difficult to distinguish between the materials' and the macroscopic corrugations' contributions to the hydrophobicity. The computational contact angle determination is also difficult as the standard "geometric" technique of liquid nanodroplet deposition is prone to many artifacts. Here, we characterise ZIF-8 hydrophobicity via: (i) the "geometric" approach and (ii) the "energetic" method, utilising the Young-Dupré formula and computationally determining the liquid-solid adhesion energy. Both approaches reveal nanoscale Wenzel-like bathing of the corrugated surface. Moreover, we illustrate the importance of surface linker termination in ZIF-8 hydrophobicity, which reduces when varied from sp3 N to sp2 N termination. We also consider halogenated analogues of the methyl-imidazole linker, which promote the transition from nanoWenzel-like to nanoCassie-Baxter-like states, further enhancing surface hydrophobicity. Present results reveal the complex interface physics and chemistry between water and complex porous, molecular crystalline surfaces, providing a hint to tune their hydrophobicity.

Functionalized carbon nanocones performance in water harvesting

In this work, we investigate the water capture process for functionalized carbon nanocones (CNCs) through molecular dynamic simulations in the following three scenarios: a single CNC in contact with a reservoir containing liquid water, a single CNC in contact with a water vapor reservoir, and a combination of more than one CNC in contact with vapor. We found that water flows through the nanocones when in contact with the liquid reservoir if the nanocone tip presents hydrophilic functionalization. In contact with steam, we observed the formation of droplets at the base of the nanocone only when hydrophilic functionalization is present. Then, water flows through in a linear manner, a process that is more efficient than that in the liquid reservoir regime. The scalability of the process is tested by analyzing the water flow through more than one nanocone. The results suggest that the distance between the nanocones is a fundamental ingredient for the efficiency of water harvesting.

A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water

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.

Optimization of the wetting-drying characteristics of hydrophobic metal organic frameworks via crystallite size: The role of hydrogen bonding between intruded and bulk liquid

HYPOTHESIS

The behavior of Heterogeneous Lyophobic Systems (HLSs) comprised of a lyophobic porous material and a corresponding non-wetting liquid is affected by a variety of different structural parameters of the porous material. Dependence on exogenic properties such as crystallite size is desirable for system tuning as they are much more facilely modified. We explore the dependence of intrusion pressure and intruded volume on crystallite size, testing the hypothesis that the connection between internal cavities and bulk water facilitates intrusion via hydrogen bonding, a phenomenon that is magnified in smaller crystallites with a larger surface/volume ratio.

EXPERIMENTS

Water intrusion/extrusion pressures and intrusion volume were experimentally measured for ZIF-8 samples of various crystallite sizes and compared to previously reported values. Alongside the practical research, molecular dynamics simulations and stochastic modeling were performed to illustrate the effect of crystallite size on the properties of the HLSs and uncover the important role of hydrogen bonding within this phenomenon.

FINDINGS

A reduction in crystallite size led to a significant decrease of intrusion and extrusion pressures below 100 nm. Simulations indicate that this behavior is due to a greater number of cages being in proximity to bulk water for smaller crystallites, allowing cross-cage hydrogen bonds to stabilize the intruded state and lower the threshold pressure of intrusion and extrusion. This is accompanied by a reduction in the overall intruded volume. Simulations demonstrate that this phenomenon is linked to ZIF-8 surface half-cages exposed to water being occupied by water due to non-trivial termination of the crystallites, even at atmospheric pressure.

Transferability of vibrational spectroscopic map from TIP4P to TIP4P-like water models

We computed the IR, Raman, and sum frequency generation spectra of water in the OH-stretch region by employing the quantum/classical mixed approach that consists of a vibrational spectroscopic map and molecular dynamics (MD) simulation. We carried out the MD simulation with the TIP4P, TIP4P/2005, and TIP4P/Ice models and applied the map designed for TIP4P by Skinner et al. to each MD trajectory. Although the map is not tuned for TIP4P-like models, TIP4P/2005 and TIP4P/Ice provide the best reproduction of the experimental vibrational spectra of liquid water and crystalline ice, respectively. This result demonstrates the transferability of the map from TIP4P to TIP4P/2005 and TIP4P/Ice, meaning that one can choose an appropriate TIP4P-like model to calculate the vibrational spectra of an aqueous system without rebuilding the map.

The water cavitation line as predicted by the TIP4P/2005 model

The formation of vapor bubbles in a metastable liquid, cavitation, is an activated process due to the free energy cost of having both phases at contact. Such an energetic penalty enables the existence of the liquid beyond its thermodynamic borders. Establishing the stability limits of a liquid as ubiquitous as water has important practical implications and has thereby attracted a lot of attention. Different experimental strategies and theoretical analyses have been employed to measure and predict the cavitation line, or the pressure-temperature kinetic stability border of liquid water. Understanding the location of the cavitation line requires knowing the cavitation rate dependence on pressure and temperature. Such dependency is difficult to obtain in experiments, and we use molecular simulations with the TIP4P/2005 model to fill this gap. By deeply overstretching liquid water below the saturation pressure, we are able to observe and quantify spontaneous cavitation. To deal with a lower overstretching regime, we resort to the Seeding technique, which consists of analyzing simulations of a liquid containing a vapor bubble under the theoretical framework of Classical Nucleation Theory. Combining spontaneous cavitation with Seeding, we get a wide overview of the cavitation rate. We study two different temperatures (450 and 550 K) and complement our perspective with the results previously obtained at 296.4 K [Menzl et al., Proc. Natl. Acad. Sci. 113, 13582 (2016)] to establish a broad simulation-experiment comparison. We find a good agreement between simulations and both isobaric heating and isochoric cooling experiments using quartz inclusions. We are, however, unable to reconcile simulations with other experimental techniques. Our results predict a decrease in the solid-liquid interfacial free energy as the liquid becomes increasingly overstretched with a temperature independent Tolman length of 0.1 nm. Therefore, the capillarity approximation underestimates the nucleation rate. Nonetheless, it provides a fair indication of the location of the cavitation line given the steep rate vs pressure dependence. Overall, our work provides a comprehensive view of the water cavitation phenomenon and sets an efficient strategy to investigate it with molecular simulations.

Effects of molecular size and orientation on the interfacial properties and wetting behavior of water/n-alkane systems: a molecular-dynamics study

Molecular dynamics simulations (MD) are performed to study the interfacial structure/tension and wetting behavior of water/n-alkane systems (water/nC5 to water/nC16 where nCx = CxH(2x + 2)). In particular, we study complete-to-partial wetting transitions by changing the n-alkane chain length (N_C) at a constant temperature, T = 295 K. Simulations are carried out with a united-atom TraPPE model for n-alkanes and the TIP4P-2005 model of water. Simulation results are in excellent agreement with the initial spreading coefficients and contact angles calculated using experimental values of the surface and interfacial tensions. In addition, it has been determined that water/(nC5-nC7) and water/(nC8-nC16), respectively, exhibit complete and partial initial wetting modes. Simulations show that the interfacial structures of water/(nC5-nC7) are different from water/(nC8-nC16) systems. In the latter, water preferentially orients near the interface to increase the number of hydrogen bonds and the charge and mass densities. Moreover, the orientation of n-alkane molecules at water/(nC8-nC16) interfaces has a long-range persistence, resulting in layered structures that increase with N_C. In addition, simulation results of the orientational order parameter S_z show alignment behavior of the n-alkane molecules with respect to the interfaces. Simulations predict that the central segments of n-alkane are strongly packed in the interfaces while the end segments (methyl groups) form smaller peaks in the outer edge of the layer. This observation confirms the "horseshoe" or "C-shaped" structure of n-alkane molecules in the water/n-alkane interfaces. At constant temperature, the interface widths of both water and the n-alkanes decrease with increasing n-alkane molecular length. These results suggest that increasing the n-alkane chain length affects the water/n-alkane interfacial properties in a manner similar to that of cooling.

Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water

In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na⁺ and Cl^- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.

Atmospheric water harvesting using functionalized carbon nanocones

In this work, we propose a method to harvest liquid water from water vapor using carbon nanocones. The condensation occurs due to the presence of hydrophilic sites at the nanocone entrance. The functionalization, together with the high mobility of water inside nanostructures, leads to a fast water flow through the nanostructure. We show using molecular dynamics simulations that this device is able to collect water if the surface functionalization is properly selected.

The lung surfactant activity probed with molecular dynamics simulations

The surface of pulmonary alveolar subphase is covered with a mixture of lipids and proteins. This lung surfactant plays a crucial role in lung functioning. It shows a complex phase behavior which can be altered by the interaction with third molecules such as drugs or pollutants. For studying multicomponent biological systems, it is of interest to couple experimental approach with computational modelling yielding atomic-scale information. Simple two, three, or four-component model systems showed to be useful for getting more insight in the interaction between lipids, lipids and proteins or lipids and proteins with drugs and impurities. These systems were studied theoretically using molecular dynamic simulations and experimentally by means of the Langmuir technique. A better understanding of the structure and behavior of lung surfactants obtained from this research is relevant for developing new synthetic surfactants for efficient therapies, and may contribute to public health protection.

Can molecular simulations reliably compare homogeneous and heterogeneous ice nucleation?

In principle, the answer to the posed titular question is undoubtedly "yes." But in practice, requisite reference data for homogeneous systems have been obtained with a treatment of intermolecular interactions that is different from that typically employed for heterogeneous systems. In this article, we assess the impact of the choice of truncation scheme when comparing water in homogeneous and inhomogeneous environments. Specifically, we use explicit free energy calculations and a simple mean field analysis to demonstrate that using the "cut-and-shift" version of the Lennard-Jones potential (common to most simple point charge models of water) results in a systematic increase in the melting temperature of ice I_h. In addition, by drawing an analogy between a change in cutoff and a change in pressure, we use existing literature data for homogeneous ice nucleation at negative pressures to suggest that enhancements due to heterogeneous nucleation may have been overestimated by several orders of magnitude.

Flap Dynamics in Pepsin-Like Aspartic Proteases: A Computational Perspective Using Plasmepsin-II and BACE-1 as Model Systems

The flexibility of β hairpin structure known as the flap plays a key role in catalytic activity and substrate intake in pepsin-like aspartic proteases. Most of these enzymes share structural and sequential similarity. In this study, we have used apo Plm-II and BACE-1 as model systems. In the apo form of the proteases, a conserved tyrosine residue in the flap region remains in a dynamic equilibrium between the normal and flipped states through rotation of the χ1 and χ2 angles. Independent MD simulations of Plm-II and BACE-1 remained stuck either in the normal or flipped state. Metadynamics simulations using side-chain torsion angles (χ1 and χ2 of tyrosine) as collective variables sampled the transition between the normal and flipped states. Qualitatively, the two states were predicted to be equally populated. The normal and flipped states were stabilized by H-bond interactions to a tryptophan residue and to the catalytic aspartate, respectively. Further, mutation of tyrosine to an amino-acid with smaller side-chain, such as alanine, reduced the flexibility of the flap and resulted in a flap collapse (flap loses flexibility and remains stuck in a particular state). This is in accordance with previous experimental studies, which showed that mutation to alanine resulted in loss of activity in pepsin-like aspartic proteases. Our results suggest that the ring flipping associated with the tyrosine side-chain is the key order parameter that governs flap dynamics and opening of the binding pocket in most pepsin-like aspartic proteases.

Understanding the Stabilization of a Bulk Nanobubble: A Molecular Dynamics Analysis

Bulk nanobubbles (NBs) have received considerable attention because of their extensive potential applications, such as in ultrasound imaging and water management. Although multiple types of experimental evidence have supported the existence and stabilization of bulk NBs, the underlying mechanism remains unclear. This study numerically investigates the bulk NB stabilization with molecular dynamics (MD) methods: the all-atom (AA) MD simulation is used for NBs of several nanometers diameter; the coarse-grained (CG) MD simulation is for the NBs of about 100 nm. The NB properties are statistically obtained and analyzed, including the inner density, inner pressure, surface charge, interfacial hydrogen bond (HB), and gaseous diffusion. The results show that the gas inside an NB has ultrahigh density (tens of kilograms per cubic meter). A double-layer surface charge exists on the NB. The inner/outer layer is positively/negatively charged, and the electrostatic stress can counteract part of the surface tension. In addition, the interfacial HB is weakened by the interaction between gas and water molecules, causing less surface tension. The above features are beneficial to NB stabilization. The NB equilibrium radii solved by the interfacial mechanical equilibrium equation agree with the MD results, indicating that this equation can describe the force balance of an NB as small as several nanometers. Besides, supersaturation appears to be necessary for the NB thermodynamic equilibrium. Based on Henry's law and the ideal gas law, the theoretical analysis suggests that the stability of the NB thermodynamic equilibrium is conditional: the number of gas molecules in NBs should be more than half that dissolved in liquid. This study unravels a stabilized bulk NB's properties and discusses the NB equilibrium and stabilization mechanism, which will advance the understanding and application of bulk NBs.

Thermophysical properties of water using reactive force fields

The widescale importance and rich phenomenology of water continue to motivate the development of computational models. ReaxFF force fields incorporate many characteristics desirable for modeling aqueous systems: molecular flexibility, polarization, and chemical reactivity (bond formation and breaking). However, their ability to model the general properties of water has not been evaluated in detail. We present comprehensive benchmarks of the thermophysical properties of water for two ReaxFF models, the water-2017 and CHON-2017_weak force fields. These include structural, electrostatic, vibrational, thermodynamic, coexistence, and transport properties at ambient conditions (300 K and 0.997 g cm^-3) and along the standard pressure (1 bar) isobar. Overall, CHON-2017_weak predicts more accurate thermophysical properties than the water-2017 force field. Based on our results, we recommend potential avenues for improvement: the dipole moment to quadrupole moment ratio, the self-diffusion coefficient, especially for water-2017, and the gas phase vibrational frequencies with the aim to improve the vibrational properties of liquid water.

Molecular dynamics simulations of interfacial properties of the CO2-water and CO2-CH4-water systems

Molecular dynamics simulations were performed to study the interfacial behavior of the pure carbon dioxide-water system and a binary 40:60 mol. % gas mixture of (carbon dioxide + methane)-water at the temperatures of 275.15 K and 298.15 K and pressures near 4 MPa for CO₂ and up to 10 MPa for methane. The simulations are used to study the dynamic equilibrium of the gases at the water-gas interface, to determine the z-density profiles for the gases and water, and calculate the interfacial tension γ under the different temperature/pressure conditions close to those of the formation of clathrate hydrates of these gases. At the same hydrostatic gas phase pressure, the CO₂-water interface has a lower interfacial tension than the CH₄-water interface. A greater number of CO₂ molecules, as much as three times more than methane at the same pressure, were adsorbed at the interfacial layer, which reflects the stronger electrostatic quadrupolar and van der Waals interactions between CO₂ and water molecules at the interface. The water surfaces are covered by less than a monolayer of gas even when the pressure of the system goes near the saturation pressure of CO₂. The surface adsorbed molecules are in dynamic equilibrium with the bulk gas and with exchange between the gas and interface regions occurring repeatedly within the timescale of the simulations. The effects of the changes in the CO₂-water interfacial tension with external temperature and pressure conditions on the formation of the clathrate hydrates and other CO₂ capture and sequestration processes are discussed.

Molecular modeling of aqueous electrolytes at interfaces: Effects of long-range dispersion forces and of ionic charge rescaling

Molecular dynamics simulations of aqueous electrolytes generally rely on empirical force fields, combining dispersion interactions-described by a truncated Lennard-Jones (LJ) potential-and electrostatic interactions-described by a Coulomb potential computed with a long-range solver. Recently, force fields using rescaled ionic charges [electronic continuum correction (ECC)], possibly complemented with rescaling of LJ parameters [ECC rescaled (ECCR)], have shown promising results in bulk, but their performance at interfaces has been less explored. Here, we started by exploring the impact of the LJ potential truncation on the surface tension of a sodium chloride aqueous solution. We show a discrepancy between the numerical predictions for truncated LJ interactions with a large cutoff and for untruncated LJ interactions computed with a long-range solver, which can bias comparison of force field predictions with experiments. Using a long-range solver for LJ interactions, we then show that an ionic charge rescaling factor chosen to correct long-range electrostatic interactions in bulk accurately describes image charge repulsion at the liquid-vapor interface, and the rescaling of LJ parameters in ECCR models-aimed at capturing local ion-ion and ion-water interactions in bulk- describes well the formation of an ionic double layer at the liquid-vapor interface. Overall, these results suggest that the molecular modeling of aqueous electrolytes at interfaces would benefit from using long-range solvers for dispersion forces and from using ECCR models, where the charge rescaling factor should be chosen to correct long-range electrostatic interactions.

Liquid exfoliation of multilayer graphene in sheared solvents: A molecular dynamics investigation

Liquid-phase exfoliation, the use of a sheared liquid to delaminate graphite into few-layer graphene, is a promising technique for the large-scale production of graphene. However, the microscale and nanoscale fluid-structure processes controlling the exfoliation are not fully understood. Here, we perform non-equilibrium molecular dynamics simulations of a defect-free graphite nanoplatelet suspended in a shear flow and measure the critical shear rate γ̇_c needed for the exfoliation to occur. We compare γ̇_c for different solvents, including water and N-methyl-pyrrolidone, and nanoplatelets of different lengths. Using a theoretical model based on a balance between the work done by viscous shearing forces and the change in interfacial energies upon layer sliding, we are able to predict the critical shear rates γ̇_c measured in simulations. We find that an accurate prediction of the exfoliation of short graphite nanoplatelets is possible only if both hydrodynamic slip and the fluid forces on the graphene edges are considered and if an accurate value of the solid-liquid surface energy is used. The commonly used "geometric-mean" approximation for the solid-liquid energy leads to grossly incorrect predictions.

Charge transfer as a ubiquitous mechanism in determining the negative charge at hydrophobic interfaces

The origin of the apparent negative charge at hydrophobic–water interfaces has fueled debates in the physical chemistry community for decades. The most common interpretation given to explain this observation is that negatively charged hydroxide ions (OH^–) bind strongly to the interfaces. Using first principles calculations of extended air–water and oil–water interfaces, we unravel a mechanism that does not require the presence of OH^–. Small amounts of charge transfer along hydrogen bonds and asymmetries in the hydrogen bond network due to topological defects can lead to the accumulation of negative surface charge at both interfaces. For water near oil, some spillage of electron density into the oil phase is also observed. The computed surface charge densities at both interfaces is approximately \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-0.015\ {\rm{e}}/{{\rm{nm}}}^{2}$$\end{document}−0.015e∕nm2 in agreement with electrophoretic experiments. We also show, using an energy decomposition analysis, that the electronic origin of this phenomena is rooted in a collective polarization/charge transfer effect.

The accumulation of negative charge at hydrophobic–water interfaces has been a source of debate for a long time. Here the authors use ab initio calculations to show that the charge accumulation at air–water and oil–water interfaces is caused by subtle charge transfer processes.

Leveraging local MP2 to reduce basis set superposition errors: An efficient first-principles based force-field for carbon dioxide

Pairwise additive model potentials for CO₂ were developed by fitting to gradients computed with the local second order Møller Plesset Perturbation theory (LMP2) method, with and without consideration of 3-body dispersion using adaptive force matching. Without fitting to experiments, all models gave good predictions of properties of CO₂, such as the density-temperature diagram, diffusion constants, and radial distribution functions. For the prediction of vibrational spectra, the inclusion of a bond-bond coupling term has been shown to be important. The CO₂ models developed only have pairwise additive terms, thus allowing microsecond time scale simulations to be performed with practical computational cost. LMP2 performed significantly better than second order Møller Plesset Perturbation theory (MP2) for the development of the CO₂ model. This is attributed to the appreciable reduction in the basis set superposition error when the localized method was used. It is argued that LMP2 is a more appropriate method than MP2 for force matching for systems where the basis set superposition error is large.

Ultrafast Propulsion of Water Nanodroplets on Patterned Graphene

The directed transport of liquids at the nanoscale is of great importance for nanotechnology applications ranging from water filtration to the cooling of electronics and precision medicine. Here we demonstrate such unidirectional, pumpless transport of water nanodroplets on graphene sheets patterned with hydrophilic/phobic areas inspired by natural systems. We find that spatially varying patterning of the graphene surfaces can lead to water transport at ultrafast velocities, far exceeding macroscale estimates. We perform extensive molecular dynamics simulations to show that such high transport velocities ( O(10² m/s)) are due to differences of the advancing and receding contact angles of the moving droplet. This contact angle hysteresis and the ensuing transport depend on the surface pattern and the droplet size. We present a scaling law for the driving capillary and resisting friction forces on the water droplet and use it to predict nanodroplet trajectories on a wedge-patterned graphene sheet. The present results demonstrate that graphene with spatially variable wettability is a potent material for fast and precise transport of nanodroplets with significant potential for directed nanoscale liquid transport and precision drug delivery.

Cavitation inception of water with solid nanoparticles: A molecular dynamics study

Cavitation in liquid with impurities is important in heterogeneous nucleation applications. One of the most widely existing kinds of impurities is solid particles, which can be found in natural water from rivers and specially prepared water such as nanofluids. Understanding the effects caused by the existence of nanoparticles on cavitation in water is vital to the rapidly developed nanotechnologies and medical researches. In this study, cavitation in water with nanoparticles is investigated through molecular dynamics simulations. The effects by nanoparticle materials and sizes on cavitation are discussed by using SiO₂ and polyethylene spherical nanoparticles with different diameters. The nucleation rate and the formation of critical bubbles in cavitation are studied via the Voronoi tessellation and the mean first passage time methods. The hydrogen bond network in water is also analyzed. Results reveal that SiO⁠₂ and polyethylene nanoparticles may destabilize the hydrogen bond network in water. With the same particle size, cavitation in water with polyethylene nanoparticles is promoted to a greater extent than that with SiO₂ nanoparticles. With the same nanoparticle material, cavitation is promoted with the increase in particle size in a range spanning half to ten times the critical bubble radius. Beyond this range, particle size has little influence on cavitation. Reasons for those effects on cavitation due to the presence of solid nanoparticles are discussed by analysing the changes of hydrogen bonds network in water.

Unexpected finite size effects in interfacial systems: Why bigger is not always better-Increase in uncertainty of surface tension with bulk phase width

We present an unexpected finite size effect affecting interfacial molecular simulations that is proportional to the width-to-surface-area ratio of the bulk phase L_l/A. This finite size effect has a significant impact on the variance of surface tension values calculated using the virial summation method. A theoretical derivation of the origin of the effect is proposed, giving a new insight into the importance of optimising system dimensions in interfacial simulations. We demonstrate the consequences of this finite size effect via a new way to estimate the surface energetic and entropic properties of simulated air-liquid interfaces. Our method is based on macroscopic thermodynamic theory and involves comparing the internal energies of systems with varying dimensions. We present the testing of these methods using simulations of the TIP4P/2005 water forcefield and a Lennard-Jones fluid model of argon. Finally, we provide suggestions of additional situations, in which this finite size effect is expected to be significant, as well as possible ways to avoid its impact.

Atomistic Model for Nearly Quantitative Simulations of Langmuir Monolayers

Lung surfactant and a tear film lipid layer are examples of biologically relevant macromolecular structures found at the air-water interface. Because of their complexity, they are often studied in terms of simplified lipid layers, the simplest example being a Langmuir monolayer. Given the profound biological significance of these lipid assemblies, there is a need to understand their structure and dynamics on the nanoscale, yet there are not many techniques able to provide this information. Atomistic molecular dynamics simulations would be a tool fit for this purpose; however, the simulation models suggested until now have been qualitative instead of quantitative. This limitation has mainly stemmed from the challenge to correctly describe the surface tension of water with simulation parameters compatible with other biomolecules. In this work, we show that this limitation can be overcome by using the recently introduced four-point OPC water model, whose surface tension for water is demonstrated to be quantitatively consistent with experimental data and which is also shown to be compatible with the commonly employed lipid models. We further establish that the approach of combining the OPC four-point water model with the CHARMM36 lipid force field provides nearly quantitative agreement with experiments for the surface pressure-area isotherm for POPC and DPPC monolayers, also including the experimentally observed phase coexistence in a DPPC monolayer. The simulation models reported in this work pave the way for nearly quantitative atomistic studies of lipid-rich biological structures at air-water interfaces.

Multiple Staggered Mesh Ewald: Boosting the Accuracy of the Smooth Particle Mesh Ewald Method

The smooth particle mesh Ewald (SPME) method is the standard method for computing the electrostatic interactions in the molecular simulations. In this work, we develop the multiple staggered mesh Ewald (MSME) method, which averages the SPME forces computed on, for example, M, staggered meshes. We prove, from a theoretical perspective, that the MSME is as accurate as the SPME, but uses M² times less mesh points in a certain parameter range. In the complementary parameter range, the MSME is as accurate as the SPME with twice the interpolation order. The theoretical conclusions are numerically validated both by a uniform and uncorrelated charge system, and by a three-point-charge water system that is widely used as solvent for the biomacromolecules.

The liquid-vapor equilibria of TIP4P/2005 and BLYPSP-4F water models determined through direct simulations of the liquid-vapor interface

In this paper, the surface tension and critical properties for the TIP4P/2005 and BLYPSP-4F models are reported. A clear dependence of surface tension on the van der Waals cutoff radius (rvdw) is shown when van der Waals interactions are modeled with a simple cutoff scheme. A linear extrapolation formula is proposed that can be used to determine the infinite rvdw surface tension through a few simulations with finite rvdw. A procedure for determining liquid and vapor densities is proposed that does not require fitting to a profile function. Although the critical temperature of water is also found to depend on the choice of rvdw, the dependence is weaker. We argue that a rvdw of 1.75 nm is a good compromise for water simulations when long-range van der Waals correction is not applied. Since the majority of computational programs do not support rigorous treatment of long-range dispersion, the establishment of a minimal acceptable rvdw is important for the simulation of a variety of inhomogeneous systems, such as water bubbles, and water in confined environments. The BLYPSP-4F model predicts room temperature surface tension marginally better than TIP4P/2005 but overestimates the critical temperature. This is expected since only liquid configurations were fit during the development of the BLYPSP-4F potential. The potential is expected to underestimate the stability of vapor and thus overestimate the region of stability for the liquid.

Coherent description of transport across the water interface: From nanodroplets to climate models

Transport of mass and energy across the vapor-liquid interface of water is of central importance in a variety of contexts such as climate models, weather forecasts, and power plants. We provide a complete description of the transport properties of the vapor-liquid interface of water with the framework of nonequilibrium thermodynamics. Transport across the planar interface is then described by 3 interface transfer coefficients where 9 more coefficients extend the description to curved interfaces. We obtain all coefficients in the range 260-560 K by taking advantage of water evaporation experiments at low temperatures, nonequilibrium molecular dynamics with the TIP4P/2005 rigid-water-molecule model at high temperatures, and square gradient theory to represent the whole range. Square gradient theory is used to link the region where experiments are possible (low vapor pressures) to the region where nonequilibrium molecular dynamics can be done (high vapor pressures). This enables a description of transport across the planar water interface, interfaces of bubbles, and droplets, as well as interfaces of water structures with complex geometries. The results are likely to improve the description of evaporation and condensation of water at widely different scales; they open a route to improve the understanding of nanodroplets on a small scale and the precision of climate models on a large scale.

Computer modelling of the surface tension of the gas–liquid and liquid–liquid interface

This review presents the state of the art in molecular simulations of interfacial systems and of the calculation of the surface tension from the underlying intermolecular potential.

Wetting and evaporation of salt-water nanodroplets: A molecular dynamics investigation

We employ molecular dynamics simulations to study the wetting and evaporation of salt-water nanodroplets on platinum surfaces. Our results show that the contact angle of the droplets increases with the salt concentration. To verify this, a second simulation system of a thin salt-water film on a platinum surface is used to calculate the various surface tensions. We find that both the solid-liquid and liquid-vapor surface tensions increase with salt concentration and as a result these cause an increase in the contact angle. However, the evaporation rate of salt-water droplets decreases as the salt concentration increases, due to the hydration of salt ions. When the water molecules have all evaporated from the droplet, two forms of salt crystals are deposited, clump and ringlike, depending on the solid-liquid interaction strength and the evaporation rate. To form salt crystals in a ring, it is crucial that there is a pinned stage in the evaporation process, during which salt ions can move from the center to the rim of the droplets. With a stronger solid-liquid interaction strength, a slower evaporation rate, and a higher salt concentration, a complete salt crystal ring can be deposited on the surface.