par Ryckaert, Jean-Paul ;Ciccotti, Giovanni
Référence The Journal of Chemical Physics, 78, 12, page (7368-7374)
Publication Publié, 1982
Référence The Journal of Chemical Physics, 78, 12, page (7368-7374)
Publication Publié, 1982
Article révisé par les pairs
Résumé : | The method of constant pressure molecular dynamics (MD), developed by Andersen for monoatomic fluids is extended to the MD, in Cartesian coordinates, of molecular systems with constraints. Andersen's proof is easily generalized after decoupling internal degrees of freedom from the center-of-mass. Only these last degrees are directly affected by Andersen's transformation (demon). The Cartesian equations of motion of individual atoms are derived from a generalized Andersen's Lagrangian. The equations are quite similar to those of the usual MD simulation at constant volume apart from an additional term coupling the molecular center-of-mass and the volume of the sample. The volume appears now as a dynamical variable evolving from the imbalance between imposed external pressure and instantaneous values of the molecular stress tensor. Some numerical aspects are discussed and the technique is briefly illustrated for the case of rigid diatomic molecules. © 1983 American Institute of Physics. |