implemented OpenMP, parallel Mersenne Twister RNG, phonon trap avoidance