Advances in Condensed Matter Physics
Volume 2021 (2021), Article ID 3928308, 8 pages
https://doi.org/10.1155/2021/3928308

Research Article

Theoretical Study of Excitonic Complexes in GaAs/AlGaAs Quantum Dots Grown by Filling of Nanoholes

Mohamed Omri,¹ Amor Sayari,2,3 and Larbi Sfaxi4,5

¹DeanShip of Scientific Research (DSR), King Abdulaziz University, Jeddah 21589, Saudi Arabia
²Department of Physics, Faculty of Science, University of Jeddah, P.O. Box 80327, Jeddah 21589, Saudi Arabia
³Université de Tunis El Manar, Faculté des Sciences de Tunis, Unité de Recherche Spectroscopie Raman UR13ES31, Tunis 2092, Tunisia
⁴High School of Sciences and Technology of Hammam Sousse, Université de Sousse, Sousse, Tunisia
⁵Laboratoire de Micro-Optoélectroniques et Nanostructures de la Faculté des Sciences de Monastir, Université de Monastir, Monastir, Tunisia

Correspondence should be addressed to Amor Sayari; amor.sayari@laposte.net

Received 19 June 2021; Accepted 11 August 2021; Published 30 August 2021

Academic Editor: Da-Ren Hang

Copyright © 2021 Mohamed Omri et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In this work, a theoretical study of the electronic and the optical properties of a new family of strain-free GaAs/AlGaAs quantum dots (QDs) obtained by AlGaAs nanohole filling is presented. The considered model consists of solving the three-dimensional effective-mass Schrödinger equation, thus providing a complete description of the neutral and charged complex excitons’ fine structure. The QD size effect on carrier confinement energies, wave functions, and s-p splitting is studied. The direct Coulomb interaction impact on the calculated s and p states’ transition energies is investigated. The behaviour of the binding energy of neutral and charged excitons (X⁻ and X⁺) and biexciton XX versus QD height is studied. The addition of the correlation effect allows to explain the nature of biexcitons often observed experimentally.

1. Introduction

Optimizing the optical qualities of quantum dots (QDs) remains a challenge for researchers. For this purpose, new GaAs/AlGaAs semiconductor QDs have been grown using a new technique based on filling self-assembled nanoholes obtained by local etching of droplets [1, 2]. This method provides QDs with particular structural properties such as unstrained, uniform, not very dense, with high symmetry, and of different shapes with reference to conventional strained QDs. Due to these characteristics, the optical properties of unstrained QDs are interesting for various increasing high-performance applications in optoelectronic devices, namely, lasers and solar cells [3, 4], in addition to quantum cryptography [5].

In this context, we provide a modelling of the GaAs/AlGaAs QDs’ shape reported by Heyn et al. [2], compatible with available atomic force microscopy (AFM) images. The theoretical study in this paper will focus on the dependence of the confinement energy of the carriers, for the ground state s and state p, on the QD height (hQD). The calculations are performed by adapting a simple configuration based on BenDaniel and Duke, Hamiltonian [6], which requires much computation time. The obtained results will be compared to the literature experimental results for validation. In addition, our theoretical approach will be compared with Graf et al.’s [1] model. Thereby, we will then estimate the effect of the QDs’ size on the binding energy of the excitonic complexes in terms of direct Coulomb interaction and correlation effects. Therefore, the final alignment of excitonic states will be explained.

2. Theoretical Model

The electronic states of a QD strongly depend on the shape chosen for the dot and its symmetry. In the literature, unstrained self-assembled QDs have been modelled by different shapes such as cone [7, 8], pyramid with a square base [9, 10], and lens [11, 12]. The electronic structure of QDs has been calculated by adopting various approaches such as the pseudo-potential model proposed by Williamson and Zunger [13, 14], the strong bond model suggested by Lee et al. [15], and the formalism of the envelope function at several bands by Stier et al. [16].

In this paper, calculations are limited to the formalism of the envelope function with one band. This method consists in choosing a large quantization box, containing the studied physical system (QD), and having all the possible symmetries of the system (Figure 1). This quantization box is a cylinder of radius R and height H. These dimensions are chosen such as and to avoid side effects. The parameters and are, respectively, the radius and height of the QD (Figure 1). The Hamiltonian of each electron (e) or hole (h) is written in the form

