-
Potential Energy Landscape of a Flexible Water Model: Equation-of-State, Configurational Entropy, and Adam-Gibbs Relationship
Authors:
Ali Eltareb,
Gustavo E. Lopez,
Nicolas Giovambattista
Abstract:
The potential energy landscape (PEL) formalism is a tool within statistical mechanics that has been used in the past to calculate the equation of states (EOS) of classical rigid model liquids at low temperatures, where computer simulations may be challenging. In this work, we use classical molecular dynamics (MD) simulations and the PEL formalism to calculate the EOS of the flexible q-TIP4P/F wate…
▽ More
The potential energy landscape (PEL) formalism is a tool within statistical mechanics that has been used in the past to calculate the equation of states (EOS) of classical rigid model liquids at low temperatures, where computer simulations may be challenging. In this work, we use classical molecular dynamics (MD) simulations and the PEL formalism to calculate the EOS of the flexible q-TIP4P/F water model. This model exhibits a liquid-liquid critical point (LLCP) in the supercooled regime, at ($P_c = 150$ MPa, $T_c = 190$ K, $ρ_c = 1.04$ g/cm$^3$) [using the reaction field technique]. The PEL-EOS of q-TIP4P/F water, and the corresponding location of the LLCP, are in very good agreement with the MD simulations. We show that the PEL of q-TIP4P/F water is Gaussian which allows us to calculate the configurational entropy of the system, $S_{conf}$. The $S_{conf}$ of q-TIP4P/F water is surprisingly similar to that reported previously for rigid water models, suggesting that intramolecular flexibility does not necessarily add roughness to the PEL. We also show that the Adam-Gibbs relation, which relates the diffusion coefficient $D$ with $S_{conf}$, holds for the flexible q-TIP4P/F water model. Overall, our results indicate that the PEL formalism can be used to study molecular systems that include molecular flexibility, the common case in standard force fields. This is not trivial since the introduction of large bending/stretching mode frequencies is problematic in classical statistical mechanics. For example, as shown previously, we find that such high-frequencies lead to an unphysical (negative) entropy for q-TIP4P/F water (yet the PEL formalism can be applied successfully).
△ Less
Submitted 17 January, 2024;
originally announced January 2024.
-
Searching for crystal-ice domains in amorphous ices
Authors:
Fausto Martelli,
Nicolas Giovambattista,
Salvatore Torquato,
Roberto Car
Abstract:
We employ classical molecular dynamics simulations to investigate the molecular-level structure of water during the isothermal compression of hexagonal ice (I$h$) and low-density amorphous (LDA) ice at low temperatures. In both cases, the system transforms to high-density amorphous ice (HDA) via a first-order-like phase transition. We employ a sensitive local order metric (LOM) [Martelli et. al.,…
▽ More
We employ classical molecular dynamics simulations to investigate the molecular-level structure of water during the isothermal compression of hexagonal ice (I$h$) and low-density amorphous (LDA) ice at low temperatures. In both cases, the system transforms to high-density amorphous ice (HDA) via a first-order-like phase transition. We employ a sensitive local order metric (LOM) [Martelli et. al., Phys. Rev. B, 97, 064105 (2018)], that can discriminate among different crystalline and non crystalline ice structures and is based on the positions of the oxygen atoms in the first and/or second hydration shell. Our results confirm that LDA and HDA are indeed amorphous, i.e., they lack of polydispersed ice domains. Interestingly, HDA contains a small number of domains that are reminiscent of the unit cell of ice IV, although the hydrogen-bond network (HBN) of these domains differ from the HBN of ice IV. The presence of ice IV-like domains provides some support to the hypothesis that HDA could be the result of a detour on the HBN rearrangement along the I$h$-to-ice IV pressure induced transformation. Both nonequilibrium LDA-to-HDA and I$h$-to-HDA transformations are two-steps processes where a small distortion of the HBN first occurs at low pressures and then, a sudden, extensive re-arrangement of hydrogen bonds at the corresponding transformation pressure follows. Interestingly, the I$h$-to-HDA and LDA-to-HDA transformations occur when LDA and I$h$ have similar local order, as quantified by the site-averaged LOMs. Since I$h$ has a perfect tetrahedral HBN, while LDA does not, it follows that higher pressures are needed to transform I$h$ into HDA than that for the conversion of LDA to HDA. In correspondence with both first-order-like phase transitions, the samples are composed of a large HDA cluster that percolates within the I$h$/LDA samples.
△ Less
Submitted 7 May, 2018;
originally announced May 2018.
-
Effect of pressure on the phase behavior and structure of water confined between nanoscale hydrophobic and hydrophilic plates
Authors:
Nicolas Giovambattista,
Peter J. Rossky,
Pablo G. Debenedetti
Abstract:
We perform systematic molecular dynamics simulations of water confined between two nanoscale plates at T=300 K. We investigate the effect of pressure (-0.15 GPa <= P <= 0.2 GPa) and plate separation (0.4 nm <= d <= 1.6 nm) on the phase behavior of water when the plates are either hydrophobic or hydrophilic. When water is confined between hydrophobic plates, capillary evaporation occurs between t…
▽ More
We perform systematic molecular dynamics simulations of water confined between two nanoscale plates at T=300 K. We investigate the effect of pressure (-0.15 GPa <= P <= 0.2 GPa) and plate separation (0.4 nm <= d <= 1.6 nm) on the phase behavior of water when the plates are either hydrophobic or hydrophilic. When water is confined between hydrophobic plates, capillary evaporation occurs between the plates at low enough P. The threshold value of d at which this transition occurs decreases with P (e.g., 1.6 nm at P ~ -0.05 GPa; 0.5 nm at P \~ 0.1 GPa), until, at high P, no capillary evaporation occurs. For d ~ 0.6 nm and P >= 0.1 GPa, the system crystallizes into a bilayer ice.
A P-d phase diagram showing the vapor, liquid, and bilayer ice phases is proposed. When water is confined by hydrophilic (hydroxylated silica) plates, it remains in the liquid phase at all P and d studied. Interestingly, we observe, for this case, that even at the P at which bulk water cavitates, the confined water remains in the liquid state. We also study systematically the state of hydration at different P for both kinds of plates. For the range of conditions studied here, we find that in the presence of hydrophobic plates the effect of P is to enhance water structure and to push water molecules toward the plates. The average orientation of water molecules next to the hydrophobic plates does not change upon pressurization. In contrast, in the presence of hydrophilic plates, water structure is insensitive to P. Hence, our results suggest that upon pressurization, hydrophobic plates behave as `soft' surfaces (in the sense of accommodating pressure-dependent changes in water structure) while hydrophilic walls behave as `hard' surfaces.
△ Less
Submitted 1 March, 2006;
originally announced March 2006.