[gmx-users] question regarding high temperature simulation in Gromacs

[Mark Abraham]

Hi,

I don't know, sorry. This is a critical part of your study design, and I
can't do it for you. Look in the literature of its use for any at high
temperature that was actually validated. Or see if e.g. the density and
diffusion rate of your solvent is represented accurately.

Mark

On Fri, Jul 6, 2018 at 1:35 PM Dr. Puspita Halder <puspita at cse.iitkgp.ac.in>
wrote:

> Hi Mark,
>
> I am using OPLS-AA force field for simulating my protein at high
> temperature under nvt condition. Does this work good for high temperature
> system under nvt? Can you suggest which force field I should use then?
>
> Thanks
>
> Puspita
>
> ----- Original Message -----
> From: "Mark Abraham" <mark.j.abraham at gmail.com>
> To: "gmx-users" <gmx-users at gromacs.org>
> Cc: "gromacs org gmx-users" <gromacs.org_gmx-users at maillist.sys.kth.se>
> Sent: Friday, July 6, 2018 1:40:10 PM
> Subject: Re: [gmx-users] question regarding high temperature simulation
> in      Gromacs
>
> Hi,
>
> You should expect that when you heat something up under constant volume
> that its pressure increases. Average velocity has gone up, so the momentum
> transferred in collisions is also up.
>
> The more important question is what evidence you have that this force field
> provides a useful simulation at this point of V,T?
>
> Mark
>
> On Fri, Jul 6, 2018 at 9:43 AM Dr. Puspita Halder <
> puspita at cse.iitkgp.ac.in>
> wrote:
>
> > Hi All,
> >
> > Thanks for your helpful suggestions. I'd like to share the method that I
> > am following now for simulating my protein at high temperature. First, I
> > carried out energy minimization run for my protein system using steepest
> > decent method. Then I performed nvt equilibration run for 100-200 ps at
> > 300K followed by npt equilibration at the same temperature and at 1atm
> > pressure for 1-2 ns with protein position restrain. Then I used simulated
> > annealing protocol (without position restrain and under nvt condition)
> for
> > reaching the temperature of the system to 500K starting from the
> structure
> > of 300 K npt simulation with 50K temperature increment each time followed
> > by 2-4 ns of equilibration run at respective temperature. Finally I did
> the
> > production run at 500 K for 30-40 ns under nvt condition. Do you think
> that
> > the protocol I am using is ok ? The problem is I am getting high value of
> > pressure average (~3000-4000) for the simulated annealing run. Please
> share
> > your comments or suggestion
> >  regarding this.
> >
> > Thanks for your help.
> >
> > Puspita Halder
> >
> > ----- Original Message -----
> > From: "Dr. Puspita Halder" <puspita at cse.iitkgp.ac.in>
> > To: "gromacs org gmx-users" <gromacs.org_gmx-users at maillist.sys.kth.se>
> > Sent: Wednesday, July 4, 2018 12:41:56 PM
> > Subject: question regarding high temperature simulation in Gromacs
> >
> > Hi All,
> >
> > I have been recently using Gromacs 5.1.4 version for simulating my
> protein
> > systems (mainly prion protein and some of its mutants). I'd like to
> perform
> > high temperature (400K or 500K) simulations with those. Now my question
> is
> > what should be the values of the simulation parameters e.g.,
> > compressibility, tau_p or tau_t at 400K or 500K? I used compressibility =
> > 4.5e-5 and tau_p = 2.0 and tau_t = 0.1 for pressure and temperature
> > coupling for simulations at 300K. What other parameters should I change
> for
> > high temperature simulations? Should I always use simulated annealing
> > protocol for simulations at higher temperature or I can directly heat the
> > system to such higher temperatures? Any comments or suggestions in this
> > regard will be highly appreciated.
> >
> > Thanks for your help in advance.
> >
> > Regards
> > Puspita Halder
> > --
> > Gromacs Users mailing list
> >
> > * Please search the archive at
> > http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> > posting!
> >
> > * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> >
> > * For (un)subscribe requests visit
> > https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> > send a mail to gmx-users-request at gromacs.org.
> >
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-request at gromacs.org.
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-request at gromacs.org.
>