Phase diagram of water under an applied electric field

Molecular dynamics simulations of the effect of static electric field on progressive ice formation

Ice accumulation under static electric fields presents a significant hazard to transmission lines and power grids. Contemporary computational studies of electrofreezing predominantly probed excessive electric fields (109 V/m) that are significantly higher than those typically encountered in proximity to transmission lines. To elucidate the influence of realistic electric fields (105 V/m) on ice crystallization, we run extensive molecular dynamics (MD) simulations across dual ice-water coexistence systems. Three aspects of work were accordingly examined. First, we investigated the influence of the effect of static electric fields, with a strength of 105 V/m, along three orthogonal axes on the phase transition during the encountered freezing and melting processes. Second, we established the mechanism of how the direction of an electric field, the initial ice crystallography, and the adjacent crystal planes influence the solidification process. Third, the results of our MD simulations were further post-processed to determine the dipole moment, radial distribution, and angle distribution resulting from the static electric field. Our results indicate that while weak electric fields do not cause complete polarization of liquid water molecules, they can induce a transition to a more structured ice-like geometry of the water molecules at the ice-water interphase region, particularly when applied perpendicular to the ice-water interphase. Notably, the interface adjacent to cubic ice exhibits a greater response to the electric fields than that adjacent to hexagonal ice. This is attributable to the intrinsic differences in their original hydrogen bonding networks.

The Magnetic Field Freezes the Mercedes-Benz Water Model

In this study, we investigate the impact of magnetic fields on the structural and thermodynamic properties of water. To accomplish this, we employed the Mercedes-Benz (MB) model, a two-dimensional representation of water using Lennard-Jones disks with angle-dependent interactions that closely mimic hydrogen bond formation. We extended the MB model by introducing two charges to enable interaction with the magnetic field. Employing molecular dynamics simulations, we thoroughly explored the thermodynamic properties concerning various magnetic flux intensities. As a result, we observed that under a weak magnetic flux, the property of water remained unaltered, while a stronger flux astonishingly led to the freezing of water molecules. Furthermore, our study revealed that once a specific flux magnitude was reached, the density anomaly disappeared, and an increase in flux caused the MB particles to form a glassy state.

Electro-Freezing of Supercooled Water Is Induced by Hydrated Al3+ and Mg2+ Ions: Experimental and Theoretical Studies

This work reports that the octahedral hydrated Al3+ and Mg2+ ions operate within electrolytic cells as kosmotropic (long-range order-making) "ice makers" of supercooled water (SCW). 10-5 M solutions of hydrated Al3+ and Mg2+ ions each trigger, near the cathode (-20 ± 5 V), electro-freezing of SCW at -4 °C. The hydrated Al3+ ions do so with 100% efficiency, whereas the Mg2+ ions induce icing with 40% efficiency. In contrast, hydrated Na+ ions, under the same experimental conditions, do not induce icing differently than pure water. As such, our study shows that the role played by Al3+ and Mg2+ ions in water electro-freezing is impacted by two synchronous effects: (1) a geometric effect due to the octahedral packing of the coordinated water molecules around the metallic ions, and (2) the degree of polarization which these two ions induce and thereby acidify the coordinated water molecules, which in turn imparts them with an ice-like structure. Long-duration molecular dynamics (MD) simulations of the Al3+ and Mg2+ indeed reveal the formation of "ice-like" hexagons in the vicinity of these ions. Furthermore, the MD shows that these hexagons and the electric fields of the coordinate water molecules give rise to ultimate icing. As such, the MD simulations provide a rational explanation for the order-making properties of these ions during electro-freezing.

Molecular dynamics simulation of water condensation with nucleus under electromagnetic wave irradiation

The condensation process of water with different nuclei under electromagnetic wave irradiation was studied by molecular dynamics simulation. It was found that there is a different electric-field effect when the condensation nucleus was a small (NH₄)₂SO₄ cluster than a CaCO3 nucleus. Through the analysis of the hydrogen-bond number, energy change, and dynamic behavior, we found that the effect of external electric field on the condensation process mainly comes from the change of potential energy caused by the dielectric response and there is a competition effect between the dielectric response and the dissolution in the system with (NH₄)₂SO₄.

Electric-field induced entropic effects in liquid water

