[gmx-developers] Velocity Verlet Algorithm
Berk Hess
hess at kth.se
Wed Nov 9 14:46:41 CET 2016
Hi,
There is nothing special with the VV integrator for any parameters. The
only thing is that it calculates the kinetic energy from the full step
instead of the half step kinetic energy (and it does an extra constraint
call for the velocities).
Cheers,
Berk
On 09/11/16 10:05 , Elena della Valle wrote:
> Dear All,
>
> I'm writing to ask you about velocity verlet integrator.
>
> I need to use md-vv because I modified in the file "update.c" the
> do_update_vv_vel and do_update_vv_pos in order to introduce the larmor
> frequency therm for simulating a magnetic field. I'm doing simulation
> with a ion in vacuum and a box with water molecules and one CL ion. I
> hope that I am clear.
>
> When I need to use as integrator the "md-vv", in the .mdp file there
> are boundary conditions that I have to take into account? If I'm doing
> pme silmulations I need to change some parameters?
>
> Thanks in advance
>
> Any suggestions will be really appreciated
>
> Best Regards
>
> Elena della Valle
>
>
>
More information about the gromacs.org_gmx-developers
mailing list