Phase diagrams and surface properties of modified water models

A Neural-Network-Optimized Hydrogen Peroxide Pairwise Additive Model for Classical Simulations

We have developed an all-atom pairwise additive model for hydrogen peroxide using an optimization procedure based on artificial neural networks (ANNs). The model is based on experimental molecular geometry and includes a dihedral potential that hinders the cis-type configuration and allows for crossing the trans one, defined between the planes that have the two oxygen atoms and each hydrogen. The model's parametrization is achieved by training simple ANNs to minimize a target function that measures the differences between various thermodynamic and transport properties and the corresponding experimental values. Finally, we evaluated a range of properties for the optimized model and its mixtures with SPC/E water, including bulk-liquid properties (density, thermal expansion coefficient, adiabatic compressibility, etc.) and properties of systems at equilibrium (vapor and liquid density, vapor pressure and composition, surface tension, etc.). Overall, we obtained good agreement with experimental data.

Surface tension of liquids and binary mixtures from molecular dynamics simulations

In this work we assess and extend strategies for calculating surface tension of complex liquids from molecular dynamics simulations: the mechanical route and the instantaneous liquid interface (ILI) approach. The former employs the connection between stress tensor and surface tension, whereas the latter involves computation of instantaneous density field. Whereas the mechanical route is general, the ILI method involves system-dependent parameters restricting its original application to liquid water only. Here we generalize the approach to complex molecular liquids using atomic van der Waals radii. The performance of the approaches is evaluated on two liquid systems: acetonitrile and water-methanol mixture. In addition, we compare the effect of the computational models for interaction potentials based on semi-empirical electronic structure theory and classical force fields on the estimate of the surface tension within both stress tensor and ILI approaches.

Coupled Monte Carlo and Molecular Dynamics Simulations on Interfacial Properties of Antifouling Polymer Membranes

We use molecular simulation to study the wetting behavior of antifouling polymer-tethered membranes. We obtain the interfacial properties (e.g., contact angle) of water at various temperatures for five polymer membranes, including a base polysulfone (PSF) membrane and four other PSF membranes grafted with antifouling polymers (two poly(ethylene glycol) (PEG) tethers and two zwitterionic tethers). We implement a coupled Monte Carlo (MC)/molecular dynamics (MD) approach to determine the interface potentials of water on the membrane surfaces in an efficient manner. Within this method, short MC and MD simulations are performed in cycles to collect the surface excess free energy of a thin water film on polymer membrane surfaces. Simulation results show that the grafting of zwitterionic tethers provides a more significant enhancement in the hydrophilicity of the PSF membrane than that of the PEG tethers. Water completely wets the surface of zwitterionic polymer membranes.

Insight into the molecular mechanism of water evaporation via the finite temperature string method

The process of water's evaporation at its liquid/air interface has proven challenging to study experimentally and, because it constitutes a rare event on molecular time scales, presents a challenge for computer simulations as well. In this work, we simulated water's evaporation using the classical extended simple point charge model water model, and identified a minimum free energy path for this process in terms of 10 descriptive order parameters. The measured free energy change was 7.4 kcal/mol at 298 K, in reasonable agreement with the experimental value of 6.3 kcal/mol, and the mean first-passage time was 1375 ns for a single molecule, corresponding to an evaporation coefficient of 0.25. In the observed minimum free energy process, the water molecule diffuses to the surface, and tends to rotate so that its dipole and one O-H bond are oriented outward as it crosses the Gibbs dividing surface. As the water molecule moves further outward through the interfacial region, its local density is higher than the time-averaged density, indicating a local solvation shell that protrudes from the interface. The water molecule loses donor and acceptor hydrogen bonds, and then, with its dipole nearly normal to the interface, stops donating its remaining hydrogen bond. At that point, when the final, accepted hydrogen bond is broken, the water molecule is free. We also analyzed which order parameters are most important in the process and in reactive trajectories, and found that the relative orientation of water molecules near the evaporating molecule, and the number of accepted hydrogen bonds, were important variables in reactive trajectories and in kinetic descriptions of the process.

Predicting cavity formation free energy: how far is the Gaussian approximation valid?

We examine the range of validity of the Gaussian model for various water-like liquids whose intermolecular potentials differ from SPC/E water, to provide insight into the temperature dependence of the hydrophobic effect for small hard sphere solutes. We find that low compressibility liquids that have more close-packed network structures show much larger deviations from Gaussian fluctuations for low or zero occupancies relative to more compressible fluids with more open networks. Water appears to be a unique molecular fluid in possessing equilibrium density fluctuations that are faithfully described by the Gaussian theory. We ascribe this success to the fact, shown here, that the orientational correlations near a small hard sphere solute involve remarkably little reorganization from the bulk, which is a consequence of water's low solvent reorganization enthalpy and entropy.

COMPUTER SIMULATION OF LIQUID–VAPOR INTERFACIAL TENSION: LENNARD-JONES FLUID AND WATER REVISITED

We review several commonly used simulation methods for computing liquid–vapor surface tension and associated theoretical treatments of the long-range correction for inhomogeneous systems. Prototype model systems considered in this review are the Lennard-Jones (LJ) fluid and the SPC/E model water. In addition, we examine a variety of factors that can affect calculation of the surface tension γ via the mechanical approach (i.e. using either KB or IK method). It is found that for the LJ fluid, the size of simulation box and the number of particles in the system can have notable effects on the computed surface tension. For SPC/E water, the Ewald parameters can influence computed surface tensions (γ) as well, e.g., very small Ewald parameters tend to overestimate γ. It is also found that the IK method consistently gives γ that are 0.6 - 0.9 mN/m greater than γ computed based on the KB method. When computing the first reciprocal–space contribution to the surface tension, the Ghoufi's strategy gives rise to more sensible profile of pressure difference PN(z)-PT(z) than the Alejandre's strategy although both strategies result in nearly the same average surface tension through the integration of PN(z)-PT(z).

Water and aqueous solutions: simple non-speculative model approach

Different ways of molecular modeling of water are analyzed and their similarities and differences identified. An up-to-date summary of achievements of a general approach to common rigid site-site interaction models of molecular fluids applied to water and aqueous solutions is then presented and discussed. The method is based on considering only a short-range part of a total realistic potential (such as SPC/E or TIPxP) which determines the structure of water (and fluids in general). A simplification of the interactions at short intermolecular separations leads then to simple models, called primitive models. Quite accurate results in an analytic form for the thermodynamic properties of the models are obtained using the thermodynamic perturbation theory. It is shown that the properly constructed primitive models reproduce, qualitatively, anomalies of pure water and basic characteristics of hydrophobic hydration. The concept of an extended excluded volume, based on pseudo-hard bodies, is introduced and exemplified by considering the partial molar volume of apolar solutes. Potential future development towards a theory of water based on the primitive models as a reference with the long-range contributions added as a perturbation is discussed.

Towards understanding water: simulation of modified water models

The unusual properties of liquid water have been attributed to factors including hydrogen bonds, the presence of a local tetrahedral network, or the high number density. Most of the properties of water can be reproduced in simulations using simple potential models with site-site pairwise interactions. Although in real water the factors listed above cannot be separated, one may construct liquids from modified water models ('not-waters') in which these factors are altered independently in a controlled way. With various collaborators we have studied the changes in liquid properties, anomalies and solvation which result from various modifications.