
512 Handbook of Self Assembled Semiconductor Nanostructures for Novel Devices in Photonics and Electronics
[58, 60] that such a simpler approach is more suitable for the description of fundamental transi-
tions in QDs, as those under consideration here. For a comprehensive review of theoretical work
on QD modelling, we refer the reader to [61] .
We considered QDs with truncated conical shapes and with base diameters and heights given
by the most frequent values in the respective distributions derived by AFM characterization of
uncapped structures; values of the ratio between base and top diameters were taken similar
to those reported in the literature (2.5–3.0) [62, 63] . The composition of QDs was considered
uniform and equal to the nominal one; this assumption is justifi ed since, as discussed in section
16.3.1, the relatively low growth temperatures of QDs and UCLs are expected to limit segregation
and interdiffusion of elements.
The band discontinuities between QDs and CLs were calculated according to the model solid
theory (MST) [64] , inclusive of strain effects. The MST parameters of InGaAs were deduced from
those of the binary constituents as proposed by Cardona and Christensen [65] , while the effective-
mass values (in units of free electron mass) for InAs QDs were taken as 0.04 for electrons, to
include the effect of strain, and 0.59 for heavy holes [57] . The values of other parameters such
as elastic constants and electronic band structure ones were taken from the recent literature
[66] . The binding energy of excitons was assumed to be 20 meV [67] .
The calculation of the hydrostatic and biaxial components of the QD strain, determining
the band-edge profi les of QDs and CLs, were carried over following the analytical approach of
Andreev et al. [47] , where values for InAs QDs embedded in GaAs are given as linear functions
of the QD–CL mismatch. While the hydrostatic strain is constant inside the QDs, the biaxial
one is position dependent, therefore we have taken an average value derived from the results of
[47, 68] ; in particular, we assumed the value of 0.05 for InAs QDs in GaAs confi ning layers and
scaled the value for InAs QDs in InGaAs CLs. Finally, the biaxial component of strain in the CLs
induced by the presence of QDs was disregarded.
The Schrödinger equation for electrons and heavy holes was solved considering both the QD
and a 1.6 ML thick WL as included in a cylindrical box larger than the QD; the envelope func-
tions of carriers were written in terms of a complete orthonormal set of eigenfunctions of the
box [56] . Then, the Schrödinger equation was transformed in a system of linear algebric equa-
tions, which were solved with well-known computational techniques.
The experimental data of Figs 16.2–16.5 proved that model calculations are indeed accurate
in predicting the dependence of emission energy and activation energy of carrier thermal escape
on the two design parameters, the CL composition x and the QD–CL mismatch f . Thanks to this
experimental validation, a deeper analysis of model calculations can provide a more quantitative
insight on the fundamental features of QD strain engineering and may allow for a rationalization
of our and literature experimental data.
350
300
250
200
Activation energy (meV)
150
100
50
0
0.045 0.050 0.055 0.060
QD–LCL mismatch
InAs/In
x
Ga
1⫺x
As QDs
metamorphic LCL
0.065 0.070 0.075
x ⫽ 0.15
x ⫽ 0.28
x ⫽ 0.35
Figure 16.5 Activation energies for PL thermal quenching as a function of the QD–CL mismatch. The dotted lines
are guides for the eye, while continuous, dashed and dotted–dashed lines represent the results of model calculations.
Reprinted with permission from [38]. Copyright 2007, American Institute of Physics.
The upper surface in Fig. 16.6 represents the calculated 10 K emission energy of QDs as a func-
tion of x and f [38] , while the lower one represents the sum of the energy differences between
CH016-I046325.indd 512CH016-I046325.indd 512 6/25/2008 3:16:19 PM6/25/2008 3:16:19 PM