Computing the phonon dispersion and DOS: Difference between revisions
| Line 9: | Line 9: | ||
To plot the phonon dispersion the tag {{TAG | LPHON_DISPERSION}}=.TRUE. must be set in the {{FILE|INCAR}} file. | To plot the phonon dispersion the tag {{TAG | LPHON_DISPERSION}}=.TRUE. must be set in the {{FILE|INCAR}} file. | ||
Additionally, the {{FILE | QPOINTS}} file must be present containing a '''q''' points path for plotting the phonon dispersion or a mesh to compute the phonon density-of-states. | |||
If the material is polar (i.e. has more than one different atom type) the long-range dipole-dipole interaction has to be treated by an Ewald summation. This is achieved by setting {{TAG | LPHON_POLAR}}=.TRUE., supplying the static dielectric tensor ({{TAG | PHON_DIELECTRIC}}) and the Born-effective charges ({{TAG | PHON_BORN_CHARGES}}). The ({{TAG | PHON_BORN_CHARGES}}) have to be obtained from a separate VASP run using the {{TAG|LEPSILON}} tag. Additionally, the user has to set a | If the material is polar (i.e. has more than one different atom type) the long-range dipole-dipole interaction has to be treated by an Ewald summation. This is achieved by setting {{TAG | LPHON_POLAR}}=.TRUE., supplying the static dielectric tensor ({{TAG | PHON_DIELECTRIC}}) and the Born-effective charges ({{TAG | PHON_BORN_CHARGES}}). The ({{TAG | PHON_BORN_CHARGES}}) have to be obtained from a separate VASP run using the {{TAG|LEPSILON}} or {{TAG|LCALCEPS}} tag. | ||
Additionally, the user has to set a reciprocal space cutoff radius ({{TAG | PHON_G_CUTOFF}}) for the Ewald summation. | |||
The amount of information written to the {{FILE | OUTCAR}} file can be tuned using the ({{TAG|PHON_NWRITE}} tag). | The amount of information written to the {{FILE | OUTCAR}} file can be tuned using the ({{TAG|PHON_NWRITE}} tag). | ||
Revision as of 12:37, 11 August 2022
After computing the force constants using the finite differences or density functional perturbation theory approaches it is possible to compute the phonon dispersion using the primitive cell as well as the phonon density of states.
Computing force constants
There are two possible approaches for computing the force constants and then building the dynamical matrix:
- Using finite differences with (IBRION=5, 6)
- Using density functional perturbation theory with (IBRION=7, 8)
Computing the phonon dispersion
To plot the phonon dispersion the tag LPHON_DISPERSION=.TRUE. must be set in the INCAR file. Additionally, the QPOINTS file must be present containing a q points path for plotting the phonon dispersion or a mesh to compute the phonon density-of-states.
If the material is polar (i.e. has more than one different atom type) the long-range dipole-dipole interaction has to be treated by an Ewald summation. This is achieved by setting LPHON_POLAR=.TRUE., supplying the static dielectric tensor (PHON_DIELECTRIC) and the Born-effective charges (PHON_BORN_CHARGES). The (PHON_BORN_CHARGES) have to be obtained from a separate VASP run using the LEPSILON or LCALCEPS tag. Additionally, the user has to set a reciprocal space cutoff radius (PHON_G_CUTOFF) for the Ewald summation. The amount of information written to the OUTCAR file can be tuned using the (PHON_NWRITE tag).
Computing the density of states
To compute the phonon density of states a sufficiently dense q-point mesh should be specified in the QPOINTS file. The density of states is computed between with and the lowest and highest phonon frequency and the broadening PHON_SIGMA.
![{\displaystyle [\omega _{\text{min}}-5\sigma ,\omega _{\text{max}}+5\sigma ]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2ad6fcb4675ac638da1a625bc8cd0423c791c735)
The number of energy points in this energy range is specified by the PHON_NEDOS tag. To use a gaussian smearing method for the computation of the DOS set PHON_DOS=1 to use the tetrahedron method set PHON_DOS=2.
Related tags and articles
QPOINTS, LPHON_DISPERSION, PHON_NWRITE, LPHON_POLAR, PHON_DIELECTRIC, PHON_BORN_CHARGES, PHON_G_CUTOFF