[gmx-users] Nonbonded energy of 1-4 interactions
zhonghua.xia at helmholtz-muenchen.de
Fri Sep 21 14:25:35 CEST 2018
Thank you so much for your answer!
So, if I want to calculate the nonbonded energy of solute-solute and solute-solvent, I should set the parameter as below:
energygrps = MOL MOL MOL SOL
And If I want to customise FudgeLJ and FudgeQQ, is it feasible to directly modify their values in the forcefield.itp file?
On 9/20/18 5:23 AM, Zhonghua Xia wrote:
> Dear all,
> I want to get the nonbonded energy (1-4 electrostatic and 1-4 LJ) of solute-solute and solute-solvent, respectively. How can I make it? What?s more, how can I set the scaling-factors of them?
These energies are in LJ-14 and Coul-14 in the .edr file. If you want
them decomposed for different groups, you need to specify energygrps in
the .mdp file, re-create a .tpr file with those groups, and post-process
your trajectory. Note that 1-4 interactions are intramolecular, so there
is no such thing as a solute-solvent 1-4 energy. The term also has no
physical meaning as it is completely dependent upon force field
convention and dihedral terms.
In GROMACS, the 1-4 scaling factors are set with FudgeLJ and FudgeQQ in
the [defaults] directive of forcefield.itp.
More information about the gromacs.org_gmx-users