[gmx-users] Barostat Selection

Justin Lemkul jalemkul at vt.edu
Wed Sep 30 03:10:28 CEST 2015



On 9/29/15 9:06 PM, Parker de Waal wrote:
> Thank you for the prompt response Justin.
>
> I was unaware that the relaxation time had no effect on simulation performance. Would you expect a difference between tau_p = 1.0 and 5.0 for production runs using Parrinello-Rahman?
>

Expect a difference with respect to what?  All you're doing is changing how 
quickly the barostat responds to changes in pressure.  Either value is probably 
fine for most purposes.

-Justin

> -----Original Message-----
> From: gromacs.org_gmx-users-bounces at maillist.sys.kth.se [mailto:gromacs.org_gmx-users-bounces at maillist.sys.kth.se] On Behalf Of Justin Lemkul
> Sent: Tuesday, September 29, 2015 9:00 PM
> To: gmx-users at gromacs.org
> Subject: Re: [gmx-users] Barostat Selection
>
>
>
> On 9/29/15 8:57 PM, Parker de Waal wrote:
>> Hi everyone,
>>
>> I recently got into an interesting discussion regarding barostat selection for NPT equilibration and later NPT production runs and was wondering what the community thought on this issue.
>>
>> Specifically, for the initial NPT equilibrium, would it be acceptable
>> to use the berendsen barostat (tau_p = 1.0 ) and then later switch to
>> the Parrinello-Rahman barostat (tau_p = 5) for production, or would it
>> be best to stick with berendsen (keep tau_p = 1.0 or increase to 5.0
>> to increase ns/day?)
>>
>
> Barostat relaxation time should have no effect on performance.  The Berendsen weak coupling algorithms do not produce correct distributions and they should not be used for production simulations.  The algorithm is good for equilibration, but shouldn't be used for much else.
>
> -Justin
>
> --
> ==================================================
>
> Justin A. Lemkul, Ph.D.
> Ruth L. Kirschstein NRSA Postdoctoral Fellow
>
> Department of Pharmaceutical Sciences
> School of Pharmacy
> Health Sciences Facility II, Room 629
> University of Maryland, Baltimore
> 20 Penn St.
> Baltimore, MD 21201
>
> jalemkul at outerbanks.umaryland.edu | (410) 706-7441 http://secure-web.cisco.com/1ojLfbMlRYuefZ_1olbKjXPZp-h4B9atVOLOfNXQ7pczhcSv3VkmPlIgK3dmbOVGiGwMhCTCoZzv_jrbSSl3rJC35FBF34au6FR1lWqNwjUmLk1i-p8rHhXfjXJ8RiLsvoVuIxKtlum35oYjAa7eqmHBr2xfQ0a2cH-V3FI4YDUtttitIzcF1c177lFZnkxF9NEeJigdBkiDZVBjuINkIF7gVdNGETE1r8Eu6XjM24e4/http%3A%2F%2Fmackerell.umaryland.edu%2F%7Ejalemkul
>
> ==================================================
> --
> Gromacs Users mailing list
>
> * Please search the archive at http://secure-web.cisco.com/1MQxYm8ywdBbJUNmGzQcwWfpEQQxdnDq282IW0mysWDnGAQ0I8EMMJWQZMUVJ_74e6dqtE8O_YjjWgQzqOhP4V1isGBx8iJMiWzNMyxP84aObzeqLlyE78aJsMOl0z8fWj2pqya_FEuaToFHEcxZ3FJ9_U25-Aoag24zDNoTIllYrX6geCNO9PdaK1dLWKA2sMLZMZ58auzoI1qPy7uhviD25PX_PNrrI_XMZ0u2MBjc/http%3A%2F%2Fwww.gromacs.org%2FSupport%2FMailing_Lists%2FGMX-Users_List before posting!
>
> * Can't post? Read http://secure-web.cisco.com/1yt_VioDuTK5YDl9Q_jRUvWaeaLxVnxvdzuVuu9PFL0lNtvesrLx2OSpmohlMSWw_fb9QKrUp6QxYH-LoV5VtUO5caLkKdWpihSLYQbK5FmyRVVzsVYJhiT7OYx6wTjeYzKUw2Z0SpWtDnYsbGqql11QL7qT39jmxhd6xHd1bIhD7yIThKVXHcI4cJrOjqBTuMEDWLm7Xd5eLzGWHr_n58K9BQchkE3L4ymXEfoUa4c0/http%3A%2F%2Fwww.gromacs.org%2FSupport%2FMailing_Lists
>
> * For (un)subscribe requests visit
> https://secure-web.cisco.com/1KMOCw8Kr5dUhkR-1S8WcbMa7DVZxUYsntzcAAfExbboYW-CEhCU2LLYhQsa4mzTaQ8FWv5KSJAAWxXI80xIwMUu8GDN-5InZw5iphCzhDvwjNOskv_FMTQuhJygoheWDld5UCSfx2SDnuSy-JhDkJZxqu-NfrSvk7SHwV9SCoicC62ET08z1LuzXERDyQ6D32dBCcWywTloaOKZwL_c5dR1JfmhryH_VVGr7GZi9O5E/https%3A%2F%2Fmaillist.sys.kth.se%2Fmailman%2Flistinfo%2Fgromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.
>

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

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalemkul at outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==================================================


More information about the gromacs.org_gmx-users mailing list