angular distribution function - periodic boundary conditions

[bananenpampe]

I am computing the angular distribution function for a system with and without periodic boundary conditions. I am under the impression, that when periodic boundary conditions are enforced, there are artifacts introduced for small angles.

frames = ase.io.read("test.xyz",index="::5")
for frame in frames:
    frame.wrap(eps=1e-12)

mol = MultiMolecule.from_ase(frames)
adf = mol.init_adf(r_max="3.73", weight=None)

# or  adf = mol.init_adf(r_max="3.73", weight=None, periodic="xyz")