Can the AMOEBA forcefield be used for high pressure simulations? The extreme case of methane and water.

We have performed classical molecular dynamics simulations using the fully polarizable Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) forcefield implemented within the Tinker package to determine whether a more adequate treatment of electrostatics is sufficient to correctly describe the mixing of methane with water under high pressure conditions. We found a significant difference between the ability of AMOEBA and other classical, computationally cheaper forcefields, such as TIP3P, simple point charge-extended, TIP4P, and optimized potentials for liquid simulations-all atom. While the latter models fail to detect any effect of pressure on the miscibility of methane in water, AMOEBA qualitatively captures the experimental observation of the increased solubility of methane in water with pressure. At higher temperatures, the solubility of water in methane also increases; this seems to be associated with the breakdown of the fourfold hydrogen-bonded water network structure: bonding in water is weaker, so the energy cost of solution is lowered.

The re-entrant transition from the molecular to atomic phases of dense fluids: The case of hydrogen.

A simple phenomenological thermodynamic model is developed to describe the chemical bonding and unbonding in homonuclear diatomic systems. This model describes the entire phase diagram of dimer-forming systems and shows a transition from monomers to dimers, with monomers favored at both very low and very high pressures, as well as at high temperatures. In the context of hydrogen, the former region corresponds to hydrogen present in most interstellar gas clouds, while the latter is associated with the long sought-after fluid metallic phase. The model predicts a molecular to atomic fluid transition in dense deuterium, which is in agreement with recently reported experimental measurements.

From Atoms to Colloids: Does the Frenkel Line Exist in Discontinuous Potentials?

The Frenkel line has been proposed as a crossover in the fluid region of phase diagrams between a "nonrigid" and a "rigid" fluid. It is generally described as a crossover in the dynamical properties of a material and as such has been described theoretically using a very different set of markers from those with which is it investigated experimentally. In this study, we have performed extensive calculations using two simple yet fundamentally different model systems: hard spheres and square-well potentials. The former has only hardcore repulsion, while the latter also includes a simple model of attraction. We computed and analyzed a series of physical properties used previously in simulations and experimental measurements and discuss critically their correlations and validity as to being able to uniquely and coherently locate the Frenkel line in discontinuous potentials.

Krypton and the Fundamental Flaw of the Lennard-Jones Potential.

We have performed a series of neutron scattering experiments on supercritical krypton. Our data and analysis allow us to characterize the Frenkel line crossover in this model monatomic fluid. The data from our measurements was analyzed using Empirical Potential Structure Refinement to determine the short- and medium-range structure of the fluids. We find evidence for several shells of neighbors which form approximately concentric rings of density about each atom. The ratio of second to first shell radius is significantly larger than in any crystal structure. Modeling krypton using a Lennard-Jones potential is shown to give significant errors, notably that the liquid is overstructured. The true potential appears to be longer ranged and with a softer core than the 6-12 powerlaws permit.

How to determine solubility in binary mixtures from neutron scattering data: The case of methane and water.

It has recently been discovered that, when subjected to moderate amounts of pressure, methane dissolves in water to form binary mixtures of up to 40% molar methane. No significant solubility of water in methane is known. In these mixtures, the water hydrogen-bond network is largely complete and surrounds the methane molecules. The discovery of this dense mixture has once again highlighted the technical difficulties involved in accurately describing and sampling mixing phenomena both computationally and experimentally. Here, we present a systematic and critical study of the methods employed to characterize binary mixtures and their robustness. This study highlights the requirements needed to develop a quantitative understanding, and it proposes new and more accessible measures of miscibility to investigators, particularly for in silico analysis.

Frenkel Line in Nitrogen Terminates at the Triple Point.

Recent studies on supercritical nitrogen revealed clear changes in structural markers and dynamical properties when the coordination number approaches its maximum value. The line in P and T space where these changes occur is referred to as the Frenkel line. Here, we qualitatively reproduce such changes in the supercritical regime using the popular "optimized potential for liquid simulation" (OPLS) classical force field for molecular dynamics. Unfortunately, at 160 K, OPLS nitrogen predicts sublimation rather than producing a liquid phase; therefore, we developed our own force field to achieve quantitative agreement with experimental data. We confirm the asymptotic behavior of the coordination number on crossing the Frenkel line and note an associated change in the diffusion constant, consistent with the non-rigid to rigid liquid-like character of the "transition". The simulations allow us to track the Frenkel line to subcritical temperatures and demonstrate that it terminates at the triple point. This establishes the experimentally measurable changes, which could unequivocally determine the Frenkel line in other systems.

