Hello all,

can someone give me a hint, please?
Do you need more information?
Has anyone had a similar problem.

I really need to figure that out, because I will simulate other systems
which can be run only with user tables.

Anja



-
> ----- Original Message -----
> ----- Original Message -----
> From Anja Kuhnhold <anja.kuhnh...@physik.uni-halle.de>
> Date Wed, 06 Jun 2012 15:07:40 +0200
> To gmx-users@gromacs.org
> Subject [gmx-users] energy conservation: shift vs shifted user potential
> Dear gmx-users,
> 
> I have a problem concerning energy conservation when using user 
> potentials (tables).
> I'm using gromacs 4.5.4
> I simulate a simple Lennard-Jones(6-12) +Fene polymer melt (1600 
> chains a 10 beads in a 26x26x26 periodic box).
> 
> I tried different vdwtypes (cutoff always 3.24):
> The cut-off version does not conserve energy -- okay.
> The shift and switch versions conserve energy -- fine.
> 
> Now I wanted to do the same with user tables:
> Simple Lennard-Jones table gives really the same results as the 
> cut-off version -- good.
> 
> But if I use a table with shifted Lennard-Jones potential it is not 
> comparable to the shift version
> and the energy is not conserved -- ?
> 
> I use a shift function as written in the manual (chapter 4.1.5) -- 
> there must be a factor alpha added in the constants A and B --
> (r1=0).
> 
> The parameters are the same for shift version and shifted user version.
> 
> Has someone an idea why the shifted user potential doesn't work in 
> this way?
> 
> Here is the mdp:
> 
> 
> integrator            = md-vv
> dt                    = 0.0035
> nsteps                        = 1000
> nstxout                       = 1
> nstvout                       = 1
> nstfout                       = 1
> nstlog                        = 1
> ns_type                       = grid
> pbc                   = xyz
> rvdw                  = 3.24
> rlist                 = 3.6
> tcoupl                        = no
> tc-grps                       = System
> tau_t                 = 2.0
> ref_t                 = 127.2717
> vdwtype                       = user;Shift
> rcoulomb              = 3.6;2.24;1.12
> coulombtype           = Cut-off
> rvdw-switch           = 0.0
> 
> energygrps            = bead
> energygrp_table               = bead bead
> 
> 
> Thanks in advance
> Anja
> 

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