MD-kinetics

This is a short tutorial on how to use MD-kinetics program.


General description

In order to run the program please execute the following command in the terminal.

Where input.inp  is the input file.

Program produces the following files:

  • cluster.xyz – XYZ trajectory file with atomic positions on each time step.
  • kinetics.dat – file with the evolution of populations of atomic states and free agents in time.
  • restart.dat – file with cluster and free agents after propagation. It can be used to run propagation on later timescale.

Moreover, the program prints input and reaction analysis to the output stream.

Available releases

  • v1.0 are available upon request.

Input

In order to run the MD-kinetics program one needs to provide it with the input file. The following section should be specified in the input.

  • propagation section (required)

This section specifies propagation parameters of the integrators.

  • potentials section (optional)

This section specifies interaction potentials between atomic states of the cluster.

  • reactions section (optional)

This section specifies chemical reactions between states of atoms. Involved atoms are specified by the following keyword  $atom{A:4s} indicating that 4s state of atom A participate in the reaction. When reaction occurs between different atoms, one needs to specify another atom by different id (for example B). Reaction also can include free agents  $agent{e_7.6188} which do not belong to any atom. Last number in reaction line is the rate constant.

  • cluster section (required)

This section specifies atomic cluster. Each line represents atom which has mass, position, velocity and some atomic states with corresponding populations.

  • agents section (optional)

This section specifies initial concentrations of free reaction agents which do not belong to any atom.