All cores in a parallel run with particle-matter-interaction use same sequence of random numbers
This is similar to the case we had in the Distribution class where the sequences of random numbers were the same. There we could alleviate this problem using two different approaches. The first approach is to discard parts of the sequence on all cores except of one (see here). This doesn't scale well with increasing number of cores and particles. The other approach was to use a different seed on each core (see here). This scales well but yields different results when different number of cores are used.
In particle-matter-interaction the first approach doesn't seem feasible. The second approach is however easy to implement.