How Water's Properties Are Encoded in Its Molecular Structure and Energies.

How are water's material properties encoded within the structure of the water molecule? This is pertinent to understanding Earth's living systems, its materials, its geochemistry and geophysics, and a broad spectrum of its industrial chemistry. Water has distinctive liquid and solid properties: It is highly cohesive. It has volumetric anomalies-water's solid (ice) floats on its liquid; pressure can melt the solid rather than freezing the liquid; heating can shrink the liquid. It has more solid phases than other materials. Its supercooled liquid has divergent thermodynamic response functions. Its glassy state is neither fragile nor strong. Its component ions-hydroxide and protons-diffuse much faster than other ions. Aqueous solvation of ions or oils entails large entropies and heat capacities. We review how these properties are encoded within water's molecular structure and energies, as understood from theories, simulations, and experiments. Like simpler liquids, water molecules are nearly spherical and interact with each other through van der Waals forces. Unlike simpler liquids, water's orientation-dependent hydrogen bonding leads to open tetrahedral cage-like structuring that contributes to its remarkable volumetric and thermal properties

Force Field Optimization, Advanced Sampling, And Free Energy Methods With Gpu-Optimized Monte Carlo (gomc) Software

In this work, to address the sampling problem for systems at high densities and low temperatures, a generalized identity exchange algorithm is developed for grand canonical Monte Carlo simulations. The algorithm, referred to as Molecular Exchange Monte Carlo (MEMC), is implemented in the GPU-Optimized Monte Carlo (GOMC) software and may be applied to multicomponent systems of arbitrary molecular topology, and provides significant enhancements in the sampling of phase space over a wide range of compositions and temperatures. Three different approaches are presented for the insertion/deletion of the large molecules, and the pros and cons of each method are discussed. Next, the MEMC method is extended to Gibbs ensemble Monte Carlo (GEMC). The utility of the MEMC method is demonstrated through the calculation of the free energies of transfer of n-alkanes from vapor into liquid 1-octanol, n-hexadecane, and 2,2,4-trimethylpentane, using isobaric-isothermal GEMC simulations. Alternatively, for system with strong inter-molecular interaction (e.g. hydrogen bonds), it’s more efficient to calculate the free energies of transfer, using standard thermodynamic integration (TI) and free energy perturbation (FEP) methods. The TI and FEP free energy calculation methods are implemented in GOMC and utility of these methods are demonstrated by calculating the hydration and solvation free energies of fluorinated 1-octanol, to understand the role of fluorination on the interactions and partitioning of alcohols in aqueous and organic environments. Additionally, using GOMC, a transferable united-atom (UA) force field, based on Mie potentials, is optimized for alkynes to accurately reproduce experimental phase equilibrium properties. The performance of the optimized Mie potential parameters is assessed for 1-alkynes and 2-alkynes using grand canonical histogram-reweighting Monte Carlo simulations. For each compound, vapor-liquid coexistence curves, vapor pressures, heats of vaporization, critical properties, and normal boiling points are predicted and compared to experiment

Adapting the semi-explicit assembly solvation model for estimating water-cyclohexane partitioning with the SAMPL5 molecules

We describe here some tests we made in the SAMPL5 communal event of 'Semi-Explicit Assembly' (SEA), a recent method for computing solvation free energies. We combined the prospective tests of SAMPL5 with followup retrospective calculations, to improve two technical aspects of the field variant of SEA. First, SEA uses an approximate analytical surface around the solute on which a water potential is computed. We have improved and simplified the mathematical model of that surface. Second, some of the solutes in SAMPL5 were large enough to need a way to treat solvating waters interacting with 'buried atoms', i.e. interior atoms of the solute. We improved SEA with a buried-atom correction. We also compare SEA to Thermodynamic Integration molecular dynamics simulations, so that we can sort out force field errors

Adapting the semi-explicit assembly solvation model for estimating water-cyclohexane partitioning with the SAMPL5 molecules

We describe here some tests we made in the SAMPL5 communal event of `Semi-Explicit Assembly' (SEA), a recent method for computing solvation free energies. We combined the prospective tests of SAMPL5 with followup retrospective calculations, to improve two technical aspects of the field variant of SEA. First, SEA uses an approximate analytical surface around the solute on which a water potential is computed. We have improved and simplified the mathematical model of that surface. Second, some of the solutes in SAMPL5 were large enough to need a way to treat solvating waters interacting with `buried atoms', i.e. interior atoms of the solute. We improved SEA with a buried-atom correction. We also compare SEA to Thermodynamic Integration molecular dynamics simulations, so that we can sort out force field errors