This article is part of the series Online First articles in Volume 6 (2011). These articles were published as Online First on SpringerLink in 2010. They should be cited with a 2011 publication year.

Open Access Nano Express

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

L Villegas-Lelovsky12*, MD Teodoro13, V Lopez-Richard1, C Calseverino14, A Malachias4, E Marega35, BL Liang37, Yu I Mazur3, GE Marques1, C Trallero-Giner6 and GJ Salamo3

Author Affiliations

1 Departamento de Física, Universidade Federal de São Carlos, São Carlos, SP 13565-905, Brazil

2 Instituto de Física, Universidade Federal de Uberlândia, Uberlândia, Minas Gerais 38400-902, Brazil

3 Arkansas Institute for Nanoscale Materials Science and Engineering, University of Arkansas, Fayetteville, AR 72701, USA

4 Laboratório Nacional de Luz Síncrotron, Campinas, Brazil

5 Instituto de Física de São Carlos, Universidade de São Paulo, São Carlos, SP 13560-970, Brazil

6 Faculty of Physics, Havana University, 10400, Havana, Cuba

7 Department of Electrical Engineer, University of California, Los Angeles, CA 90095, USA

For all author emails, please log on.

Nanoscale Res Lett 2011, 6:56  doi:10.1007/s11671-010-9786-8


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/6/1/56


Received:6 July 2010
Accepted:9 September 2010
Published:1 October 2010

© 2010 Villegas-Lelovsky et al.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

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.

Keywords:
Molecular beam epitaxy; Self-assembled quantum dots; Inter-dot coupling; Anisotropic effects; Linear polarized photoluminescence emission; Grazing-incidence X-ray diffraction synchrotron; Optoelectronic

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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a>, 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.

