[gmx-users] system cooling down when runing NVE
Erik Marklund
erikm at xray.bmc.uu.se
Wed May 9 11:04:20 CEST 2012
In general, NVE is tricky business, since small integration errors build up and cause energy drift in the long run. Do you see energy drift too? In my experience one needs double precision and a shorter timestep than what you employ. But that is of course system dependent.
Erik
9 maj 2012 kl. 03.17 skrev Peter C. Lai:
> you have position restraints on, which I expect would damp collisions between solvent and solute. temp drops towards some sort of equilibrium, which doesn't necessarily match your starting temp even though energy of the system is conserved ..sounds like expected behavior to me?
> --
> Sent from my Android phone with K-9 Mail. Please excuse my brevity.
>
> Thanh Binh NGUYEN <nguyentb at bii.a-star.edu.sg> wrote:
> Dear Gromacs experts,
> > I'm just a newbaby in Gromacs, and hence I have a lot of problem
> > when running this program. I try to run NVE simulation of a
> > protein. First, I run an NTP ensemble, follow by NVT and finally,
> > NVE. In NPT and NVT, T remains as constant, however, in NVE
> > simulation, the temperature drops gradually instead of maintain
> > approximately constant. Could you give me any advice to solve this
> > problem?
> > Thank you very much.
> > Regards,
> > Nguyen
> >
> > P.S: below is my mdp file for three simulations:
> > NPT ensemble:
> > title = DEN2pro_MD
> > define = -DFLEXIBLE
> > constraints = all-bonds
> > integrator = md
> > dt = 0.002 ; ps !
> > nsteps = 1000000 ; total 2000 ps.
> > nstcomm = 1
> > nstxout = 5000 ; output
> coordinates every 10 ps
> > nstvout = 5000 ; output velocities every 10 ps
> > nstfout = 0
> > nstlog = 10
> > nstenergy = 10
> > nstlist = 10
> > ns_type = grid
> > rlist = 0.9
> > coulombtype = PME
> > rcoulomb = 0.9
> > rvdw = 1.0
> > fourierspacing = 0.12
> > fourier_nx = 0
> > fourier_ny = 0
> > fourier_nz = 0
> > pme_order = 6
> > ewald_rtol = 1e-5
> > optimize_fft = yes
> > ; Berendsen temperature coupling is on in four groups
> > Tcoupl = berendsen
> > tau_t = 0.1 0.1
> > tc_grps = protein non-protein
> > ref_t = 300 300
> > ; Pressure coupling is on
> > Pcoupl = berendsen
> > pcoupltype = isotropic
> > tau_p = 1.0
> > compressibility = 4.5e-5
> > ref_p = 1.0
> > ; Generate velocites is on at 300 K.
> > gen_vel = yes
> > gen_temp = 300.0
> > gen_seed = 173529
> >
> ><
> br
> />>
> > NVT ensemble
> > title = DEN2pro_MD NVT equilibration
> > define = -DPOSRES ; position restrain the protein and ligand
> > ; Run parameters
> > integrator = md ; leap-frog integrator
> > nsteps = 100000 ; 2 * 100000 = 200 ps
> > dt = 0.002 ; 2 fs
> > ; Output control
> > nstxout = 5000 ; save coordinates every 10 ps
> > nstvout = 5000 ; save velocities every 10 ps
> > nstenergy = 100 ; save energies every 0.2 ps
> > nstlog = 100 ; update log file every 0.2 ps
> > energygrps = Protein
> > ; Bond parameters
> > continuation = no ; first dynamics run
> > constraint_algorithm = lincs ; holonomic constraints
> > constraints = all-bonds ; all bonds (even heavy atom-H
> > bonds) constrained
> > lincs_iter = 1 ; accuracy of LINCS
> &
> gt;
> lincs_order = 4 ; also related to accuracy
> > ; Neighborsearching
> > ns_type = grid ; search neighboring grid cells
> > nstlist = 5 ; 10 fs
> > rlist = 0.9 ; short-range neighborlist cutoff (in nm)
> > rcoulomb = 0.9 ; short-range electrostatic cutoff (in nm)
> > rvdw = 1.4 ; short-range van der Waals cutoff (in nm)
> > ; Electrostatics
> > coulombtype = PME ; Particle Mesh Ewald for long-range
> > electrostatics
> > pme_order = 4 ; cubic interpolation
> > fourierspacing = 0.12 ; grid spacing for FFT
> > ; Temperature coupling is on
> > tcoupl = V-rescale ; modified Berendsen thermostat
> > tc-grps = Protein Non-protein ; two coupling groups - more accurate
> > tau_t = 0.1 0.1 ; time constant, in ps
> > ref_t =
> 300
> 300 ; reference temperature,
> > one for each group, in K
> > ; Pressure coupling is off
> > pcoupl = no ; no pressure coupling in NVT
> > ; Periodic boundary conditions
> > pbc = xyz ; 3-D PBC
> > ; Dispersion correction
> > DispCorr = EnerPres ; account for cut-off vdW scheme
> > ; Velocity generation
> > gen_vel = yes ; assign velocities from Maxwell distribution
> > gen_temp = 300 ; temperature for Maxwell distribution
> > gen_seed = -1 ; generate a random seed
> >
> >
> > NVE ensemble:
> > title = DEN2pro_MD NVE equilibration
> > define = -DPOSRES ; position restrain the protein and ligand
> > ; Run parameters
> > integrator = md-vv ; leap-frog integrator
> > nsteps = 1000000 ; 2 * 1000000 = 2000 ps
> > dt = 0.002 ; 2 fs
> > ; Out
> put
> control
> > nstxout = 5000 ; save coordinates every 10 ps
> > nstvout = 5000 ; save velocities every 10 ps
> > nstenergy = 100 ; save energies every 0.2 ps
> > nstlog = 100 ; update log file every 0.2 ps
> > energygrps = Protein
> > ; Bond parameters
> > continuation = no ; first dynamics run
> > constraint_algorithm = lincs ; holonomic constraints
> > constraints = all-bonds ; all bonds (even heavy atom-H
> > bonds) constrained
> > lincs_iter = 1 ; accuracy of LINCS
> > lincs_order = 4 ; also related to accuracy
> > ; Neighborsearching
> > ns_type = grid ; search neighboring grid cells
> > nstlist = 5 ; 10 fs
> > rlist = 1.6 ; short-range neighborlist cutoff (in nm)
> > rcoulomb = 0.9 ; short-range electrostatic cutoff (in nm)
> >
> rlistlong = 1.6
> > ;treatment of van der waals interactions
> > vdwtype = Shift
> > rvdw = 0.95 ; short-range van der Waals cutoff (in nm)
> > rvdw-switch = 0.9
> > ; Electrostatics
> > coulombtype = PME-Switch; Particle Mesh Ewald for long-range
> > electrostatics
> > pme_order = 4 ; cubic interpolation
> > fourierspacing = 0.12 ; grid spacing for FFT
> > ; Temperature coupling is on
> > tcoupl = no ; no temperature coupling in NVE
> > ; Pressure coupling is off
> > pcoupl = no ; no pressure coupling in NVE
> > ; Periodic boundary conditions
> > pbc = xyz ; 3-D PBC
> > ; Dispersion correction
> > DispCorr = EnerPres ; account for cut-off vdW scheme
> > ; Velocity generation
> > gen_vel = yes ; assign velocities from Maxwell distribution
> > gen_temp = 3
> 00
> ; temperature for Maxwell distribution
> > gen_seed = -1 ; generate a random seed
> >
> >
>
>
>
> --
> 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
-----------------------------------------------
Erik Marklund, PhD
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596, 75124 Uppsala, Sweden
phone: +46 18 471 6688 fax: +46 18 511 755
erikm at xray.bmc.uu.se
http://www2.icm.uu.se/molbio/elflab/index.html
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20120509/a1aa830f/attachment.html>
More information about the gromacs.org_gmx-users
mailing list