On 28/01/2011 4:25 AM, WU Yanbin wrote:
Dear GMXers,
I would like to reproduce the water droplet contact angle on graphite
surface, as in Werder, T.; Walther, J. H.; Halicioglu, T.;
Koumoutsakos, P. /J. Phys. Chem. B/ *2003*, /107/, 1345-1352.
The box size is 20nm by 20nm by 30nm. Graphite is represented by a
two-layer carbon sheet. After equilibration, the water droplet is
around 3.5nm in radius. A large part of the box is vacuum. Periodic
boundary condition (PBC) is applied. Graphite is fixed throughout the
simulation.
I'm not sure what PBC is gaining for you. I think you want vacuum
boundary conditions, which may as well be non-periodic.
From the manual, the "comm_mode=Linear" option should be applied,
which means the translational motion of the center of mass of the
water droplet should be removed. My first question is: If the time
step is 2fs, what's the recommended value for nstcomm? (nstlist is
updated very 3 time step)
10 is probably fine.
Also in the manual and this mail list, it's mentioned that for
non-PBC, isolated system, "comm_mode=Angular" should be applied. Does
this also applied to my partial filled PBC box?
No, this is referring to simulations like a protein in vacuum.
Mark
Do let me know if the question is not clear.
Thank you.
Best,
Yanbin
--
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