thumbnailFigure 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, jz = ±3/2) and the light-hole (lh) (J = 3/2, jz = ±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

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M2">View MathML</a>

(1)

where,

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M3">View MathML</a>

(2)

with the Luttinger parameters γi (i = 1, 2, 3), and the momentum operators <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M4">View MathML</a>.

The Hamiltonian of the hole in the quantum dot system is

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M5">View MathML</a>

(3)

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; ℋBP 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

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M6">View MathML</a>

(4)

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, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M7">View MathML</a> are the weight coefficients in the basis set of envelope wavefunctions, Fj(z)fn,m(ρ, φ), at a position (ρ, φ) inside the dot. The solutions for the in-plane motion, fn,m(ρ, φ), are given by [24]

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M8">View MathML</a>

(5)

for semi-cylindrical confinement (Figure 1c). In these expressions, μn,m is the mth zero of the Bessel function of order n, Jn(x), whereas the form of function Fj(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 ΔEv = 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M9">View MathML</a>. Taking advantage of this fact, the Hilbert space for the hole wavefunctions Ψv(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

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M10">View MathML</a>

(6)

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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M11">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M12">View MathML</a> involved in the off-diagonal terms Eq. 2 of the Hamiltonian ℋKL 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M13">View MathML</a>, 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M14">View MathML</a>, 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 (<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M15">View MathML</a>) 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> direction [25,26].

thumbnailFigure 2. Oscillator strength contours <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M16">View MathML</a> for correlated dots with inter-dot distance d = 160 Å and strain order factor ε|| = –0.3%.

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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> and parallel to [011] directions (see Figure 1) can be calculated as <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M17">View MathML</a>. For uncorrelated dot arrays, showing mostly random distribution (case (i)), they are given by

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M18">View MathML</a>

(7)

where P = ⟨s|px|x⟩ = ⟨s|py|y⟩ = ⟨s|pz|z⟩ is the isotropic conduction-valence-band momentum matrix element between functions at the Γ-point, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M19">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M7">View MathML</a>, 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,

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M20">View MathML</a>

(8)

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 - aw)/aw, where the lateral lattice constant (a) can be calculated as [27]

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M21">View MathML</a>

(9)

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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> favors the confinement of a carrier with higher in-plane effective mass, which leads to the exchange of the ground-state character, since <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M22">View MathML</a>.

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

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M23">View MathML</a>

(10)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M24">View MathML</a>

(11)

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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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: <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M25">View MathML</a> while <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M26">View MathML</a>. 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a>. Under weak confinement regime, the total energy determining the level position is mainly inversely proportional to the effective mass, as

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M27">View MathML</a>

(12)

and

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M28">View MathML</a>

(13)

where ⟨D[011]⟩ and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M29">View MathML</a> denote mean confining lengths. Consequently, by tuning the confinement anisotropically, the condition Elh < Ehh can be attained due to the mass anisotropy of carriers. As a result, the corresponding envelope functions must be more extended in one direction than the other. Thus, the corresponding 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M30">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M31">View MathML</a>, produces smaller oscillator strength for optical transition polarized along the inter-dot coupling direction <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a>.

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, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M13">View MathML</a> (see Figure 1c). Here, we have estimated the strain strength to hold with the uncorrelated dot array condition <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M32">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M33">View MathML</a> 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.

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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M29">View MathML</a>. 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, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a>. 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.

thumbnailFigure 3. Oscillator strength contours fulfilling <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M34">View MathML</a> for correlated dots with inter-dot distances d = 160 Å (solid line) and 330 Å (dashed line) and strain order factors ε|| = –0.3% (red), ε|| = –0.4% (green) and ε|| = –0.5% (blue).

For the limiting cases (see Figure 4), <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M35">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M36">View MathML</a>, the oscillator strengths for polarized emissions attain the conditions <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M37">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M38">View MathML</a>, 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M39">View MathML</a> > 400 Å and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M39">View MathML</a> < 200 Å where the crossing points coincide were found out, respectively, for various strain strengths and for different inter-dot distances.

thumbnailFigure 4. Calculated oscillator strengths for crossed linear optical polarizations along the directions <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> (red line) and [011] (blue line) for two coupled QD's (Figure 1c) with semicylindrical shape and axis in the <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> direction with (a) <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M13">View MathML</a> = 280 Å, d = 160 Å upon strain order of ε|| = - 0.3% (solid line) and -0.9% (dashed line). (b) <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M13">View MathML</a> = 350 Å, d = 330 Å and ε|| = - 0.1% (solid line), -0.2% (dashed line) and -0.3% (dashed-dotted line). The crossing point stands for isotropic optical emission.

thumbnailFigure 5. Calculated oscillator strengths for crossed linear optical polarizations for a strained system of two coupled QDs with different inter-dot distances d = 160 Å (dashed line), 330 Å (dashed-dotted line) and infinite (solid line) corresponding to a single (isolated) QD. 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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.

thumbnailFigure 6. One layer AFM 1 × 1 μm image of In0.4Ga0.6As QDs in samples grown under different conditions. Sample A (left) shows 1D chain-like ordering along the <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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

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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M40">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M41">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M42">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M40">View MathML</a> scans is not as pronounced as observed for sample A, indicating a less anisotropic relaxation.

thumbnailFigure 7. Radial scans at the vicinity of the GaAs (022) (solid dots) and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M40">View MathML</a> (open dots) reflections for sample A (a) and B (b). 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 qa = (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 Δqa of transversal scans, L = 2π/Δqa [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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M42">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M43">View MathML</a>, 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 aGaAs 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

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M44">View MathML</a>

(14)

where x and y are the in-plane coordinates along the [011] and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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.

thumbnailFigure 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.

From Figure 8a, one clearly observes that iso-strain contour lines from one QD of sample A almost reach the neighbor QDs along the <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> chain direction. An asymmetric ratio <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M43">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> direction. In sample B (Figure 8b), the asymmetric shape of iso-strain regions is still observed, but with a ratio <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M43">View MathML</a> of 1.35. Although an elongation is observed along the <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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.

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]

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M45">View MathML</a>

(15)

where εxy is the in-plane shear strain and εxx, εyy 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.

thumbnailFigure 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> directions. Similar cuts for the average island of sample B are seen on (c), (e) and (g).

The maps generated by in-plane cuts in the simulation of QDs in sample A clearly exhibit elongated contours along the <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> and along [011].

thumbnailFigure 10. PL spectra for crossed linear polarizations, taken at T = 10 K with excitation wave length λ = 532 nm along <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> and [011] directions for samples (a) A and (b) B. The degree of linear polarization: <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M46">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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 onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M47">View MathML</a>, 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M16">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M1">View MathML</a> 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 ∂z ln Fj(z|z=l = β∂z ln Fj(z)|z=l and ∂z ln Fj(z|z=l+d = β-1z ln Fj(z)|z=l+d, we are able to obtain the transcendental equation

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M48">View MathML</a>

(16)

where <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M49">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M50">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M51">View MathML</a> is the rate between the hole effective masses in the barrier and the well region, respectively, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M52">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M53">View MathML</a>, 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 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M54">View MathML</a> with the corresponding wavefunctions of the symmetric (+) and antisymmetric (-) hole states,

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M55">View MathML</a>

(17)

where

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M56">View MathML</a>

Appendix 2: Matrix Elements

The matrix elements of the momentum operators are necessary in order to build the Hamiltonian matrix form of the ℋKL. In polar coordinates the operators <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M57">View MathML</a> are written as,

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M58">View MathML</a>

Projecting on the wavefunctions (5) <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M59">View MathML</a> with t = (n, m) and t' = (n', m'), it is straightforward to show that

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M60">View MathML</a>

(18)

with

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M61">View MathML</a>

and

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M62">View MathML</a>

(19)

where <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M63">View MathML</a> are numbers ruled by

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M64">View MathML</a>

(20)

In the particular case where t = t' Eq. 20 can be reduced to <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M65">View MathML</a>

Also, it follows the relation

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M66">View MathML</a>

(21)

The other matrix elements for the high-order operators <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M67">View MathML</a> are evaluated numerically using Eqs. 18–21 and the matrix identity <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M68">View MathML</a>.

It is worth to show that the element matrix of the diagonal terms in Dhh(lh) accomplished

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M69">View MathML</a>

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

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/56/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/56/mathml/M70">View MathML</a>

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

Acknowledgements

Authors acknowledge financial support from the agencies FAPESP and CNPq (GEM, VL-R), CONACYT/Mexico and FAPEMIG (LV-L) and LNLS-MCT (AM) and ICTP/Trieste (CT-G) and the National Science Foundation of the U.S. trough Grant DMR-0520550 (BLL, YuIM). LV-L thanks E. Gomez for technical assistance.

References

  1. Ohshima T:

    Phys Rev A. 2000, 62:062316. Publisher Full Text OpenURL

  2. Li SS, Xia JB, Liu JL, Yang FH, Niu ZC, Feng SL, Zheng HZ:

    J Appl Phys. 2001, 90:6151. Publisher Full Text OpenURL

  3. Li SS, Long GL, Bai FS, Feng SL, Zheng HZ:

    Proc Natl Acad Sci USA. 2001, 98(21):11847. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Prado SJ, Trallero-Giner C, Alcalde AM, Lopez-Richard V, Marques GE:

    Phys Rev B. 2004, 69:201310. Publisher Full Text OpenURL

  5. Lopez-Richard V, Alcalde AM, Prado SJ, Marques GE, Trallero-Giner C:

    Appl Phys Lett. 2005, 87:231101. Publisher Full Text OpenURL

  6. Lopez-Richard V, Prado SJ, Marques GE, Trallero-Giner C, Alcalde AM:

    Appl Phys Lett. 2006, 88:052101. Publisher Full Text OpenURL

  7. Mlinar V, Tadić M, Peeters FM:

    Phys Rev B. 2006, 73:235336. Publisher Full Text OpenURL

  8. Wang ZM, Holmes K, Mazur YI, Salamo GJ:

    Appl Phys Lett. 2004, 84:1931. Publisher Full Text OpenURL

  9. Karlsson KF, Troncale V, Oberli DY, Malko A, Pelucchi E, Rudra A, Kapon E:

    Appl Phys Lett. 2006, 89:251113. Publisher Full Text OpenURL

  10. Troncale V, Karlsson KF, Oberli DY, Byszewski M, Malko A, Pelucchi E, Rudra A, Kapon E:

    J Appl Phys. 2007, 101:081703. Publisher Full Text OpenURL

  11. Švrček V:

    Nano Micro Lett. 2009, 1:40. OpenURL

  12. Botsoa J, Lysenko V, Géloën A, Marty O, Bluet JM, Guillot G:

    Appl Phys Lett. 2008, 92:173902. Publisher Full Text OpenURL

  13. Sheng W, Xu SJ:

    Phys Rev B. 2008, 77:113305. Publisher Full Text OpenURL

  14. Sheng W:

    Appl Phys Lett. 2006, 89:173129. Publisher Full Text OpenURL

  15. Favero I, Cassabois G, Jankovic A, Ferreira R, Darson D, Voisin C, Delalande C, Roussignol P, Badolato A, Petroff PM, Gerard JM:

    Appl Phys Lett. 2005, 86:041904. Publisher Full Text OpenURL

  16. Cortez S, Krebs O, Voisin P, Gerard JM:

    Phys Rev B. 2001, 63:233306. Publisher Full Text OpenURL

  17. Mlinar V, Tadić M, Partoens B, Peeters FM:

    Phys Rev B. 2005, 71:235336. Publisher Full Text OpenURL

  18. Margapoti E, Worschech L, Mahapatra S, Brunner K, Forchel A, Alves FM, Lopez-Richard V, Marques GE, Bougerol C:

    Phys Rev B. 2008, 77:073308. Publisher Full Text OpenURL

  19. Marega E Jr, Waar ZA, Hussein M, Salamo GJ:

    Mater Res Soc Symp Proc. 2007., 959

    0959-M17-16

    OpenURL

  20. Luttinger JM, Kohn W:

    Phys Rev. 1955, 97:869. Publisher Full Text OpenURL

  21. Fishman G:

    Phys Rev B. 1995, 52:11 132. Publisher Full Text OpenURL

  22. Tadić M, Peeters FM, Janssens KL:

    Phys Rev B. 2002, 65:165333. Publisher Full Text OpenURL

  23. Cesar DF, Teodoro MD, Tsuzuki H, Lopez-Richard V, Marques GE, Rino JP, Lourenço SA, Marega E Jr, Dias IFL, Duarte JL, González-Borrero PP, Salamo GJ:

    Phys Rev B. 2010, 81:233301. Publisher Full Text OpenURL

  24. Trallero-Herrero C, Trallero-Giner C, Ulloa S, Perez-Alvarez R:

    Phys Rev E. 2001, 64:056237. Publisher Full Text OpenURL

  25. Xia JB:

    Phys Rev B. 1991, 43:9856. Publisher Full Text OpenURL

  26. Fishman G:

    Phys Rev B. 1995, 52:11132. Publisher Full Text OpenURL

  27. Mathieu H, Allegre J, Chatt A, Lefebvre P, Faurie JP:

    Phys Rev B. 1988, 38:7740. Publisher Full Text OpenURL

  28. Saito H, Nishi K, Sugou S, Sugimoto Y:

    Appl Phys Lett. 1997, 71:590. Publisher Full Text OpenURL

  29. Noda S, Abe T, Tamura M:

    Phys Rev B. 1998, 58:7181. Publisher Full Text OpenURL

  30. Malachias A, Magalhães-Paniago R, Neves BRA, Rodrigues WN, Moreira MVB, Pfannes H-D, de Oliveira AG, Kycia S, Metzger TH:

    Appl Phys Lett. 2001, 79:4342. Publisher Full Text OpenURL

  31. Roch T, Holý V, Hesse A, Stangl J, Fromherz T, Bauer G, Metzger TH, Ferrer S:

    Phys Rev B. 2002, 65:245324. Publisher Full Text OpenURL

  32. Kegel I, Metzger TH, Fratz P, Peisl J, Lorke A, Garcia JM, Petroff PM:

    Europhys Lett. 1999, 45:222. Publisher Full Text OpenURL

  33. Magalhães-Paniago R, et al.:

    Phys Rev B. 2002, 66:245312. OpenURL

  34. Morkoc D, Sverdlov B, Gao G-B:

    Proc IEEE. 1993, 81:493. Publisher Full Text OpenURL

  35. Mashita M, Hiyama Y, Arai K, Koo B-H, Yao T:

    Jpn J Appl Phys. 2000, 39:4435. Publisher Full Text OpenURL

  36. Chou PC, Pagano NJ: Elasticity: Tensor, Dyadic, and Engineering Approaches. Dover Publications, New York; 1992. OpenURL

  37. Granados D, García JM:

    Appl Phys Lett. 2303, 82:2401. Publisher Full Text OpenURL

  38. Lorke A, Blossey R, García JM, Bichler M, Abstreiter G:

    Mater Sci Eng B. 2002, 88:225. Publisher Full Text OpenURL

  39. Teodoro MD, Campo VL Jr, Lopez-Richard V, Marega E Jr, Marques GE, Galvao Gobato Y, Iikawa F, Brasil MJSP, AbuWaar ZY, Dorogan VG, Mazur YI, Benamara M, Salamo GJ:

    Phys Rev Lett. 2010, 104:086401. PubMed Abstract | Publisher Full Text OpenURL