I'm the culprit here. I came up with that when we thought we needed
parameters off the usual combining matrix. What the numbers stand for is
the Lennard-Jones sigma and epsilon of the first atom type when it
interacts with anything under the normal rules, then the sigma and epsilon
of the second atom type when it interacts with the first. So all of those
OW's in the parm14ipq.dat file have the sigma and epsilon of TIP4P-Ew water
oxygen (the OW parameters are recast from the TIP3P which I sloppily left
in the list above--TIP4P-Ew is the water model of choice for ff14ipq).
Then, when type ND interacts with OW, it does so as if the ND has sigma and
epsilon 2^(5/6) * 1.960000 A and 0.170000 kcal/mol. So the interaction of
ND and OW takes place according to
sigma = 0.5 * (1.775931 + 1.960000) * 2^(5/6) = 3.3283A
epsilon = sqrt(0.162750 * 0.170000) = 0.16634 kcal/mol
That's the way LJEDIT works. I'll see about getting this into the manual,
but it's rare that someone has wanted the feature.
Dave
On Mon, Jul 23, 2018 at 2:16 PM Matias Machado <mmachado.pasteur.edu.uy>
wrote:
> Dear AMBER experts,
>
> I'm trying to set some Lennard-Jones interactions out-off the combination
> rule. I was able to do it by using Parmed, however I want to set them in a
> frcmod file so they can be applied by leap without further modification of
> the topology.
>
> I know the LJEDIT section (available from amber14 on) would allow me to do
> so, however its syntax is poorly documented, indeed I couldn't find any
> mention to that keyword in the Amber manuals, nor in the file formats at
> the AMBER web page [http://ambermd.org/FileFormats.php]
>
> The only reference I found was this post in the mailing list [
> http://dev-archive.ambermd.org/201301/0018.html], in which its
> implementation is discussed. From those comments, I would expect the
> following syntax:
>
> LJEDIT
> <Type 1> <Type 2> <R* pair> <Eps pair>
> ...
> END
>
> However, looking at the implementation of parm14ipq.dat [
> https://www.ncbi.nlm.nih.gov/pubmed/25328495] and CUFIX
> (frcmod.ff99cufix) [http://bionano.physics.illinois.edu/CUFIX]
> I found 4 terms instead of 2 for defining the LJ interaction, my guess is
> that 1-4 interactions can also be set separately but I'm not sure about
> that...
>
> In any case... what does each set of parameters stand for?
>
> Best regards,
>
> MatÃas
>
> ------------------------------------
> PhD.
> Researcher at Biomolecular Simulations Lab.
> Institut Pasteur de Montevideo | Uruguay
> [http://pasteur.uy/en/laboratorios-eng/lsbm]
> [http://www.sirahff.com]
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Jul 23 2018 - 12:30:03 PDT
