[gmx-users] Coulomb and vdW energies

Shreyas Kaptan shreyaskaptan at gmail.com
Sun Dec 1 11:29:02 CET 2019


Hi.

I would not say that it is completely unrealistic/unphysical. Especially,
the Electrostatic part. If within two subgroups, you could show that there
is systematic repulsive or attractive potential, that could lead to some
hypothesis about how they interact (e.g. pushing/pulling). Of course, this
is only true of ratios of decomposed energies that are large proportions of
the total energy, as smaller portions might not be really
significant indicators of strong interactions.

Shreyas

On Sat, Nov 30, 2019 at 6:47 PM Justin Lemkul <jalemkul at vt.edu> wrote:

>
>
> On 11/30/19 3:14 AM, elham wrote:
> > Dear all, I want to obtain the coulomb and vdW energies between two
> atoms of separate molecules, However  the .edr output calculate coulomb
> and vdW energies of whole system. There is not any option(like -n .ndx) to
> define the the special atoms in the gmx energy command. How can I calculate
> the coulomb and vdW energies between the two atoms which are related to two
> separate molecules?Best
>
> Set energygrps in the .mdp file, generate a new .tpr, and use mdrun
> -rerun to re-evaluate the energies from an existing trajectory. Note
> that decomposition of energy in this way has no physical meaning. Total
> interaction energy may be physically realistic if the force field was
> specifically parametrized against e.g. QM target data but decomposing
> the energy in this way is not meaningful.
>
> -Justin
>
> --
> ==================================================
>
> Justin A. Lemkul, Ph.D.
> Assistant Professor
> Office: 301 Fralin Hall
> Lab: 303 Engel Hall
>
> Virginia Tech Department of Biochemistry
> 340 West Campus Dr.
> Blacksburg, VA 24061
>
> jalemkul at vt.edu | (540) 231-3129
> http://www.thelemkullab.com
>
> ==================================================
>
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-request at gromacs.org.



-- 
Shreyas Sanjay Kaptan


More information about the gromacs.org_gmx-users mailing list