where and are, respectively, the effective mass of the electron or hole and the potential for conduction and valence offset. and are the position vectors of electrons and holes inside the cylinder. The offset potential depends on the geometry of the QD; thus, we consider that the confinement potential has a Gaussian shape, and we describe it as follows:

with and .

Figure 1: Schematic illustration of the shape modelling of the GaAs QD obtained by filling the AlGaAs nanohole.

To access the energy levels of QDs, we used the matrix method initially proposed by Marzin and Bastard [8]. Due to the cylindrical symmetry of the QDs, the wave function is expressed as a Fourier–Bessel series written as

where are the coefficients to be determined and are periodic and orthogonal functions which are given by

where n, m, and are integers and is the Bessel function of order m and is its nth root.

The matrix element is defined as

with being the cylinder volume. The energy levels and wave functions are calculated using the parameters listed in Table 1.

Table 1: Numerical values of the input parameters of the Al0.36Ga0.44As barrier layer used in our calculations.

Electron
Heavy hole

0.067
0.093
313.20
0.51
0.60
205.20

mw and mb are the effective masses, respectively, in the quantum well material (GaAs) and in the barrier material (AlGaAs) given in free electron mass unit (me) and the offset of the conduction (valence) band Ve(h). The energy band gap of GaAs at low temperature is (GaAs) = 1519 meV.

The AlGaAs nanoholes are modelled by holes of the same shape and dimensions as those of [1]. Nanoholes of radius and depth have been considered in our calculations.

The excitonic energies of ground state s and excited state p are, respectively, denoted by and . Their expressions are given by

where is the GaAs QD band gap energy. Ei and HHi are the confinement energies of the electron and hole (i = 1 for ground state s and i = 2, 3, 4, ..., for exited states p, d, f, …, respectively). is the electron-hole Coulomb energy, in the s or p state, treated in the perturbative approach, and it is written as

with εr being the relative dielectric constant of GaAs [17] and and representing the volume of the cylinder associated with electrons and holes, respectively.

The recombination energies of the neutral exciton , charged excitons , and biexciton are calculated as follows [18]:

The binding energies of the trions , and the biexciton , without correlation effects, are estimated within the framework of the Hartree–Fock approximation and defined by

The new binding energies of the trions , and the biexciton , with correlation effects, are defined as follows [18]:

with , , , and being the correlation energies which are written in terms of in the following form:

with

where and and are the energies of a single particle (i or j) of the ground state. and are the energies of the excited states defined by the indices ni and mj.

3. Modelling Results and Discussion

3.1. Confinement Energy and Wave Function

In Figure 2, we have presented the first five calculated squared electron wave functions for a GaAs QD of 12.9 nm height. The plots show that the electron and hole wave functions are highly localized inside the GaAs QD. In Figure 3, we have represented the variation of electron (En) and heavy hole (HHn) confinement energies for the states s (n = 1) and p (n = 2) as functions of the QD height . From Figure 3, the confinement energies are very sensitive to the QD height. We observe a progressive decrease of these energies when increases. This behaviour is also observed for high excited states (n > 2) of the GaAs QD (Table 2). The electrons’ confinement energies are higher than those of heavy holes due to their lower effective mass compared to that of holes.

Figure 2: Top and side views of the probability densities of single-electron states in a GaAs QD with 12.9 nm height.

Figure 3: Electron (En) and hole (HHn) confinement energies for the fundamental s (n = 1) and the first excited p (n = 2) states as a function of the QD height .

Table 2: Energies of the ground state E₁ (HH₁) and the two excited states E₂ (HH₂) and E₃ (HH₃) for electrons (holes).

E₁ (meV)
HH₁ (meV)
E₂ (meV)
HH₂ (meV)
E₃ (meV)
HH₃ (meV)

4.3
150.8
71.5
167.0
77.5
184.4
83.9

6.4
87.7
30.9
102.3
34.4
119.0
39.4

8.6
56.8
12.9
70.0
15.9
85.1
20.1

10.7
42.7
7.9
54.3
10.2
67.4
13.1

12.9
33.6
5.6
44.0
7.3
56.0
9.3

3.2. Excitonic Energies