Externally applied electric fields in liquid water can induce a plethora of effects with wide implications in electrochemistry and hydrogen-based technologies. Although some effort has been made to elucidate the thermodynamics associated with the application of electric fields in aqueous systems, to the best of our knowledge, field-induced effects on the total and local entropy of bulk water have never been presented so far. Here, we report on classical TIP4P/2005 and ab initio molecular dynamics simulations measuring entropic contributions carried by diverse field intensities in liquid water at room temperature. We find that strong fields are capable of aligning large fractions of molecular dipoles. Nevertheless, the order-maker action of the field leads to quite modest entropy reductions in classical simulations. Albeit more significant variations are recorded during first-principles simulations, the associated entropy modifications are small compared to the entropy change involved in the freezing phenomenon, even at intense fields slightly beneath the molecular dissociation threshold. This finding further corroborates the idea that electrofreezing (i.e., the electric-field-induced crystallization) cannot take place in bulk water at room temperature. In addition, here, we propose a molecular-dynamics-based analysis (3D-2PT) that spatially resolves the local entropy and the number density of bulk water under an electric field, which enables us to map their field-induced changes in the environment of reference H2O molecules. By returning detailed spatial maps of the local order, the proposed approach is capable of establishing a link between entropic and structural modifications with atomistic resolution.

Simple rose model of water in constant electric field

A simple two-dimensional statistical mechanical water model, called the rose model, was used in this work. We studied how a homogeneous constant electric field affects the properties of water. The rose model is a very simple model that helps explain the anomalous properties of water. Rose water molecules are represented as two-dimensional Lennard-Jones disks with potentials for orientation-dependent pairwise interactions mimicking formations of hydrogen bonds. The original model is modified by addition of charges for interaction with the electric field. We studied what kind of influence the electric field strength has on the model's properties. To determine the structure and thermodynamics of the rose model under the influence of the electric field we used Monte Carlo simulations. Under the influence of a weak electric field the anomalous properties and phase transitions of the water do not change. On the other hand, the strong fields shift the phase transition points as well as the position of the density maximum.

Incorporating Polarization and Charge Transfer into a Point-Charge Model for Water Using Machine Learning

Rigid nonpolarizable water models with fixed point charges have been widely employed in molecular dynamics simulations due to their efficiency and reasonable accuracy for the potential energy surface. However, the dipole moment surface of water is not necessarily well-described by the same fixed charges, leading to failure in reproducing dipole-related properties. Here, we developed a machine-learning model trained against electronic structure data to assign point charges for water, and the resulting dipole moment surface significantly improved the predictions of the dielectric constant and the low-frequency IR spectrum of liquid water. Our analysis reveals that within our atom-centered point-charge description of the dipole moment surface, the intermolecular charge transfer is the major source of the peak intensity at 200 cm^-1, whereas the intramolecular polarization controls the enhancement of the dielectric constant. The effects of exact Hartree-Fock exchange in the hybrid density functional on these properties are also discussed.

The electric field changes the anomalous properties of the Mercedes Benz water model

The influence of a homogeneous constant electric field on water properties was assessed. We used a simple two-dimensional statistical mechanical model called the Mercedes-Benz (MB) model of water in the study. The MB water molecules are two-dimensional disks with Gaussian arms that mimic the formation of hydrogen bonds. The model is modified with added charges for interaction with the electric field. The influence of the strength of the electric field on the water's properties was studied using Monte Carlo simulations. The structure and thermodynamics of the water were determined as a function of the strength of the electric field. We observed that the properties and phase transitions of the water in the low strength electric field does not change. In contrast, the high strength electric field shifts boiling and melting points as well as the position of the density maxima. After further increasing the strength of the electric field the density anomaly disappears.

Gate Alignment of Liquid Water Molecules in Electric Double Layer

The behavior of liquid water molecules near an electrified interface is important to many disciplines of science and engineering. In this study, we applied an external gate potential to the silica/water interface via an electrolyte-insulator-semiconductor (EIS) junction to control the surface charging state. Without varying the ionic composition in water, the electrical gating allowed an efficient tuning of the interfacial charge density and field. Using the sum-frequency vibrational spectroscopy, we found a drastic enhancement of interfacial OH vibrational signals at high potential in weakly acidic water, which exceeded that from conventional bulk-silica/water interfaces even in strong basic solutions. Analysis of the spectra indicated that it was due to the alignment of liquid water molecules through the electric double layer, where the screening was weak because of the low ion density. Such a combination of strong field and weak screening demonstrates the unique tuning capability of the EIS scheme, and would allow us to investigate a wealth of phenomena at charged oxide/water interfaces.