Structural Markers of the Frenkel Line in the Proximity of Widom Lines.

We have performed a neutron scattering experiment on supercritical fluid nitrogen at 160 K (1.27 TC) over a wide pressure range (7.8 MPa/0.260 g/mL-125 MPa/0.805 g/mL). This has enabled us to study the process by which nitrogen changes from a fluid that exhibits gaslike behavior to one that exhibits rigid liquidlike behavior at a temperature close to, but above, the critical temperature by crossing the Widom lines followed by the Frenkel line on pressure (density) increase. We find that the Frenkel line transition is indicated by a transition to a regime of rigid liquidlike behavior in which the coordination number remains constant within error, in agreement with our previous work at 300 K. The Frenkel line transition takes place at approximately the same density at 160 and 300 K. The data do not conclusively show an additional transition at the location of the known Widom lines. We find that behavior remains gaslike until the Frenkel line is crossed and our data support the hypothesis that Widom line transitions are density increase-driven.

Squeezing Oil into Water under Pressure: Inverting the Hydrophobic Effect.

The molecular structure of dense homogeneous fluid water-methane mixtures has been determined for the first time using high-pressure neutron-scattering techniques at 1.7 and 2.2 GPa. A mixed state with a fully H-bonded water network is revealed. The hydration shell of the methane molecules is, however, revealed to be pressure-dependent with an increase in the water coordination between 1.7 and 2.2 GPa. In parallel, ab initio molecular dynamics simulations have been performed to provide insight into the microscopic mechanisms associated with the phenomenon of mixing. These calculations reproduce the observed phase change from phase separation to mixing with increasing pressure. The calculations also reproduce the experimentally observed structural properties. Unexpectedly, the simulations show mixing is accompanied by a subtle enhancement of the polarization of methane. Our results highlight the key role played by fine electronic effects on miscibility and the need to readjust our fundamental understanding of hydrophobicity to account for these.

Pressure-Induced Miscibility Increase of CH4 in H2O: A Computational Study Using Classical Potentials.

Methane and water demix under normal (ambient) pressure and temperature conditions because of the polar nature of water and the apolar nature of methane. Recent experimental work has shown, though, that increasing the pressure to values between 1 and 2 GPa (10-20 kbar) leads to a marked increase of methane solubility in water, for temperatures which are well below the critical temperature for water. Here, we perform molecular dynamics simulations based on classical force fields-which are well-used and have been validated at ambient conditions-for different values of pressure and temperature. We find the expected increase in miscibility for mixtures of methane and supercritical water; however, our model fails to reproduce the experimentally observed increase in methane solubility at large pressures and below the critical temperature of water. This points to the need to develop more accurate force fields for methane and methane-water mixtures under pressure.

When immiscible becomes miscible-Methane in water at high pressures.

At low pressures, the solubility of gases in liquids is governed by Henry's law, which states that the saturated solubility of a gas in a liquid is proportional to the partial pressure of the gas. As the pressure increases, most gases depart from this ideal behavior in a sublinear fashion, leveling off at pressures in the 1- to 5-kbar (0.1 to 0.5 GPa) range with solubilities of less than 1 mole percent (mol %). This contrasts strikingly with the well-known marked increase in solubility of simple gases in water at high temperature associated with the critical point (647 K and 212 bar). The solubility of the smallest hydrocarbon, the simple gas methane, in water under a range of pressure and temperature is of widespread importance, because it is a paradigmatic hydrophobe and occurs widely in terrestrial and extraterrestrial geology. We report measurements up to 3.5 GPa of the pressure dependence of the solubility of methane in water at 100°C-well below the latter's critical temperature. Our results reveal a marked increase in solubility between 1 and 2 GPa, leading to a state above 2 GPa where the maximum solubility of methane in water exceeds 35 mol %.