Dear, amber10 users,
i am trying to constrain a torsion angle which is formed by the centers of
mass of four different atom groups.
But no matter what group members i were choosing, the initial dihedral angel
in "log" file is the same wrong
value.
The different machine and compiler have been tried and came out with the
same wrong value.
I tried to constrain angle and distance which are formed by the centers of
mass of atom groups,
and the results are correct. Anyone has an idea what might cause this group
torsion constrain goes wrong?
Thanks a lot!!!
Best,
Han Zheng
Here are my IN and RST file (tor.RST):
----------------------------------------------------------------------------------------
minization of DNA in water
&cntrl
imin = 1,
maxcyc = 50,drms=0.01,
ncyc = 50,
ntb = 1,
igb = 0,
cut = 12,
NTR=0,
iwrap=1
nmropt=1,
/
&wt type=DUMPFREQ, istep1=1, /
&wt type=END, /
DISANG=tor.RST
DUMPAVE=minrestraint.log
LISTIN=POUT
LISTOUT=POUT
--------------------------------------------------------------------------------------------
&rst iat = -1, -1, -1, -1, 0,
r1=-180., r2=-51.659, r3=-51.659, r4=180., rk2 =1000.0, rk3=1000.0,
igr1 = 513, 514,
igr2 = 511, 512,
igr3 = 478, 479,
igr4 = 197, 198,
&end
--------------------------------------------------------------------------------------------
In the output file:
Here are comments from the DISANG input file:
******
COM ( -1)-COM ( -1)-COM ( -1)-COM ( -1) NSTEP1= 0
NSTEP2= 0
Atom ranges defining Center of Mass Group in the first position:
513 -> 514/
Atom ranges defining Center of Mass Group in the second position:
511 -> 512/
Atom ranges defining the Center of Mass Group in the third position:
478 -> 479/
Atom ranges defining the Center of Mass Group in the 4th position:
197 -> 198/
R1 =-180.000 R2 = -51.659 R3 = -51.659 R4 = 180.000 RK2 =1000.000 RK3 =
1000.000
Rcurr: 85.008 Rcurr-(R2+R3)/2: 136.667 MIN(Rcurr-R2,Rcurr-R3):
136.667
Number of restraints read = 1
It seems the group constrain has been added successfully, and the initial
torsion angle is calculated correctly as 85.008.
But at the end of the output file:
------------------------------------------------------------------------------
First atom Last atom curr. value target deviation penalty
------------------------------------------------------------------------------
*H1' DT 17 -- *H1' DG 16: -51.708 85.000 136.708 5693.016 t
Total torsion penalty: 5693.016
------------------------------------------------------------------------------
>From the minrestraint.log file:
0 -51.659
1 -51.659
2 -51.660
3 -51.661
4 -51.662
5 -51.663
6 -51.665
7 -51.667
8 -51.670
9 -51.673
....
49 -51.708
*
No matter what atom groups i chose for this dihedral angle, the initial
value in log file is always -51.659.
It seems the program recognizes the center of mass dihedral angel, but
during the minimization/dynamics,
it is calculating something else. *
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Oct 26 2009 - 11:30:03 PDT
