When we simulate a bulk molecular system using molecular dynamics in (say) NVE ensemble (microcanonical), we apply some periodic boundary conditions (PBC) to avoid surface effects. What effect does this have on quantities like Pair Distribution function, since it so happens that part of a molecule is near one edge of the box and the other is in a different edge of the box. What effect does this have on PDF, MSD and other such thermodynamic quantities ?
PS : I am currently simulating SPC/E water molecules using cubic box, in some frames I find water molecules which are split apart, but connected by PBC.