# Donor Impurity-Related Optical Absorption in GaAs Elliptic-Shaped Quantum Dots

Authors: M. A. Londoño; R. L. Restrepo; J. H. Ojeda; Huynh Vinh Phuc; M. E. Mora-Ramos; E. Kasapoglu; A. L. Morales; C. A. Duque
Journal: Journal of Nanomaterials
 (2017)

---

## Abstract

The conduction band and electron-donor impurity states in elliptic-shaped GaAs quantum dots under the effect of an externally applied electric field are calculated within the effective mass and adiabatic approximations using two different numerical approaches: a spectral scheme and the finite element method. The resulting energies and wave functions become the basic information needed to evaluate the interstate optical absorption in the system, which is reported as a function of the geometry, the electric field strength, and the temperature.


---

## Body

## 1. Introduction

Semiconductor elliptical quantum dots (QDs) have been the subject of investigation for a number of years due to their prospective applications in optoelectronics. Recent works on electronics and optical properties in this kind of nanosystems can be referred to in [1–10]. Among them it is worth highlighting, for instance, the practical realization of elliptical QDs in the InGaN/GaN system which can be used as single photon sources, allowing predefining photon states according to the QD orientation [9]. This points at the prospective applications of this kind of nanostructures in fields such as quantum information.The calculation of charge carrier states in 3D confined systems heavily depends on the geometry of the structure. In the particular case of elliptic-like QDs the analytical solution of effective mass equations is not possible in general. Furthermore, the inclusion of the effect of an externally applied electric field makes this possibility be unreachable. Therefore, numerical ways of solution for conduction and valence band states are required. In this sense, different approaches appear reported: finite- and boundary-element calculations [6, 7], variational solutions [8], and meshless schemes for the numerical solution of the effective mass equation [11, 12], to mention but a few.It is well known that even in the case of high quality and high purity samples, semiconductor compounds contain atoms of external elements acting as impurities, not to mention the intentional doping aimed at obtaining a desired carrier concentration in the material. In consequence, the investigation of the influence of impurity atoms on the spectrum of carrier states in semiconductor-based systems is always of significance. On the other hand, the optical response associated with quantum confined carrier states is an important element for both the understanding of the very energy spectrum of the systems and the design of their practical applications. In the case of elliptic-shaped systems a numerical calculation of optical transitions in InAs/GaAs QDs is reported, for example, in [6]. By using the adiabatic approximation, Gusev et al. have reported the electronic structure in low dimensional systems, with parabolic and rectangular potential, including impurity and exciton states [13, 14]. Their results are developed for quantum wells, wires, and dots, with particular spheroidal-shapes [14]. Their findings are in excellent agreement with variational calculations.In this article we are going to present the results of a study on the electron and electron-impurity states in elliptic-shaped GaAs-based QDs, including the influence of external static electric fields. Two kinds of structures are considered, for which Figure1 depicts their corresponding pictorial views. The calculation is carried out within the framework of the effective mass and parabolic band approximations. With the outcome information we shall discuss the light absorption related with transitions between a number of the lowest confined electron and electron-impurity states. Accordingly, the organization of the paper is as follows: Section 2 is devoted to presenting in detail the theoretical framework. Section 3 contains the obtained results and the corresponding discussions, and, finally, in Section 4 we give the conclusions.Figure 1
(a, b) Pictorial view of the two kinds of quantum dots considered in the present work: with constant height (a) and with variable height (b).
(a)
(b)


## 2. Theoretical Model

