Anisotropic Confinement, Electronic Coupling and Strain Induced Effects Detected by Valence-Band Anisotropy in Self-Assembled Quantum Dots.

A method to determine the effects of the geometry and lateral ordering on the electronic properties of an array of one-dimensional self-assembled quantum dots is discussed. A model that takes into account the valence-band anisotropic effective masses and strain effects must be used to describe the behavior of the photoluminescence emission, proposed as a clean tool for the characterization of dot anisotropy and/or inter-dot coupling. Under special growth conditions, such as substrate temperature and Arsenic background, 1D chains of In0.4Ga0.6 As quantum dots were grown by molecular beam epitaxy. Grazing-incidence X-ray diffraction measurements directly evidence the strong strain anisotropy due to the formation of quantum dot chains, probed by polarization-resolved low-temperature photoluminescence. The results are in fair good agreement with the proposed model.

Introduction

Recent attention has been given to the study of coupled quantum dot (QD) arrays for their potential application in quantum information processing [1-3]. The self-assembling process and its control become essential concerns in the search for new proposals of optoelectronic and quantum computing devices. Also, the spinor states in quasi-zero dimensional systems and their electronics have become features of renewed interest [4-7]. High uniformity of size, shape and distribution control of dot arrays are required in many application proposals like detectors, low-threshold lasers and photonic crystals. The lack of control over the self-assembly process of formation of these QDs leads to inhomogeneous broadening in size and/or shape that may degrade the quality of a device application. Therefore, the need for probing size, shape and effective inter-dot coupling has become an important area of research in recent years [8-12].

The anisotropy observed in linearly polarized PL-emissions from self-assembled QDs has been studied in recent years, and several works have detected some correlation with the anisotropic shape of the QD array [13-16]. There is also an agreement about the complexity of valence-band effects in QDs as a relevant issue when dealing with optical response from transitions between these completely localized states [7,17,18].

In the present work, we addressed mechanisms of testing simultaneously one-dimensional (1D) lateral ordering of dots, inter-dot coupling and 2D anisotropy of self-assembled QDs from studies of grazing-incidence X-ray diffraction (GID) and polarized photoluminescence (PL) emissions under different excitation power. This work has been motivated by the plausibility of controlled self-assembling growth of 1D dot arrays (QD chains) [19] and their potential use for testing important quantum effects such as correlation of information and optical coupling between dots where the relevant aspects of effects associated with inter-dot coupling and QD shape, size and distribution deserve special attention. It is also discussed the interplay between shape and strain fields with the inter-dot correlation that is revealed in the GID measurements and PL-emission spectra from QD arrays. Two sets of samples are investigated: one shows chain-like 1D correlation between neighboring dots and the other exhibits a mostly random island distribution. Two different QD shape models are used in order to calculate and test the polarized optical emission spectra dependence with spatial dot correlation and local geometry. The experimental confirmation included in this work highlights and supports the importance of probing correlated distribution in QD arrays for the characterization and improving of the growth-controlled processes.

Theoretical Model

A multi-band k · p model based on the standard Kohn–Luttinger [20] and parabolic Hamiltonians to probe the electronic structure of holes and electrons, respectively, in dots grown along the [100] direction was developed. Due to strong valence-band admixture, such a procedure provides straightforward information on the relaxation of the inter-band optical transition selection rules, using lower computational efforts than in tight-binding calculation model, for example [13,14]. The built-in strain field distribution, which lead to the formation of self-assembled QD arrays, has been considered within the Bir–Pikus deformation potential model [21]. Uniform strain tensors are assumed, a model that neglects effects caused by variations at the QD interfaces [22,23]. This approximation works reasonably well for the study of ground-state properties of medium (~150 Å) and large (>250 Å) size dots.