It is known that exciton formation in a QD is different than in a bulk crystal. By comparing the QD size to the Bohr radius, it is possible to define three regimes: strong confinement regime, weak confinement regime, and intermediate confinement regime. In our calculations, we will employ the strong confinement approximation, where the electron-hole Coulomb interaction is considered as a small perturbation against the single-particle terms in the Hamiltonian [9]. In our case, the nanoholes have a lateral shape with a Gaussian profile, their depth is typically 16 nm, and their radius is around 30 nm. The GaAs QDs have a lateral size similar to that of the nanoholes, but their height can vary between 3 and 12 nm. Therefore, in view of the 15 nm exciton Bohr radius in GaAs, the confinement along the growth axis (z-axis) is stronger than the lateral confinement in the (xy) plane. Yet, the wave function spatial dependence in the growth direction and lateral direction is only slight, but it is significantly affected by the Coulomb interaction. Thus, can be considered as a disturbance [9, 19].

In Table 3, the calculated values of are given as a function of QD height . We note that, for , increases as the height of the QDs decreases. This behaviour is similar to that observed by Abbarchi et al. [20]. The impact of the direct Coulomb interaction is more pronounced for state s than for state p because of the small spatial extension of the wave function of state s compared to that of state p. From Table 3, we attribute the decrease of , for , to the overlapping integrals of the wave functions associated with the electrons and holes of states s and p. The minimum recovery is obtained for . Similar behaviour has been observed in self-assembled InAs/GaAs QDs [20, 21].

Table 3: Variation of the Coulomb interaction term and the overlapping integrals and as a function of the QD height .

4.4
−14.94
−12.26
0.947
0.896

6.4
−15.58
−12.35
0.928
0.864

8.7
−14.34
−11.27
0.930
0.874

10.5
−12.86
−10.43
0.956
0.935

12.9
−11.54
−9.78
0.979
0.965

To investigate the accuracy of the numerical approach, the emission energies of states s and p are compared to those calculated by Graf et al. [1] and experimental data reported by Heyn et al. [2]. In Figure 4, we have represented the theoretical and experimental variations of the neutral exciton transition energies for the ground and first excited states as a function of . In our approach, we neglected the N body effects. This choice is justified based on the results of Heyn et al. [2], who showed that photoluminescence (PL) peaks shift slightly towards red by about 2 meV when the excitation power was increased. This value remains very low compared to the confinement and Coulomb energies. Graf et al. [1] theoretically studied the optical properties of GaAs QDs. Their approach is based on the eight-band k.p model, considering N body effects. Although their model is sophisticated, they neglected the band curvature of the top surface of the QD and considered it flat. This has a direct impact on the confinement energies of the carriers. However, an agreement is obtained between our theoretical results and experimental data reported by Heyn et al. as shown in Figure 4. This agreement demonstrates the suitability of our modelling approach.

(a)

(b)

Figure 4: Comparison of transition energies of the neutral exciton obtained by our model and those calculated by Graf et al. [1] and measured by Heyn et al. [2]: (a) for the ground state and (b) for the first excited state .

3.3. Excitonic Complexes

The binding energies of the excitonic complexes in this type of QDs have been calculated. Thereby, more particular interest is given to the neutral exciton X, the biexciton XX, and the positively X⁺ and negatively X⁻ charged excitons. These excitons have been the subject of several theoretical and experimental studies for different types of self-assembled QDs such as InGaAs/GaAs [21, 22] and InAs/InP [23]. In these high-confinement systems, the bond energies can be determined via micro-PL spectroscopy on single dots to study the effect of size and shape on the bond energies. Unlike self-assembled InAs QDs, the addition of correlation effect, in strain-free GaAs QDs, helps to explain the binding nature of the biexciton, often observed.

Figure 5 shows the variations of the correlation energies and the bond energies of the complex excitons, with and without correlation effect, as a function of QD height . The degree of correlation is specific to excitonic complexes and is sensitive to the variation of . The correlation effect for a biexciton is larger due to the higher number of charge carriers involved. We also obtain , in agreement with the prediction established by Schliwa et al. [22]. By comparing the bond energies without and with correlation effect, we observe the formation of bound excitonic states. Indeed, we underline a transition from an unbound state to a bound state for the biexciton and for the negatively charged exciton when the correlation effect is introduced. The binding energies of the biexciton obtained via our model are between −2.78 and −5.01 meV. These values are higher in magnitude than the experimental values (−1.3 to −2 meV) available in the literature [6, 24] for GaAs QDs based on nanoholes’ shape and size. However, our results agree with those of self-assembled III–V QDs, where it is well known that the binding energy of the biexciton varies between 1 and 6 meV [25, 26]. The comparison between experimental and theoretical data is complex given the lack of atomic force microscopy data for the studied GaAs QDs. Our study provides a comprehensive understanding of the correlation and size effects on the final alignment of excitonic states in strain-free GaAs QDs.

