The positional correction potential to eliminate density artifacts is constructed along the unitless relative distance a for the corresponding system geometry (Eqs. ( ) and ( )). The distance coordinate is discretized using a one-dimensional grid of points k = 1...K. Each grid point is the origin of a kernel function G (Eq. 15). The corresponding potential energy function U ρ is applied to the reference point R R R I of each solvent molecule I and reads