[gmx-users] Regular vs. CHARMM TIP3P water model
Berk Hess
gmx3 at hotmail.com
Tue Sep 14 13:27:43 CEST 2010
> Date: Tue, 14 Sep 2010 13:11:47 +0200
> Subject: RE: [gmx-users] Regular vs. CHARMM TIP3P water model
> From: nicolas.sapay at cermav.cnrs.fr
> To: gmx-users at gromacs.org
>
> >
> > Hi,
> >
> > I don't understand what exactly you want to reproduce.
> > Standard tip3p and Charmm tip3p are different models, so the density does
> > not have to be identical.
> > The Gromacs Charmm FF implementation paper:
> > http://pubs.acs.org/doi/full/10.1021/ct900549r
> > gives 1002 for tip3p and 1015 for charmm tip3p (this is with LJ switched
> > from 1 to 1.2 nm).
> > I have in an old paper 986 for tip3p, with 0.85 nm cut-off plus dispersion
> > correction.
> > These numbers are quite different. I'll do some checking.
> >
> > Berk
>
> Actually, TIP3P and CHARMM TIP3P are supposed to give (almost) identical
> results as they are not exactly different models but different
> implementations of a same model. That is why I'm surprised of the rather
> large difference in the density values observed in my simulations or in
> the JCTC paper.
>
> In CHARMM, the additionnal Lennard-Jones parameters for the TIP3P hydrogen
> are not here to better reproduce water properties but for practical
> purpose - normally, the hydrogens are almost entirely included in the
> oxygen VDW radius. As far as I understand, those parameters are here to
> improve the stability of simulations by limiting artificial close contacts
> between a large negative charge and the point charge on the hydrogen
> atoms.
Note that even if the H's are nearly included in the O's LJ, their dispersion interaction
still works over larger distances, which would indeed lower the density.
Berk
>
> In fact, I was wondering if the Settle algorithm is appropriate to CHARMM
> TIP3P. In the original CHARMM RTF file, there is a SHAKE constraint on the
> H-H distance which is not present in the GROMACS tips3p.itp file. So, a
> better use of the CHARMM TIP3P model in GROMACS would be to add this
> constraint in tips3p.itp and use Shake instead of Settle, isn't it?
>
> >
> >> Date: Tue, 14 Sep 2010 11:00:55 +0200
> >> From: spoel at xray.bmc.uu.se
> >> To: gmx-users at gromacs.org
> >> Subject: Re: [gmx-users] Regular vs. CHARMM TIP3P water model
> >>
> >> On 2010-09-14 10.23, Nicolas SAPAY wrote:
> >> > Hello everybody,
> >> >
> >> > I have many difficulties to reproduce TIP3P simulation results with
> >> CHARMM
> >> > TIP3P. Regular TIP3P gives systematically a lower density than its
> >> CHARMM
> >> > counterpart, independantly from the cutoff for non-bonded
> >> interactions,
> >> > the version of GROMACS (4.5.1 or 4.0.7) or the double precision.
> >> > Regular 961,067 +/-0,756 g/L
> >> > CHARMM 980,860 +/-0,492 g/L
> >> >
> >> > The Enthalpy of vaporization follows a similar scheme:
> >> > regular -39,992 +/-0,021 kJ/mol
> >> > CHARMM -40,665 +/-0,009 kJ/mol
> >>
> >> How about the dispersion correction?
> >> If that is not turned on densities will be too low (in both cases).
> >>
> >> >
> >> > In fact, CHARMM gives results closer to what I should obtain at 300K
> >> and 1
> >> > bar for TIP3P. I suspect an issue with the bond constraints, but I
> >> can't
> >> > locate precisely where it is. Settle parameters are exactly the same
> >> (as
> >> > well as constraints for flexible water). Did someone ever face a
> >> similar
> >> > problem?
> >> >
> >> > Rgards,
> >> > Nicolas
> >> >
> >> > P.S. I My system is just a box of 1728 water molecules
> >> pre-equilibrated
> >> > for 500 ps at 300 K and 1 bar.
> >> >
> >> > P.S. II I'm using the following simulation parameters:
> >> > constraints = hbonds
> >> > constraint_algorithm = LINCS
> >> > continuation = no
> >> > lincs-order = 4
> >> > lincs-iter = 1
> >> > lincs-warnangle = 30
> >> > morse = no
> >> > (LINCS or SHAKE should not make any difference since TIP3P is normally
> >> > treated by SETTLE).
> >> >
> >> > Other parameters are:
> >> > ; COUPLING
> >> > ; Temperature coupling
> >> > tcoupl = Berendsen
> >> > nsttcouple = -1
> >> > nh-chain-length = 10
> >> > tc_grps = System
> >> > tau_t = 0.1
> >> > ref_t = 300
> >> > ; Pressure coupling
> >> > pcoupl = Berendsen
> >> > pcoupltype = isotropic
> >> > nstpcouple = -1
> >> > tau_p = 1
> >> > compressibility = 4.5e-5
> >> > ref_p = 1
> >> >
> >> > ; NEIGHBORSEARCHING PARAMETERS
> >> > ; nblist update frequency
> >> > nstlist = 10
> >> > ns_type = grid
> >> > pbc = xyz
> >> > rlist = 1.15
> >> >
> >> > ; OPTIONS FOR ELECTROSTATICS AND VDW
> >> > ; Method for doing electrostatics
> >> > coulombtype = PME
> >> > rcoulomb-switch = 0
> >> > rcoulomb = 1.15
> >> > vdwtype = Cutoff
> >> > rvdw = 1.0
> >> > fourierspacing = 0.10
> >> > pme_order = 6
> >> > ewald_rtol = 1.0e-6
> >> > ewald_geometry = 3d
> >> > epsilon_surface = 0
> >> > optimize_fft = yes
> >> >
> >> > ; RUN CONTROL PARAMETERS
> >> > integrator = md
> >> > dt = 0.002
> >> > nsteps = 200000
> >> > comm_mode = Linear
> >> > nstcomm = 1
> >> > comm_grps = System
> >> >
> >> >
> >>
> >>
> >> --
> >> David van der Spoel, Ph.D., Professor of Biology
> >> Dept. of Cell & Molec. Biol., Uppsala University.
> >> Box 596, 75124 Uppsala, Sweden. Phone: +46184714205.
> >> spoel at xray.bmc.uu.se http://folding.bmc.uu.se
> >> --
> >> gmx-users mailing list gmx-users at 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-request at gromacs.org.
> >> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> > --
> > gmx-users mailing list gmx-users at 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-request at gromacs.org.
> > Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
>
> --
> [ Nicolas Sapay - Post-Doctoral Fellow ]
> CERMAV - www.cermav.cnrs.fr
> BP53, 38041 Grenoble cedex 9, France
> Phone: +33 (0)4 76 03 76 44/53
>
> --
> gmx-users mailing list gmx-users at 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-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20100914/c603facd/attachment.html>
More information about the gromacs.org_gmx-users
mailing list