AMBER: Periodic boundary simulation

From: Soon <protege.snu.ac.kr>
Date: Tue, 1 Jun 2004 18:08:31 +0900 (KST)

Hello all,

I am using Amber8 and I made octahedron simulation box by

xx=loadpdb xx.pdb
Solvateoct xx TIP3PBOX 4

on leap with ff03 paramters.


   After minimization of the water solvated system, I performed NVT molecular dynamics simulation to equilibrate the system for about 50.0 ps.

   Now, my question is why on earth the water molecules are coming out of the octahedron simulation box during this NVT molecular dynamics run ? As I know the water molecules are suppose to move around inside the simulation box during the MD simulation. Not out of the box.
The system is not cluster system after all !

Here is my MD input.

 test
 &cntrl
    ntb=1, ig=33333,
    tempi=300.0, temp0=300.0, ntt=1,
    irest=0, ntx=1,
    nstlim=40000,
    dt=0.0017, ntc=2, cut=9.0,
    ntwe=200, ntwx=200, ntpr=200,
 /


Am I missing something ? Any comments will be appreciated. Thank you.

Soon

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Tue Jun 01 2004 - 10:53:00 PDT
Custom Search