Spin-orbit splitting of semiconductor band structures¶
The spin-orbit coupling (SOC) is a relativistic effect that causes splitting of the electronic bands in many materials, including semiconductors. This effect is not included in standard calculations.
In this tutorial, you will learn how to take spin-orbit coupling into account in density functional theory calculations (ATK-DFT) and also in semi-empirical Slater-Koster calculations (ATK-SE) for semiconductors.
- We start off with a short introduction to relativistic effetcs in electronic structure theory, including spin-orbit coupling.
- You will then use ATK-DFT to study how the SOC splits the electronic bands of Silicon around the \(\Gamma\)-point, leading to the so-called “split-off” valence band and bands with “heavy” and “light” holes.
- You will then employ a SO+MGGA method to get not only the correct band splits but also a reliable estimate of the indirect band gap.
- Lastly, you will consider the gallium-arsenide crystal, and use the Slater–Koster method with SOC parameters to compute the GaAs band structure.
Relavistic effects in Kohn-Sham DFT¶
The standard time-independent Kohn-Sham Hamiltonian describes non-relativistic electrons moving in an external field set up by atomic nuclei (and possibly any other external time-independent fields). Relavistic effects are therefore completely disregarded. This is usually a very good approximation for valence electrons outside the atomic cores, but for heavy elements, like gold and lead, relativistic contributions to the electronic structure can be crucial. Moreover, spin-orbit coupling, which cannot be captured in a strictly non-relativistic description, will often break degeneracy in solid band dispersions, and thereby lead to experimentally observed band splittings.
Electronic-structure software employing pseudopotentials, like the ATK-DFT engine does, will usually incorporate relativistic effects on core electrons in an approximate manner by using scalar-relativistic pseudopotentials. This is a computationally efficient and quite reliable approach.
However, calculations including spin-orbit coupling require fully relativistic pseudopotentials and a noncollinear representation of the electronic spin degrees of freedom. This can be computationally intense, but is mandatory if relativistic effects on the electronic ground state are to be taken fully into account.
Silicon band splitting with ATK-DFT¶
You will now use ATK-DFT to calculate the LDA bandstructure of silicon with and without spin-orbit coupling.
Spin-orbit DFT calculations are more expensive than standard LDA or GGA methods and will often require more SCF iterations to converge. Usually, it is therefore more efficient to first do a spin-polarized ground state calculation, and then use this state as an initial guess for the spin-orbit DFT calculation. This procedure reduces the number of spin-orbit SCF steps and significantly reduces the overall computational time.
Not all pseudopotentials can be used for spin-orbit DFT calculations. The OMX potentials supplied with QuantumATK contain the required SO terms, while the FHI potentials do not.
Please remember that OMX pseudopotentials have in general more valence electrons in the cores than the FHI potentials do. Therefore, the OMX potentials usually require good basis sets and higher density-mesh cutoff energies (e.g. 200 Hartree).
- Use the Builder to set up the silicon crystal with the default (experimental) lattice constant.
- Send the bulk configuration to the Script Generator and add two sets of the following three blocks: New Calculator, Initial State, and Bandstructure.
LSDA initial guess for the ground state¶
- Select the OMX pseudopotential with 4 valence electrons, and use the “Medium” accuracy basis set:
- In the Basic calculator settings,
increase the density-mesh cutoff to 150 Hartree,
select “Polarized” spin,
set up a 7x7x7 k-point grid,
lsda.ncas output file:
- Click OK to close the widget.
- Change Initial state type to “User spin”. Then click OK.
This defines a non-zero initial guess for the silicon polarization, even though we know that silicon is not magnetic. Indeed, we will find that it is not, because the self-consistent polarization will be close to zero. However, a self-consist LSDA calculation is usually a good starting point for noncollinear calculation, with or withour spin-orbit.
That fully defines the LSDA ground state calculation followed by a band structure analysis.
Next, you will set up a similar workflow for SOLDA calculations, but use the LSDA ground state as an initial guess for the SOLDA ground state.
Remember to also choose the correct basis set and pseudopotential.
If you use a non-polarized calculation as initial state for the spin-orbit calculation, you should set the initial spin to zero.
The ATK-DFT script is now done. It can also be downloaded from here:
You should be able to monitor the progress of the job while it is running:
Once the results have appeared on the LabFloor, highlight the two bandstructure objects and use the Compare Data plugin to directly visualize the LSDA and SO+LDA band structures in the same plot:
Several observations can be made:
- As expected, the LSDA spin-up and spin-down bands are identical (blue lines).
- At first sight, SO+LDA and LSDA bands appear indistinguishable. In particular, they both yield and indirect band gap of around 0.5 eV.
- However, if you zoom in around the highest valence bands at the \(\Gamma\)-point (see below), the SO coupling has partly lifted the degeneracy of those bands by splitting them:
- The top band has split into two bands, corresponding to “heavy” and “light” holes.
- The third band is now separated from the two others by the so-called “split-off” energy, experimentally measured to be 42.6 meV [dYHS89]. With SO+LDA we find 44 meV.
SO+MGGA band gap¶
The band splitting obtained above with the SO+LDA method is in very good agreement with the experimental splitting (42.6 meV [dYHS89]) However, the band gap is less than half of the experimental value (1.12 eV).
This is largely remedied by appliyng the TB09LDA method, which employs a metaGGA exchange-correlation potential and gives remarkably accurate band gaps for semiconductors [dTB09]. We should therefore expect the SO+MGGA method to accurately predict band splittings as well as band gaps.
- Go back to the Script Generator and delete the first set of three blocks from the script.
- Change the exchange correlation functional to “SOMGGA”
- Make sure that the Initial State block reads in the ground state from
lsda.nc, and that the Bandstructure analysis block saves its output in
- Save the script as
mgga.py. You can also download it if needed:
- Run the job – it may take 5-10 minutes to finish if run in serial. This can be reduced to 3 minutes if the job is run in parallel with 4 MPI processes.
The scf convergence can for small bulk systems be speed-up by increasing the damping factor to 0.5 in the iteration control parameters setup.
Using again the Compare Data plugin, you can now compare the SO+LDA band structure to the SO+MGGA one.
In the figure above it is clear that SO+MGGA has increased the distance from valence to conduction bands. Try to use the Bandstructure Analyzer plugin to measure the SO+MGGA indirect band gap. You should find ~1.3 eV.
The figure below zooms in on the split valence bands. They have not changed much from the SO+LDA results, apart from a vertical shift.
GaAs band structure with ATK-SE and SO coupling¶
The ATK-SE engine offers tight-binding models for electronic structure calculations as a parametrized alternative to DFT. The Slater–Koster (SK) model is one such method, and includes spin-orbit coupling through a speciel SOC parameter, which is implemented in the ATK-SE engine.
Choose the ATK-SE: Slater-Koster calculator engine and the Bassani.sp3d5s basis set, and choose Noncollinear Spin-Orbit for the spin.
In the Bandstructure block, increase the Points per segment to 101 in order to finely sample the band structure between each high-symmetry point, and trim the default Brillouin zone route such that it reads “G, X, W, L, G, K”.
Run the calculation – it will only take a few seconds – and then use the Bandstructure Analyzer to plot the results. The SO induced spin-splitting of the GaAs band structure is clearly visible. If you zoom in a bit you can measure a direct band gap of 1.51 eV at the \(\Gamma\)-point, which is virtually identical to the experimental gap at zero Kelvin.
|[dTB09]||F. Tran and P. Blaha. Accurate band gaps of semiconductors and insulators with a semilocal exchange-correlation potential. Phys. Rev. Lett., 102:226401, 2009. doi:10.1103/PhysRevLett.102.226401.|
|[dYHS89]||(1, 2) Zhiyi Yu, YX Huang, and SC Shen. Spin-orbit splitting of the valence bands in silicon determined by means of high-resolution photoconductive spectroscopy. Physical Review B, 39(9):6287, 1989.|