Molecular dynamics simulations at constant pressure and/or temperature

Journal of Chemical Physics - Tập 72 Số 4 - Trang 2384-2393 - 1980
Hans Christian Andersen1
1Department of Chemistry, Stanford University, Stanford, California 94305

Tóm tắt

In the molecular dynamics simulation method for fluids, the equations of motion for a collection of particles in a fixed volume are solved numerically. The energy, volume, and number of particles are constant for a particular simulation, and it is assumed that time averages of properties of the simulated fluid are equal to microcanonical ensemble averages of the same properties. In some situations, it is desirable to perform simulations of a fluid for particular values of temperature and/or pressure or under conditions in which the energy and volume of the fluid can fluctuate. This paper proposes and discusses three methods for performing molecular dynamics simulations under conditions of constant temperature and/or pressure, rather than constant energy and volume. For these three methods, it is shown that time averages of properties of the simulated fluid are equal to averages over the isoenthalpic–isobaric, canonical, and isothermal–isobaric ensembles. Each method is a way of describing the dynamics of a certain number of particles in a volume element of a fluid while taking into account the influence of surrounding particles in changing the energy and/or density of the simulated volume element. The influence of the surroundings is taken into account without introducing unwanted surface effects. Examples of situations where these methods may be useful are discussed.

Từ khóa


Tài liệu tham khảo

1971, Annu. Rev. Phys. Chem., 22, 563, 10.1146/annurev.pc.22.100171.003023

1973, Ann. Rev. Phys. Chem., 24, 159, 10.1146/annurev.pc.24.100173.001111

1975, Ann. Rev. Phys. Chem., 26, 145, 10.1146/annurev.pc.26.100175.001045

1976, Rev. Mod. Phys., 48, 587, 10.1103/RevModPhys.48.587

1977, Annu. Rev. Phys. Chem., 28, 373, 10.1146/annurev.pc.28.100177.002105

1978, Annu. Rev. Phys. Chem., 29, 441, 10.1146/annurev.pc.29.100178.002301

1967, Phys. Rev. Lett., 18, 988, 10.1103/PhysRevLett.18.988

1968, J. Chem. Phys., 49, 3609, 10.1063/1.1670641

1976, J. Chem. Soc. Faraday II, 72, 1667, 10.1039/f29767201667

1964, Phys. Rev., 136, A405, 10.1103/PhysRev.136.A405

1967, Phys. Rev., 159, 98, 10.1103/PhysRev.159.98

1968, Phys. Rev., 165, 201, 10.1103/PhysRev.165.201

1976, J. Chem. Phys., 64, 1564, 10.1063/1.432380

1978, Phys. Rev. Lett., 41, 1244, 10.1103/PhysRevLett.41.1244

1974, J. Chem. Phys., 60, 1545, 10.1063/1.1681229

1977, J. Am. Chem. Soc., 99, 7413, 10.1021/ja00465a002

1973, J. Chem. Phys., 59, 1323

1979, J. Am. Chem. Soc., 101, 1913, 10.1021/ja00502a001

1976, Proc. R. Soc. London A, 348, 485, 10.1098/rspa.1976.0051

1953, J. Chem. Phys., 21, 1087, 10.1063/1.1699114