Simulations of vapor water clusters at vapor–liquid equilibrium

Liquid-Vapor Coexistence and Spontaneous Evaporation at Atmospheric Pressure of Common Rigid Three-Point Water Models in Molecular Simulations

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.

On the Different Faces of the Supercritical Phase of Water at a Near-Critical Temperature: Pressure-Induced Structural Transitions Ranging from a Gaslike Fluid to a Plastic Crystal Polymorph

The present study reports a systematic analysis of a wide variety of structural, thermodynamic, and dynamic properties of supercritical water along the near-critical isotherm of T = 1.03T_c and up to extreme pressures, using molecular dynamics and Monte Carlo simulations. The methodology employed provides solid evidence about the existence of a structural transition from a liquidlike fluid to a compressed, tightly packed liquid, in the density and pressure region around 3.4ρ_c and 1.17 GPa, introducing an alternative approach to locate the crossing of the Frenkel line. Around 8.5 GPa another transition to a face-centered-cubic plastic crystal polymorph with density 5.178ρ_c is also observed, further confirmed by Gibbs free energy calculations using the two-phase thermodynamic model. The isobaric heat capacity maximum, closely related to the crossing of the Widom line, has also been observed around 0.8ρ_c, where the local density augmentation is also maximized. Another structural transition has been observed at 0.2ρ_c, related to the transformation of the fluid to a dilute gas at lower densities. These findings indicate that a near-critical isotherm can be divided into different domains where supercritical water exhibits distinct behavior, ranging from a gaslike one to a plastic crystal one.

Collisional cross-section of water molecules in vapour studied by means of 1H relaxation in NMR

In gas phase, collisions that affect the rotational angular momentum lead to the return of the magnetization to its equilibrium (relaxation) in Nuclear Magnetic Resonance (NMR). To the best of our knowledge, the longitudinal relaxation rates R₁ = 1/T₁ of protons in H₂O and HDO have never been measured in gas phase. We report R₁ in gas phase in a field of 18.8 T, i.e., at a proton Larmor frequency ν₀ = 800 MHz, at temperatures between 353 and 373 K and pressures between 9 and 101 kPa. By assuming that spin rotation is the dominant relaxation mechanism, we estimated the effective cross-section σ_J for the transfer of angular momentum due to H₂O-H₂O and HDO-D₂O collisions. Our results allow one to test theoretical predictions of the intermolecular potential of water in gas phase.

The upside-down water collection system of Syntrichia caninervis

Desert plants possess highly evolved water conservation and transport systems, from the root structures that maximize absorption of scarce ground water(1-5), to the minimization of leaf surface area(6) to enhance water retention. Recent attention has focused on leaf structures that are adapted to collect water and promote nucleation from humid air(7-9). Syntrichia caninervis Mitt. (Pottiaceae) is one of the most abundant desert mosses in the world and thrives in an extreme environment with multiple but limited water resources (such as dew, fog, snow and rain), yet the mechanisms for water collection and transport have never been completely revealed. S. caninervis has a unique adaptation: it uses a tiny hair (awn) on the end of each leaf to collect water, in addition to that collected by the leaves themselves. Here we show that the unique multiscale structures of the hair are equipped to collect and transport water in four modes: nucleation of water droplets and films on the leaf hair from humid atmospheres; collection of fog droplets on leaf hairs; collection of splash water from raindrops; and transportation of the acquired water to the leaf itself. Fluid nucleation is accomplished in nanostructures, whereas fog droplets are gathered in areas where a high density of small barbs are present and then quickly transported to the leaf at the base of the hair. Our observations reveal nature's optimization of water collection by coupling relevant multiscale physical plant structures with multiscale sources of water.

Thermodynamic curvature from the critical point to the triple point

