On Thu, Jul 22, 2004, Johnson Agbo wrote:
> I am running minimization and dynamics on cubic ice with periodic
> boundary conditions,
Amber has limitations in simulating small unit cells. Basically, the
cutoff + "skin" has to be less than half of the unit cell size.
>
> Ewald PARAMETER RANGE CHECKING:
> parameter skinnb+cutoffnb: (nonbond list cut) has value 0.10000E+02
> This is outside the legal range
> Lower limit: 0.00000E+00 Upper limit: 0.95370E+01
So, you need to use a smaller cutoff + skin (probably by reducing skinnb
from its default value of 2 to a value of 1), or create a bigger unit cell.
> I went into the program file ( ew_legal.h ) and changed the
> coefficients to the values shown below,
> ( _REAL_ ew_coefflo,ew_coeffhi
> parameter(ew_coefflo=0.0d0,ew_coeffhi=17.d0)
Your problem has really nothing to do with the ewald coefficeint....
....dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Jul 23 2004 - 10:53:00 PDT