I’ve been trying to simulate an interfacial system with the COMB3 potential, but have had some problems that I can’t seem to diagnose. For simplicity I recreated the problem with a small bulk alpha-alumina (Al2O3) system, with the initial data file and two input scripts attached.
The input in.COMB3.corundum_newCOMB3 is an input file I am using for a newly compiled version of LAMMPS (29 Mar 2019) which includes an updated version of the COMB3 forcefield files and fix qeq/comb source code from the group of Susan Sinnott https://research.matse.psu.edu/sinnott/software.
in.COMB3.corundum was used with LAMMPS(22 Aug 2018).
Regardless of whether I use the newly compiled version or the 22 Aug 2018 version I have the same issue.
Running an NVE simulation with COMB3 and fix qeq/comb the total energy decays extremely rapidly (~0 K in < 20 ps). While I know that due to the charge dynamics, the energy isn’t conserved, this seems to be something beyond a small drain on the energy due to the charge dynamics.
If I instead use
fix 2 all nvt temp 1043.0 1043.0 0.01
I instead see something different. For an NVT simulation, there is a buildup in the center-of-mass momentum, while the thermal motions go to zero. The thermostat keeps the temperature at 1043 K, but all of the kinetic energy is in center of mass translation.
If I use the same simulation setup for bulk aluminum instead of alumina, I don’t see this behavior. For an NVE simulation of bulk Al, there is a much slower decay of the total energy, more consistent with the loss of energy conservation due to the charge dynamics. In an NVT simulation, There doesn’t appear to be the buildup of total momentum, even over 200 ps or more.
As a final note, if I track the global scalar for the new version of fix qeq/comb from the Sinnott group, it fluctuates around zero.
An example thermodynamic output file is attached as comb3-3728228.
Any ideas for what might be going on here would be greatly appreciated.
Department of Chemistry
The University of Kansas
data.small_corundum (10 KB)
in.COMB3.29Mar19.corundum (781 Bytes)
in.COMB3.22Aug18.corundum (785 Bytes)
comb3-3728228 (947 KB)