The double quantum dots structures under investigation are schematically illustrated in Figure 1. According to realistic dimensions the dots are assumed to have semi-cylindrical shape with radius ρ, laterally separated by an inter-dot barrier of thickness d. Since the main focus is concentrated in the tunneling along the lateral direction , the confining potential is defined as V(, z) = V() + V(z), where the infinite barrier model have been used, as represented in Figure 1b (Figure 1c), at top (left) and bottom (right) interfaces, whereas the finite barrier model at the internal interfaces have been adopted, as represented in Figure 1c, in order to account for inter-dot coupling effects.

Figure 1

(a) Schematic modeling of QD size and inter-dot coupling used in this study of self-assembled dots formed along the indicated crystalline directions. (b) Confinement model for random distribution dots in the (100) plane. (c) Confinement model for testing anisotropic size and plausible inter-dot electronic coupling.

For the GaInAs alloys in consideration, at the center of the Brillouin zone, the split-off band is energetically well separated from the topmost valence subbands. In the limit of decoupled split-off band, the four-band Hamiltonian provides a good description of low-lying hole states by considering the coupling between the heavy-hole (hh) (J = 3/2, j = ±3/2) and the light-hole (lh) (J = 3/2, j = ±1/2). In the effective-mass approximation, when spanned in this basis, the kinetic energy of the hole is described by the 4 × 4 Kohn–Luttinger Hamiltonian

where,

with the Luttinger parameters γ (i = 1, 2, 3), and the momentum operators .

The Hamiltonian of the hole in the quantum dot system is

where V() is an infinite barrier outside of the semi-cylindrical cross-section, and V(z) is a double quantum well potential with infinite high outside walls, whose finite barrier is due to the offset between the band edges in the well and barrier materials; ℋ is the Bir–Pikus Hamiltonian [21].

By exploring cylindrical symmetry in the Kohn–Luttinger model, the wave function of a hole state can be written in the form

The indexes (j, n, m) label the quantization along the z-direction (j) and in-plane (n, m) quantum numbers, respectively, α denotes the spin-up (| ↑ >) and spin-down (| ↓ >) periodic Bloch function character, namely: |hh ↑⟩, |lh ↑⟩, |hh ↓⟩ or |lh ↓⟩ and, finally, are the weight coefficients in the basis set of envelope wavefunctions, F(z)f(ρ, φ), at a position (, φ) inside the dot. The solutions for the in-plane motion, f(ρ, φ), are given by [24]

for semi-cylindrical confinement (Figure 1c). In these expressions, μ is the mth zero of the Bessel function of order n, J(x), whereas the form of function F(z) depends on the profile potential along z-direction between the dots. The depth of the quantum well is determined by the offset between the valence-band edges in the dot and the barrier materials. For the GaAs/In0.4Ga0.6As interface, the valence-band offset can be estimated as ΔE = 214 meV. By analytically solving the Schrödinger equation for holes and regarding the mismatch between the Luttinger parameters in the GaAs/In0.4Ga0.6As interfaces, the transcendental equation is derived, which determines all subband energies (j) and the corresponding wavefunctions (see Appendix 1). The signal (±) in the Eq. 16 provides them, respectively, with symmetric or asymmetric character . Taking advantage of this fact, the Hilbert space for the hole wavefunctions Ψ(r) can be split into two orthogonal subspaces, labeled I and II, that are classified according to the parity of the quantum number j. As a result, the Hilbert subspace I(II) gathers spinor states with spin-up (spin-down) components having odd j-values (even j-values) that are coupled with states with spin-down (spin-up) and odd j-values (even j-values). Hence, the eigenvalue problem for the Hamiltonian in Eq. 3 can be solved independently for each class of states I and II. The hole state wavefunction (4) for a given subspace can then be written as

The hole states of the semi-cylindrical QDs system are calculated by exact diagonalization of the Hamiltonian ℋ, on a finite basis set expansion given by Eq. 6 using a standard numerical diagonalization technique. The matrix elements of the momentum operators and involved in the off-diagonal terms Eq. 2 of the Hamiltonian ℋ are given in Appendix 2.

