Re: AMBER: Distance Restraints Problem

From: David A. Case <case.scripps.edu>
Date: Tue, 2 Jan 2007 15:37:18 -0800

On Tue, Jan 02, 2007, Esther Brugger wrote:
>
> (1) Is it correct if I put ntc=2, ntf=2 in this distance restraints MD
> simulation? or I shouldn't use it in this restraint calcuation.

Should be OK, assuming you are not trying to restrain a C-H (or N-H, O-H) bond
distance.


> Here is the input file:
> NPT MD 100ps with Distance Restrains
> &cntrl
> nmropt=1, pencut = -0.001,
> ntb = 2, pres0 = 1, ntp = 1,
> taup = 2.0, cut = 12,
> ntc = 2, ntf = 2,
> tempi = 300.0, temp0 = 300.0,
> ntt = 3, gamma_ln = 1.0,
> nstlim = 50000, dt = 0.002, vlimit=20,
> ntpr = 100, ntwx = 100, ntwr = 1000
> /
> & wt type = 'END'/
> DISANG=DIST.RST
>
>
> (3) Here is the file of my DISANG file, I really wish someone can help me
> to write it in a simple way.

I don't understand the queestion. What is wrong with the DISANG file (looked
OK to me at a brief glance....)? Looks like you are restraining the various
H-H distances to be between 2 and 5 Ang.

[I don't think there is anyway to make it "simpler": DISANG files are meant to
be read by a computer; they can also be read and edited by humans, but I don't
see why you need to do that here.]

...regards...dac

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Jan 03 2007 - 06:07:30 PST
Custom Search