From: Francesco Pietra (chiendarret_at_gmail.com)
Date: Tue Jul 24 2012 - 03:23:19 CDT
Hi:
I was running MD without problems for a protein with namd 2.9-cuda,
amber PARM7 ff.
Now, I have inserted a short lipid between two aa, using GAFF params
for the lipid. Now, minimization could not be carried out: "Error: the
periodicity of dihedral #249 is zero!"
Here the relevant part of prmtop file, where the said zero is at line
3 from the bottom:
%FLAG DIHEDRAL_PERIODICITY
%FORMAT(5E16.8)
1.00000000E+00 2.00000000E+00 2.00000000E+00 2.00000000E+00 3.00000000E+00
4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 3.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00 3.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 1.00000000E+00 3.00000000E+00
1.00000000E+00 3.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00
4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 1.00000000E+00 2.00000000E+00
3.00000000E+00 4.00000000E+00 3.00000000E+00 2.00000000E+00 4.00000000E+00
3.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00
2.00000000E+00 3.00000000E+00 4.00000000E+00 2.00000000E+00 4.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00
2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00
3.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00
4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00
2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00
3.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00
4.00000000E+00 2.00000000E+00 2.00000000E+00 1.00000000E+00 2.00000000E+00
3.00000000E+00 4.00000000E+00 1.00000000E+00 3.00000000E+00 2.00000000E+00
4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00 2.00000000E+00
3.00000000E+00 4.00000000E+00 2.00000000E+00 1.00000000E+00 2.00000000E+00
3.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00
4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00
2.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00
4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 1.00000000E+00 2.00000000E+00
3.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00
4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00
2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00
3.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00
4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00
3.00000000E+00 4.00000000E+00 2.00000000E+00 1.00000000E+00 2.00000000E+00
2.00000000E+00 2.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00
3.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00
2.00000000E+00 3.00000000E+00 4.00000000E+00 2.00000000E+00 2.00000000E+00
2.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00
2.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00
2.00000000E+00 3.00000000E+00 4.00000000E+00 3.00000000E+00 1.00000000E+00
2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00
3.00000000E+00 4.00000000E+00 2.00000000E+00 2.00000000E+00 2.00000000E+00
2.00000000E+00 2.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00
1.00000000E+00 2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00
2.00000000E+00 3.00000000E+00 4.00000000E+00 1.00000000E+00 2.00000000E+00
3.00000000E+00 4.00000000E+00 3.00000000E+00 0.00000000E+00 3.00000000E+00
1.00000000E+00 2.00000000E+00 1.00000000E+00 2.00000000E+00 2.00000000E+00
2.00000000E+00
I am aware that CHARMM27 uses n= 1, 2, 3, 4 in
V DA = V n [ 1 + cos (n φ2 − φ 0 )]
while I don't know if the above zero is allowed in PARM7, or if it is
an error anyway.
Frankly, I don't know from where to beging to fix the problem, and
which other data should I produce to get help.
Thanks for your understanding
francesco pietra
This archive was generated by hypermail 2.1.6 : Tue Dec 31 2013 - 23:22:18 CST