As shown schematically in Figure 1, effects associated with isotropic and anisotropic spatial confinements are simulated in the calculation by changing the lateral sizes, D011 and , in the (100) plane as well as the inter-dot distance (d). Two geometry cases will be studied: (i) Uncorrelated dots, which consider isotropic spatial confinement in the (100) plane, with , without inter-dot coupling. The spin quantization axis (z-axis) is chosen along direction [100] (Figure 1a) and 2D dot distribution is random; (ii) Correlated dots, which consider anisotropic spatial confinement () and include inter-dot coupling (Figure 1b) that leads to a chain-like 1D dot alignment. Here, the spin quantization axis (z-axis) must be set along the direction [25,26].

These two models were tested and compared in order to search for the main qualitative differences between optical emission probabilities for light polarized along and perpendicular to the z-axis, respectively. This modeling tests the different behavior of optical emissions associated fundamentally with the difference between heavy-hole (hh) and light-hole (lh) longitudinal and transversal ellipsoidal effective masses as well as the effects originated from the strain fields on these hole energy levels.

The oscillator strength for optical electric fields linearly polarized along the and parallel to [011] directions (see Figure 1) can be calculated as . For uncorrelated dot arrays, showing mostly random distribution (case (i)), they are given by

where P = ⟨s|p|x⟩ = ⟨s|p|y⟩ = ⟨s|p|z⟩ is the isotropic conduction-valence-band momentum matrix element between functions at the Γ-point, is the overlap between jth electron and hole envelope functions along z-axis, and the factor 2 is due to double spin degeneracy.

All coefficients , shown in Eq. 7, are real when calculated for cylindrical uncorrelated dot array case, and using the expansion set in Eq. 3. This result leads to identical oscillator strengths and, consequently, equal PL intensities for both optical linear polarizations. More specifically,

according to Eq. 4, and this identity is independent of QD size. Besides, neither hydrostatic nor axial strain contributions would induce changes to Eq. 8 in this symmetric case (unless anisotropic strains are applied). Therefore, a distribution of cylindrical uncorrelated dots over the (100) plane would lead to identical linear PL-emission intensities polarized along and perpendicular to the z-axis.

In correlated arrays showing preferential dot diffusion, the compressive strain can be relaxed by forming 1D arrangement, as occurring for strain distribution in free-standing superlattices. In this case, the in-plane strain is defined by ε|| = εxx = εyy = (a⊥ - a)/a, where the lateral lattice constant (a⊥) can be calculated as [27]

Here, Sα = (S11 + S12)α is the sum of elastic compliance constants, Lα (aα) is the width (bulk lattice constant) of the corresponding layers regions α = w (well) or b (barrier). In this way, a 3% strain can be relaxed to a value near 1%. Although shear strain contribution, which affects the separation between hh and lh subbands, becomes relaxed, the hydrostatic strain component leads to the effective reduction of the inter-dot potential barrier, which enhances the inter-dot coupling and tunneling. The envelope function spreading along the direction favors the confinement of a carrier with higher in-plane effective mass, which leads to the exchange of the ground-state character, since .

The effects associated with the anisotropic confinement, within the inter-dot coupled model and simulated by a semi-cylindrical dot shape (see Figure 1c), uses only the subset of the expansion functions in Eq. 5 that complies with null boundary conditions at the flat part of the semi-cylinder. The corresponding linear crossed polarized optical matrix elements, for this correlated dot array model (case (ii)), are given by

Here, the factor 2 occurs due to the summation over subbands j = 1,2 since these states are nearly degenerate for large inter-dot separation, d. It is clear that the identity in Eq. 8 has changed and no longer holds for all values of the inter-dot distance and QD sizes. We will be showing below that mass anisotropy of hole ground-state might be hold responsible for these anisotropic optical emission intensities once the dot confinement strength becomes relaxed in certain directions, whether by dot size anisotropy and/or by inter-dot coupling tuned by the strain fields.