Here we shall consider the motion of conduction band electrons in the 3D elliptically shaped QD under the effect of in-plane applied electric field[F→=(Fx,Fy,0)] and considering the presence of a shallow-donor impurity. The Hamiltonian of the system, within the framework of the effective mass and parabolic band approximations, is given by(1)H=-ħ22m∥∗∇ρ2-ħ22m⊥∗∂2∂z2+Vρ,z+eF→·r→-e2εrei.The quantities |e|, m∥∗ (m⊥∗), and ε are the absolute value of the electron charge, the in-plane (growth direction) electron effective mass, and the GaAs static dielectric constant, respectively. ∇ρ2 is the in-plane Laplacian operator, r→=(x,y,z) is the electron position, and rei is the electron-impurity distance, with the donor impurity localized at (xi,yi,0). Considering that the height of the structure is very small comparing with the dimensions of the base, we can use the adiabatic approximation where the 3D-confinement potential can be written as the sum between two potentials: the first one depending on the in-plane xy-dimensions and the second one depending on the z-dimension; that is, V(ρ,z)=V1(ρ)+V2(z;x,y), where the z-dependent confinement potential (V2(z;x,y)) depends of the in-plane coordinates (x,y) (in this work this dependence is associated with the variations of the height of the QD). V1(ρ) is zero into the region delimited by S1 and infinite outside the region delimited by S1, where S1 is given by x2/Rx2+y2/Ry2=1 with z=0. Note that the frame of reference is localized at the gravity center of the elliptic region. Also, due to the restriction of heights with respect to the base of the dot, we take the approximation rei≃|(x,y,0)-(xi,yi,0)| for the electron-impurity distance. This approximation does not mean changes in the calculated Coulomb correlation more than 2% with respect to the corresponding exact value.According to the schematic view for the shape of the 3D elliptically shaped QD depicted in Figure1, two kinds of configurations will be studied: elliptical base with constant value of the height (a) and elliptical base with elliptical dependence of the height (b).The 3D problem in (1) does not have analytical solutions. To find the eigenfunctions and eigenvalues of the Schrödinger equation we will use the adiabatic approximation. In this paper we deal with heterostructures where the vertical dimension of the dot is small as compared with the in-plane dimensions. Under such conditions the electron has a fast-movement along the z-direction and slow-movement in the x-y plane. This is the spirit of the adiabatic approximation used in this paper where we decouple the two movements in (1). In this case we propose that the eigenfunctions of the Hamiltonian in (1) can be written as the product Ψ(x,y,z)=Nϕ(x,y)χ(z;x,y), where N is the normalization constant and ϕ(x,y) and χ(z;x,y) are, respectively, the eigenfunctions of the following two differential equations: (2)Hx,yϕx,y=-ħ22m∥∗∇ρ2+V1ρ+eF→·r→-e2εrei+E1h;x,yϕx,y=Eϕx,y,(3)-ħ22m⊥∗∂2∂z2+V2z;x,yχz;x,y=E1h;x,yχz;x,y,where E1(h;x,y) is the adiabatic potential energy and corresponds to the ground state of (3) which is a function of the in-plane coordinates. It is clear that the eigenvalues and eigenfunctions of (3) are the same as those for a single quantum well with width h(x,y). Three possibilities have been considered for the upper and lower surfaces (separated by the distance h(x,y)) in Figure 1: (i) a QW with two rigid potential barriers, (ii) a symmetric QW with finite potential barriers (V0), and (iii) an asymmetric structure with infinite and finite (V0) barriers. In the first case the adiabatic potential is given by(4)E1h;x,y=ħ2π22m⊥∗hx,y2whereas in the second and third cases, respectively, it is obtained by the first solution of the transcendental equations(5)tan⁡hx,yβ2-ξβ=0,cot⁡hx,yβ+ξβ=0,where ξ=2m⊥,B∗/ħ2[V0-E1(h;x,y)], β=2m⊥,W∗/ħ2E1(h;x,y), and m⊥,W∗ (m⊥,B∗) refers to the effective mass inside (outside) the dot region.The energyV0 is equal to the barrier potential for electrons in GaAs surrounded by a Ga1-xAlxAs material. V0=QΔEg, where Q is the band offset and ΔEg is the gap difference between the Ga1-xAlxAs and GaAs materials; here we take Q=0.6 and ΔEg=(1155x+370x2) meV and x=0.3 for the aluminum concentration.The eigenfunctions and eigenvalues of (2) can be obtained by several methods. In this work, particularly, we have used two different methods to confirm our findings: (i) expansion of the wave function in a set of complete wave functions associated with a confined electron in a square rectangular box (with dimensions Lx×Ly) with infinite confinement potential [18–21] and (ii) the finite elements method (FEM). In the first case, we have chosen to approach the solution of the in-plane problem in (2) by means of the following 2D expansion [18–21]: (6)ϕx,y=2LxLy∑n,mCn,msin⁡nπx+Lx/2Lxsin⁡mπy+Ly/2Ly.This means that we are writing the eigenfunctions of the in-plane problem as a linear combination of the independent solutions of a rectangular quantum box with infinite potential barriers. For a more realistic approximation, the dimensions Lx and Ly must be chosen appropriately larger than the typical in-plane ellipse-size. The calculation procedure then implies the construction of a Hamiltonian matrix from Hx,y and its subsequent diagonalization in order to obtain the eigenenergies Ej and the set of expansion coefficients Cn,m(j) that describe the allowed quantum states.In this study we have used a basis of sine functions in a regionLx=Ly=5a0 with quantum numbers n and m both running from 1 to 50, which gives us a total of 2500 functions. The resulting set of states is sorted out in increasing order of energy and after that we proceeded to choose the first N=200 states to diagonalize the Hamiltonian. With that choice we obtained a convergence of the lowest ten energy states up to 0.05R0 for the highest fifteen states (where R0=ħ2/2m∥∗a02 is the effective unit for energy, a0=ħ2ε/m∥∗e2 being the effective unit for length). Calculations were made also for Lx=Ly=10a0 with N=300 and for Lx=Ly=7a0 with N=700 and the obtained set of the lowest fifteen energies did not change more than 0.05R0. Also the condition Lx=Ly=10a0 with N=700 was studied with the same convergence.With the aim of solving the eigenvalues problem given by (2) by means of FEM, we introduce theSobolev SpaceH01(S1), which is formed by functions u, such that u and ∇u are defined and square-integrable functions on S1. We consider the variational form for (2) in S1 joined with a homogeneous Dirichlet boundary condition, which arises from considering infinity confinement outside S1: Find ϕ∈H01(S1) and E∈R such that, for any test function ψ∈H01(S1), (7) is satisfied. Numerical solutions of this problem were obtained using solvers implemented inFreeFem++ [22].(7)∫S1ħ22m∥∗∇ϕ·∇ψ+V1ρ+eF→·r→+E1h;x,yϕx,yψx,ydxdy=E∫S1ϕx,yψx,ydxdy.To use the FEM we have defined three contour lines: (i) S1-boundary, with 300 points, which corresponds to the xy-boundary region of the QD, (ii) an auxiliary circular boundary to define the maximum variation of the Coulomb interaction, centered at the impurity position with radius 0.2a0 and 30 points, and (iii) a second auxiliary circular boundary, also centered at the impurity position, enclosing an empty space with radius 10-4a0, to prevent the singularity associated with the Coulomb interaction. Note that these two auxiliary boundaries were chosen so as to converge to the 2D-hydrogenic binding energy (4R0) for large enough dimensions of the elliptical region.The knowledge of the wave functionsΨj(x,y,z) and their corresponding energy values is crucial to the study of the optical response related with interstate electron transitions in the QD. We have chosen to investigate the optical absorption coefficient, which is derived from the imaginary part of the dielectric susceptibility that characterizes the electric polarization of the electron system under the influence of the electric field of an electromagnetic signal of frequency ω [23, 24]. The working expression for the optical absorption coefficient is(8)αω=∑l>j=14ωμ0εϵ0Mjl2e2ρjlħΓjlEjl-ħω2+ħΓjl2,where the quantity Ejl=El-Ej is the energy difference associated with the transition between an initial state with energy Ej and a final state El. Mjl is the electric dipole moment off-diagonal matrix element involving those very states (in this work calculations are for circular xy-polarization of the incident radiation), whereas Γjl are the corresponding transition damping rates. The quantities μ0 and ϵ0 are the free space magnetic permeability and dielectric permittivity. Finally, ρjl=2(fj-fl)/V is the corresponding three-dimensional electron concentration [23], where(9)fr=11+expEr-EF/kBTis the Fermi-Dirac occupation number (EF is the Fermi level energy, V is the system’s volume, and T is the temperature).The input parameters for a prototypical GaAs QD are as follows:ε=12.35, c=3×108m/s, ϵ0=8.85×10-12F/m, μ0=1.256×10-6Tm/A, e=1.6×10-19 C, and me∗=0.067m0, where m0 is the free electron mass. In GaAs we deal with spherical-symmetric conduction effective mass; that is, m∥∗=m⊥∗=m∗. Calculations are for m∗=0.067m0 in the GaAs region and m∗=0.086m0 in the Ga0.7Al0.3As barriers when the finite confinement is considered. Here m0 is the free electron mass. For the length and energy effective units the considered values are, respectively, a0=97.98 Å and R0=9.95 meV. For the broadening we took the value Ωjl=ħΓjl=3meV [25, 26], in agreement with the corresponding value determined experimentally and supported by theoretical calculations [27] for a structure similar to ours. The broadening obviously depends on the temperature, dimensions, roughness, doping concentration, alloy disorder, and nonparabolicity of subbands, as well as the influence of elastic and inelastic scattering mechanisms. We did not take into account these dependencies in our calculations and, for that reason, we justify our choice for the broadening parameter due to the very good resolution obtained for the main structures in the absorption coefficient.


