Water is routinely exposed to external electric fields. Whether, for example, at physiological conditions, in contact with biological systems, or at the interface of polar surfaces in countless technological settings, water responds to fields on the order of a few V Å⁻¹ in a manner that is under intense investigation. Dating back to the 19th century, the possibility of solidifying water upon applying electric fields – a process known as electrofreezing – is an alluring promise that has canalized major efforts since, with uncertain outcomes. Here, we perform long (up to 500 ps per field strength) ab initio molecular dynamics simulations of water at ambient conditions under external electric fields. We show that fields of 0.10 − 0.15 V Å⁻¹ induce electrofreezing to a ferroelectric amorphous phase which we term f-GW (ferroelectric glassy water). The transition occurs after ~ 150 ps for a field of 0.15 V Å⁻¹ and after ~ 200 ps for a field of 0.10 V Å⁻¹ and is signaled by a structural and dynamic arrest and the suppression of the fluctuations of the hydrogen bond network. Our work reports evidence of electrofreezing of bulk liquid water at ambient conditions and therefore impacts several fields, from fundamental chemical physics to biology and catalysis.

We investigate the large-scale structure of amorphous ices and transitions between their different forms by quantifying their large-scale density fluctuations. Specifically, we simulate the isothermal compression of low-density amorphous ice (LDA) and hexagonal ice to produce high-density amorphous ice (HDA). Both HDA and LDA are nearly hyperuniform; i.e., they are characterized by an anomalous suppression of large-scale density fluctuations. By contrast, in correspondence with the nonequilibrium phase transitions to HDA, the presence of structural heterogeneities strongly suppresses the hyperuniformity and the system becomes hyposurficial (devoid of “surface-area fluctuations”). Our investigation challenges the largely accepted “frozen-liquid” picture, which views glasses as structurally arrested liquids. Beyond implications for water, our findings enrich our understanding of pressure-induced structural transformations in glasses.

We investigate the microscopic origin of water’s anomalies by inspecting the hydrogen bond network (HBN) and the spatial organization of low-density-liquid (LDL) like and high-density-liquid (HDL) like environments. Specifically, we simulate—via classical molecular dynamics simulations—the isobaric cooling of a sample composed of 512 water molecules from ambient to deeply undercooled conditions at three pressures, namely, 1 bar, 400 bars, and 1000 bars. In correspondence with the Widom line (WL), (i) the HDL-like dominating cluster undergoes fragmentation caused by the percolation of LDL-like aggregates following a spinodal-like kinetics; (ii) such fragmentation always occurs at a “critical” concentration of ∼20%–30% in LDL; (iii) the HBN within LDL-like environments is characterized by an equal number of pentagonal and hexagonal rings that create a state of maximal frustration between a configuration that promotes crystallization (hexagonal ring) and a configuration that hinders it (pentagonal ring); (iv) the spatial organization of HDL-like environments shows a marked variation. Moreover, the inspection of the global symmetry shows that the intermediate-range order decreases in correspondence with the WL and such a decrease becomes more pronounced upon increasing the pressure, hence supporting the hypothesis of a liquid-liquid critical point. Our results reveal and rationalize the complex microscopic origin of water’s anomalies as the cooperative effect of several factors acting synergistically. Beyond implications for water, our findings may be extended to other materials displaying anomalous behaviours.