Dear Lammps users,
I have a query regrading the thermo command. I have been trying to learn how to use the “heat” “fix” command to perform various non-equilibrium simulations. In my latest test however I have noticed that the thermo command does not appear to be working when used in conjunction with this “fix”. I have attached my input for the mentioned case.
To expand on the specifics of the simulation, I create a box of LJ particles on an FCC lattice which is 18x18x7 lattice units and then proceed to equilibriate it at T=1.0 (LJ units) in the NVT ensemble. I then switch to NVE and add and removed heat each time step from 4 spherical regions which I define earlier in the script. When I attempt to run this part of the simulation the thermo command appears to have no effect as the various information it usually prints to standard output only appears for the zeroth step and no subsequent steps.
Many thanks in advance, Jeff.
in.spheres (905 Bytes)