27 #include "../base/base_uses.f90"
32 CHARACTER(len=*),
PARAMETER,
PRIVATE :: moduleN =
'velocity_verlet_control'
47 TYPE(md_environment_type),
POINTER :: md_env
48 TYPE(global_environment_type),
POINTER :: globenv
50 CHARACTER(LEN=*),
PARAMETER :: routinen =
'velocity_verlet'
53 TYPE(force_env_type),
POINTER :: force_env
54 TYPE(simpar_type),
POINTER :: simpar
56 CALL timeset(routinen, handle)
59 CALL get_md_env(md_env, force_env=force_env, simpar=simpar)
62 IF (simpar%do_respa .AND.
nve_ensemble .NE. simpar%ensemble)
THEN
63 cpabort(
"RESPA integrator not implemented for this ensemble")
67 SELECT CASE (simpar%ensemble)
69 cpabort(
"Integrator not implemented")
71 IF (simpar%do_respa)
THEN
74 CALL nve(md_env, globenv)
77 CALL nvt(md_env, globenv)
83 CALL npt_i(md_env, globenv)
85 CALL npt_f(md_env, globenv)
95 CALL npt_f(md_env, globenv)
Interface for the force calculations.
Define type storing the global information of a run. Keep the amount of stored data small....
Provides integrator routines (velocity verlet) for all the ensemble types.
subroutine, public nvt(md_env, globenv)
nvt integrator for particle positions & momenta
subroutine, public isokin(md_env)
simplest version of the isokinetic gaussian thermostat
subroutine, public reftraj(md_env)
uses coordinates in a file and generates frame after frame of these
subroutine, public nph_uniaxial(md_env)
nph_uniaxial integrator (non-Hamiltonian version) for particle positions & momenta undergoing uniaxia...
subroutine, public nph_uniaxial_damped(md_env)
nph_uniaxial integrator (non-Hamiltonian version) for particle positions & momenta undergoing uniaxia...
subroutine, public langevin(md_env)
Langevin integrator for particle positions & momenta (Brownian dynamics)
subroutine, public npt_f(md_env, globenv)
Velocity Verlet integrator for the NPT ensemble with fully flexible cell.
subroutine, public nve_respa(md_env)
RESPA integrator for nve ensemble for particle positions & momenta.
subroutine, public nvt_adiabatic(md_env, globenv)
nvt adiabatic integrator for particle positions & momenta
subroutine, public nve(md_env, globenv)
nve integrator for particle positions & momenta
subroutine, public npt_i(md_env, globenv)
npt_i integrator for particle positions & momenta isotropic box changes
subroutine, public get_md_env(md_env, itimes, constant, used_time, cell, simpar, npt, force_env, para_env, reftraj, t, init, first_time, fe_env, thermostats, barostat, thermostat_coeff, thermostat_part, thermostat_shell, thermostat_baro, thermostat_fast, thermostat_slow, md_ener, averages, thermal_regions, ehrenfest_md)
get components of MD environment type
Type for storing MD parameters.
Provides an interface to the velocity-verlet based integrator routines for all ensembles.
subroutine, public velocity_verlet(md_env, globenv)
...