I’m studying, and trying to reproduce as best I can the study in Parrinello, and Rahman’s, “Crystal Structure and Pair Potentials: A Molecular Dynamics Study”, PRL, (45) no. 14, 1980, with lammps.
Among the numerous choices for deforming the simulation box, including the howto section of the doc on triclinic boxes,
http://lammps.sandia.gov/doc/Section_howto.html#howto_12
I’ve decided on using fix npt as the simplest choice in this self-training exercise.
How can I correctly implement the quantity in the paper, “W*”, which the authors call a mass parameter in my npt lammps simulation? (Parrinello & Rahman used a value of W* = 20)
Am I also using the npt fix correctly to indepently vary the x, y, and z dimensions of the simulation box correctly?
I output the box sides as thermo during the runs and they all look approximately equal throughout the entire run.
Any other guidance/comments appreciated.
respectfully,
test.in (1.09 KB)