## 3. Results and Discussion

Figure2 shows the functional dependence of the energy of lowest confined states with regard to the change in the horizontal semilength for the case of the elliptical quantum dot with constant height depicted in Figure 1(a). The results are for zero applied field and correspond to two different values of the semilength Ry. The upper row corresponds to conduction band states [graphics (a) and (b)], and the lower row contains the energies of electrons coupled to the donor impurity located at the elliptical QD center. As can be expected, the increase of the QD size implies the reduction of the allowed confined energy values due to the loss in carrier localization. The inclusion of the Coulombic center (donor impurity) leads to an additional decrement of the energies as a result of the attractive electron-impurity interaction. Although the calculation was performed using the FEM in all cases, the spectral scheme that involves the diagonalization of the Hamiltonian matrix resulting from expansion (6) was also applied and the corresponding results appear superimposed to the FEM ones in Figure 2(a), showing a remarkably good coincidence. Note that in both cases of the numerical calculations the confinement potential around the elliptical region is infinite.Figure 2
(a, b) Energy of the first ten confined-electron states in an elliptical GaAs quantum dot as a function of the horizontal semilength (Rx) of the structure for two different values of the vertical length Ry: 1.0a0 (a) and 1.5a0 (b). In (c) and (d) the corresponding results including the effects of on-center donor impurity are depicted. The results are for zero electric field. A constant dot-height of 0.3a0 has been considered. In (a) the open symbols correspond to calculations using a diagonalization method. In (c) the full symbols correspond to theoretical calculations considering the 3D-character of the electron-impurity distance and using a variational procedure [15–17].
(a)
(b)
(c)
(d)In Figure2(c) we have also included the electron-impurity correlated energy for the ground state calculated via a variational technique and considering the 3D-shape of the structure (it means considering the 3D-character of the electron-impurity distance) [15–17]. Calculations were done with a one-variational parameter trial wave function. As can be seen, the FEM-results and the variational calculations agree with differences less than 2%, as was established above after (1).In Figure3 we are presenting the lowest electron-impurity energies, depending on the position of the donor impurity atom along the y-axis [cases (a) and (c)] and the x-axis [cases (b) and (d)] for the QD system of Figure 1(a). The static electric field is absent in cases (a) and (b). At the same time, a static electric field of 100kV/cm is assumed to be applied along the y-direction [case (c)] and along the x-direction [case (d)]. The zero-field configurations show a symmetric dependence of the coupled electron-impurity energies which, depending on the symmetries of the unconnected conduction band electron states, will show maxima and minima at different donor impurity positions. For all depicted level is coincident if the donor atom is placed at one of the QD borders the effect of the electron-impurity attraction is weaker and therefore the energy levels tend to go up (this is due to the an infinite barrier effect). It is worth noticing the presence of level anticrossings, typical of nonisotropic 2D or 3D systems.Figure 3
(a, b) Energy of the first ten confined-electron states in an elliptical GaAs quantum dot as a function of the impurity position along they-direction (a, c) and x-direction (b, d). The results are for Rx=1.5a0 and Ry=1.0a0 for zero applied electric field (a, b), F=(0,100,0) kV/cm (c), and F=(100,0,0) kV/cm (d). A constant dot-height of 0.3a0 has been considered.
(a)
(b)
(c)
(d)Including the influence of the static electric field modifies the above-mentioned picture. This time the electron density of probability becomes displaced inside the QD region. As a consequence, depending on the localization of the donor impurity, the energy values will be shifted upwards or downwards. This is a consequence of the rise or the fall of the expected electron-impurity distance, respectively, implying the weakening or the strengthening of the attractive Coulombic interaction. The field effect is more pronounced when it is applied along the horizontal direction of the elliptical cross section because in that case the field-induced electron delocalization is stronger. It can be seen that the electric field causes the appearance of level anticrossings in cases not appearing whenF=0. This is another consequence of the field-induced additional symmetry loss.These features can be confirmed by observing the evolution of the lowest electron confined energy levels as functions of the electric field strength shown in Figure4. Again, we are dealing with the constant height elliptic-shaped disc geometry and the dimensions Rx and Ry are fixed. Upper row [cases (a) and (b)] corresponds to the conduction band states in absence of impurity and the graphics in lower row [cases (c) and (d)] contain the results that include the donor center atom in the structure. The decrease of the energy level position as a result of the increment in the field intensity (greater loss in electron localization) is readily apparent in all cases. The above discussed situation where the presence of the attractive electron-impurity interaction implies the larger rate of decrement is also clearly present. Furthermore, the more noticeable effect of the x-oriented electric field is observed as well.Figure 4
(a, b) Energy of the first ten confined-electron states in an elliptical GaAs quantum dot as a function of the applied electric field along they-direction (a, c) and x-direction (b, d). The results are for Rx=1.5a0 and Ry=1.0a0 without impurity (a, b) and with on-center impurity (c, d). A constant dot-height of 0.3a0 has been considered.
(a)
(b)
(c)
(d)The results for the variation of the energy levels as functions ofRx in the case of the variable height QD [see Figure 1(b)] is shown in Figure 5. There, the zero-field situation is explored for two different geometrical configurations. In one of them, the QD height is assumed to vary in the form h=Rz1-x2/Rx2-y2/Ry2 [cases (a) and (c)]. In the second one, the height changes according to h=h0+Rz1-x2/Rx2-y2/Ry2 [cases (b) and (d)]. This means that close to the QD bottom there is a strip of constant height, resembling the presence of a possible wetting layer. In addition, the presence of the donor impurity center is included in the results appearing in Figures 5(c) and 5(d).Figure 5
(a, b) The results are as in Figure2, but here we are interested in an quantum dot with variable height. In (a, c) the height of the dot is given by h=Rz1-x2/Rx2-y2/Ry2, whereas in (b, d) the results are for h=h0+Rz1-x2/Rx2-y2/Ry2. Calculations are for Ry=1.0a0, Rz=0.3a0, and h0=0.1a0.
(a)
(b)
(c)
(d)The reduction in the energy values as a consequence of the increase in the horizontal dot size, as discussed above in the constant height case, is observed as well. It is interesting to note that the inclusion of a nonzero-height strip at the base of the QD that amounts33% of increase in the top height causes a significant redshift of the confined energies. For instance, in the case of the ground state, the energy decreases approximately 40% with respect to the unstriped geometry, both with and without the presence of the impurity atoms. The presence of anticrossings is preserved.In Section2 we have commented the possibility of the numerical calculation of the confined states in the QD which has the elliptical-shaped configuration with the finite confinement potential via FEM. Accordingly, Figure 6 contains the zero electric field results for the variation of the confined states in the case of nonconstant height geometry [analogous to that appearing in Figure 1(b)] as functions of Rx. Graphics (a) and (b) correspond to conduction band electron states without impurity effect, whereas graphics (c) and (d) show the allowed electron-donor impurity energies. The configurations chosen for the finite potential energy confinement profiles are as follows: (i) With QD height given by h=Rz1-x2/Rx2-y2/Ry2 [cases (a) and (c)], the quantum well along the z-direction is taken as symmetric finite barrier one with depth provided by the corresponding band offset of the Al0.3Ga0.7As/GaAs system. (ii) With QD height given by h=h0+Rz1-x2/Rx2-y2/Ry2, with h0=0.25a0 [cases (b) and (d)], the quantum well along the z-direction has a bottom finite barrier and an infinite barrier at the top that could simulate an upper free-standing configuration.Figure 6
(a, b) The results are as in Figure2, but here we are dealing with a variable height of the dot. In (a, c) the height of the dot is given by h=Rz1-x2/Rx2-y2/Ry2 with symmetric finite confinement potentials along the z-direction, whereas in (b, d) the results are for h=h0+Rz1-x2/Rx2-y2/Ry2 with asymmetric finite-infinite-confinement potentials along the z-direction. Calculations are for Ry=1.0a0, Rz=0.3a0, and h0=0.25a0.
(a)
(b)
(c)
(d)In first place one may readily noticed that, as it should be expected, the values of the obtained energies are significantly smaller compared with the previously discussed infinite potential barrier cases. However, the behavior of the energies as functions of the horizontal semilength are qualitatively similar, with the same arguments related with the electron spatial localization and the expected electron-impurity distance as the key elements for their explanation. It should be kept in mind that the confining configuration (ii) implies a larger size of the width of the quantum well along the vertical direction. As a consequence, in spite of the presence of a top infinite barrier, the obtained energies are lower in both the uncoupled and Coulombic-coupled situations.In order to proceed with the investigation of the optical absorption coefficient in the elliptic-shaped QDs we need to evaluate the off-diagonal electric dipole matrix elements associated with the interstate transitions that will be considered. As indicated in (8), only the contributions coming from the transitions that involve the lowest four confined states shall be included. Following the same geometrical, compositional, and interaction setups that led to the energy level results in Figures 2–6, the corresponding variations of the Mij are displayed in Figures 7–11, respectively. Note that all the results presented in Figures 7–11 are for the low temperature regime (4 K).Figure 7
Dipole matrix elements for circular polarization of the incident radiation considering combinations between the first four confined states. The results are obtained under the same geometrical, compositional, and interaction configurations leading to Figure2.
(a)
(b)
(c)
(d)Figure 8
Dipole matrix elements for circular polarization of the incident radiation considering combinations between the first four confined states. The results are obtained under the same geometrical, compositional, and interaction configurations leading to Figure3.
(a)
(b)
(c)
(d)Figure 9
Dipole matrix elements for circular polarization of the incident radiation considering combinations between the first four confined states. The results are obtained under the same geometrical, compositional, and interaction configurations leading to Figure4.
(a)
(b)
(c)
(d)Figure 10
Dipole matrix elements for circular polarization of the incident radiation considering combinations between the first four confined states. The results are obtained under the same geometrical, compositional, and interaction configurations leading to Figure5.
(a)
(b)
(c)
(d)Figure 11
Dipole matrix elements for circular polarization of the incident radiation considering combinations between the first four confined states. The results are obtained under the same geometrical, compositional, and interaction configurations leading to Figure6.
(a)
(b)
(c)
(d)The different, and sometimes jumbled, behaviors of these quantities are governed by two main elements: the spatial extension and the symmetry of the involved wave functions. They, combined with the polarization of the incident light (here considered to be circular), determine whether the dipole matrix element vanishes or remains finite as well as its functional dependence with respect to the varying quantity in the system. For instance, the effect of the anticrossings, which imply the sudden change in the symmetry of certain wave function and, therefore, the imposition of a particular selection rule for the transition under specific circumstances, can be clearly noticed.The total light absorption coefficient that includes the contribution of the above transitions commented on is shown in Figure12 as a function of the incident photon energy for two different values of the temperature T=4K [(a), (c)] and T=300 [(b), (d)]. The geometry is that with constant height, depicted in Figure 1(a), with fixed dimensions. The electric field, oriented along the y-direction, intensity is varied from 8kV/cm to 200kV/cm. Graphics (a) and (b) correspond to the case without the presence of the impurity atom, and in the graphics (c) and (d) the system with on-center donor impurity atom is taken into account. In the case without impurity we can identify the resonant peaks associated with transitions having the ground state as the initial one. The shift of the peaks as a result of the increment in F can be identified by observing the behavior of the involved state energies as functions of the electric field strength in Figure 4. The evolution of the resonant peak amplitudes can also be noticed which follow the dependence generated by the combined effects of the variation of the resonant peak position and the electric dipole moment matrix element through the product Eij|Mij|2/ħ.Figure 12
Total absorption coefficient as a function of the incident photon energy for several values of the electric field applied in they-direction and varying in steps of 8 kV/cm from 0 to 200 kV/cm. Rx=1.5a0, Ry=1.0a0, and constant dot-height of 0.3a0. The curves are shifted 0.8×107 m-1. In (a) and (b) the results are without donor impurity whereas in (c) and (d) they are for on-center impurity. Two values of temperature have been considered as indicated in the figures. The green lines correspond to the calculated energy transitions from Figure 4(a).
(a)
T
=
4
K
(b)
T
=
300
K
(c)
T
=
4
K
(d)
T
=
300
KIn our calculation, the effect of the temperature is supposed to fall on the population of the transition states. We are assuming that the influence ofT on the transition energy differences and state wave functions is not significant and can be neglected. Note that in the case of finite quantum wells in the z-direction, the temperature coefficient of the gap for the well and barrier materials is the same and consequently there is no change with temperature of the height of the confinement potential. In the case of the effects of the temperature on the electron effective mass and static dielectric constant, they are global and involve almost rigid modifications on the whole spectra without changes in the transition energies. Accordingly, the results depicted in Figures 12(b) and 12(d) reveal the consequence of the change in the level populations, with the appearance of absorption peaks for T=300K, additional to those observed at low temperature. At the same time, the variation in ρij conditions the change in the resonant peak amplitudes as can be noticed.Finally, Figure13 shows the calculated total optical absorption coefficient in the symmetric finite confinement configuration of the elliptic-shaped QD with nonconstant height. The light absorption coefficient is presented as a function of the incident photon energy and curves are drawn taking the horizontal elliptic size Rx as a parameter, with Ry and Rz fixed. Again, two values of the temperature are taken into account. Graphics of the results without the presence of any impurity are labeled (a) and (b), whilst (c) and (d) correspond to the case with an on-center ionized donor atom.Figure 13
Total absorption coefficient as a function of the incident photon energy for several values of the horizontal length varying in steps of0.04a0 from 0.5a0 to 1.5a0. The curves are shifted 3.2×107 m-1. The height of the dot is given by h=Rz1-x2/Rx2-y2/Ry2 with symmetric finite confinement potentials along the z-direction. Calculations are for Ry=1.0a0 and Rz=0.3a0. In (a) and (b) the results are without donor impurity whereas in (c) and (d) they are for on-center impurity. Two values of temperature have been considered as indicated in the figures.
(a)
T
=
4
K
(b)
T
=
300
K
(c)
T
=
4
K
(d)
T
=
300
KThe shifts of the resonant absorption peak positions can be explained as above, that is, by observing the behavior of the confined state energies appearing in Figures6(a) and 6(c) and the variation of the differences between them. The change in their amplitudes is, once again, justified by the behavior of the involved transition energies and the electric dipole matrix elements shown in Figures 11(a) and 11(c). Also, as we previously commented, the influence of a higher temperature, affecting the electron concentration ρ, allows for the resolution of other transition peaks, highlighting the thermal effect in this case.