First of all, let us analyze the effect of the spatial confinement in the case of a single dot with the semi-cylindrical shape, namely: the limiting case d → ∞ shown in Figure 1c. As the strength of the spatial confinement is relaxed along the direction by the QD size increase, the topmost valence band becomes occupied by a state with a stronger lh-character and reduced hh-contribution [4,28,29]. This effect is caused by the strong hole mass anisotropy, namely: while . It can be noted, from simple arguments, that hh- or lh-mass character of the valence-band ground-state can be interchanged by weakening the spatial confinement strength in the direction . Under weak confinement regime, the total energy determining the level position is mainly inversely proportional to the effective mass, as

and

where ⟨D[011]⟩ and denote mean confining lengths. Consequently, by tuning the confinement anisotropically, the condition E PL transitions allowed for certain light polarization can probe the anisotropic character of the Bloch functions that, in the multi-band calculations, are determined by the values of the expansion coefficients in Eq. 4. It is noted, from Eq. 10, that a state having small hh-character and, consequently, small values of coefficients and , produces smaller oscillator strength for optical transition polarized along the inter-dot coupling direction .

Figure 2 shows the oscillator strength values calculated for two coupled semi-cylindrical QDs with two values of the transverse diameter, D[011], as a function of the axial length, (see Figure 1c). Here, we have estimated the strain strength to hold with the uncorrelated dot array condition and confirm that the bigger the transverse size of dot array is (Figure 1b) the smaller must the strain order factor be. Furthermore, for compressive strain ε|| > ε⊥, the crossing point can be shifted toward the dotted line. For dilation strain, with ε|| < ε⊥, the crossing point is shifted away from the uncorrelated dot condition, and this condition can be attained in self-assembled QDs grown along the [100] direction. Certainly, shear strain field distribution is able to tune the equal oscillator strength condition for these mutually perpendicular polarized emissions in isolated anisotropic QDs.

Figure 2

Oscillator strength contours .

Analogously to the exchange of ground-state character induced by anisotropic confinement and shear strain fields, this effect can be also produced by electronic coupling between nearest-neighboring QDs, an effect that leads to the enhancement of the effective value . The interchange of ground-state character is highly favored in coupled dots by increasing the inter-dot tunneling, as can be seen in Figure 3, which leads to the envelope function spreading along the coupling direction, . In order to show this effect, we have used the combination of dots with finite inter-dot separation, d. Note, in Figure 3, that coupled dots will show a left-shifted crossing point for equal oscillator strength, when compared to the uncorrelated dot case. As discussed before, this shift can be further modified by shear strain fields.

Figure 3

Oscillator strength contours fulfilling .

For the limiting cases (see Figure 4), and , the oscillator strengths for polarized emissions attain the conditions and , respectively, and these results are attributed to the anisotropy of hole effective masses. The crossing point where the polarized emissions have equal intensities can be shifted by the shear strain contribution to hh- and lh-energy level positions. In Figure 5, it can be observed that the crossing points are shifted to the right as the strain order factor and/or inter-dot distance are increased. Furthermore, two asymptotic limits > 400 Å and < 200 Å where the crossing points coincide were found out, respectively, for various strain strengths and for different inter-dot distances.

Figure 4

Calculated oscillator strengths for crossed linear optical polarizations along the directions . The crossing point stands for isotropic optical emission.

Figure 5

Calculated oscillator strengths for crossed linear optical polarizations for a strained system of two coupled QDs with different inter-dot distances . Here was taken a lateral size D[011] = 350 Å and a strain order factor ε|| = - 0.2%.

Experimental Confirmation of the Purposed Theory

Experiments that confirm this modeling were performed using In0.4Ga0.6As QDs grown by molecular beam epitaxy on semi-insulating (100)GaAs. The QDs were obtained using the Stranski–Krastanov growth mode. Two set of samples were prepared for the experiments: (A) QDs with strong anisotropy in shape along direction and with partial ordering along that; (B) QDs with weak or no anisotropy on the (100) surface and large separation in both in-plane directions. The shape and the distribution of QDs were controlled by the Arsenic background. The use of As2 or As4 background during the growth allows the control of group III element diffusion on GaAs (100) surfaces, providing choices for different dot samples with the same composition but different shapes and distribution. Details of growth mechanisms and the processes involved in diffusion controlling by the background Arsenic environment are described in Ref. [19].

