It is evident from the big-system C-O RDF (bottom plot, black line) that the majority of the correlations have died out already at 14 Å. Thus, the absence of convergence to 1 for the small-system RDFs is thus not due to the solvent shell-structure extending far beyond the simulation box size of the small system, but rather, due to the issue of normalizing to the density using the full box volume instead of the effective volume, as discussed above.