Constant temperature molecular dynamics calculation on Lennard-Jones fluid and its application to water

Hideki Tanaka, Koichiro Nakanishi, Nobuatsu Watanabe

Research output: Contribution to journalArticle

48 Citations (Scopus)

Abstract

Constant temperature molecular dynamics calculation has been carried out on Lennard-Jones liquid simulating argon. This method, proposed recently by Andersen, intends to transform the system from microcanonical to canonical ensemble and to keep the temperature of the system a constant value by the generation of random velocities when molecules exchange their energies with heat reservoir with a certain probability. A simple scheme is given for the estimation of collision probability and the effect of introducing this probability on dynamic behavior is examined in detail for 108 argon atoms as a test simulation. It is established that a collision probability of 0.01 is sufficient to realize the constancy of temperature reasonably well with no appreciable disturbance in the dynamic behavior. A model of pure water with Matsuoka-Clementi-Yoshimine (MCY) potential has also been simulated in the same manner. In the case of 2/16 MCY water at 298.15 K, the temperature difference is only 0.73 K with a collision probability of 0.005. Various static properties of MCY water have been calculated with reasonable agreement with those by previous Monte Carlo calculation, and moreover, the dynamic behavior calculated for MCY water gives a satisfactory picture on both translational and rotational motions in water, including reasonable agreement of diffusion coefficient with experimental datum.

Original languageEnglish
Pages (from-to)2626-2634
Number of pages9
JournalThe Journal of Chemical Physics
Volume78
Issue number5
Publication statusPublished - 1983
Externally publishedYes

Fingerprint

Molecular dynamics
molecular dynamics
Fluids
Water
fluids
water
Argon
collisions
Temperature
temperature
argon
translational motion
temperature gradients
disturbances
diffusion coefficient
energy transfer
heat
Atoms
Molecules
Liquids

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

Constant temperature molecular dynamics calculation on Lennard-Jones fluid and its application to water. / Tanaka, Hideki; Nakanishi, Koichiro; Watanabe, Nobuatsu.

In: The Journal of Chemical Physics, Vol. 78, No. 5, 1983, p. 2626-2634.

Research output: Contribution to journalArticle

@article{c2922709521f47328045e535442273b6,
title = "Constant temperature molecular dynamics calculation on Lennard-Jones fluid and its application to water",
abstract = "Constant temperature molecular dynamics calculation has been carried out on Lennard-Jones liquid simulating argon. This method, proposed recently by Andersen, intends to transform the system from microcanonical to canonical ensemble and to keep the temperature of the system a constant value by the generation of random velocities when molecules exchange their energies with heat reservoir with a certain probability. A simple scheme is given for the estimation of collision probability and the effect of introducing this probability on dynamic behavior is examined in detail for 108 argon atoms as a test simulation. It is established that a collision probability of 0.01 is sufficient to realize the constancy of temperature reasonably well with no appreciable disturbance in the dynamic behavior. A model of pure water with Matsuoka-Clementi-Yoshimine (MCY) potential has also been simulated in the same manner. In the case of 2/16 MCY water at 298.15 K, the temperature difference is only 0.73 K with a collision probability of 0.005. Various static properties of MCY water have been calculated with reasonable agreement with those by previous Monte Carlo calculation, and moreover, the dynamic behavior calculated for MCY water gives a satisfactory picture on both translational and rotational motions in water, including reasonable agreement of diffusion coefficient with experimental datum.",
author = "Hideki Tanaka and Koichiro Nakanishi and Nobuatsu Watanabe",
year = "1983",
language = "English",
volume = "78",
pages = "2626--2634",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "5",

}

TY - JOUR

T1 - Constant temperature molecular dynamics calculation on Lennard-Jones fluid and its application to water

AU - Tanaka, Hideki

AU - Nakanishi, Koichiro

AU - Watanabe, Nobuatsu

PY - 1983

Y1 - 1983

N2 - Constant temperature molecular dynamics calculation has been carried out on Lennard-Jones liquid simulating argon. This method, proposed recently by Andersen, intends to transform the system from microcanonical to canonical ensemble and to keep the temperature of the system a constant value by the generation of random velocities when molecules exchange their energies with heat reservoir with a certain probability. A simple scheme is given for the estimation of collision probability and the effect of introducing this probability on dynamic behavior is examined in detail for 108 argon atoms as a test simulation. It is established that a collision probability of 0.01 is sufficient to realize the constancy of temperature reasonably well with no appreciable disturbance in the dynamic behavior. A model of pure water with Matsuoka-Clementi-Yoshimine (MCY) potential has also been simulated in the same manner. In the case of 2/16 MCY water at 298.15 K, the temperature difference is only 0.73 K with a collision probability of 0.005. Various static properties of MCY water have been calculated with reasonable agreement with those by previous Monte Carlo calculation, and moreover, the dynamic behavior calculated for MCY water gives a satisfactory picture on both translational and rotational motions in water, including reasonable agreement of diffusion coefficient with experimental datum.

AB - Constant temperature molecular dynamics calculation has been carried out on Lennard-Jones liquid simulating argon. This method, proposed recently by Andersen, intends to transform the system from microcanonical to canonical ensemble and to keep the temperature of the system a constant value by the generation of random velocities when molecules exchange their energies with heat reservoir with a certain probability. A simple scheme is given for the estimation of collision probability and the effect of introducing this probability on dynamic behavior is examined in detail for 108 argon atoms as a test simulation. It is established that a collision probability of 0.01 is sufficient to realize the constancy of temperature reasonably well with no appreciable disturbance in the dynamic behavior. A model of pure water with Matsuoka-Clementi-Yoshimine (MCY) potential has also been simulated in the same manner. In the case of 2/16 MCY water at 298.15 K, the temperature difference is only 0.73 K with a collision probability of 0.005. Various static properties of MCY water have been calculated with reasonable agreement with those by previous Monte Carlo calculation, and moreover, the dynamic behavior calculated for MCY water gives a satisfactory picture on both translational and rotational motions in water, including reasonable agreement of diffusion coefficient with experimental datum.

UR - http://www.scopus.com/inward/record.url?scp=0039274975&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0039274975&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0039274975

VL - 78

SP - 2626

EP - 2634

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 5

ER -