Enhancement of the local asymmetry in the hydrogen bond network of liquid water by an ultrafast electric field pulse

Condensed phase electron decomposition analysis based on density functional theory has recently revealed an asymmetry in the hydrogen-bond network in liquid water, in the sense that a significant population of water molecules are simultaneously donating and accepting one strong hydrogen-bond and another substantially weaker one. Here we investigate this asymmetry, as well as broader structural and energetic features of water's hydrogen-bond network, following the application of an intense electric field square pulse that invokes the ultrafast reorientation of water molecules. We find that the necessary field-strength required to invoke an ultrafast alignment in a picosecond time window is on the order of 108 Vm-1. The resulting orientational anisotropy imposes an experimentally measurable signature on the structure and dynamics of the hydrogen-bond network, including its asymmetry, which is strongly enhanced. The dependence of the molecular reorientation dynamics on the field-strength can be understood by relating the magnitude of the water dipole-field interaction to the rotational kinetic energy, as well as the hydrogen-bond energy.

The Impact of the Electric Field on Surface Condensation of Water Vapor: Insight from Molecular Dynamics Simulation

In this study, molecular dynamics simulations were carried out to study the coupling effect of electric field strength and surface wettability on the condensation process of water vapor. Our results show that an electric field can rotate water molecules upward and restrict condensation. Formed clusters are stretched to become columns above the threshold strength of the field, causing the condensation rate to drop quickly. The enhancement of surface attraction force boosts the rearrangement of water molecules adjacent to the surface and exaggerates the threshold value for shape transformation. In addition, the contact area between clusters and the surface increases with increasing amounts of surface attraction force, which raises the condensation efficiency. Thus, the condensation rate of water vapor on a surface under an electric field is determined by competition between intermolecular forces from the electric field and the surface.

Viscous field-aligned water exhibits cubic-ice-like structural motifs

Strong electric fields are known to greatly accelerate the freezing of water in molecular dynamics simulations, and have also been shown to affect the thermodynamics of the phase transition. In this work, a mechanistic explanation for field-induced crystallization of water is presented. Due to the coupling between the rotational and the translational degrees of freedom of individual water molecules, an applied field can directly drive the formation of cubic-ice like local motifs in water. Analysis of the angular distributions of water molecules in TIP4P-2005 water at field strengths between 0.0 and 0.32 V Å-1 demonstrates the existence of such motifs in the field-aligned liquid phase that is observed prior to the onset of the freezing transition. The dynamic properties of this field-aligned liquid phase are also studied, and its viscosity is shown to be within a factor of two of that of regular liquid water using the Green-Kubo method as well as mean squared displacements. The choice between the NPT and the NVT ensembles is shown to have a strong impact on the evolution of molecular dynamics trajectories at field strengths close to the threshold for the freezing transition, and the importance of properly accounting for the electric field terms in the pressure virial is emphasized.

Phase boundaries, nucleation rates and speed of crystal growth of the water-to-ice transition under an electric field: a simulation study

We investigate with computer simulations the effect of applying an electric field on the water-to-ice transition. We use a combination of state-of-the-art simulation techniques to obtain phase boundaries and crystal growth rates (direct coexistence), nucleation rates (seeding) and interfacial free energies (seeding and mold integration). First, we consider ice Ih, the most stable polymorph in the absence of a field. Its normal melting temperature, speed of crystal growth and nucleation rate (for a given supercooling) diminish as the intensity of the field goes up. Then, we study polarised cubic ice, or ice Icf, the most stable solid phase under a strong electric field. Its normal melting point goes up with the field and, for a given supercooling, under the studied field (0.3 V nm^-1) ice Icf nucleates and grows at a similar rate as Ih with no field. The net effect of the field would then be that ice nucleates at warmer temperatures, but in the form of ice Icf. The main conclusion of this work is that reasonable electric fields (not strong enough to break water molecules apart) are not relevant in the context of homogeneous ice nucleation at 1 bar.

Polarizing Oxygen Vacancies in Insulating Metal Oxides under a High Electric Field

