A nonequilibrium molecular dynamics method for calculating the thermal conductivity is presented. It reverses the usual cause and effect picture. The “effect,” the heat flux, is imposed on the system and the “cause,” the temperature gradient is obtained from the simulation. Besides being very simple to implement, the scheme offers several advantages such as compatibility with periodic boundary conditions, conservation of total energy and total linear momentum, and the sampling of a rapidly converging quantity (temperature gradient) rather than a slowly converging one (heat flux). The scheme is tested on the Lennard-Jones fluid.
REFERENCES
1.
D. J. Evans and G. P. Morriss, Statistical Mechanics of Nonequilibrium Liquids (Academic, London, 1990).
2.
H. J. C. Berendsen, in Computer Simulation in Materials Science, edited by M. Meyer and V. Pontikis (Kluwer, Dordrecht, 1991).
3.
R.
Vogelsang
, C.
Hoheisel
, and M.
Luckas
, Mol. Phys.
64
, 1203
(1988
).4.
M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Oxford Science, Oxford, 1987).
5.
W. R. P. Scott (private communication).
6.
M. E.
Tuckerman
, B. J.
Berne
, and G. J.
Martyna
, J. Chem. Phys.
94
, 6811
(1991
).7.
H. J. C.
Berendsen
, J. P. M.
Postma
, W. F.
van Gunsteren
, A.
DiNola
, and J. R.
Haak
, J. Chem. Phys.
81
, 3684
(1984
).8.
9.
R. Sok, Internal Report, Groningen University, 1989;
H. J. C. Berendsen, in Computer Simulation in Materials Science, edited by M. Meyer and V. Pontikis (Kluwer, Dordrecht, 1991), p. 139.
10.
Handbook of Chemistry and Physics, 74th ed., edited by D. R. Lide (Chemical Rubber, Boca Raton, 1993).
11.
W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, Numerical Recipes (Cambridge University, New York, 1986).
12.
This content is only available via PDF.
© 1997 American Institute of Physics.
1997
American Institute of Physics
You do not currently have access to this content.
