RMPCDMD
Functions/Subroutines
single_body.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

program single_body
 Simulate a single colloidal rigid-body particle. More...
 
subroutine flag_particles
 
subroutine change_species
 
subroutine refuel
 
subroutine read_links (filename, n_links, links, links_d)
 
double precision function, dimension(3) get_unit_r ()
 
subroutine compute_cell_wise_max_v
 

Function/Subroutine Documentation

subroutine single_body::change_species ( )

Definition at line 852 of file single_body.f90.

subroutine single_body::compute_cell_wise_max_v ( )

Definition at line 946 of file single_body.f90.

subroutine single_body::flag_particles ( )

Definition at line 825 of file single_body.f90.

double precision function, dimension(3) single_body::get_unit_r ( )

Definition at line 928 of file single_body.f90.

subroutine single_body::read_links ( character(len=*), intent(in)  filename,
integer, intent(out)  n_links,
integer, dimension(:,:), intent(out), allocatable  links,
double precision, dimension(:), intent(out), allocatable  links_d 
)

Definition at line 897 of file single_body.f90.

subroutine single_body::refuel ( )

Definition at line 874 of file single_body.f90.

program single_body ( )

Simulate a single colloidal rigid-body particle.

This simulation models a chemically active colloid particle in either a periodic simulation box or with confinement in the y direction.

The coordinates of the colloid particle's beads must be provided in a H5MD file, as a "fixed-in-time" dataset. The body of the particle can operate as a rigid-body (with either RATTLE or quaternion dynamics) or as an elastic network.

Parameters
Llength of simulation box in the 3 dimensions
rhofluid number density
TTemperature. Used for setting initial velocities and (if enabled) bulk thermostatting.
tauMPCD collision time
alphaangle of collision
probabilityprobability to change A to B upon collision
bulk_raterate of B->A reaction
N_MDnumber MD steps occuring in tau
N_loopnumber of MPCD timesteps
colloid_samplinginterval (in MD steps) of sampling the colloid position and velocity
do_solvent_ioif true (T), a snapshot of the solvent in the final step is dump to the datafile
equilibration_loopsnumber of MPCD steps for equilibration
epsilon_Cinteraction parameter of C sphere with both solvent species (2 elements)
epsilon_Ninteraction parameter of N sphere with both solvent species (2 elements)
data_filenamefilename for input Janus coordinates
data_groupparticles group in the input file
epsilon_colloidinteraction parameter for colloid-colloid interactions
reaction_radiusradius for the reaction around the Janus particle
link_tresholddistance criterion for finding rigid-body links
do_read_linksread link information from a file
links_filefilename for the links data
do_rattleperform RATTLE
do_lennard_jonescompute colloid-colloid Lennard-Jones forces
do_elasticcompute colloid-colloid elastic network forces
elastic_kelastic constant for the elastic network
rattle_pos_toleranceabsolute tolerance for Rattle (position part)
rattle_vel_toleranceabsolute tolerance for Rattle (velocity part)
do_quaternionperform quaternion velocity Verlet
quaternion_tresholdtreshold for the iterative procedure for the quaternion integrator
sigmaradius of the colloidal beads for colloid-solvent interactions
sigma_colloidradius of the colloidal beads for colloid-colloid interactions
polar_r_maxmaximal radius for the polar fields
do_ywalluse a confining potential in the y direction, 9-3 Lennard-Jones
wall_sigmawall LJ sigma
wall_epsilonwall LJ epsilon
wall_shiftwall shift
fluid_wallboundary condition for the fluid

Definition at line 52 of file single_body.f90.