I am modeling a simple 2D Graphitic Carbon Nitride (g-C3N4) nanosheet with 54 atoms. I used several force fields with both LAMMPS and GROMACS and the potential energies of the system in npt at 1K and 1 atmosphere are:

  1. LAMMPS using ffield.reax.CHOSFCl = -6,024 Kcal/mol
  2. LAMMPS using tersoff = -7,148 Kcal/mol
  3. LAMMPS using ATB server - gromos54a7 forcefield = -270 Kcal/mol
  4. GROMACS: using ATB server - gromos54a7 forcefield = -228 Kcal/mol

We know that the zero for potential energy is not the same for all the force fields but is it usual to be an order of magnitude difference between REAX and GROMOS force fields? P.S. All the force fields were recommended and used by the literature for the g-C3N4.

  • 1
    $\begingroup$ I would be less interested in the absolute energies and more interested in a relative energy. Do you have any relative energies benchmarked against previous computational studies? $\endgroup$ – jezzo Aug 9 '20 at 0:47
  • $\begingroup$ You might get a better answer asking this in mattermodeling.stackexchange.com $\endgroup$ – Tristan Maxson Aug 9 '20 at 5:23

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.