I'm studying a bio-molecular system using different ensembles(NVE, NVT, and NPT). What method(algorithm) for each ensemble was used for LAMMPS? I think there are a lot of variations even for Nose-Hoover thermostat for NVT. Or would you let me know the references you(for developers) used to develop the codes, esp. for NPT simulation?
My another question is related to MSD in the NPT ensemble under periodic boundary condition which I asked a few weeks ago. I guess I've created a code which calculates MSD from LAMMPS output. But, I think there are two components(?) of MSD, one from the forces by the thermostat and one from the use of barostat. I think the displacement by the barostat should be neglected for true MSD calculation because it comes from some artificial(?) movement(dilation or contraction of the simulation box). I see that fix_msd can be combined with fix_npt to calculate MSD in LAMMPS. In that case how does it work for correct MSD calculation? It seems to me that the code uses the true displacement from current box size(including the displacement by the expansion or contraction of the box), which is contradictory to my understanding. Am I right or wrong? Or how can I correctly calculate MSD in the NPT ensemble?
Thanks in advance,
Dept. of Chem. Eng.
TAMU 3122 Texas A&M University
College Station, TX 77840