<div dir="ltr"><div>Dear all,</div><div>I was running a simulation with bulk water and CNT, with the protocol -&gt; 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&#39;t figure it out.<br></div><div><br></div><div>GROMACS: (in kJ/mol)<br></div><div><br></div><div>Coulomb (SR) Coul. recip. <br> -3.59079e+05 8.41910e+02 <br></div><div><br></div><div>in kcal/mol this translates to <br></div><div>Potential Kinetic En. Total Energy Temperature Pres. DC <br>-85821.94072 201.221319</div><div><br></div><div>LAMMPS: (in kcal/mol)<br></div><div><br></div><div>E_coul = 177903.2534 E_long = -270193.3273 <br></div><div><br></div><div><br></div></div>