> Thanks Justin for your reply. To avoid compressing down the cell I thought
> using semiisotropic option with 0 compressibility in Z would be
> appropriate.
>
>

You must use NVT only. Otherwise, the cell will compress in XY directions
to compensate its inability to compress in Z direction.




> Here are the steps:
>
> 1- First I did a 10 ns NPT to equilibrate the box and used the last frame
> gro and cpt file as input for the next step
> 2- I extended the box in Z to more than twice the initial box size and
> issued the following:
>
> grompp -f md.mdp -c extendedZ.gro -t .cpt -p .top -o extendedZ.tpr
> mpirun -np 4 mdrun_mpi -deffnm .extendedZ -s -o -c -g -e -x -v
>
> My questions are:
>
> 1) Can I get the proper solvent-air interface to obtain density profile by
> extending the Z direction (last line of the gro file obtained from NPT at
> each pressure)? Is that what yo mean by " then build a new unit cell and
> run under NVT" ?
>



If you want solvent/air interface, you should add $air.

If you want solvent/vacuum interface, it is enough to extend the box.

If you want liquid/vapor interface, you need to re-equilibrate the system
in NVT with a space for vapor available in the box.


Dr. Vitaly Chaban
-- 
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to