The first term captures the rotational response of solvent molecules with density N 0 (r) ( Ñ0 is its Fourier transform) at temperature T , with ρ lm mol being the charge density distribution of a solvent molecule projected to each spherical harmonic component indexed by l, m. The second term captures the electronic polarization response of solvent molecules from atomic sites with densities N 0 α (⃗ r), each with susceptibilities of strength χ α and non-local weight function w α (r). All the above parameters are extracted from DFT calculations of solvent molecules, while the correlation strength functions C rot (ω) and C pol (ω) are constrained by the experimental rotational and electronic dielectric response of the solvent. See Ref. 31 for the derivation and complete parameterization of the SaLSA solvation model.