[AMBER] Issues with 12-6-4

From: Nell Karpinski via AMBER <amber.ambermd.org>
Date: Fri, 14 Feb 2025 14:29:03 -0500

Hello,

I am following the "Modeling a Magnesium-DNA System using the 12-6-4
LJ-Type Nonbonded Model, Example A" tutorial to use LJ_1264 for Hg2+. I am
attempting to use 'add_12_6_4 :HG2', however, I am getting all zeros for
my LJ CCoeff.

%FLAG LENNARD_JONES_CCOEF
%COMMENT For 12-6-4 potential used for divalent metal ions
%FORMAT(5E16.8)
  0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
 0.00000000E+00
  0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
 0.00000000E+00
  0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
 0.00000000E+00
  0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
 0.00000000E+00
  0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
 0.00000000E+00
  0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
 0.00000000E+00
  0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
 0.00000000E+00
  0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
 0.00000000E+00
  0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
 0.00000000E+00

I have tried to use parmed with Amber/18.0, Amber/22.0, and Amber/24.0, but
I obtain the same result.

My parmed.in file contents are:

loadRestrt hg_final_mango_ligand.inpcrd

setOverwrite True

add12_6_4 :HG2 watermodel SPCE

outparm hg_final_mango_ligand_1264.prmtop hg_final_mango_ligand_1264.inpcrd

Am I doing something wrong here? Please let me know, thank you in advance!

All the best,
Nell
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Feb 14 2025 - 11:30:02 PST
Custom Search