LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Reaxff c/h/o energy conservation and temperature drift in NVE
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Reaxff c/h/o energy conservation and temperature drift in NVE


From: Benjamin Jensen <bdjensen@...1754...>
Date: Thu, 14 Dec 2017 07:51:03 -0500

The Chenoweth CHO parameters have some pretty steep gradients, and usually requires smaller than typical time steps for energy conservation. I wouldn't be surprised if 0.5 fs was still to large. I'd rerun with decreasing time steps and look for convergence.
-Ben

On Wed, Dec 13, 2017 at 5:08 PM, Ray Shan <rshan@...1795...> wrote:
The ReaxFF combustion force field was fitted to describe the combustion of hydrocarbons, not graphene oxides.  You might just be using a force field that is not appropriate for your simulation.

Ray

On 12/13/17, 1:58 PM, "Thomas Zhang" <tk.zhang@...3421...> wrote:

Dear LAMMPS community,


This question is regarding the Reaxff c/h/o combustion force field (Chenoweth, Van duin, and Goddard, 2008). I am simulating Graphene Oxide (GO), and using the NVE ensemble with NEMD to calculate thermal conductivity. There are issues with energy conservation in the NVE ensemble. For example, one system is GO with 50% oxygen (hydroxyl) concentration. The system is first equilibrated using NPT at 300K for 4ns, and a timestep of 0.5fs. Then it is transitioned from NVT to NVE for 25ps each. When the calculation is done in NVE, and a heat flux is applied (using the M-P algorithm), the temperature drifts from 300K - ~370K in the span of 300ps, as shown in the first attached file. Most likely the equilibration and timestep used are sufficient. This issue is less pronounced when the temperature gradient is removed.

I have tested the same methodology (NEMD) with Silicon using the SW potential, and the temperature does not drift at all, as shown by the second attached file.


Please let me know if you can be of any help,


Thanks!


Thomas


------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
lammps-users mailing list
lammps-users@...6297....sourceforge.net
https://lists.sourceforge.net/lists/listinfo/lammps-users