David:
OK, thanks.
Now, suppose I did create a supercell with good box dimensions, but the atom coordinates came from my own program (Not xleap).
How do I manually set periodic boundary conditions to hexagonal crystal lattice with alpha = 90, beta = 90, gamma = 120 (NOT box with alpha = 90, beta = 90, gamma = 90) to run MD in sander? I used the setbox command to set initial PBC close to actual lattice dimensions prior to changing it in last line of *crd (coordinate) file.
What about the *top (topography) file BOX_DIMENSION parameter such as:
%FLAG BOX_DIMENSIONS
%FORMAT(5E16.8)
9.00000000E+01 1.15728620E+01 1.17719310E+01 1.31228620E+01
How does sander use these parameters as oppose to PBC dimensions in last line of *.crd file
Does sander recognize my atom coordinates as opposed to the ones created by xleap SaveamberParm in applying the my hex lattice parameters?
Any info would be greatly appreciated.
Thanks again in advance
-keith
-----Original Message-----
From: owner-amber.scripps.edu [mailto:owner-amber.scripps.edu] On Behalf Of David A. Case
Sent: Wednesday, November 19, 2008 3:36 PM
To: amber.scripps.edu
Subject: Re: AMBER: How to manually adjust periodic boundary conditions (PBC)
On Wed, Nov 19, 2008, Chin, Keith B wrote:
>
> Suppose rather than a box lattice, I wish to apply PBC for hex lattice,
> where a=4.516, b= 4.516, c = 7.435, alpha = 90, beta = 90, gamma = 120.
Amber cannot do this small a unit cell. Basically, the box size must be at
least twice the cutoff + "skin" (usually 1 to 2 Ang.) Your only real
alternative is to simulate several unit cells as one "supercell". A cell size
of 18 Ang. is about the practical minimum for Amber.
>
> Ewald PARAMETER RANGE CHECKING:
> parameter skinnb+cutoffnb: (nonbond list cut) has value 0.12000E+02
> This is outside the legal range
> Lower limit: 0.00000E+00 Upper limit: 0.58665E+01
> Check ew_legal.h
This is basically saying (in jargon) what I said above. It looks like you are
setting the cutoff to 10, which is too large -- you should probably stick to
the default of 8.
...good luck...dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Fri Dec 05 2008 - 16:05:19 PST