(a)

(b)

(c)

Figure 5: Variation of correlation energies (δ) of X, X⁺, and X⁻ excitons and biexciton XX as a function of (a). Bond energies of excitons X, X⁺, and X⁻ and biexciton XX as a function of (b) without correlation effect and (c) with correlation effect.

4. Conclusion

We have theoretically studied GaAs/AlGaAs QDs obtained by filling AlGaAs nanoholes with shape and profile provided by AFM measurements. The effect of size on carrier confinement, wave functions, and s-p splitting has been studied. The confinement energies of electrons and holes are very sensitive to the quantum dot height . The excitonic energies of the s and p states are calculated. The impact of the direct Coulomb interaction is more pronounced for state s than for p due to the small spatial extension of the wave function of state s compared to that of state p. An agreement is obtained between our theoretical results and experiment data of the literature, which indicates the suitability of our modelling approach. The behaviour of the binding energy of charged excitons (X⁻ and X⁺) and biexciton XX has been studied. The addition of the correlation effects allowed the formation of bound excitonic states as expected by experiments, and they are sensitive to the QD height. Our study provides a comprehensive understanding of the correlation and size effects on the final alignment of excitonic states in strain-free GaAs QDs.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This project was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah (Grant no. G:1715-305-1440). The authors, therefore, acknowledge with thanks the DSR for the technical and financial support.

References

A. Graf, D. Sonnenberg, V. Paulava, A. Schliwa, C. Heyn, and W. Hansen, “Excitonic states in GaAs quantum dots fabricated by local droplet etching,” Physical Review B, vol. 89, no. 11, Article ID 115314, 2014.

C. Heyn, A. Stemmann, and W. Hansen, “Dynamics of self-assembled droplet etching,” Applied Physics Letters, vol. 95, no. 17, Article ID 173110, 2009.

W. Rouis, M. Haggui, S. Rekaya, L. Sfaxi et al., “Local photocurrent mapping of InAs/InGaAs/GaAs intermediate-band solar cells using scanning near-field optical microscopy,” Journal of Solar Energy Materials and Solar Cells, vol. 144, pp. 324–330, 2016.

A. Sayari, M. Ezzidini, B. Azeza, S. Rekaya et al., “Improvement of performance of GaAs solar cells by inserting self-organized InAs/InGaAs quantum dot superlattices,” Journal of Solar Energy Materials and Solar Cells, vol. 113, pp. 1–6, 2013.

A. K. Ekert, “Quantum cryptography based on bell’s theorem,” Physical Review Letters, vol. 67, no. 6, pp. 661–663, 1991.

D. J. BenDaniel and C. B. Duke, “Space-charge effects on electron tunneling,” Physical Review, vol. 152, no. 2, pp. 683–692, 1966.

P. H. Lelong, O. Heller, and G. Bastard, “Multiparticle states and coulomb blockade in InAs/GaAs quantum dots,” Physica E:Low-Dimensional Systems and Nanostructures, vol. 2, no. 1-4, pp. 678–681, 1998.

J. Y. Marzin and G. Bastard, “Calculation of the energy levels in InAsGaAs quantum dots,” Solid State Communications, vol. 92, no. 5, pp. 437–442, 1994.

M. Grundmann, O. Stier, and D. Bimberg, “InAs/GaAs pyramidal quantum dots: strain distribution, optical phonons, and electronic structure,” Physical Review B, vol. 52, no. 16, pp. 11969–11981, 1995.

L. R. C. Fonseca, J. L. Jimenez, J. P. Leburton, and R. M. Martin, “Self-consistent calculation of the electronic structure and electron-electron interaction in self-assembled InAs-GaAs quantum dot structures,” Physical Review B, vol. 57, no. 7, pp. 4017–4026, 1998.

