I am carrying out a hugoniostat simulation on Al. The script I have used
is the one available on LAMMPS website for a LJ fcc solid. I have
converted this script into metal units and want to perform a hugoniostat
simulation for Pzz=20 GPa for ALuminium FCC crystal. I have first
minimized the system using ‘box/relax’ and ‘miminize’ commands. Further,
compressed the system using ‘hugoniostat’ as given in the script.
When the run corresponding to the compression starts, all the
parameters (temp, kinetic energy , pressure etc) starts diverging with time as can
be seen in the output file attached. I could observe that the potential
Energy (4th column, E_pot in the output file) being calculated by LAMMPS
is 0 which may be a source of this divergence during the evolution in
the system. Can anyone tell why this script is giving zero potential
energy throughout the evolution. Input file is also attached.
in.txt (5.01 KB)
Al.eam.alloy (576 KB)
out.out (115 KB)