Two sets of samples A were grown under As4 background. In one set, the layer of dots was left uncapped for morphology analysis, and in the other, the QDs were buried with GaAs for low-temperature PL analysis. The other two sets samples B were grown under the same conditions as the sets A, except that under As2 background. Surface morphologies of the two uncapped samples were performed by using atomic force microscopy (AFM), as shown in Figure 6, imaged by Nanoscope IV in the tapping mode and using a high-resolution Silicon tip. The (1 × 1) μm AFM images show the morphologies of the In0.4Ga0.6As uncapped dot samples. The mean dot size and the center-to-center distance along the direction of both sets are displayed in Table 1. The AFM pictures show clearly the effect of different Arsenic background both on dot formation and distribution. The predominantly anisotropic dot shape and distribution obtained along direction is for samples grown under As4 environment. Finally, these sets of samples enable us to use sample (B) as the reference for uncoupled QD arrays with mostly isotropic distribution on the (100) plane.

Figure 6

One layer AFM 1 × 1 μm image of . Sample A (left) shows 1D chain-like ordering along the direction. Sample B (right) shows mostly isotropic or randomized dot distribution in the (001) plane.

Table 1

Average QD parameters with dispersion obtained from a Gaussian fit of the AFM data

Sample	Density (cm-2)	D_[01¯1](Å)	D[011] (Å)	h (Å)	d (Å)
A	3.9 × 1010	280 ± 12	220 ± 10	63 ± 12	160 ± 30
B	1.9 × 1010	350 ± 15	350 ± 15	90 ± 15	330 ± 75

: Dot width along ; D[011]: dot width along [011]; h: QD height; d: inter-dot distance along the "chain-direction"

Grazing-incidence X-ray diffraction (GID) measurements were performed in both samples at the XRD2 beamline of the Brazilian Synchrotron Light Laboratory (LNLS), using a 4 + 2 axis diffractometer. The X-ray photon energy was fixed to 10 keV. Since both samples were capped by a GaAs 50 nm layer, the incident angle was fixed at 0.28°, slightly above the GaAs critical angle, maximizing the signal from the buried quantum dots. The diffracted signal was measured by integrating the exit angle from 0 to 1.2° [30].

Figure 7a and 7b show longitudinal θ - 2θ scans in the vicinity of the in-plane (022) and reflections for samples A and B, respectively. Such scans are sensitive to the strain relaxation inside the In0.4Ga0.6As QDs and GaAs surrounding lattice. For all scans, diffuse intensity is observed surrounding the narrow and intense GaAs Bragg peak, located at |H| = |K| = 2. For sample A, the longitudinal scan performed at the vicinity of the GaAs (022) reflection exhibits a much broader profile than the scan measured with the sample rotated by 90°, close to the reflection. Such a behavior indicates that a more effective strain relaxation for the islands may take place along the [022] direction, while a more strained lattice profile is found along the direction. The intensity distribution in both profiles of Figure 7a is almost symmetric with respect to the GaAs peak position, denoting the existence of compressively strained InGaAs inside the QDs, as well as on the GaAs matrix surrounding the QDs [31]. Similar diffraction profiles are observed in the longitudinal scans performed on sample B (Figure 7b). For this sample, the difference of widths of diffuse intensity on (022) and scans is not as pronounced as observed for sample A, indicating a less anisotropic relaxation.

Figure 7

Radial scans at the vicinity of the . Lateral size from iso-strain regions in samples A (c) and B (d) obtained from the width of transversal scans.

