Document Type


Publication Date



ThermoPhonon is a stand-alone code, which can be integrated into other software packages. Typically, it is used together with a density functional theory (DFT) code (such as VASP, Wien2k, AbInit, SIESTA) and a phonon code (such as Phonopy or Phon). The workflow is the following.

  • Molecular dynamics (MD) in a supercell at a given temperature T is performed using another code. After sufficient equilibration, the output in the form of atomic positions and forces for a large number of selected MD steps is recorded into a file. If needed, one can modify this file by applying additional constraints, such as enforced crystal symmetry or subtracted motion of the center of mass.
  • ThermoPhonon reads the file with atomic positions and forces and writes a new file with the force constants.
  • Force constants can be used by another code (such as Phonopy or Phon) to produce phonon spectrum for plotting, in the assumption of known equilibrium atomic positions provided in a separate file.

Separate utilities for file conversion are not included into the ThermoPhonon code and not covered by its license and copyright. Compiling instructions are in the source code. The attachment consists of 4 files:

  • ThermoPhonon source code (Fortran),
  • ThermoPhonon manual (PDF),
  • a truncated sample input file “force.txt” (text),
  • a sample output file “FORCE_CONSTANTS” (text).

Copyright Owner

Nikolai Zarkevich, the Ames Laboratory, Iowa State University


Author Affiliations: Ames Laboratory