<div dir='auto'><div>Hi,<div dir="auto"><br></div><div dir="auto">This is a question for the gmx-users list, but I'll answer here.</div><div dir="auto"><br></div><div dir="auto">In general decomposition of Coulomb PME terms is meaningless, because there is no unique definition of the components. Also you need to have all parameters identical. Having said that, I think the main difference is because GROMACS adds the self correction to the real space term and LAMMPS seems to add it to the the reciprocal space.</div><div dir="auto"><br></div><div dir="auto">Cheers,</div><div dir="auto"><br></div><div dir="auto">Berk</div><br><div class="gmail_extra"><br><div class="gmail_quote">On Oct 24, 2019 00:11, Jeet Majumdar &lt;mail2majumdar@gmail.com&gt; wrote:<br type="attribution"><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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'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>&nbsp;-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>
</blockquote></div><br></div></div></div>