Liquid water is a ubiquitous and essential component of life, playing a critical role in a wide variety of chemical and biological processes . A comprehensive understanding of water at the atomic scale is vital for advancing research in diverse domains such as aqueous chemistry , biochemistry , atmospheric science , and environmental engineering . Furthermore, unraveling the intricate behavior of water molecules enables deeper insights into solvation dynamics , water-materials in-scale phenomena in water through experiments like NMR , IR , X-ray or neutron spectroscopy for which measurement interpretations often rely on theoretical models. Although a variety of computationally efficient and relatively accurate empirical force fields have been developed , those remain intrinsically incapable of describing bond breaking in chemical reactions. Therefore, the quantitative understanding of condensed phase water, in particular its reactivity, and role as a universal solvent can only fully emerge from the development of accurate ab initio molecular dynamics (AIMD) simulations . These simulations need to faithfully represent both electronic reorganization and nuclear quantum effects (NQEs) associated with hydrogen bonding but, at present, such a comprehensive predictive quantum picture at ambient conditions remains quite elusive. In addition, the cost of most accurate wavefunction-based approaches such as post-Hartree-Fock (e.g., MP2 or RPA ), coupled cluster (CCSD(T)) , or configuration interaction (CI) hinders their potential application across the entire water phase diagram.