In order to quantify the strain relaxation inside QDs in both samples, transversal scans were performed at several positions along the longitudinal profiles shown in Figure 1a and 1b. These scans (not shown here) are measured by fixing the θ - 2θ condition and varying the sample rotation angle θ solely. In momentum transfer space, the angular momentum transfer q = (4π/λ)sin(2θ/2)sing(Δθ) is varied, where Δθ = (θ/2θ)/2. Such a procedure allows to obtain the average lateral size L of regions inside the QDs with constant strain status by evaluating the width Δq of transversal scans, L = 2π/Δq[30,32]. Values obtained for the local lateral size of iso-strain regions as a function of the in-plane strain status for samples A and B in the [022] and directions are shown in Figure 7c and 7d, respectively. For both samples, the lateral size of iso-strain regions along the QDs chain direction is larger than along the axis parallel to the chains. The ratio , which is a quantitative indicator of the anisotropic lattice relaxation inside QDs, is larger for sample A than for sample B, corroborating the qualitative information inferred from the widths of longitudinal scans.

Some considerations must be drawn before extending the analysis of the data shown in Figure 7c and 7d. For uncapped QDs, the relaxation of lattice parameter is monotonic from their base to their apex [32]. Capped QDs, in contrast, exhibit a non-monotonic gradient, with lateral and vertical strain variations. This condition generally implies in the existence of similar in-plane strain status on the island base and apex, both in contact with the GaAs matrix. It is therefore impossible to resolve vertically the position of iso-strain areas for the capped QDs with our GID measurements. Nevertheless, the lateral sizes observed represent a good approximation of the in-plane area of iso-strain regions projected on the substrate surface plane. Such approach allows for a visualization of the anisotropic strain relaxation. Since the diffraction signal observed at (|H|, |K|) > 2 is related to the existence of compressively strained GaAs surrounding the QDs [31], maps with the projected view of iso-strain areas were extracted from the experimental data by taking into account the L values for (|H|, |K|) < 2, which are directly related to compressively strained In0.4Ga0.6As from the islands (Tensile strained GaAs at the bottom and at the apex of the island also contribute to the diffracted intensity at |H| = |K| < 2. However, the total volume of material with local lattice parameter larger than a outside the island is considerably smaller than the amount of material contained inside the islands. For a discussion on tensile strained substrate material see [33]). These projection maps for QDs from samples A and B are shown in Figure 8a and 8b, respectively. The iso-strain projection areas were drawn following the condition that they are contained on curves delimited by

Figure 8

In-plane projection of iso-strain regions for a field of view with several islands for samples A (a) and B (b). The in-plane strain represented in the color scale is relative to the GaAs bulk lattice.

where x and y are the in-plane coordinates along the [011] and directions, respectively, considering the plane origin at the central QD position. Figure 8 shows the iso-strain areas for an in-plane region of approximately 1,100 Å × 700 Å, which contains 9 QDs for sample A and 4 QDs for sample B (see Table 1). The color scale in these maps refers to the in-plane strain with respect to bulk GaAs.

From Figure 8a, one clearly observes that iso-strain contour lines from one QD of sample A almost reach the neighbor QDs along the chain direction. An asymmetric ratio of 1.7 is found for the broader iso-strain contour lines of QDs in this sample, pointing out again to a more pronounced strain relaxation along the [011] axis. The physical presence of very close QDs along the chains may therefore induce a modulation of the strain field that allows for a gentle strain relaxation in the direction. In sample B (Figure 8b), the asymmetric shape of iso-strain regions is still observed, but with a ratio of 1.35. Although an elongation is observed along the direction, the QDs are too apart from each other and do not strongly influence the strain field of the neighbor QDs in this direction.

Since the GID measurements do not reveal directly the height above the substrate of each iso-strain region finite element method, simulations were performed using a commercial software package to provide complementary information on the strain configuration of capped islands. In our simulations, a three-dimensional box containing a single GaAs capped In0.4Ga0.6As island was created for each sample, with periodic contour conditions at all lateral edges in order to take into account the symmetry of QD chains and the possible interaction with the strain field from neighbor QDs. A 15 Å thick wetting layer of nominal concentration was inserted between the islands and the substrate, following Ref. [34]. The island profiles used in this simulation were extracted from the AFM measurements in uncapped islands (Figure 6) that resulted in the dimensions from Table 1. The nominal composition was kept, assuming thus a negligible deviation of island stoichiometry from the nominal values (Anomalous grazing-incidence diffraction measurements performed at the Ga - K edge do not point out to deviations (within an error bar of 7%) from the nominal In/Ga content inside QDs.). Such assumptions consider that islands do not undergo dramatic changes in morphology or composition under capping, which is a valid approximation for the growth temperatures used here and the reduced strain with respect to pure InAs islands [35]. Finally, a 500 Å-thick cap layer was added to the simulation, as represented in Figure 9a.

