This is a testing run. I first used “minimize” with “box/deform” to reduce pxx, pyy, pxy and fmax to almost zero. Then I ran nve. I found that the pxx, pyy, pxy and fmax in the last line of the thermo output for “minimize” is slightly different from those in the first line of the thermo output for “run”. Though the difference is very small, I was wondering why that was happening. I expected them to be identical.
In addition, the thermo output is still changing slightly during “run” even when I set all force components and velocities to zero.
The in file and log file are as follows. Thank you very much!
Dan
In file:
units metal
dimension 2
atom_style atomic
boundary p p p
This is a testing run. I first used "minimize" with "box/deform" to reduce
pxx, pyy, pxy and fmax to almost zero. Then I ran nve. I found that the pxx,
pyy, pxy and fmax in the last line of the thermo output for "minimize" is
slightly different from those in the first line of the thermo output for
"run". Though the difference is very small, I was wondering why that was
happening. I expected them to be identical.
have you ever heard some stories about the limitations of
floating-point math? and what happens, if you compute numbers, that
depend on summing up values of equal magnitude but opposite sign.
In addition, the thermo output is still changing slightly during "run" even
when I set all force components and velocities to zero.
unless you do a run with "pre no", you have a reneighboring happen
during setup and that can result in a reordering of the atoms and thus
affect your numbers. without knowing the exact system setup and how
large your individual force contributions are (before they get summed
up), it is difficult to say, whether you have a genuine problem or are
just overinterpreting your data. the fact that you provide epsilon and
sigma with excessive precision and yet use an extremely short cutoff
hints at the latter.
Thank you, Axel. I understand why the thermo at the beginning of run differs from the thermo at the end of minimize. But I do not understand very well why the thermo still changes during the first run step, i.e.from Step 456 to Step 457 since the setup and reneighboring have finished before the first run step. When the system evolves during the step, no interatomic forces exist, so the atoms should stay right at the same positions. Also, no reneighboring/reordering is taking place. So I would expect that the thermos remain unchanged. I was wondering if you could explain a little more about this and point out which part of my reasoning is wrong.
Thank you, Axel. I understand why the thermo at the beginning of run differs
from the thermo at the end of minimize. But I do not understand very well
why the thermo still changes during the first run step, i.e.from Step 456
to Step 457 since the setup and reneighboring have finished before the first
run step. When the system evolves during the step, no interatomic forces
exist, so the atoms should stay right at the same positions. Also, no
reneighboring/reordering is taking place. So I would expect that the thermos
remain unchanged. I was wondering if you could explain a little more about
this and point out which part of my reasoning is wrong.
i will try provide some more explanations, if you answer *my* question
first and request for the missing data to reproduce your test. it is a
waste of time to discuss in this hypothetical way without being able
to reproduce what you are seeing and making additional checks and
tests to support my assessments.
if you believe there is a bug in LAMMPS, it is your task to prove it
is there, not my task to prove it isn't.