The region not occupied by the solute, i.e., Ω C = R 3 \ Ω, is considered to be occupied by the bulk solvent. Note that without many changes, the method can also be applied if x j do not necessarily represent nuclear coordinates or just a subset of those (such as is needed in coarse-grained models). The charge density of the solute is denoted by ρ.