mqwindex - Script command

result = mqwindex(stack_properties, simulation_parameters, config);

stack_properties: struct with fields that define MQW stack geometry and material properties.

simulation_parameters: struct with fields that define simulation parameters for which the output will be calculated.

config: struct with fields that configure the behaviour of the simulation.

result: struct or a cell of structs in case of multiple partitions, where each struct contains 5 datasets: spatial band diagram, band structure in (E,k) space, spatial wave functions for each (E,k) state, emission coefficients, and complex index.

result = mqwindex(stack_properties, simulation_parameters);

same as above, but using all default values for the fields in config struct.

banddiagram

dataset

Conduction and valence band edge including strain, but not including quantum confinement effects.

bandstructure

dataset

(E,kt) band diagram for conduction and valence bands. Attributes are: conduction_band, valence_band_lo, valence_band_up. See below for the meaning of lo and up. Parameters are kt and subband.

wavefunction

dataset

Spatial wavefunction for each (E,kt) point. Attributes are conduction_band_1, valence_band_lo_1, valence_band_lo_2, valence_band_up_1, valence_band_up_2, where the 4x4 k.p basis in the valence band is split into two (lo for lower and up for upper) 2x2 bases (vectors in each 2x2 basis designated with 1 and 2). For more information look at references [1] and [2]. Parameters are coordinate, kt, and subband.

emission

dataset

Gain and spontaneous emission coefficients. Attributes are: spontaneous_TE, spontaneous_TM, stimulated_TE, stimulated_TM, where TE and TM stand for electromagnetic modes. Parameters are: frequency/energy/wavelength and ndensity (charge density).

Emission coefficients are calculated for the total stack thickness, including barriers. If only quantum well thickness is of interest, excluding barriers, these coefficients should be scaled by multiplying with (total length)/(total qw length). It is important to ensure that emission coefficients apply only to the thickness used for the calculation of the mode overlap with the gain region. When using partitioning, there will be overlapping barriers between different partitions, e.g. simulation_parameters.stackpartitions = [1,3;3,5], where 1, 3, and 5 are barriers. In that case emission coefficients for each partition again apply to the total thickness of that partition, meaning there may be some double counting with respect to the mode overlap region thickness. To avoid this, emission coefficients can be scaled to apply to quantum wells only, or to apply to a portion of partition that does not overlap with adjacent partitions.

When the exciton model is turned on the attributes become: absorption_TE, absorption_TM. These represent the absorption coefficients (negative gain) in the units of [1/m]. The spontaneous emission is not calculated due to the assumption of the depleted carrier density in the quantum wells.

index

dataset

Complex index \(\bar{n}=\Delta n+i\kappa\). Attributes are: index_TE, index_TM, where TE and TM stand for electromagnetic modes. Parameters are: frequency/energy/wavelength and ndensity (charge density).

Same comments regarding index scaling for correct thicknes as for emission dataset.

ex

dataset

Exciton energies Ex. Exciton energies are a function of exciton orbital quantum number (orbital) and angular momentum quantum number (angularMomentum).

phix

dataset

Exciton wavefunctions PhiX in the momentum (in-plane wavevector) space. The wave function coefficients are parametrized in terms of conduction band subband index (cSubband), valence band subband index (vSubband), transverse wave vector (kt), angular momentum quantum number (angularMomentum), and orbital quantum number (orbital).