KPROJ: A Program for Unfolding Electronic and Phononic Bands

Jiaxin Chen,Mingxing Chen

2024-10-14

Abstract:We introduce a program named KPROJ that unfolds the electronic and phononic band structure of materials modeled by supercells. The program is based on the $\textit{k}$-projection method, which projects the wavefunction of the supercell onto the ${\textbf{k}}$-points in the Brillouin zone of the artificial primitive cell. It allows for obtaining an effective "local" band structure by performing partial integration over the wavefunctions, e.g., the unfolded band structure with layer-projection for interfaces and the weighted band structure in the vacuum for slabs. The layer projection is accelerated by a scheme that combines the Fast Fourier Transform (FFT) and the inverse FFT algorithms. It is now interfaced with a few first-principles codes based on plane waves such as VASP, Quantum Espresso, and ABINIT. In addition, it also has interfaces with ABACUS, a first-principles simulation package based on numerical atomic basis sets, and PHONOPY, a program for phonon calculations.

Materials Science,Computational Physics

What problem does this paper attempt to address?

### Problems the Paper Attempts to Solve

This paper aims to address the issue of band folding that occurs when modeling materials using supercells. Specifically:

1. **Impact of Band Folding**:
- Band folding can obscure the intrinsic features of a material's band structure. For example, an indirect band gap might appear as a direct band gap if the valence band maximum (VBM) and conduction band minimum (CBM) at different k-points fold into the same k-point.
- As the supercell size increases, the calculated band structure may differ significantly from that of the original unit cell, complicating the understanding of the effects of chemical doping and interfaces on the material's band structure.
- Folded bands can also make it difficult to interpret angle-resolved photoemission spectroscopy (ARPES) experimental results.

2. **Limitations of Existing Methods**:
- Although various band unfolding techniques exist, they still have shortcomings when dealing with complex systems such as doped systems and interfaces.

### Solution

To eliminate the effects of band folding, the authors developed a program called KPROJ, which uses the k-projection method to calculate the unfolded band structure of supercell systems. The specific methods include:

1. **k-Projection Method**:
- By projecting the wave functions of the supercell onto k-points in the artificially defined Brillouin zone (BZ) of the unit cell, the band structure is unfolded.
- Projection operators are constructed using translation operators and their irreducible representations.
- Partial integration of wave functions within a specific spatial window is performed to obtain layer-projected band structures, with fast Fourier transform (FFT) and inverse FFT algorithms used to accelerate calculations.

2. **Program Interface**:
- The KPROJ program has been interfaced with several plane-wave-based first-principles software packages (such as VASP, Quantum Espresso, ABINIT) as well as the numerical atomic basis set first-principles simulation package ABACUS and the phonon calculation program PHONOPY.

### Application Examples

The paper demonstrates the application of the KPROJ program in the following areas:

1. **Graphene Supercell**:
- In a √7×√7 graphene supercell, the band structure obtained from standard DFT calculations is similar to that of the original unit cell near the K-point but differs significantly near the Γ and M points. The unfolded band structure matches that of the original unit cell.

2. **Defective Graphene**:
- By projecting the wave functions of defective graphene onto k-points in the BZ of the corresponding unit cell, non-dispersive bands are eliminated, revealing a band gap at the Dirac cone.

3. **Bilayer Graphene/6H-SiC Interface**:
- Layer-projected calculations yield the "effective" band structure of bilayer graphene, showing a band structure very similar to that of a free-standing system, with only a slight band gap in the electronic bands.

4. **Surface Bulk Projection**:
- By projecting the wave functions of the Ag(111) surface onto k-points in the BZ of bulk Ag, a band structure similar to that of the bulk is obtained, revealing the effect of surface states on band broadening near the Fermi level.

5. **Phonons of Defective Graphene and Moiré Superlattices**:
- The KPROJ program is used to unfold the phonon dispersion relations of defective graphene and a 9° twisted graphene moiré superlattice, revealing the impact of carbon defects on phonon dispersion.

### Conclusion

The KPROJ program effectively unfolds the electronic and phonon band structures of supercell systems and is compatible with various first-principles software packages. It is particularly suitable for studying interface systems and provides an efficient layer-projection algorithm. Additionally, it offers several practical tools to help users prepare input files and post-process data. Through these applications, the KPROJ program helps to reveal the effects of defects and interfaces on the electronic and phonon dispersion of materials.