On 9/28/12 3:18 AM, Felipe Pineda, PhD wrote:
Hi again,

I posted yesterday the query below, but have not received any feedback up to
now. I would like to put my question in another way: in Justin's "CAP-15 in
DPPC" tutorial he uses during NVT equilibration Bussi's thermostat (V-rescale)
together with three separate coupling groups:

; Temperature coupling is on
tcoupl          = V-rescale                 ; modified Berendsen thermostat
tc-grps         = Protein DPPC SOL_CL   ; three coupling groups - more accurate
tau_t           = 0.1   0.1     0.1             ; time constant, in ps
ref_t           = 323   323     323             ; reference temperature, one for
each group, in K),

Is this separate coupling needed at all taking into account that this
temperature coupling method produces the correct coupling? Is it probably better
to use a single coupling constant (system) and so avoid problems using this
artificial correction to deal with the "hot solvent / cold solute" artifact?


Technically the correct way to do it is to use one global thermostat rather than individual ones, but multiple thermostats are used for the reason described. I think it is actually more to do with electrostatics approximations than the thermostat itself (i.e. plain cutoffs leading to poor energy conservation and thus heating of the solvent, which diffuses rapidly) than anything else, so using PME is more robust. Someone please correct me if I'm remembering wrong. I don't recall having seen a test of an inhomogeneous system with a single thermostat, so multiple thermostats remain common practice. There are numerous, extensive discussions in the list archive about these topics so I would encourage you to search around a bit.

Could somebody kindly tell me how large could the coupling constant, tau_t, be?
Is it OK to switch, e.g., from 0.1 to 0.3 ps after T equilibration is reached,
i.e., during the production phase?


This issue is discussed in the reference for the V-rescale thermostat (Bussi et 
al).

-Justin

--
========================================

Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
--
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