k⋠p calculations for double p-type δ-doped quantum wells in GaAs

The hole subband structure in double p-type δ-doped quantum wells in GaAs is computed with the use of the 4×4  k⋅p Hamiltonian. The Thomas-Fermi-Dirac approach is implemented for the description of the valence band bending, and the hole states at the Brillouin zone center are calculated along its lines, within the effective mass approximation. The zone center eigenstates obtained are then used to diagonalize the k⋅p Hamiltonian for non-zero k. The hole subband structure is analyzed as a function of the impurity density and the distance between δ wells.