Logo of PhononDirect method

PHONON Software


An ab initio program optimizes the structure of the crystalline supercell within constraints imposed by a crystallographic space group. In optimized configuration the forces acting on all atoms of the supercell vanish. Displacing one atom from its equilibium position, one generates non-zero Hellmann-Feynman forces acting on all atoms of the supercell.

The lattice dynamics requires to know the force constants, which are the second derivatives of the potential energy. In the direct method one calculates the force constants from the Hellmann-Feynman forces. The Hellmann-Feynman forces, in turn, are found when a single atom is displaced. In this way all force constants of the supercell can be found.

To obtain the reliable phonon dispersion relations the supercell diameter should be about 8 - 10 angsterms. The magnitude of the force constants beyond that distnce is usually negligible, and in such a case the phonon dispersion curve can be calculated exactly. But even if the supercell size is small, exact phonon frequencies can be obtained for the wave vectors commensurate with the supercell size. Then, the phonon dispersion curves are the symmetry controlled interpolations between the exact points.

For high-symmetry supercells treated by an ab initio program, Phonon calculates the phonon dispersion curves within minutes.
Within the direct method approach and according to the symmetry imposed by the crystal space group Phonon calculates the following items:
See: