AMBER Archive (2008)

Subject: AMBER: How to manually adjust periodic boundary conditions (PBC)

From: Chin, Keith B (keith.b.chin_at_jpl.nasa.gov)
Date: Wed Nov 19 2008 - 13:58:41 CST

Hi,

I was wandering if anyone know how to adjust periodic boundary conditions in sander *crd and *top files for MD.

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.

I tried using the xleap SetBox cmd to assign PBC parameter close to my hex parameters, and then manually adjust the gamma angle to 120 instead of 90 in the last line of the *.crd file. However, this gave me an error of when executing the MD, which indicates that my PBC parameters is still too small for my lattice.

Here's the error in my md output file:

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

Only when I expand my PBC parameters way beyond my original lattice size does the MD executes completely, but obviously it'll provide bad results, especially for density.

Any advise on how to do this OR where I can find more info on resolving this issue.

Thanks in advance.
-------------------------------------------------------
Keith Chin
Electrochemical Technologies Group
JPL/NASA
Office: (818) 354 - 4051

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
      to majordomo_at_scripps.edu