Excitons

The exciton states are computed using the approach from [ChuangOpto1995]. The explanation below only covers the most important aspects of the model, for detailed derivation please refer to the [ChuangOpto1995] book. The exciton computation is only available for 1D systems.

Model

Assuming the effective mass approximation for the electron and hole, the Schrödinger equation for the exciton can be written as:

(2.2.41)(He(r¯e)+Hh(r¯h)e24πε|r¯er¯h|)Φ(re,rh)=EΦ(re,rh)

where He and He are the Hamiltonians for the electron and hole, respectively, r¯e and r¯h are the coordinates of the electron and hole, respectively, ε is the dielectric constant, and E is the energy of the exciton.

The wavefunciton of the exciton, formed by electron n and hole m, will be calculated in the form:

(2.2.42)Φ(re,rh)=exp(iKt¯Rt¯)F(ρ,ze,zh)=exp(iKt¯Rt¯)ϕnm(ρ)fn(ze)gm(zh)

where Kt¯ is the in-plane wavevector of the exciton, Rt¯ is the in-plane coordinate of the exciton, F(ρ,ze,zh) is the exciton envelope function, fn(ze) and gm(zh) are the single particle envelope wavefuncitons of electron and hole in the growth direction. Then, the equation for the unknown ϕnm(ρ) is given by:

(2.2.43)(22mrρ2Vnm(ρ))ϕnm(ρ)=Ebindingϕnm(ρ)

where mr is the reduced mass of the exciton, Ebinding is the binding energy of the exciton, and Vnm(ρ) is expressed as:

(2.2.44)Vnm=dze|fn(ze)|2dzh|gm(zh)|2e24πεs(ρ2+|zezh|2)

The solution of the equation for ϕnm(ρ) can be found variationally by minimizing the binding energy of the exciton. The form of the solution is assumed to be similar to 1S state of 2D hydrogen atom:

(2.2.45)ϕ(ρ)=2π1λexp(ρ/λ)

where λ is the variational parameter, which has an interpreatation of exciton inplane Bohr radius. The variational parameter λ is determined by minimizing the binding energy of the exciton from equation (2.2.43).

Averaging model parameters

The model depends on dielectric constant ε, effective masses of the electron and hole me and mh, which are not constant in heterostructures. If not given in the input file, the volume averaged values of these parameters are used. For effective masses, density weighted average is also possible.

Excitons in multiband Hamiltonians

The computation of the exciton in the case of 8-band kp Hamiltonian is complicated by the fact, that the electron and hole Hamltonians are no longer separabable. In that case the equations derived from effective mass Hamiltonians are used, using wavefunctions computed with the 8-band Hamiltonian. As the effective masses are not longer parameters of the Hamiltonian, the effective masses used are computed from the paramters for the 8-band Hamiltonian: L,M,N,EP,S,Egap. The same approach is used for 6-band kp Hamiltonian, where the effective masses are computed from the parameters L,M,N.