Structural Origin of Water's Anomalous Heat Capacity ✛ ✦
Context
Heat capacity is one of the key properties of matter, defining its ability to store energy. For most liquids, it is low and weakly dependent on temperature, yet water presents a striking exception: at room conditions, its specific heat capacity is almost double that of ethanol or hexane, and upon deep supercooling (−40°C), a sharp maximum appears. This anomaly is fundamental for climate, geophysics, and biochemistry: oceans act as a giant thermal buffer, and cellular processes are tuned to the unique thermal properties of the aqueous environment. From a physics standpoint, the anomalous behavior points to hidden degrees of freedom associated with rearrangements of the hydrogen-bond network and possibly with the existence of two distinct liquid phases. Classical theories long failed to quantitatively explain the observed picture, especially the role of quantum effects in nuclear motion (quantum thermodynamics) and their connection to entropic fluctuations.
Methods
To separate contributions, molecular dynamics simulations were used in both classical and quantum versions (the path-integral method, dating back to Richard Feynman). Interatomic interactions were described by neural network potentials of a new generation (NEP), trained on ultra-accurate quantum-chemical calculations at the MB-pol level and the hybrid functional revPBE0-D3. Such potentials combine first-principles accuracy with the speed of empirical force fields, enabling simulations of thousands of H₂O molecules (comprising hydrogen and oxygen) over times up to tens of nanoseconds. Nuclear quantum effects were accounted for by "smearing" atoms into ring polymers of 32 beads, which is equivalent to incorporating the uncertainty principle and the thermal de Broglie wavelength. To estimate heat capacity from first principles, the vibrational density of states was computed and the Planck–Einstein quantum reweighing of occupancies was applied, subtracting the classical contribution governed by the Boltzmann distribution.
Results
Classical molecular dynamics overestimates heat capacity by 4–5 gas constants (R) compared to experiment, whereas quantum corrections (path-integral or spectral reweighing) almost completely eliminate the discrepancy. Importantly, both classical and quantum approaches reproduce the heat capacity peak near 230 K—the temperature of the CP maximum in experiments for water in nanopores. A key role is played by the structural order parameter SSSI(3), which counts "intruders" into the second coordination shell. Its temperature derivative practically coincides with CP, and its variance sharply increases in the peak region, indicating intense fluctuations between low-density (LDL) and high-density (HDL) local structures. A two-level model yields an effective transition enthalpy of 2.8–3.3 kJ/mol (classical) and 3.6–3.9 kJ/mol (with Planck quanta), which is comparable to the energy difference between amorphous ices. Decomposing heat capacity into baseline and structural components shows that it is precisely the LDL↔HDL interconversion that creates both the anomalously high plateau at room temperature and the sharp maximum upon supercooling.
Implications
The work provides, for the first time, a microscopic explanation of water's anomalous heat capacity, linking it to interconversions between local structures of ordered and disordered types. This supports the scenario of a liquid–liquid critical point in supercooled water and indicates that quantum effects are not the direct cause of the anomaly but merely set the baseline level of heat capacity. Such a separation of contributions paves the way for more accurate water models for climate modeling and materials science.
Future development
Future plans include investigating isotopic effects (replacing hydrogen with deuterium), which noticeably shift the temperature of the CP maximum, clarifying the role of quantum effects. It is also promising to apply the approach to other anomalous liquids (silica, liquid phosphorus), where a liquid–liquid phase transition is suspected.
Impact
The results have direct significance for climatology (ocean heat uptake), cryobiology (cell resistance to supercooling), and the chemical physics of aqueous solutions.
Next steps
The next step will involve simulations on a finer temperature grid to accurately determine critical parameters, as well as experiments with ultrafast calorimetry to test the model's predictions.
Key open problems
The research is directly connected to one of the major unsolved problems in condensed matter physics—the nature of the liquid–liquid phase transition in water. The hypothesis of a second critical point in deeply supercooled water at negative pressures remains unconfirmed experimentally, and this work provides computational evidence for its reality.
The knowledge machine's view
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
Using molecular dynamics with neural network potentials and accounting for quantum nuclear effects, the work shows that water's anomalous heat capacity is determined by structural fluctuations between configurations of low and high local density, and the heat capacity peak upon supercooling is linked to a sharp increase in the variance of the order parameter SSSI(3).
Where its strength is
The classical and quantum contributions to heat capacity were quantitatively separated, and a direct connection was established between the macroscopic response and microscopic rearrangements of the second coordination shell, providing a structural order parameter as a convenient tool for predicting thermodynamics.
Where it could go next
- It might be interesting to compare the observed fluctuations between LDL- and HDL-like configurations with the dynamics of hydrogen bonds, which, as shown in a neighboring work, has a characteristic collective rhythm and could serve as a temporal marker for structural transitions.based on: Water dances: the hidden rhythm of hydrogen bonds arXiv:2607.06780
- Nearby lies a study on how quantum chaos in many-body systems leads to thermalization; the formalism proposed there could help more rigorously justify the link between quantum nuclear fluctuations and the growth of SSSI(3) variance near the heat capacity peak.based on: How quantum chaos gives birth to thermal equilibrium arXiv:2606.28660
Worth trying
- Test whether the lifetime of hydrogen bonds correlates with the relaxation time of SSSI(3) fluctuations in the temperature peak region.
- Assess whether the temperature evolution of the SSSI(3) variance can be described through the eigenstate statistics of an effective bond network Hamiltonian, using the approach from the work on quantum chaos.
Why this may matter
If the structural order parameter is indeed universally linked to heat capacity, this opens the way to simple engineering models for the thermophysical properties of aqueous solutions and to understanding the role of hydrogen bonds in thermoregulation of biological systems.
Nearby works in our archive
- How quantum chaos gives birth to thermal equilibriumsimilarity 0.65
- Water dances: the hidden rhythm of hydrogen bondssimilarity 0.63
- Thermal Split: How Uneven Heating Divides a Mixturesimilarity 0.62
- Living Things Glow with Heat, Like the Sunsimilarity 0.62
- Waves That Remember Left from Rightsimilarity 0.61
- Topology Holds Memorysimilarity 0.61
The closest work in the archive sits at 0.65 by meaning — the usual level for our corpus.
Original abstract
as written by the authors arXiv:2607.08957
Water isobaric heat capacity is anomalously large under ambient conditions and exhibits a sharp maximum upon supercooling. Using classical and path-integral molecular dynamics with accurate machine-learning interatomic potentials, we show that nuclear quantum effects primarily act by suppressing high-frequency vibrations, while the anomalous temperature dependence of the isobaric heat capacity originates from structural fluctuations, quantified by the second-solvent-shell intruder order parameter. A simple two-state mapping reveals an effective enthalpy scale of about 4 kJ/mol associated with the interconversion of low- and high-density-like local structures, providing a microscopic link between their population changes and the excess heat capacity from supercooled to ambient conditions.
Water has one of the highest specific heat capacities among all known liquids—4.18 J/(g·K) at 25°C. Thanks to this property, the oceans have absorbed more than 90% of the excess heat from global warming over the past half-century.
Key numbers
- specific heat capacity of water: 4.18 J/(g·K) at 25°C
- specific heat capacity of methanol (for comparison): 2.56 J/(g·K) at 25°C
- temperature of CP maximum: ~230 K (−43°C)
- effective enthalpy of LDL→HDL transition: 3–4 kJ/mol
- number of beads in PIMD quantum simulation: 32