From: Peter Freddolino (petefred_at_ks.uiuc.edu)
Date: Fri Dec 01 2006 - 11:13:04 CST
Hi Aytug,
the order you're going to be wanting is minimize-NPT-NVT. Since the only
function of the NPT phase is to get the right cell size/pressure, you
should only run it until the cell size levels off. You should then
equilibrate in the NVT ensemble normally.
Peter
AYTUG TUNCEL wrote:
> Well I guess I got the answer. Somehow my parameters are similar to your example. Here are the parameters,
>
> # Periodic Boundary Conditions
> cellBasisVector1 106. 0. 0.
> cellBasisVector2 0. 117. 0.
> cellBasisVector3 0. 0. 90.
> cellOrigin 32.0 10.0 50.0
> wrapAll on
>
> # PME (for full-system periodic electrostatics)
> PME yes
> PMEGridSizeX 108
> PMEGridSizeY 120
> PMEGridSizeZ 90
>
> The differences between the vectors and gridsizes seem to be in an acceptable range. If that is OK, then NPT water equilibration (before the actual run) should be responsible for the deformation, since I did not do in my simulation. But the question is when (before or after minimization) and how long should it be done? I mean is there anyway to understand whether the water molecules reached equilibrium so that I can continue with NVT run. thanks for the reply.
>
> aytug
>
> -----Original Message-----
> From: "Ana Celia Araujo Vila Verde" <avilaverde_at_engr.psu.edu>
> To: <namd-l_at_ks.uiuc.edu>
> Date: Thu, 30 Nov 2006 18:15:35 -0500
> Subject: RE: namd-l: water box deformation
>
> Hi,
>
> I didn't quite understand what you said related to PMEGRidSize, but you should definitely equilibrate your molecule in NPT before running in NVT, otherwise the number of water molecules for the volume you're simulating will most likely be wrong (if you want water to have a density around 1 g/cm3) and you'll start observing voids in the water.
>
> The PMEgridSize values may be a little smaller or bigger than the actual values you use for the cellBasisVector but they should always be factorizable (if there is such a word) in 2s, 3s or (less good) 5s to speedup the calculations. Eg:
> For
> cellBasisVector1 99.9 0. 0.
> cellBasisVector2 0. 99.9 0.
> cellBasisVector3 0. 0. 63.
>
> I would use
> PME yes
> PMEGridSizeX 96
> PMEGridSizeY 96
> PMEGridSizeZ 64
>
> because 64=2**6 and 96=3 * 2**5
>
> Hope it helps.
> _________________________________
> Ana Célia Araújo Vila Verde
> Penn State University
> Department of Chemical Engineering
> Fenske Laboratory
> University Park, PA 16802
> USA
>
>
> Phone: +(1) (814) 863-2879
> Fax: +(1) (814) 865-7846
> avilaverde_at_engr.psu.edu
> http://www.che.psu.edu/Faculty/JMaranas/group/avv.html
> http://www.gfct.fisica.uminho.pt/Members/acvv
> _________________________________
>
>
>
> -----Original Message-----
> From: owner-namd-l_at_ks.uiuc.edu [mailto:owner-namd-l_at_ks.uiuc.edu] On Behalf Of AYTUG TUNCEL
> Sent: Thursday, November 30, 2006 5:59 PM
> To: namd-l_at_ks.uiuc.edu
> Subject: namd-l: water box deformation
>
>
> Dear all,
>
> I made explicit solvent simulation (NVT) for my protein. While solvating the molecule I made the box size distance 10 ang. away from the protein. But in the early steps of the run after minimization some deformations occured in the water box. When I load the dcd symmetrical collapse of the water molecules at the sides of the box can easily be seen. I did some mining in the archive and figured that it may related to PMEGridSize values.For x and y axis they are greater than the cellBasisVectors by 2 and 3 ang. respectively. Might this be the only reason for the deformation or are there other factors. Another question is, how can I overcome this problem, should I run a short NPT simulation before the actual run to establish a uniform water distribution (density) in the box? thanks in advance.
>
> Aytug
>
>
>
>
>
>
>
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:44:13 CST