I evaluate the thermodynamic curvature R for fourteen pure fluids along their liquid-vapor coexistence curves, from the critical point to the triple point, using thermodynamic input from the NIST Chemistry WebBook. In this broad overview, R is evaluated in both the coexisting liquid and vapor phases. R is an invariant whose magnitude |R| is a measure of the size of mesoscopic organized structures in a fluid, and whose sign specifies whether intermolecular interactions are effectively attractive (R<0) or repulsive (R>0). I discuss five principles for R in pure fluids: (1) Near the critical point, the attractive part of the interactions forms loose structures of size |R| proportional to the correlation volume ξ(3), and the sign of R is negative. (2) In the vapor phase, there are instances of compact clusters of size |R| formed by the attractive part of the interactions and prevented from collapse by the repulsive part of the interactions, and the sign of R is positive. (3) In the asymptotic critical point regime, the R's in the coexisting liquid and vapor phases are equal to each other, i.e., commensurate. (4) Outside the asymptotic critical-point regime incommensurate R's may be associated with metastability. (5) The compact liquid phase has |R| on the order of the volume of a molecule, with the sign of R being negative for a liquidlike state held together by attractive interactions and the sign of R being positive for a solidlike state held up by repulsive interactions. These considerations amplify and extend the application of thermodynamic curvature in pure fluids.

Water absorption in polyethylene under external electric fields

Monte Carlo simulations of the solubility and structure of water in polyethylene in thermodynamic equilibrium with liquid water were performed in external fields ranging from 2 x 10(5) to 4 x 10(9) V/m. For a given equilibrium temperature and pressure, the water solubility decreases at higher fields. This occurs since it is energetically favorable for water molecules to be in the pure water phase than in the polyethylene matrix at high field strengths, and results in an increased density in the water phase. However, fields relevant to high voltage conduction (in the absence of defects that can lead to large local field strengths) do not change the solubility. In addition, at large fields the number of water clusters decreases for all cluster sizes. The rate of decrease is highest for large clusters, and a larger fraction of water molecules exist as monomers in the polyethylene matrix at high fields. Large fields also cause alignment of the water molecules, which leads to more clusters with linear topologies and hence an increase in the cluster radius of gyration.

Monte Carlo simulations of equilibrium solubilities and structure of water in n-alkanes and polyethylene

Gibbs ensemble Monte Carlo methods based on a force field that combines the simple point charge [Berendsen et al., in Intermolecular Forces, edited by Pullman (Reidel, Dordrecht, 1981), p. 331] and transferable potentials for phase equilibria [Martin and Siepmann, J. Phys. Chem. B 102, 2659 (1998)] models were used to study the equilibrium properties of binary systems consisting of water and n-alkanes with chain lengths from hexane to hexadecane. In addition, systems where extended linear alkane chains (up to 300 carbon units long) were used to represent amorphous polyethylene were simulated in the presence of water using a connectivity altering osmotic Gibbs ensemble. In these simulations the equilibrium between a liquid water phase and a polymer phase into which water was inserted was studied. The predicted solubilities, which were determined between 350 and 550 K, are in good agreement with experiment, where experimental results are available, and the density of water molecules in the hydrocarbons is approximately 63% as high as in saturated water vapor under the same conditions. At the lower temperatures most of the water exists as monomers; increasing the temperature leads to an increase in the density of water in the alkane phase and hence in the fraction of molecules that participate in clusters. Dimers are the most prevalent clusters in all hydrocarbons and at all temperatures studied, and the fraction of clusters of given size decrease with increasing cluster size. A large fraction of trimers, tetramers, and pentamers, which are the cluster sizes for which topologies have been studied, are cyclic at low temperatures, but at higher temperatures linear structures predominate. The same properties are observed for pure water vapor clusters in equilibrium with the liquid phase, showing that the cluster topologies are not significantly affected by the surrounding hydrocarbon.

Local Density Inhomogeneities and Dynamics in Supercritical Water: A Molecular Dynamics Simulation Approach

Molecular dynamics atomistic simulations in the canonical ensemble (NVT-MD) have been used to investigate the "Local Density Inhomogeneities and their Dynamics" in pure supercritical water. The simulations were carried out along a near-critical isotherm (Tr = T/Tc = 1.03) and for a wide range of densities below and above the critical one (0.2 rho(c) - 2.0 rho(c)). The results obtained reveal the existence of significant local density augmentation effects, which are found to be sufficiently larger in comparison to those reported for nonassociated fluids. The time evolution of the local density distribution around each molecule was studied in terms of the appropriate time correlation functions C(Delta)rhol(t). It is found that the shape of these functions changes significantly by increasing the density of the fluid. Finally, the local density reorganization times for the first and second coordination shell derived from these correlations exhibit a decreasing behavior by increasing the density of the system, signifying the density effect upon the dynamics of the local environment around each molecule.