We demonstrate a thermodynamic formulation to quantify defect formation energetics in an insulator under a high electric field. As a model system, we analyzed neutral oxygen vacancies (color centers) in alkaline-earth-metal binary oxides using density functional theory, Berry phase calculations, and maximally localized Wannier functions. The work of polarization lowers the field-dependent electric Gibbs energy of formation of this defect. This is attributed mainly to the ease of polarizing the two electrons trapped in the vacant site, and secondarily to the defect induced reduction in bond stiffness and softening of phonon modes. The formulation and analysis have implications for understanding the behavior of insulating oxides in electronic, magnetic, catalytic, and electrocaloric devices under a high electric field.

Macroscopic electric field inside water-filled biological nanopores

Multi-drug resistance bacteria are a challenging problem of contemporary medicine. This is particularly critical for Gram-negative bacteria, where antibiotics are hindered by the outer membrane to reach internal targets. Here more polar antibiotics make use of nanometric water-filled channels to permeate inside. We present in this work a computational all-atom approach, using water as a probe, for the calculation of the macroscopic electric field inside water-filled channels. The method allows one to compare not only different systems but also the same system under different conditions, such as pH and ion concentration. This provides a detailed picture of electrostatics in biological nanopores shedding more light on how the charged residues of proteins determine the electric field inside, and also how medium can tune it. These details are central to unveil the filtering mechanism behind the permeation of small polar molecules through nanometric water-filled channels.

A Theoretical Study of the Hydration of Methane, from the Aqueous Solution to the sI Hydrate-Liquid Water-Gas Coexistence

Monte Carlo and molecular dynamics simulations were done with three recent water models TIP4P/2005 (Transferable Intermolecular Potential with 4 Points/2005), TIP4P/Ice (Transferable Intermolecular Potential with 4 Points/ Ice) and TIP4Q (Transferable Intermolecular Potential with 4 charges) combined with two models for methane: an all-atom one OPLS-AA (Optimal Parametrization for the Liquid State) and a united-atom one (UA); a correction for the C–O interaction was applied to the latter and used in a third set of simulations. The models were validated by comparison to experimental values of the free energy of hydration at 280, 300, 330 and 370 K, all under a pressure of 1 bar, and to the experimental radial distribution functions at 277, 283 and 291 K, under a pressure of 145 bar. Regardless of the combination rules used for σ_C,O, good agreement was found, except when the correction to the UA model was applied. Thus, further simulations of the sI hydrate were performed with the united-atom model to compare the thermal expansivity to the experiment. A final set of simulations was done with the UA methane model and the three water models, to study the sI hydrate-liquid water-gas coexistence at 80, 230 and 400 bar. The melting temperatures were compared to the experimental values. The results show the need to perform simulations with various different models to attain a reliable and robust molecular image of the systems of interest.

Perspectives on external electric fields in molecular simulation: progress, prospects and challenges

In this review, the application of a wide variety of external electric fields in molecular simulation shall be discussed, including time-varying and electromagnetic, as well as the utility and potential impact and prospects for exploitation of such simulations for real-world and industrial end use. In particular, non-equilibrium molecular dynamics will be discussed, as well as challenges in addressing adequate thermostatting and scaling field amplitudes to more experimentally relevant levels. Attention shall be devoted to recent progress and advances in external fields in ab initio molecular simulation and dynamics, as well as elusive challenges thereof (and, to some extent, for molecular dynamics from empirical potentials), such as timescales required to observe low-frequency and intensity field effects. The challenge of deterministic molecular dynamics in external fields in sampling phase space shall be discussed, along with prospects for application of fields in enhanced-sampling simulations. Finally, the application of external electric fields to a wide variety of aqueous, nanoscale and biological systems will be discussed, often motivated by the possibility of exploitation in real-world applications, which serve to underpin our molecular-level understanding of field effects in terms of microscopic mechanisms, and possibly with a view to control thereof.

Communication: Kinetic and pairing contributions in the dielectric spectra of electrolyte solutions

