Dissipative particle dynamics with energy conservation: Modelling of heat flow

Abstract

The stochastic differential equations for a model of dissipative particle dynamics, with both total energy and total momentum conservation at every time step, are presented. The algorithm satisfies detailed balance as well as the fluctuation–dissipation theorems that ensure that the proper thermodynamic equilibrium can be reached. The Fokker–Planck equation for the evolution of the probability distribution for the system is given. Macroscopic equilibrium probability distributions as well as equations of state for the model are also derived, and an appropriate definition of the free energy of the system is consistently proposed. Several simulations, results of equilibrium as well as transport properties, including heat transport and thermal convection in a box, are shown as proof of the internal consistency of the model.

References

1. P. J. Hoogerbrugge and J. M. V. A. Koelman, Europhys. Lett., 1992, 19, 155.
2. J. J. Monhagan, Comput. Phys. Rep., 1985, 3, 71.
3. E. S. Boek, P. V. Coveney, H. N. W. Lekkerkerker and P. van der Schoot, Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top., 1997, 55, 3124.
4. R. D. Groot and T. J. Madden, J. Chem. Phys., 1998, 108, 8713.
5. U. Frisch, B. Hasslacher and Y. Pomeau, Phys. Rev. Lett., 1986, 56, 1505.
6. A. J. C. Ladd, J. Fluid. Mech., 1994, 271, 285.
7. C. Marsh, G. Backx and M. H. Ernst, Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top., 1997, 56, 1676.
8. J. B. Avalos and A. D. Mackie, Europhys. Lett., 1997, 40, 141.
9. P. Español and P. Warren, Europhys. Lett., 1995, 30, 191.
10. H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, John Wiley & Sons, New York, 1985.
11. S. R. de Groot and P. Mazur, in Non-equilibrium Thermodynamics, Dover, New York, 1984.
12. W. van Saarloos, D. Bedeaux and P. Mazur, Physica A (Amsterdam), 1982, 110, 147.
13. J. B. Avalos and I. Pagonabarraga, Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top., 1995, 52, 5881.
14. N. G. van Kampen, Stochastic Processes in Physics and Chemistry, North Holland, Amsterdam, 2nd edn. 1992.
15. L. D. Landau and E. Lifshitz, Fluid Mechanics, Pergamon, New York, 1984.
16. P. Mazur and D. Bedeaux, Physica A (Amsterdam), 1991, 173, 155.
17. J. B. Avalos, Ph.D. thesis, University of Barcelona, 1992; I. Pagonabarraga and J. B. Avalos, unpublished results..
18. I. Pagonabarraga, M. H. J. Hagen and D. Frenkel, Europhys. Lett., 1998, 42, 377.
19. The relation between the equations of motion for a given system and the detailed balance condition can be found, for instance, in Section V.6 of ref. 16.
20. J.-P. Hansen and I. R. McDonald, Theory of Simple Liquids, Academic Press, London, 1986.
21. The requirement that the heat capacity of the particle must be larger than k is due to the fact that the mesoscopic random heat flow is modelled as a continuum process subjet to Gaussian statistics. Thus, the probability that a given particle, i say, at a given instant of time randomly loses more energy than the actual value u_i is exponentially reduced if ø > k. This is due to the fact that the characteristic particle internal energy is of the order of øT while fluctuations in the same variable are of the order of kT. Therefore, in general, the model is only able to deal with “classical” or “high temperature” situations, where the DPD particles have an average internal energy substantially larger than the thermal energy kT.