N A N O E X P R E S S
Impacts of Coulomb Interactions on the Magnetic Responses
of Excitonic Complexes in Single Semiconductor Nanostructures
Wen-Hao Chang• Chia-Hsien Lin• Ying-Jhe Fu•Ta-Chun Lin•Hsuan Lin•Shuen-Jen Cheng•
Sheng-Di Lin•Chien-Ping Lee
Received: 27 November 2009 / Accepted: 5 January 2010 / Published online: 21 January 2010 Ó The Author(s) 2010. This article is published with open access at Springerlink.com
Abstract We report on the diamagnetic responses of different exciton complexes in single InAs/GaAs self-assembled quantum dots (QDs) and quantum rings (QRs). For QDs, the imbalanced magnetic responses of inter-par-ticle Coulomb interactions play a crucial role in the dia-magnetic shifts of excitons (X), biexcitons (XX), and positive trions (X-). For negative trions (X-) in QDs, anomalous magnetic responses are observed, which cannot be described by the conventional quadratic energy shift with the magnetic field. The anomalous behavior is attributed to the apparent change in the electron wave function extent after photon emission due to the strong Coulomb attraction by the hole in its initial state. In QRs, the diamagnetic responses of X and XX also show different behaviors. Unlike QDs, the diamagnetic shift of XX in QRs is considerably larger than that of X. The inherent structural asymmetry combined with the inter-particle Coulomb interactions makes the wave function distribution of XX very different from that of X in QRs. Our results suggest that the phase coherence of XX in QRs may survive from the wave function localization due to the structural asym-metry or imperfections.
Keywords Quantum dots Quantum rings Magnetophotoluminescence Diamagnetic shift
Introduction
Single semiconductor nanostructures, such as quantum dots (QDs) and quantum rings (QRs), have attracted much attention due not only to their fundamental interest, but also to their potential applications in the prospective quantum information technology [1–8]. In semiconductor nanostructures, excitonic effects play a central role in their optical properties [9–15]. Electrons and holes form a variety of exciton complexes in QDs and QRs due to the enhanced Coulomb interactions by the spatial confine-ments. The spatial extents of excitonic wave functions thus reflect the both the combined effects of spatial confine-ments and the inter-particle Coulomb interactions among constituent charged carriers [9, 10, 13, 14]. Applying a magnetic field B is an efficient way to probe the spatial extent of excitonic wave functions. In a magnetic field B, the exciton emission energy increases quadratically with B, i.e., the diamagnetic shift DE = cB2, with a diamagnetic coefficient c proportional to the area of the excitonic wave function. Depending on the confinement regime of the nanostructures, the measured c could have different phys-ical meanings [9]. In the strong confinement regime, where the single-particle energy dominates over the Coulomb energies, c is a measure of the spatial confinement of the QDs, while the magnetic responses of inter-particle Cou-lomb energies only appear as correction terms to the overall diamagnetism. On the other hand, when confined excitons are in the weak confinement regime, the Coulomb energies become prominent, such that c will be mainly determined by the magnetic response of the inter-particle Coulomb interactions. For exciton complexes confined in nanostructures [14, 15], the diamagnetic behaviors are more complicated because of the more elaborate Coulomb interactions among the constituent charged carriers.
W.-H. Chang (&) C.-H. Lin H. Lin S.-J. Cheng
Department of Electrophysics, National Chiao Tung University, Hsinchu 30010, Taiwan
e-mail: whchang@mail.nctu.edu.tw Y.-J. Fu T.-C. Lin S.-D. Lin C.-P. Lee
Department of Electronics Engineering, National Chiao Tung University, Hsinchu 30010, Taiwan
Therefore, a systematic investigation is essential into fur-ther understanding the magnetic responses of different exciton complexes in QDs or QRs.
In this work, we report on the diamagnetic response of different exciton complexes, including neutral excitons (X), biexcitons (XX), positive and negative trions (X?and X-) confined in single InAs/GaAs self-assembled QDs and QRs. In QDs, the influences of imbalanced magnetic responses of Coulomb energies on the overall diamagne-tisms of X, XX, and X?are observed and discussed. As for the magnetic response of X- in QDs, we observe an anomalous diamagnetic behavior, which cannot be descri-bed by the conventional quadratic energy shift with the applied B. Such anomalous behaviors for X- can be attributed to the apparent change in the electron wave function extent after photon emission due to the strong Coulomb attraction by the hole in its initial state. In QRs, the diamagnetic responses of X and XX also show different behavior. The impacts of inherent structural asymmetry of the self-assembled QRs, as well as the inter-particle Cou-lomb interactions, on the distribution of X and XX wave functions are discussed. Our results suggest that the phase coherence of XX in QRs may survive from the wave function localization due to the structural asymmetry or imperfections.
Experimental Details
The investigated samples were grown on GaAs (100) substrates by a molecular beam epitaxy system. For the QD sample, a layer of InAs self-assembled QDs (2.0 MLs) was grown on GaAs at 480°C without substrate rotations, yielding a gradient in area dot density ranging from 108to 1010cm-2. Atomic force microscopy (AFM) of uncapped samples reveals that the dots are lens shaped, with an average height and diameter of &2(±0.5) and &15 nm, respectively. The QDs were finally capped by a 100-nm undoped GaAs layer.
For the QR sample, a low-density InAs QD layer was first grown by depositing 2-ML InAs at 520°C as QR precursors. The substrate temperature was then lowered to 500°C for depositing a thin GaAs layer (1.7 nm) to cover the QD sidewalls, followed by a 50-s growth interruption for expelling the indium atoms from the QD center for the QR formation [16, 17]. Surface topography of uncapped QRs has also been investigated by AFM. The area density of surface QRs is about 107cm-2. Typical surface QRs have a rim diameter of 35 nm, a height of &1.3 nm, and a center dip of about 2 nm. The realistic dimension of the embedded QRs is expected to be much smaller [18]. We also found that the QR is anisotropic; the rim of the surface QR is higher along [110].
To isolate individual QDs or QRs for optical measure-ments, a 100-nm-thick Al metal mask was fabricated on the sample surface with arrays of 0.3-lm apertures using electron-beam lithography. Single QD spectroscopes were carried out at 5–8 K in a specially designed micro-photo-luminescence (l-PL) setup (see Fig.1), where the sample was mounted in a low-temperature stage and inserted into the bore of 6-T superconducting magnet for magneto-PL measurements. A He–Ne laser beam or a cw Ti:Sapphire laser was focused onto the aperture through a microscope objective (NA = 0.5). The PL signals were collected by the same objective, analyzed by a 0.75-m grating mono-chromator, and detected by a liquid-nitrogen-cooled charged-coupled device (CCD) camera, which yield a resolution-limited spectral linewidth of about 60 leV.
Results and Discussions Quantum Dots
Several apertures containing only one QD have been investigated, and all of which showed similar spectral features. Typical single QD spectra consists of four dom-inant emission lines associated with the recombination of neutral excitons (X), biexcitons (XX), positive and negative trions (X? and X-). These excitonic spectra have been unambiguously identified according to power-dependent and polarization-resolved PL measurements [14]. Mag-neto-PL measurements have been taken on a total of seven QDs with X energy distributed over the range of 1,349– 1,385 meV. Typical spectra selected from one particular QD are shown in Fig.2. When a magnetic field, B, was applied along the QD growth direction (Faraday geometry), each line splits into a doublet through the Zeeman effect. In
addition, the average energy of each Zeeman doublet increases with the increasing B, known as the diamagnetic shift. In Fig.3, the measured diamagnetic shifts for the four excitonic emission lines are plotted as a function of B2. For X, X?, and XX, the measured diamagnetic shifts show a quadratic dependence DE = cB2, with a clear trend of cX[ cXXffi cXþ. This trend holds for all investigated dots
and is a consequence of different magnetic responses of inter-particle Coulomb energies.
To explain this behavior, we consider here a 2D parabolic potential as the effective lateral confinement for both elec-trons and holes in the QDs, with a quantization energy of hxb
and a wave function extent of ‘b¼
ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi h=mbxb
p
, where b denotes e or h and mbare the effective masses [20,21]. The
use of parabolic model simplifies the problem considerably, because of the availability of analytical formula for all Coulomb interactions being parameterized by ‘b. The
inter-particle direct Coulomb energy is given by [20–22]
Vab¼ e2 4pee0 ffiffiffi p 2 r 1 ‘ where ‘ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ‘2 aþ ‘2b =2 r
. In a magnetic field B, xbcan be
replaced by the B-dependent hybridize frequency XbðBÞ ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi x2 bþ ðeB=2mbÞ 2 q
. In the weak-field limit, the
B-dependent single-particle energy hXb(B) and the
Coulomb energies Vab(B) can be expanded analytically as,
hXbðBÞ hxbþ cSPb B 2þ ð1Þ VabðBÞ Vabð0Þ þ cCoulab B 2þ ð2Þ where cSP
b ¼ e2‘2b=8mb is the single-particle diamagnetic
coefficient, and cCoul
ab ¼ kð‘6aþ ‘6bÞ=2‘3 accounts for the
magnetic response of Vab(B), with a constant k defined as
e4pffiffiffiffiffiffiffiffip=2=64h2
pee0
. If we take the Coulomb energies as perturbations to the single-particle energy in the strong confinement regime, the diamagnetic shift of X should be corrected as cX ¼ cSPe þ c
SP h
cCoul
eh . On the other hand,
the diamagnetic shift of XX will deviate from that of X by an amount of Dc¼ cX cXX¼ 2cCouleh c Coul ee c Coul hh . Because cCoul ab varies as *‘
3, a slight difference between ‘ e
and ‘hcan lead to very different values for cCoulee , c Coul eh and
cCoul
hh . Since ‘h\ ‘e and hence ‘3h ‘3e, we have cCoulee
cCoul eh c
Coul
hh that leads to Dc¼ cX cXX cCoulee . The
same arguments can be applied to the X?case, where cXþis
also reduced by a similar amount, in agreement with our experimental observations. The reason for cXþ cXX can
also be understood from the fact that the emission energies of XX and X?are usually correlated even in different QDs [23], because their energy difference EXX EXþ equals
roughly to Vee–Veh, which is value less sensitive to the
magnetic field.
It is worth to emphasize that the systematic differences among cX, cXX, and cXþcould be observed only when ‘eand
‘hexhibit a large difference. For small InAs QDs, electron
gradually lose confinement as the dot size reduces, which in turn push the electron level toward the wetting-layer continuum, resulting in a more extended electron wave function penetrating into the barrier material [24]. On the other hand, the hole wave function remains well localized in the QDs, leading to a large difference between ‘eand ‘h.
For larger In(Ga) As QDs, since ‘eand ‘hare expected to
be similar, the diamagnetic response of all exciton com-plexes will be nearly identical [13], so that the magnetic response of interparticle Coulomb energies of such strongly confined few-particle systems becomes unable to resolve experimentally.
For the case of X-in QDs, its magnetic response shows an anomalous behavior, which does not follow the con-ventional B2-dependence energy shift, as can be seen from Fig.3b. Our measurements show that the X-diamagnetic shifts in most QDs investigated are smaller and non-qua-dratic. To understand the peculiar behavior, the charge configurations of the initial and the final states of X -should be considered. As depicted in Fig. 4, the initial state of X-consists of two electrons and one hole, which will
Fig. 2 Magneto-PL spectra for a typical QD measured under an applied B = 0–6 T
Fig. 3 aThe energy shifts of X, X?, and XX in a QD as a function of B2. b The energy shifts of X- as a function of B2. A clear non-quadratic B dependence can be seen
leave one electron in its final state after recombination. The X- diamagnetic shift thus reflects the diamagnetic responses of both the initial and final states. We have performed numerical model simulations using single-band Hamiltonian of cone-shape QDs including Hartree approximations to calculate the lateral extents of the initial-state ‘e,i and the final-state ‘e,felectron wave functions of
X-. We found that the ‘e,fis always more extended than
‘e,i. This can be realized from the presence of hole in its
initial state, which contracts the electron wave function by the Coulomb attraction. As the QD sizes reduces, ‘e,f
increases rapidly, due to the penetration of electron wave function into the barrier material [24]. The very extended electron becomes sensitive to the long-range Coulomb attractive potential produced by the hole in the initial state, by which ‘e,i will be contracted and become apparently
smaller than ‘e,f. As a result, the final-state diamagnetic
shift increases, so that the overall diamagnetic shift in X-is reduced. This phenomenon is not expected for X?in QDs, since the wave function extents of the initial- and final-state holes are almost identical due to the fact that hole has a larger effective mass, such that the hole wave function remains well confined even in small QDs. Therefore, the Coulomb attractive potential produced by the weakly confined electron becomes less important.
Now, we explain why the X-diamagnetic shift exhibits a non-quadratic B dependence. To illustrate this behavior, we consider the diamagnetic shift in the carrier’s single-particle energy and expand it in powers of B as, hXbðBÞ ¼
hxbþ cSPb B 2þ jSP
b B
4þ , where the quartic coefficient
jSP
b is negative and is given by j SP b ¼ e
4‘6
b=128mbh 2. By
taking into account the difference between ‘e,i and ‘e,f, a
simple algebraic analysis for the X- diamagnetic shift, DEXðBÞ, gives the following expression:
DEXðBÞ cXB2þ jXB4þ ; ð3Þ
where cX¼ 2cSPe;iþ cSPh;i
cSP
e;f and jX¼ 2jSPe;iþ jSPh;i
jSP
e;f. Because of ‘h\ ‘e and mh me, cSPh;i and jSPh;i only
have minor influences on the overall diamagne-tism. Accordingly, we obtain cX 2cSPe;i cSPe;f ¼
e22‘2 e;i ‘2e;f
=8me and jX 2jSPe;i jSPe;f ¼ e22‘6e;i
‘6 e;f
=128meh2. Equation (3) makes clear how the
dif-ference in ‘e,i and ‘e,fcan lead to anomalous diamagnetic
shifts for X-. For a normal case of ‘e,i& ‘e,f= ‘e, we have
cX ce cXand cX jX. That is, the X-diamagnetic
shift behaves as the usual B2dependence with a coefficient similar to that of X. A very interesting case occurs whenffiffiffi
2 p
‘e; i¼ ‘e;f, which leads to 2ce,i= ce,fand cancels out the
B2term, resulting in a dominant B4-dependent energy shift. The energy shift shown in Fig.3(b) exhibits both quadratic and quartic B dependence. Therefore, the measured cur-vature can also be a sensitive probe to the change in wave function extents between the initial and the final-state electrons of X-.
Quantum Rings
Figure5a and b show the surface topography and line scans along [110] and [110] directions of uncapped QRs. Typical single-QR magneto-PL spectra are shown in Fig.5c. Two emission lines associated with the recombination from X and XX can be observed. The same measurement has been per-formed on a total of seven QRs with X emission energies ranging from 1,320 to 1,328 meV. Apart from the Zeeman splitting, both X and XX show clear diamagnetic shifts with the increasing B. The average diamagnetic coefficient of X
Fig. 4 A schematic for the charge configurations of X- in QDs before and after photon emissions. The final-state electron is more extended due the Coulomb attraction of the hole in its initial state
Fig. 5 aSurface topography and b line scans along [110] and [110] directions of uncapped QRs. c Typical single-QR magneto-PL spectra taken under different applied B. d The B-dependent energy shifts of the X and XX. Solid and dashed lines are fitting curves using quadratic B dependence
for all the investigated QRs is cX ¼ 6:8 leV=T2.
Interest-ingly, we found that XX shows a considerably larger dia-magnetic coefficient with an average value of cXX¼
14:8 leV=T2, which is more than the double of the c
Xvalue.
The larger cXX indicates a more extended XX wave
function in QRs. This is different from the case of QDs, where cXXis usually similar to, or somewhat smaller than,
cX due to the different spatial extents of the electron and
hole wave functions and their responses to the applied B. Here, we argue that the larger cXXis a consequence of the
fact that self-assembled QRs do not have perfect azimuthal symmetry. It has been reported that buried self-assembled QRs exhibit an asymmetric craterlike shape, with a diam-eter substantially smaller than the ring-shaped islands on the surface of uncapped QR structures [18]. In addition, due to the preferential out diffusion of dot material along the [110] direction, the embedded rim height is higher along the [110] direction, resulting in two separate potential valleys along the [110] direction, resembling a pair of connected QDs [19]. Based on the proposed shape for embedded QRs in Ref. [19], we have performed theoretical calculations for the wave function distributions of X and XX in QRs, which are shown in Fig.6. Although the actual shape of embedded QRs is not known exactly, the numerical model used here is still a good illustration when QRs have any structural imperfection. Using selective chemical etching to remove the GaAs capping layer may be an efficient way to reveal the structure of buried QRs [25,26].
The lack of rotational symmetry in the potential of embedded QRs has significant impacts on the diamagnetic responses of X and XX. For X confined in QRs, the height variation strongly localizes the hole in one of the potential valleys due to the large effective mass. Consequently, the electron will be bound to the same valley by the electron-hole Coulomb attraction. This means that the wave function extent of X is determined mainly by the confinement of the potential valley and the Coulomb interaction. Therefore, the diamagnetic response of X is similar to an elliptic QD.
For XX confined in QRs, the two holes may be sepa-rately localized in different valleys due to the strong hole– hole Coulomb repulsion. Due to the Coulomb attractions of the two separately localized holes, the electron wave functions of XX are more likely to spread over the two valleys and become more extended than that of X. Unlike X, the wave function extent of XX is determined mainly by the diameter of the embedded QRs. This explains why cXX
is significantly larger than cXfor QRs.
Before conclusions, we would like to comment on the implications of the larger diamagnetic shift of XX in QRs. Charged particles confined to a QR are expected to show the well-known Aharanov–Bohm (AB) effect due to the quantum interference of the carrier’s wave function in the ring-shaped geometry [6–8]. However, for an exciton confined in ring-like nanostructures, the AB effect is not expected to occur for such a charge-neutral composite unless the electron and hole can propagate coherently in different trajectories with a non-zero electric dipole moment [27]. This can occur naturally in type-II QD sys-tems, where the electron and hole are spatially separated [28,29]. However, for InGaAs self-assembled QRs where both the electron and hole are confined, it is still an open question whether the excitonic AB effect can be observed. In fact, this issue is further complicated by the inherent structural asymmetry and/or imperfections presented inevitably in self-assembled QRs [18, 19]. Our results suggest that the inherent structural asymmetry combined with the inter-particle Coulomb interactions tends to localize the X wave function and hence may smear out its phase coherence. On the other hand, due to hole–hole repulsion, the XX wave function is more extended and spread over the entire ring. Therefore, the phase coherence of XX may be more likely to survive from such wave function localizations at lower magnetic field.
Conclusions
In conclusion, the diamagnetic response of different exciton complexes, including neutral excitons (X), biexcitons (XX), positive and negative trions (X?and X-) confined in single InAs/GaAs self-assembled QDs and QRs has been
Fig. 6 a A schematic for the QR geometry with a structural anisotropy used in model calculations. b and c are the calculated electron wave functions of the lowest X and XX states. The dash circle represents the rim diameter of 14 nm
investigated. In QDs, we found a systematic trend of cX[ cXX cXþ caused by the imbalanced magnetic
responses of inter-particle Coulomb interactions, which could be observed only when the confined electron and hole wave functions exhibit a large difference in their lateral extents. Furthermore, the difference Dc¼ cX cXX was
found to scale as the cube of its single-particle wave function extent and therefore can be a sensitive probe to the electron-hole wave function asymmetry. On the other hand, the magnetic response of X-in QDs fell into a special regime where the conventional quadratic diamagnetic shift failed to describe its magnetic response. Our measurements show that the X-diamagnetic shifts in most QDs investigated are even smaller and non-quadratic. Such anomalous behaviors for X-were explained by the apparent change in the electron wave function extent after photon emission due to the strong Coulomb attraction induced by the hole in its initial state. In QRs, the diamagnetic responses of X and XX also show different behavior. Unlike single QDs, XX in QRs shows a considerably larger diamagnetic coefficient than the X. Numerical model calculations indicate that the inherent structural asymmetry and imperfection of the self-assembled QRs, combined with the inter-particle Coulomb interactions, play a crucial role in the distribution of X and XX wave functions. Our results suggest that the phase coherence of XX in QRs may survive from the wave function localization due to the structural asymmetry or imperfections.
Acknowledgment This work was supported in part by the program of MOE-ATU and the National Science Council of Taiwan under Grant Nos.: NSC-97-2112-M-009-015-MY2, NSC-97-2120-M009-004 and NSC-97-2221-E009-161. We also acknowledge the Center for Nanoscience and Technology (CNST) at National Chiao Tung University.
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which per-mits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
References
1. P. Michler, A. Kiraz, C. Becher, W.V. Schoenfeld, P.M. Petroff, L. Zhang, E. Hu, A. Imamoglu, Science 290, 2282 (2000) 2. M. Pelton, C. Santori, J. Vuckovic, B. Zhang, G.S. Solomon,
J. Plant, Y. Yamamoto, Phys. Rev. Lett. 89, 233602 (2002)
3. Z. Yuan, B.E. Kardynal, R.M. Stevenson, A.J. Shields, C.J. Lobo, K. Cooper, N.S. Beattie, D.A. Ritchie, M. Pepper, Science 295, 102 (2002)
4. W.-H. Chang, W.-Y. Chen, H.-S. Chang, T.-P. Hsieh, J.-I. Chyi, T.-M. Hsu, Phys. Rev. Lett. 96, 117401 (2006)
5. X. Li, Y. Wu, D. Steel, D. Gammon, T.H. Stievater, D.S. Katzer, D. Park, C. Piermarocchi, L.J. Sham, Science 301, 809 (2003) 6. A. Lorke, R.J. Luyken, A.O. Govorov, J.P. Kotthaus, J.M. Garcı´a,
P.M. Petroff, Phys. Rev. Lett. 84, 2223 (2000)
7. M. Bayer, M. Korkusinski, P. Hawrylak, T. Gutbrod, M. Michel, A. Forchel, Phys. Rev. Lett. 90, 186801 (2003)
8. N.A.J.M. Kleemans, I.M.A. Bominaar-Silkens, V.M. Fomin, V.N. Gladilin, D. Granados, A.G. Taboada, J.M. Garcı´a, P. Offermans, U. Zeitler, P.C.M. Christianen, J.C. Maan, J.T. Devreese, P.M. Koenraad, Phys. Rev. Lett. 99, 146808 (2007) 9. S.N. Walck, T.L. Reinecke, Phys. Rev. B 57, 9088 (1998) 10. M. Bayer, S.N. Walck, T.L. Reinecke, A. Forchel, Phys. Rev. B
57, 6584 (1998)
11. Y. Naganume, Y. Arakawa, S. Tsukamoto, M. Nishioka, S. Sasaki, N. Miura, Phys. Rev. Lett. 69, 2963 (1992)
12. T. Someya, H. Akiyama, H. Sakaki, Phys. Rev. Lett. 74, 3664 (1995)
13. C. Schulhauser, D. Haft, R.J. Warburton, K. Karrai, A.O. Govorov, A.V. Kalameitsev, A. Chaplik, W. Schoenfeld, J.M. Garcia, P.M. Petroff, Phys. Rev. B 66, 193303 (2002)
14. M.-F. Tsai, H. Lin, C.-H. Lin, S.-D. Lin, S.-Y. Wang, M.-C. Lo, S.-J. Cheng, M.-C. Lee, W.-H. Chang, Phys. Rev. Lett. 101, 267402 (2008)
15. T.-C. Lin, C.-H. Lin, H.-S. Ling, Y.-J. Fu, W.-H. Chang, S.-D. Lin, C.-P. Lee, Phys. Rev. B 80, 081304(R) (2009)
16. D. Granados, J.M. Garcı´a, Appl. Phys. Lett. 82, 2401 (2003) 17. H.S. Ling, C.P. Lee, J. Appl. Phys. 102, 024314 (2007) 18. P. Offermans, P.M. Koenraad, J.H. Wolter, D. Granados, J.M.
Garcı´a, V.M. Fomin, V.N. Gladilin, J.T. Devreese, Appl. Phys. Lett. 87, 131902 (2005)
19. V.M. Fomin, V.N. Gladilin, S.N. Klimin, J.T. Devreese, N.A.J.M. Kleemans, P.M. Koenraad, Phys. Rev. B 76, 235320 (2007) 20. S.-J. Cheng, W. Sheng, P. Hawrylak, Phys. Rev. B 68, 235330
(2003)
21. S.-J. Cheng, Phys. Rev. B 76, 075329 (2007) 22. R.J. Warburton et al., Phys. Rev. B 58, 16221 (1998)
23. S. Rodt, A. Schliwa, K. Po¨tschke, F. Guffarth, D. Bimberg, Phys. Rev. B 71, 155325 (2005)
24. K.L. Janssens, F.M. Peeters, V.A. Schweigert, Phys. Rev. B 63, 205311 (2001)
25. Z.M. Wang, L. Zhang, K. Holmes, G.J. Salamo, Appl. Phys. Lett. 86, 143106 (2005)
26. F. Ding, L. Wang, S. Kiravittaya, E. Mu¨ller, A. Rastelli, O.G. Schmidt, Appl. Phys. Lett. 90, 173104 (2007)
27. A.O. Govorov, S.E. Ulloa, K. Karrai, R.J. Warburton, Phys. Rev. B 66, 081309R (2002)
28. E. Ribeiro, A.O. Govorov, W. Carvalho Jr, G. Medeiros-Ribeiro, Phys. Rev. Lett. 92, 126402 (2004)
29. I.R. Sellers, V.R. Whiteside, I.L. Kuskovsky, A.O. Govorov, B.D. McCombe, Phys. Rev. Lett. 100, 136405 (2008)