While local and semilocal functionals fail to capture intermediate to long-range vdW forces 122 , AIMD simulations have demonstrated that GGAs enhanced with vdW representations typically lead to a softer structure of bulk water, accompanied by increased mobility that aligns more closely with experimental measurements . This improvement is achieved by incorporating dispersion-corrected atom-centered potentials (DCACPs) , empirical dispersion corrections (e.g., Grimme's D2 127 and D3 128 ), or non-local correlation terms (e.g., (r)VV10 , vdW-DF 132 , TS-vdW ). However, the performance of such corrections relies on the original GGA to which the combination may not always improve, or may even deteriorate properties . Other studies pointed out the necessity of including a fraction of exact exchange, thus resorting to rung-4 hybrids, to effectively describe hydrogen bonding but without reaching a perfect agreement with the experiment .