Hello,
I would be grateful if anyone on the list could give me the benefit of
their expertise on the following proposed simulation:
System: 2 amorphous approximately spherical silica particles of diameter
2.5 nm, each with a net charge of -9e, with about 5 nm spacing between
their centers of mass, surrounded by TIP4P water and enough Na+ ions to
neutralize the overall system (probably found near the charged silica
surfaces). The system temperature would be set at 300K. The idea is to
watch the particles move towards each other (e.g monitor particle
COM-COM distance as a function of time) due to VdW forces until stopped
by the mutual electrostatic repulsion.
To my naive thinking, there are four ways to tackle this:
1. An isolated sphere of water molecules of diameter at least 8 nm,
centered on the COM of the 2 particles, enclosing both particles, with a
Coulomb cutoff, VdW cutoff and rlist set to 8 nm.
Disadvantages: slow Coulomb and VdW calculation, possible water ordering
effects (see recent DvS paper) with large cutoff, water boiling off the
sphere.
2. Like 1. , only with the outermost water molecules in the droplet
frozen in place to form a 'vesicle', hopefully keeping the water from
boiling off.
3. Particle pair COM centered in a big periodic box (e.g. 12 x 10 x 10
nm), particle pair aligned along the x-axis, box filled with water.
Systems arranged so that each particle is closer to every part of its
immediate twin that to any periodic copies. rvdw also set to 8 nm.
Advantages: gain efficiencies of PME for Coulomb interaction
Disadvantages: Painfully large numbers of atoms, might be possible with
large parallel computing power
4. Do away with the explicit water molecules completely, using an
implicit solvent (but, presumably, keeping explicit Na+ ions and
particles). rvdw set to 8 nm. Any advice on how to do it this way - this
is unknown territory for me? Any GROMACS tutorials for implicit solvent
or GB out there?
Advantages: Small number of atoms
Disadvantages: Lose details of particle-water and ion-water interactions
(but the latter may not be so important, as polarized water models more
realistic than TIP4P for ion-water interactions and hydration shell
structure anyway?)
Insight, suggestions and additional options very gratefully received.
--
Dr. Steven R. Kirk <[EMAIL PROTECTED], [EMAIL PROTECTED]>
Dept. of Technology, Mathematics & Computer Science (P)+46 520 223215
University West (F)+46 520 223299
P.O. Box 957 Trollhattan 461 29 SWEDEN http://taconet.webhop.org
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
