[gmx-users] Queries regarding temperature gradient
s.neelima531 at gmail.com
Fri Jul 6 22:30:55 CEST 2018
Dear Dr. Lemkul,
Thank you for your answer. Can you please also advise how can I get the
temperature of a system as a function of distance using gromacs 5.1.5? gmx
energy plots the temperature of the centre of mass of energy groups against
time. gmx traj also does the same from velocity information.
On Sat, 7 Jul 2018, 00:01 Justin Lemkul, <jalemkul at vt.edu> wrote:
> On 7/6/18 1:07 PM, Sabreen Farnaz wrote:
> > Dear all,
> > This is in regard to my previous queries. The warning in NVT step can
> > possibly be ignored using -maxwarn but I am still not sure about the
> > Can someone please help?
> tau_t = -1 means temperature coupling is not applied. Whether that's
> what you want or not is up to you, but in that case, your ensembles are
> not NVT or NPT, because T is no longer conserved.
> > Thanks,
> > Neelima S.
> > On Fri, 6 Jul 2018, 01:56 Sabreen Farnaz, <s.neelima531 at gmail.com>
> >> Dear all,
> >> I am setting up a system where I want the solute to be immersed in a
> >> solvent with a temperature gradient. I have defined two different
> >> temperature groups, two portions of the solvent, and want them to be at
> >> T1 and T2 kelvins. The rest of the system will not be subjected to a
> >> thermostat (tau_t=-1).
> Realize that with PBC, you're not likely to get what you want. I don't
> know how you're defining the T1 and T2 groups, but if they're in contact
> via PBC, the "gradient" really doesn't exist - you have a hot and cold
> bath directly in contact. See previous threads on temperature gradients
> for more discussion. If the T1 and T2 groups are capable of diffusing,
> then you definitely don't have a gradient.
> >> My questions are:
> >> 1. Should NVT and NPT mdp files also have the same t_ref and tau_t? A
> >> tau_t=-1 gives a warning in nvt grompp, can it be ignored?
> >> 2. I am generating velocities in nvt run and continuing npt with
> >> velocities. For production run, I only take the equilibrated structure
> >> generate velocities again. (No continuation)
> If you're re-generating velocities, what was the point of the previous
> equilibration steps? You're just re-setting your system to a random state.
> >> Should gen_temp have the same temperature values in each of nvt and md
> >> files?
> If you intend to model the same ensemble, yes, but again re-generating
> velocities after equilibrating doesn't make much sense to me.
> >> Please help.
> >> Thank you!
> >> Regards,
> >> Neelima S.
> Justin A. Lemkul, Ph.D.
> Assistant Professor
> Virginia Tech Department of Biochemistry
> 303 Engel Hall
> 340 West Campus Dr.
> Blacksburg, VA 24061
> jalemkul at vt.edu | (540) 231-3129
> Gromacs Users mailing list
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> * 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.
More information about the gromacs.org_gmx-users