S. Sauvage, P. Boucaud, J. M. Gérard, and V. Thierry-Mieg, “In-plane polarized intraband absorption in InAs/GaAs self-assembled quantum dots,” Physical Review B, vol. 58, no. 16, pp. 10562–10567, 1998.

M. Baira, L. Sfaxi, L. Bouzaiene, H. Maaref, N. Chauvin, and C. Bru-Chevallier, “Broken symmetry in laterally coupled InAs/GaAs quantum dots molecule,” Journal of Applied Physics, vol. 104, no. 6, Article ID 064314, 2008.

A. J. Williamson and A. Zunger, “InAs quantum dots: predicted electronic structure of free-standing versus GaAs-embedded structures,” Physical Review B, vol. 59, no. 24, pp. 15819–15824, 1999.

A. J. Williamson and A. Zunger, “Pseudopotential study of electron-hole excitations in colloidal free-standing InAs quantum dots,” Physical Review B, vol. 61, no. 3, pp. 1978–1991, 2000.

S. Lee, L. Jönsson, J. W. Wilkins, G. W. Bryant, and G. Klimeck, “Electron-hole correlations in semiconductor quantum dots with tight-binding wave functions,” Physical Review B, vol. 63, no. 19, Article ID 195318, 2001.

O. Stier, M. Grundmann, and D. Bimberg, “Electronic and optical properties of strained quantum dots modeled by 8-band k⋅p theory,” Physical Review B, vol. 59, no. 8, pp. 5688–5701, 1999.

W. J. Moore and R. T. Holm, “Infrared dielectric constant of gallium arsenide,” Journal of Applied Physics, vol. 80, no. 12, pp. 6939–6942, 1996.

G. A. Narvaez, G. Bester, and A. Zunger, “Excitons, biexcitons, and trions in self-assembled (In, Ga) As/GaAs quantum dots: recombination energies, polarization, and radiative lifetimes versus dot height,” Physical Review B, vol. 72, no. 24, Article ID 245318, 2005.

R. Seguin, A. Schliwa, S. Rodt, K. Pötschke, U. W. Pohl, and D. Bimberg, “Size-dependent fine-structure splitting in self-organized InAs/GaAs quantum dots,” Physical Review Letters, vol. 95, no. 25, Article ID 257402, 2005.

M. Abbarchi, T. Kuroda, T. Mano et al., “Energy renormalization of exciton complexes in GaAs quantum dots,” Physical Review B, vol. 82, no. 20, Article ID 201301R, 2010.

R. Neffati, I. Saïdi, S. B. Radhia, K. Boujdaria, and C. Testelin, “Effect of AlGaAs composition on the type-I to type-II transition in AlInAs/AlGaAs self-assembled quantum dots,” Semiconductor Science and Technology, vol. 30, no. 8, Article ID 085008, 2015.

A. Schliwa, M. Winkelnkemper, and D. Bimberg, “Few-particle energies versus geometry and composition of InxGa1−xAs/GaAs self-organized quantum dots,” Physical Review B, vol. 79, no. 7, Article ID 075443, 2009.

M. Gong, W. Zhang, G. Can Guo, and L. He, “Atomistic pseudopotential theory of optical properties of exciton complexes in InAs/InP quantum dots,” Applied Physics Letters, vol. 99, no. 23, 231106 pages, 2011.

C. Heyn, M. Klingbeil, C. Strelow, A. Stemmann, S. Mendach, and W. Hansen, “Single-dot spectroscopy of GaAs quantum dots fabricated by filling of self-assembled nanoholes,” Nanoscale Research Letters, vol. 5, no. 10, pp. 1633–1636, 2010.

K. Brunner, G. Abstreiter, G. Böhm, G. Tränkle, and G. Weimann, “Sharp-line photoluminescence and two-photon absorption of zero-dimensional biexcitons in a GaAs/AlGaAs structure,” Physical Review Letters, vol. 73, no. 8, pp. 1138–1141, 1994.

A. Kuther, M. Bayer, A. Forchel, A. Gorbunov et al., “Zeeman splitting of excitons and biexcitons in single In0.60Ga0.40As/GaAs self-assembled quantum dots,” Physical Review B, vol. 58, no. 12, pp. R7508–R7511, 1998.