how to get the irreducible q-points for phonon calculations

Some of the most commonly heard questions in the Abinit Community, along with answers. (Read only)

how to get the irreducible q-points for phonon calculations

Postby mverstra » Wed Jul 28, 2010 5:19 pm

When doing a phonon calculation, only the irreducible perturbations need to be calculated, and abinit and anaddb will complete them by symmetry. In particular, the dynamical matrix will be completed and symmetrized at each q-point.

To get the irreducible q-points for your system, you can do a small auxiliary calculation (with low ecut and minimal bands). Make a new directory, and copy your input file for the ground state there. If you set the k-point grid (ngkpt) to be equal to the q-point grid you want for the phonons, and set kptopt to 1. Also, don't forget to set <b>shiftk</b> to 0 0 0, in order to get a non-shifted grid with the Gamma point.

From the header of the log or output file, the <b>kpt</b> points can be copied directly to the phonon input file, as q-points.
Matthieu Verstraete
University of Liege, Belgium
Posts: 637
Joined: Wed Aug 19, 2009 12:01 pm

Return to Frequently Asked Questions

Who is online

Users browsing this forum: No registered users and 1 guest