fix reaxff/bonds command

Accelerator Variants: reaxff/bonds/kk


fix ID group-ID reaxff/bonds Nevery filename
  • ID, group-ID are documented in fix command

  • reax/bonds = style name of this fix command

  • Nevery = output interval in timesteps

  • filename = name of output file


fix 1 all reaxff/bonds 100 bonds.reaxff


Write out the bond information computed by the ReaxFF potential specified by pair_style reaxff in the exact same format as the original stand-alone ReaxFF code of Adri van Duin. The bond information is written to filename on timesteps that are multiples of Nevery, including timestep 0. For time-averaged chemical species analysis, please see the fix reaxff/species command.

The specified group-ID is ignored by this fix.

The format of the output file should be reasonably self-explanatory. The meaning of the column header abbreviations is as follows:

  • id = atom id

  • type = atom type

  • nb = number of bonds

  • id_1 = atom id of first bond

  • id_nb = atom id of Nth bond

  • mol = molecule id

  • bo_1 = bond order of first bond

  • bo_nb = bond order of Nth bond

  • abo = atom bond order (sum of all bonds)

  • nlp = number of lone pairs

  • q = atomic charge

If the filename ends with “.gz”, the output file is written in gzipped format. A gzipped dump file will be about 3x smaller than the text version, but will also take longer to write.






See the Speed packages page for more instructions on how to use the accelerated styles effectively.


The fix reaxff/bonds command requires that the pair_style reaxff is invoked. This fix is part of the REAXFF package. It is only enabled if LAMMPS was built with that package. See the Build package page for more info.

To write gzipped bond files, you must compile LAMMPS with the -DLAMMPS_GZIP option.