OTF-MD

This is a short manual on how to use the On-The-Fly Molecular Dynamics program.

The program is intended to perform classical dynamics simulations of molecular systems using potential energy surface and gradients generated by external software on-the-fly.


The source code is available at the bitbucket.org. In order to download repository you should be included in the list of developers.


Input file description

The OTF-MD program uses the input file to specify parameters of calculations. Here we present an example of input file with description of parameter which can be modified by user.

In order to simplify operations with variables in the input file, one can use the following template variables which will be substituted during program run.

  • $RUN_SCRIPT – Run script file name with path. Can be useful in command string.
  • $XYZ_OUT – XYZ output file on current propagation step. This file is assumed to be used in command string in order to provide an external script with the current geometry.
  • $VEL_OUT – File with velocities on current propagation step.
  • $EG_OUT – File with energy and gradient produced by external script.
  • $TMP_DIR – Temporary directory variable.
  • $TSTEP – Current time step.

Geometry files

The OTF-MD program reads geometry from two input files: xyz_file and mol_file.

The XYZ file specifies the molecule geometry by giving the number of atoms with Cartesian coordinates that will be read on the first line, a comment on the second, and the lines of atomic coordinates in the following lines. Atomic coordinates are in ångströms.

As an example here we present the content of XYZ file for the propiolic acid molecule.

The MOL file specifies atomic numbers and atomic masses of atoms. The first column must contain the same atomic labels as specified in XYZ file. Second and third columns contain atomic numbers and atomic masses, respectively.

The VEL file has the same structure as XYZ file. However, instead of Cartesian coordinates it contains velocities written in [\angstrom/a.u. of time] units.


External script to calculate energy and gradient

The script should take molecular geometry in XYZ format. For this provided geometry, the script should evaluate energy and gradient. Output should be formatted as following.


Analysis programs

Analysis of the molecular dynamics trajectory calculated with OTF-MD program can be performed with programs available in utils/ directory.

  • calc_energies : calculates kinetic and potential energies of the system during trajectory. The third column are the total energy of the system. In case of absence of external fields the total energy of the system must be conserved.
    In order to run the program, execute the following command. The input file is the same which was used for propagation.
  • create_xyz_anim : merges XYZ files into one. The output can be used to visualize XYZ trajectory. The input file is the same which was used for propagation.
  • create_TGA_traj : creates trajectory files for analysis by TGA program.