Specifically, the fitting quality was effectively unaltered for purely electrostatic perturbations (from R 2 = 0.95 to 0.94) but drastically improved in H-bonding environments (from R 2 < 0 to ~0.9). In this fit, the previously optimized VSE and empirical H-bonding parameters remain similar to those obtained in Fig. and Fig. with |Δ𝜇 ⃗| = 0.22