I'm noticing that the code crashes pretty early in a simulation if I have multiple cores across a dimension having periodic boundary conditions. Is this a known issue and does anyone have an idea why this tends to happen? It seems that this only occurs with multisphere particles but not with regular spheres.