Figure 9

(a) Representation of the two-dimensional cuts shown in maps panels (b–g) performed on the finite element method simulations with periodic contour conditions at the substrate box edges. The color contours represent variations on the first axial principal strain, which allows a qualitative comparison with the GID data of Figure 8. Cuts on the bottom (b), middle (d) and top (f) of the average island of sample A show an elongated strain profile along the directions. Similar cuts for the average island of sample B are seen on (c), (e) and (g).

Two-dimensional cuts of the simulated data are shown in Figure 9b, d, and 9f for sample A and Figure 9c, e, and 9g for sample B. The selected cuts are schematically depicted at Figure 9a and were chosen to be at the island bottom (b) and (c), middle (d) and (e), and top (f) and (g). Since the representation used in Figure 8 cannot be directly correlated to the Cartesian in-plane strain components x and y, the maps of Figure 9b–g were drawn as a function of the axial (first) principal strain component. Such principal component analysis allows for the reduction of the dimensionality of the data set, providing a resulting representation with radial symmetry. The axial strain component is given by [36]

where ε is the in-plane shear strain and ε, ε the normal in-plane strains. For all principal strain maps, the color scale represents the deviation of the local lattice parameter with respect to the bulk local lattice parameter. Therefore, higher principal strain values are found in positions where the In0.4Ga0.6As lattice of the islands is fully strained to the GaAs lattice constant. Finite values of the axial strain component are also observed in regions surrounding the islands, in which the GaAs local lattice is affected by the proximity to the island. Selected contour level edges were marked by dark lines in all maps as a guide to the eyes.

The maps generated by in-plane cuts in the simulation of QDs in sample A clearly exhibit elongated contours along the direction, most notably for the cuts at the island basis and middle. This indicates that for lower in-plane strain conditions, the lattice surrounding the islands behave as semi-continuous wires along the direction. In the QDs of sample B, an elongation of axial strain contour levels is also observed along the chain directions for all maps. However, the anisotropic effect is much more reduced with respect to the results obtained for sample A.

The effects of different dot size distributions and inter-dot coupling have been analyzed by low-temperature linear polarized PL measurements carried out on the samples A and B buried with GaAs. The samples were placed in a closed-cycle Helium cryostat (Janis—CCS-150) and excited using a 532-nm continuous wave YAG laser (Coherent Verdi V10—10 W). The PL signal was carried out by a monochromator (SpectraPro 2500i—0.5 m focal length) and detected by a liquid-nitrogen-cooled InGaAs photodiode detector array (Princeton Instruments—model 7498-0001). Figure 10a and 10b show the PL intensity at 10 K for samples A and B, respectively, where the emission spectra, for each sample type, collected with two linear polarizations, namely: along and along [011].

Figure 10

PL spectra for crossed linear polarizations, taken at . The degree of linear polarization: has been included in these panels. (c) PL peak position as a function of the excitation intensity.

In Figure 10a, one may see a polarization degree around 6%, as might be expected due to the elongation in the quantum dots profile revealed by the AFM images (Figure 6) and strain distribution (Figure 8). As highlighted in Figures 2 and 3, the oscillator strength grows for emissions linearly polarized along the larger dot size direction. This behavior is enhanced for inter-dot separation up to d ~ 160 Å. When d is further reduced, the inter-dot tunneling probability increases considerably, and this behavior is also enhanced. The PL intensity polarized along coupling direction is also enhanced in coupled QDs by the reduction of barrier heights due to hydrostatic strain of the order of 1%. Besides, the anisotropic PL-emissions from sample A, as shown in Figure 5a, can be qualitatively reproduced by the oscillator strengths, shown in Figure 3, calculated by using the nominal values for both samples. As seen in Figure 3, an effective increase in inter-dot tunneling (distance d ~ 160) Å would lead to the relaxation of the confinement along the direction. These effects would lead to a hole ground-state character exchange from predominant-hh to -lh, and to the intensity difference between these cross-polarized emissions, experimentally confirmed by Figure 10a.

