Journal article
Molecular dynamics at constant Cauchy stress.
-
Miller RE
Department of Mechanical and Aerospace Engineering, Carleton University, Ottawa, Ontario K1S5B6, Canada.
-
Tadmor EB
Department of Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, Minnesota 55455, USA.
-
Gibson JS
Department of Mechanical and Aerospace Engineering, Carleton University, Ottawa, Ontario K1S5B6, Canada.
-
Bernstein N
Center for Materials Physics and Technology, Naval research Laboratory, Washington, DC 20375, USA.
-
Pavia F
STI IGM LAMMM, Ecole Polytechnique Federale de Lausanne, 1015 Lausanne, Switzerland.
Published in:
- The Journal of chemical physics. - 2016
English
The Parrinello-Rahman algorithm for imposing a general state of stress in periodic molecular dynamics simulations is widely used in the literature and has been implemented in many readily available molecular dynamics codes. However, what is often overlooked is that this algorithm controls the second Piola-Kirchhoff stress as opposed to the true (Cauchy) stress. This can lead to misinterpretation of simulation results because (1) the true stress that is imposed during the simulation depends on the deformation of the periodic cell, (2) the true stress is potentially very different from the imposed second Piola-Kirchhoff stress, and (3) the true stress can vary significantly during the simulation even if the imposed second Piola-Kirchhoff is constant. We propose a simple modification to the algorithm that allows the true Cauchy stress to be controlled directly. We then demonstrate the efficacy of the new algorithm with the example of martensitic phase transformations under applied stress.
-
Language
-
-
Open access status
-
closed
-
Identifiers
-
-
Persistent URL
-
https://sonar.ch/global/documents/281789
Statistics
Document views: 21
File downloads: