A Unified Statistical-Thermodynamic Model of the Density and Compressibility Anomalies of Liquid Water

This work presents an analytical statistical-thermodynamic modeling framework designed to elucidate the physical mechanisms underlying the macroscopic anomalies of liquid water. By representing water as a cooperative, non-ideal two-state statistical mixture of Low-Density (LDS) and High-Density (HDS) structures, we utilize the physical principle of shared-bond coordination (accounting for the fact that each hydrogen bond is shared by two neighboring molecules) to establish a direct, self-consistent thermodynamic link between the macroscopic transition enthalpy and the single-bond dissociation energy (ΔH_trans = 0.5 · E_eff). Through a first-order Taylor-Newton expansion of the governing transcendental equation for the Temperature of Maximum Density (TMD) locus, we derive an explicit, closed-form quadratic relation for the TMD curve: T_max(P) ≈ (3/2) · T_c(P) - B_effective · T_c(P)². This simplified formulation provides an analytical explanation for the negative slope of the TMD line and reproduces the qualitative topological features of the reentrant, "nose-shaped" TMD curve under negative pressures (tension) as proposed by Speedy (1982). Additionally, by mapping this formulation onto the Sastry thermodynamic consistency relation (1996), we show how the negative slope of the TMD line is mathematically linked to the anomalous temperature dependence of the isothermal compressibility ((∂κ_T / ∂T)_P < 0 along the TMD locus). Rather than presenting new empirical anomalies, this model offers a transparent, unified mathematical explanation showing that water's density and compressibility anomalies are co-dependent manifestations of the same underlying structural transition. Key Physical & Mathematical Highlights of the Model: Microscopic-Macroscopic Bridge: We resolve the parameter-space redundancy of two-state models by physically locking the transition enthalpy to the single hydrogen-bond breaking energy: ΔH_trans(P) = 0.5 · E_eff(P). First Closed-Form TMD Equation: Instead of relying on multi-parameter transcendental computer fits, this model analytically collapses the equation of state specifically around the structural coexistence line (T_c), where the structural fluctuation susceptibility reaches its maximum. Parameter Dimensionality Reduction: The volume coefficient B = 2·R·β_0 / (ΔV_s·ΔS_trans) is explicitly constrained by water's experimental molar volume (V ≈ 18 cm³/mol) and background thermal expansivity (α ≈ 2×10⁻⁴ K⁻¹), yielding a baseline molar expansivity of β_0 ≈ V·α ≈ 3.6×10⁻⁹ m³/(mol·K). This eliminates arbitrary parameter-fitting. The Density-Compressibility Link: Using exact Maxwell relations along the TMD line (where thermal expansion α = 0), we mathematically prove that the negative slope of the TMD line (dT_max/dP < 0) dictates the anomalous temperature dependence of the isothermal compressibility: (∂κ_T / ∂T)_P < 0 on the TMD line.