This is a question of topological insulator.

Liang Fu and C. L. Kane proposed a method to judge whether an inversion symmetric insulator is a topological insulator or not in their article(L. Fu and C.L. Kane, Phys. Rev. B 76, 045302 (2007)). The method is just to determine the parity of the occupied band eigenstates at the eight or four(in two dimensions) time-reversal invariant momenta $\Gamma_i$in the Brillouin zone. The Z2 invariant is determined by quantity $${{\delta }_{i}}=\prod\limits_{m=1}^{N}{{{\xi }_{2m}}\left( {{\Gamma }_{i}} \right)}$$Where ${{\xi _{2m}}\left( {{\Gamma _i}} \right)} $is the parity eigenvalue of 2m band at $\Gamma_i$ point.

My question is how to determine the parity of band state at these points from first principle band calculation(like Wien2K band calculation)?

  • $\begingroup$ You would normally write down the Bloch Hamiltonian $\mathcal{H}$ as a matrix in a convenient orbital and/or spin basis. Within this basis you know that $\mathcal{T}\left|\mathbf{k},\uparrow\right\rangle =\left|-\mathbf{k},\downarrow\right\rangle$. You can construct a matrix, in this basis, which satisfies that (and obviously $\left[\mathcal{H},\mathcal{T}\right]=0$). Then you can find the eigenvalues of the orbitals at the TRIMs. $\endgroup$
    – NanoPhys
    Sep 7, 2013 at 16:23
  • $\begingroup$ @NanoPhys You mean something like slater-koster method to build the Hamiltonian directly? But, Is there a way to determine the parity from the first principle data(like Wien2K) directly? From this reference"Nature Physics 5(6): 438-442",according to their description I think they can directly determine the parity . $\endgroup$
    – user15964
    Sep 8, 2013 at 2:00
  • $\begingroup$ In the literature, that I have come across so far, people construct a model Hamiltonian $\mathcal{H}(\mathbf{k})$ which respects the $\mathcal{T}$-, rotational, inversion symmetries etc. with adjustable parameters. They diagonalize the Hamiltonian and fit the spectrum to first-principles calculations or experimental data to determine the values of these adjustable parameters; in your case the former. For example, you can look at the BHZ model in eq. (6) of arxiv.org/abs/0801.0901. The values of the parameters, on page 29, are computed from more rigorous calculations. $\endgroup$
    – NanoPhys
    Sep 8, 2013 at 9:49

1 Answer 1


You should use parity analysis in wien2k by group theoretical considerations. You must first converge your material then run "x irrep -so" for spin orbit inclusion or without it. You should just notice to use case.vector file from the k-path of the band to have a comparison between case.band.agr and case.outputir[so] or case.irrep. don't forget to use compatibility relation between different k-point character tables. Acually you should select entangled band to investigate thier topological nature.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.