In the late 1970s, Hubbard and Onsager predicted that adding salt to a polar solution would result in a reduced dielectric permittivity that arises from the unexpected tendency of solvent dipoles to align opposite to the applied field. Here we develop a novel non-equilibrium molecular dynamics simulation approach to determine this decrement accurately. Using a thermodynamic consistent all-atom force field we show that for an aqueous solution containing sodium chloride around 4.8 mol/l, this effect accounts for 12% of the total dielectric permittivity. The dielectric decrement can be strikingly different if a less accurate force field for the ions is used. Using the widespread GROMOS parameters, we observe in fact an increment of the dielectric permittivity rather than a decrement, caused by ion pairing and introduced by a too low dispersion force.

Anomalies in bulk supercooled water at negative pressure

Water anomalies still defy explanation. In the supercooled liquid, many quantities, for example heat capacity and isothermal compressibility κT, show a large increase. The question arises if these quantities diverge, or if they go through a maximum. The answer is key to our understanding of water anomalies. However, it has remained elusive in experiments because crystallization always occurred before any extremum is reached. Here we report measurements of the sound velocity of water in a scarcely explored region of the phase diagram, where water is both supercooled and at negative pressure. We find several anomalies: maxima in the adiabatic compressibility and nonmonotonic density dependence of the sound velocity, in contrast with a standard extrapolation of the equation of state. This is reminiscent of the behavior of supercritical fluids. To support this interpretation, we have performed simulations with the 2005 revision of the transferable interaction potential with four points. Simulations and experiments are in near-quantitative agreement, suggesting the existence of a line of maxima in κT (LMκT). This LMκT could either be the thermodynamic consequence of the line of density maxima of water [Sastry S, Debenedetti PG, Sciortino F, Stanley HE (1996) Phys Rev E 53:6144-6154], or emanate from a critical point terminating a liquid-liquid transition [Sciortino F, Poole PH, Essmann U, Stanley HE (1997) Phys Rev E 55:727-737]. At positive pressure, the LMκT has escaped observation because it lies in the "no man's land" beyond the homogeneous crystallization line. We propose that the LMκT emerges from the no man's land at negative pressure.

Sudden switchover between the polyamorphic phase separation and the glass-to-liquid transition in glassy LiCl aqueous solutions

Lithium chloride aqueous solutions (LiClaq solutions) below 10 mol.% are vitrified by cooling from room temperature to 77 K at 0.3 GPa. We examine the solvent state of the glassy sample and its transformation by heating at 1 atm using low-temperature differential scanning calorimetry and Raman spectroscopy. This experimental study suggests strongly that the solvent state of the glassy LiClaq solution closely relates to the state of high-density amorphous ice. Moreover, we reconfirm that the separation into the low-density amorphous ice and the glassy highly concentrated LiClaq solution occurs in the glassy dilute LiClaq solution at ∼130 K, not the glass-to-liquid transition which is commonly observed in the glassy LiClaq solution above ∼10 mol.%. In order to interpret the sudden switchover between the glass-to-liquid transition and the phase separation at ∼10 mol.%, we propose a state diagram of LiClaq solution which connects with a polyamorphic state diagram of pure water and discuss a possibility that the electric field induces a polyamorphic transition of water.

Electromelting of confined monolayer ice

In sharp contrast to the prevailing view that electric fields promote water freezing, here we show by molecular dynamics simulations that monolayer ice confined between two parallel plates can melt into liquid water under a perpendicularly applied electric field. The melting temperature of the monolayer ice decreases with the increasing strength of the external field due to the field-induced disruption of the water-wall interaction induced well-ordered network of the hydrogen bond. This electromelting process should add an important new ingredient to the physics of water.

Ab initio molecular dynamics study of dissociation of water under an electric field

The behavior of liquid water under an electric field is a crucial phenomenon in science and engineering. However, its detailed description at a microscopic level is difficult to achieve experimentally. Here we report on the first ab initio molecular-dynamics study on water under an electric field. We observe that the hydrogen-bond length and the molecular orientation are significantly modified at low-to-moderate field intensities. Fields beyond a threshold of about 0.35  V/Å are able to dissociate molecules and sustain an ionic current via a series of correlated proton jumps. Upon applying even more intense fields (∼1.0  V/Å), a 15%-20% fraction of molecules are instantaneously dissociated and the resulting ionic flow yields a conductance of about 7.8  Ω-1 cm-1, in good agreement with experimental values. This result paves the way to quantum-accurate microscopic studies of the effect of electric fields on aqueous solutions and, thus, to massive applications of ab initio molecular dynamics in neurobiology, electrochemistry, and hydrogen economy.