For the isotropic case, PL-emissions occur when , a condition well-fulfilled for the cylindrical model of Figure 1b. By changing the dot shape and coupling along direction in (100) plane, the model shows that condition can be obtained for semi-cylindrical geometry only for a small combination of values that emulates uncoupled dot distribution in the (100) plane if strain effects are included into the Hamiltonian. According to the theoretical modeling, an isotropic dot distribution on the (100) plane (case (i)) accounts for isotropic crossed polarized PL-emissions, as shown in Figure 5b for sample B. However, according to Figure 10b, a small polarization degree is still present in symmetric QDs, associated with the elongation that remains, as revealed by the Figure 8b. Such feature might come from the anisotropic diffusion rate of Indium atoms during the growth, which presents a higher mobility than the Gallium atoms. Furthermore, the Indium diffusion coefficient is faster along the than along the [110] direction, and as a result, the quantum dots of sample B are not completely symmetric [37,38].

To confirm the results from X-ray measurements, Figure 10c displays the shift in peak position of the spectra as a function of the excitation intensity. Note, for sample A, a shift toward higher energies as the excitation intensity grows. Such a blue-shift for the elongated dots has been associated with the screening of the built-in electric field due to the presence of strain. On the other hand, for sample B no remarkable energy shift is observed showing that the strain is not so pronounced as in the previous case [39].

Conclusions

The control and simulation of size anisotropy and effective inter-dot tunneling effects, as described in this work, is an important issue to be addressed during the characterization of ordered sets of coupled dots. The strain fields, present during the growth process of these QDs have led to the appearance of anisotropic geometric shapes, mostly elongated along the preferential direction. These effects can be probed by polarized optical responses from different samples. In summary, we have shown that the shape, spatial distribution and the inter-dot coupling of InGaAs self-assembled QDs can be probed and characterized by using linearly polarized PL-emissions. Valence-band effects due to admixing between hole states and strong anisotropic effective masses have led to different PL intensities in samples with lateral QD ordering forming "chain-like" structures. The envelope function model used here to describe the polarized optical responses showed fairly good agreement with structural AFM and X-ray data and may be used to predict or characterize the strength of inter-dot coupling and/or anisotropic dot shape and distribution.

Appendix 1: Double Quantum Well Potential

After matching the wavefunctions fulfilling the hole-Shchrödinger equation at the interfaces using ∂ ln F(z| = β∂ ln F(z)| and ∂ ln F(z| = β-1∂ ln F(z)|, we are able to obtain the transcendental equation

where and and is the rate between the hole effective masses in the barrier and the well region, respectively, and , as well as l is the well width, d the inter-dot distance and V the barrier height. By carrying out numerical calculation, solutions of the Eq. 16 yields to the energy levels with the corresponding wavefunctions of the symmetric (+) and antisymmetric (-) hole states,

where

Appendix 2: Matrix Elements

The matrix elements of the momentum operators are necessary in order to build the Hamiltonian matrix form of the ℋ. In polar coordinates the operators are written as,

Projecting on the wavefunctions (5) with t = (n, m) and t' = (n', m'), it is straightforward to show that

with

and

where are numbers ruled by

In the particular case where t = t' Eq. 20 can be reduced to

Also, it follows the relation

The other matrix elements for the high-order operators are evaluated numerically using Eqs. 18–21 and the matrix identity .

It is worth to show that the element matrix of the diagonal terms in D accomplished

Taking into account the loss of translational invariance along the z direction by replacing the wave vector-component k by the operator -i∂ it is therefore convenient to write the resulting expression for the element matrix in a symmetrized form

where index j stands for the piecewise wavefunctions (17). The resulting integrals in z-direction are solved numerically.