Convergence behaviour of solvation shells in simulated liquids

A convenient way to analyse solvent structure around a solute is to use solvation shells, whereby solvent position around the solute is discretised by the size of a solvent molecule, leading to multiple shells around the solute. The two main ways to define multiple shells around a solute are either directly with respect to the solute, called solute-centric, or locally for both solute and solvent molecules alike. It might be assumed that both methods lead to solvation shells with similar properties. However, our analysis suggests otherwise. Solvation shells are analysed in a series of simulations of five pure liquids of differing polarity. Shells are defined locally working outwards from each molecule treated as a reference molecule using two methods: the cutoff at the first minimum in the radial distribution function and the parameter-free Relative Angular Distance method (RAD). The molecular properties studied are potential energy, coordination number and coordination radius. Rather than converging to bulk values, as might be expected for pure solvents, properties are found to deviate as a function of shell index. This behaviour occurs because molecules with larger coordination numbers and radius have more neighbours, which make them more likely to be connected to the reference molecule via fewer shells. The effect is amplified for RAD because of its more variable coordination radii and for water with its more open structure and stronger interactions. These findings indicate that locally defined shells should not be thought of as directly comparable to solute-centric shells or to distance. As well as showing how box size and cutoff affect the non-convergence, to restore convergence we propose a hybrid method by defining a new set of shells with boundaries at the uppermost distance of each locally derived shell.