## 4. Conclusions

In the present work we have addressed the calculation of the conduction band effective mass states of elliptically shaped quantum dots with finite and infinite confinement potentials, the presence of a donor impurity atom, and the influence of externally applied static electric fields. The solution of the Schrödinger-like effective mass differential equation incorporates the adiabatic approximation to uncouple the electron motion along thez and (x,y) directions, as well as the finite element method to solve the resulting differential problem. The use of the finite element method usingFreeFem++ is an important tool for solving this kind of problem. Results also were obtained via a spectral procedure which reveals a very good agreement with the finite element scheme.Our results confirm the known effects of the presence of the electron-impurity interaction and the application of static electric fields on the spectrum of carrier confined in quantum nanostructures, indeed with the particularities associated with the specific geometry of the system under study.The light absorption associated with electron transitions between the allowed quantum states in each case is studied making use of the previously calculated energies and wave functions. It is shown that the changes in the electron state energies and probability densities due to modifications in the type of confinement, the geometry, and the presence or absence of the donor impurity center, as well as the variation of the level population with temperature, are all causes for the shift of the absorption resonant peaks and/or the increment or reduction of their corresponding amplitudes.Quantum dots with the shape here discussed are, actually, experimentally realized systems, with well-identified current and prospective applications. We hope, with this work, to shed some light on the electronic and optical features of a GaAs-based structure of this type.



---

*Source: 5970540-2017-04-05.xml*