Solvation Lies Within: Simulating Condensed-Phase Properties from Local Electronic Structures

In transitions between different environmental settings, a molecular system inevitably undergoes a range of detectable changes, and the ability to accurately simulate such responses, e.g., in the form of shifts to molecular energies, remains an important challenge across physical chemistry. Based on an exact decomposition of total energies from Kohn–Sham density functional theory in a basis of spatially localized molecular orbitals, the present work outlines a robust protocol for sampling the effect of solvation within homogeneous condensed phases by focusing on perturbations to local electronic structures only. We report chemically intuitive results for binding energies of water, ethanol, and acetonitrile that all display fast convergence with respect to the bulk size. Results stay largely invariant with respect to the choice of basis set while reflecting differences in density functional approximations, and our protocol thus allows for a physically sound and efficient estimation of general effects related to bulk solvation.