VASP2KP: kp models and Lande g-factors from ab initio calculations

Sheng Zhang,Haohao Sheng,Zhi-Da Song,Chenhao Liang,Yi Jiang,Song Sun,Quansheng Wu,Hongming Weng,Zhong Fang,Xi Dai,Zhijun Wang
DOI: https://doi.org/10.1088/0256-307X/40/12/127101
2023-12-14
Abstract:The $k\cdot p$ method is significant in condensed matter physics for the compact and analytical Hamiltonian. In the presence of magnetic field, it is described by the effective Zeeman's coupling Hamiltonian with Landé $ g $-factors. Here, we develop an open-source package VASP2KP (including two parts: vasp2mat and mat2kp) to compute $k\cdot p$ parameters and Landé $g$-factors directly from the wavefunctions provided by the density functional theory (DFT) as implemented in Vienna ab initio Simulation Package (VASP). First, we develop a VASP patch vasp2mat to compute matrix representations of the generalized momentum operator $ \mathbf{\hat{\pi}}=\mathbf{\hat{p}}+\frac{1}{2mc^2}\left(\mathbf{\hat{s}}\times\nabla V(\mathbf{r})\right) $, spin operator $\mathbf{\hat{s}}$, time reversal operator $\hat{T}$ and crystalline symmetry operators $\hat{R}$ on the DFT wavefunctions. Second, we develop a python code mat2kp to obtain the unitary transformation $U$ that rotates the degenerate DFT basis towards the standard basis, and then automatically compute the $k\cdot p$ parameters and $g$-factors. The theory and the methodology behind VASP2KP are described in detail. The matrix elements of the operators are derived comprehensively and computed correctly within the projector augmented wave method. We apply this package to some materials, e.g., Bi$_2$Se$_3$, Na$_3$Bi, Te, InAs and 1H-TMD monolayers. The obtained effective model's dispersions are in good agreement with the DFT data around the specific wave vector, and the $g$-factors are consistent with experimental data. The VASP2KP package is available at <a class="link-external link-https" href="https://github.com/zjwang11/VASP2KP" rel="external noopener nofollow">this https URL</a>.
Materials Science
What problem does this paper attempt to address?
The paper mainly addresses the problem of how to directly obtain the k·p model parameters and Landé g factor from first-principles calculations (such as the VASP software). In condensed matter physics, the k·p method is used to simplify the electronic band structure, particularly in the vicinity of specific wave vectors, while the Landé g factor is related to the magnetic properties of materials in a magnetic field. Although existing methods can calculate these parameters, there are no open-source codes directly extracted from the DFT wave functions. The author developed an open-source software package called VASP2KP, which consists of two parts: vasp2mat (a patch for VASP) and mat2kp (a post-processing Python code). vasp2mat can calculate the matrix representations of the general momentum, spin, time-reversal, and crystal symmetry operators, while mat2kp is responsible for the unitary transformation from the degenerate states in DFT ground state to the k·p standard basis and automatically calculates the k·p parameters and g factor. With this software package, researchers can more conveniently construct effective models to describe the low-energy electron behavior of materials, and the obtained results are in good agreement with DFT data and experimental data. This work provides a powerful tool for understanding and designing materials with specific electronic properties, especially those exhibiting behavior under a magnetic field.