[gmx-developers] Discrepancy with LAMMPS in Coulomb (SR) Coul. recip. energy
mail2majumdar at gmail.com
Thu Oct 24 00:12:03 CEST 2019
I was running a simulation with bulk water and CNT, with the protocol ->
Minimization, NPT, NVT. Keeping as much close as possible, I did the same
simulation in LAMMPS too. I observe that the LAMMPS and GROMACS log file
shows complete disagreement in sign and value in the coulomb energy
calculation in real space and reciprocal space; in all the steps. Although
upon adding the two parts (real+reciprocal) to get total Coulomb energy,
they are matching within error limit. Can anyone please explain me why am i
getting such a weird result? Is it some in-build difference in calculation
technique, or am I going wrong somewhere? I couldn't figure it out.
GROMACS: (in kJ/mol)
Coulomb (SR) Coul. recip.
in kcal/mol this translates to
Potential Kinetic En. Total Energy Temperature Pres. DC
LAMMPS: (in kcal/mol)
E_coul = 177903.2534 E_long = -270193.3273
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the gromacs.org_gmx-developers