This work presents new parallelizable numerical schemes for the integration of Dissipative Particle Dynamics with Energy conservation (DPDE). So far, no numerical scheme introduced in the literature is able to correctly preserve the energy over long
times and give rise to small errors on average properties for moderately small timesteps, while being straightforwardly parallelizable. We present in this article two new methods, both straightforwardly parallelizable, allowing to correctly preserve the total energy of the system. We illustrate the accuracy and performance of these new schemes both on equilibrium and nonequilibrium parallel simulations.
