Most importantly, the nonlocality of the SaLSA model leads to χ lq GG ′ → 0 as G → ∞ because the nonlocal weight functions ρlm mol (G) and wα (G) → 0 as G → ∞ (Figure ). This allows the liquid susceptibility to be converged with a small plane wave cutoff, in comparison to the χ lq GG ′ ∝ G • G ′ of local models (Figure ), that require additional numerical care in constructing and inverting ε GG ′ . We implemented calculation of χ lq GG ′ (q; ω) in the plane-wave basis in the open-source JDFTx software, exported directly as a matrix in the plane-wave basis for each q and ω, which can then be included in the construction and inversion of ε GG ′ (q; ω) already implemented within the open-source BerkeleyGW software using the static subspace approximation. Although the calculations in this work only used a single k/q-point, we emphasize that this formalism is general and can be used to predict electronic structure properties of solvated surfaces as well.