Re: [AMBER] group specification for different force constants for different groups

From: Jason Swails <jason.swails.gmail.com>
Date: Wed, 29 Aug 2012 13:41:09 -0400

On Wed, Aug 29, 2012 at 1:06 PM, Sajeewa Pemasinghe <sajeewasp.gmail.com>wrote:

> Hi everyone,
>
> I want apply a force constant of 500.0 (kcalmol-1A-2) for residues 1-569,
> 400 for 580-607 and 300 for 615-992. Is my group specification below
> correct?
>
> GroupA
> 500.0
> RES 1 569
> END
> GroupB
> 400.0
> RES 580 607
> END
> GroupC
> 300.0
> RES 615 992
> END
> END
>

Appears to be.

Note that these force constants are *enormous* for MD simulations, and
could likely cause integration failures if you are running MD. For MD,
consider restraints about 2 orders of magnitude lower.

For minimizations, this difference will not have much of an effect, I don't
think.

HTH,
Jason

-- 
Jason M. Swails
Quantum Theory Project,
University of Florida
Ph.D. Candidate
352-392-4032
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Aug 29 2012 - 11:00:03 PDT
Custom Search