Water’s Hidden Switch: How Ferroelectricity May Trigger Its Mysterious Phase Transitions

Capillary waves observed at the boundary between water and glycerol in microfluidic flows. Copyright: the authors, published in Phys. Rev. Lett. 134, 054001 (2025), DOI: 10.1103/PhysRevLett.134.054001, permission granted by the American Physical Society 2025.
Capillary Waves Found Between Fully Miscible Fluids
15th December 2025
Capillary waves observed at the boundary between water and glycerol in microfluidic flows. Copyright: the authors, published in Phys. Rev. Lett. 134, 054001 (2025), DOI: 10.1103/PhysRevLett.134.054001, permission granted by the American Physical Society 2025.
Capillary Waves Found Between Fully Miscible Fluids
15th December 2025

Water’s Hidden Switch: How Ferroelectricity May Trigger Its Mysterious Phase Transitions

Water’s unique thermodynamic anomalies, manifesting even under Earth’s ambient conditions, hold the key to understanding fundamental mechanisms driving life and Earth’s processes. Finding evidence nowadays in several molecular dynamics simulation studies, a compelling hypothesis for explaining the existence of water’s equilibrium anomalies is the so-called liquid-liquid phase transition (LLPT) between high-density and low-density liquid phases (HDL/LDL) in a supercooled, metastable liquid state. Its origin has been so-far not clarified. On the other hand, water is a polar liquid and, as such, can, in principle, undergo a ferroelectric phase transition, resulting in the setting of macroscopic polarization, under proper thermodynamic conditions.

Fig. 1: Sketch of the polarization-density phase diagram in the temperature (T)-pressure (p) thermodynamic plane obtained from classical density functional theory. Copyright: 2024 the Authors. Published by PNAS under Creative Commons Attribution-NonCommercial License 4.0 (CC BY-NC-ND).
Fig. 1: Sketch of the polarization-density phase diagram in the temperature (T)-pressure (p) thermodynamic plane obtained from classical density functional theory. Copyright: 2024 the Authors. Published by PNAS under Creative Commons Attribution-NonCommercial License 4.0 (CC BY-NC-ND).

A recent study lead by a researcher from SoftComp partner Ca’ Foscari University of Venice unveiled not only a link between ferroelectric and liquid-liquid phase transitions, but also a possible role of ferroelectricity in promoting the LLPT, by analysing extensive molecular dynamics simulations and developing a classical density functional theory. The theory treats water as a polar liquid. Grounded in the characteristics of the microscopic dipolar potential interaction and the liquid’s molecules annealed positional disorder, it leads to a free energy expression that supports phase transitions both ferroelectric and liquid-liquid. This research not only characterizes but can shed light on the origin of the LLPT and thermodynamic anomalies in water.

Significant questions arise

“By demonstrating that the overlooked dipolar degrees of freedom can actually play a leading role in the liquid-liquid phase transition in water, this research seeks to open thought-provoking perspectives”, says first author Maria Grazia Izzo. “Although the expression for free energy enabling phase transitions applies universally to polar liquids, its specific coefficients, which determine whether the transitions effectively occur under certain thermodynamic conditions, depend on the microscopic details. Our research thus raises significant questions: Is water ‘merely’ a polar liquid with the microscopic characteristics suitable to allow a ferroelectric phase transition close to Earth’s environmental conditions? And what contribute to shaping its microscopic characteristics?”

Read more: M.G. Izzo, J. Russo, G. Pastore, PNAS 121 (47) e2412456121 (2024) https://doi.org/10.1073/pnas.2412456121

SoftComp partners: Ca’ Foscari University of Venice

Images Copyright: 2024 the Authors. Published by PNAS under Creative Commons Attribution-NonCommercial License 4.0 (CC BY-NC-ND).

Fig. 2: A) Polarization Binder cumulant, U L , average density, , dielectric constant, ε 0 , and isothermal compressibility, K T , at different points of the temperature (T)-pressure (p) thermodynamic plane obtained from MD simulations with TIP4P/Ice water model. Colour scale represents the quantities value, symbol size is proportional to the error, obtained via block averaging. The solid and dashed black lines are visual guides, marking the first- order LLPT line and the Widom line (WL), respectively. The diamond symbol marks the critical point. B) Same quantities as functions of T along isobaric lines intersecting the WL ( 0 bar, 1000 bar, Left) and the first-order LLPT line ( 2500 bar, Right). Copyright: 2024 the Authors. Published by PNAS under Creative Commons Attribution-NonCommercial License 4.0 (CC BY-NC-ND).
Fig. 2: A) Polarization Binder cumulant, U L , average density, <ρ>, dielectric constant, ε 0 , and isothermal compressibility, K T , at different points of the temperature (T)-pressure (p) thermodynamic plane obtained from MD simulations with TIP4P/Ice water model. Colour scale represents the quantities value, symbol size is proportional to the error, obtained via block averaging. The solid and dashed black lines are visual guides, marking the first- order LLPT line and the Widom line (WL), respectively. The diamond symbol marks the critical point. B) Same quantities as functions of T along isobaric lines intersecting the WL ( 0 bar, 1000 bar, Left) and the first-order LLPT line ( 2500 bar, Right). Copyright: 2024 the Authors. Published by PNAS under Creative Commons Attribution-NonCommercial License 4.0 (CC BY-NC-ND).
Research Gate
Research Gate