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ₜrans = 0. 5 · Eₑff). 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ₘax (P) ≈ (3/2) · Tc (P) - Bₑffective · Tc (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ₜrans (P) = 0. 5 · Eₑff (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 (Tc), where the structural fluctuation susceptibility reaches its maximum. Parameter Dimensionality Reduction: The volume coefficient B = 2·R·β₀ / (ΔVₛ·ΔSₜrans) 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 β₀ ≈ 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ₘax/dP < 0) dictates the anomalous temperature dependence of the isothermal compressibility: (∂κT / ∂T) P < 0 on the TMD line.
Vakhtang Mchedlishvili (Thu,) studied this question.