Abstract
We study the excitonic effects on the secondorder nonlinear optical properties of semispherical quantum dots considering, on the same footing, the confinement potential of the electronhole pair and the Coulomb interaction between them. The exciton is confined in a semispherical geometry by means of a threedimensional semiparabolic potential. We calculate the optical rectification and second harmonic generation coefficients for two different values of the confinement frequency based on the numerically computed energies and wavefunctions of the exciton. We present the results as a function of the incident photon energy for GaAs/AlGaAs quantum dots ranging from few nanometers to tens of nanometers. We find that the secondorder nonlinear coefficients exhibit not only a blueshift of the order of meV but also a change of intensity compared with the results obtained ignoring the Coulomb interaction in the socalled strongconfinement limit.
Introduction
Nonlinear optical properties of semiconductor quantum dots have attracted considerable interest due to their several potential applications [14]. In particular, secondorder nonlinear optical properties, such as nonlinear optical rectification (OR) and second harmonic generation (SHG), have received special theoretical [58] and experimental [9,10] attention due to their magnitudes being stronger than those of highorder ones, making them the first nonlinear optical effects experimentally observed.
The confinement of carriers provided by a quantum dot is well described by a parabolic potential when only the lowest excited states of the carriers are considered. However, selfassembled quantum dots growth in the laboratory usually exhibit asymmetrical shapes that ensure the generation of nonlinear optical effects. In order to model these asymmetries, an asymmetrical potential is required.
Recently, several authors [5,6,8] studied the effects of an exciton on the secondorder nonlinear properties in onedimensional semiparabolic quantum dots. Using analytical approximate results, they showed that the excitonic effects enhance significantly the OR and SHG coefficients. They used the socalled strongconfinement limit, ignoring in this way the Coulomb interaction between electron and hole because of the quantum dot dimensions are smaller than the effective Bohr radius, and finding that the excitonic effect reduces itself to an effectivemass model in one particle scheme.
In this study we find eigenenergies and eigenstates of an exciton in a semispherical quantum dot solving the corresponding threedimensional Schrödinger equation using a finite elements method and taking into account both the confinement and Coulomb potentials of the electronhole pair. We present the OR and SHG coefficients as a function of the incident photon energy with and without Coulomb potential. Our results show that energy and intensity of the peaks in the secondorder nonlinear optical coefficients change when Coulomb interaction is introduced.
This article is organized as follows. In "Theory" section, we present the characteristic quantities of the harmonic and Coulomb potentials, and the definitions of the weak and strongconfinement limits in terms of these parameters. In addition, we present the analytical expressions for the optical nonlinearities, such as OR and SHG, obtained by the density matrix formalism. In "Results" section, we show the OR and SHG coefficients with and without Coulomb interaction as a function of the incident photon energy for two quantum dot sizes. We also give account of the changes presented by the secondorder nonlinear coefficients focusing in the role played by the Coulomb interaction. Conclusions are summarized in final section.
Theory
The effectivemass Hamiltonian for the electronhole pair in the threedimensional quantum dot reads [11]
where and are the effective masses of the electron and the hole, respectively, ε is the background dielectric constant and V (r_{i}) is the threedimensional semiparabolic potential that we define as
The angle θ is the usual polar angle in spherical coordinates, and ω_{0 }the oscillator frequency considered in this study the same for the electron and the hole. The potential defined in Equation 2 confines the exciton in the upper half of a sphere, i.e., the quantum dot has a semispherical shape.
Hamiltonian (1) can be separated in terms of centerofmass and relative coordinates, respectively,
where is the total mass, and is the reduced mass. The centerofmass and relative position coordinates are defined as usual,
with the corresponding momenta P =  iħ∇_{R }and p =  iħ∇_{r }in terms of p_{e }and p_{h},
The explicit separability of the centerofmass and relative coordinates in Equation 3 lead to the following total envelope wave function and total energy for the system:
The centerofmass part of Hamiltonian (3) is a threedimensional semiparabolic oscillator that can be solved analytically. Therefore, the problem is now reduced to solve the relative motion Hamiltonian:
Hamiltonian (8) has been solved analytically in two limiting cases (strong and weak confinement) for onedimensional quantum dots. The eigenfunctions and eigenvalues are presented in references [5] and [8]. In onedimensional case, the confinement potential also imposes constraints to spatial coordinates, resulting in a hydrogenlike (asymmetricharmonic) reduced particle Hamiltonian for weak (strong) limit.
The harmonic potential in Equation 8 defines both the size L of the quantum dot,
and the energy quanta ħω_{0 }due to confinement, which is related to L by
On the other hand, the Coulomb potential defines the effective Bohr radius and the effective Rydberg energy of the electronhole interaction,
The strongconfinement limit is established when , or equivalently , and the weakconfinement limit when , or .
The secondorder nonlinear optical coefficients can be obtained by density matrix approach and perturbation expansion method [12,13]. The expression for the OR coefficient, within a twolevel system approach, is given by [5,6]
where e is the electron charge, σ_{s }is the density of electrons in the quantum dot, T_{1 }is the longitudinal relaxation time, T_{2 }is the transverse relaxation time, and
For the resonance condition
there is a peak intensity given by (1/T_{1 }≪ 1/T_{2 }≪ ω_{01})
The SHG coefficient in a three level system is [8]
where N is the density of carriers in the quantum dot, E_{ij }= E_{i } E_{j}, Γ_{10 }= Γ_{20 }= Γ_{0 }are the relaxtion rates, and
Under the double resonance condition, i.e., ħω ≈ E_{10 }≈ E_{20}/2, the intensity of the peak is given by
and its energy by
Results
In this study, the results are presented for a GaAs/AlGaAs structure. We have used the following parameters in the calculations: = 0.067m_{0}, = 0.09m_{0 }(m_{0 }is the mass of a free electron) [4], T_{1 }= 1 ps, T_{2 }= 0.2 ps [12], σ_{s }= 5 × 10^{24 }m ^{3 }[5], ε = 12.53, Γ_{0 }= 1/0.14ps ^{1}, N = 3 × 10^{16 }cm^{3 }[8].
In Figure 1, we plot the characteristic lengths and energies for the confined particle in a GaAs/AlGaAs quantum dot as a function of the confinement frequency ω_{0}. and are independent on ω_{0 }because they are related to the Coulomb potential. In Figure 1a, we can see that the lengths L and are of the same order of magnitude for a confinement frequency around ω_{0 }= 1 × 10^{13 }s^{1}. In Figure 1b, we observe that also ħω_{0 }and show similar values around ω_{0 }= 1 × 10^{13 }s^{1}. For this reason, we conclude that, in this frequency range, neither the strongconfinement limit nor the weak limit can be assumed because both interactions, harmonic and Coulomb, are important. Therefore, we propose a numerical technique to calculate eigenenergies and eigenstates of Hamiltonian (8), considering the harmonic and Coulomb potentials.
Figure 1. Characteristic (a) lengths and (b) energies for the confined particle in a GaAs/AlGaAs quantum dot as a function of the confinement frequency. The red (black) lines correspond to L and ħω_{0 }( and ), respectively.
With the aim of exploring the nonlinear behavior at higher frequencies, i.e., when the quantum dot size is smaller than in the previous case, we choose ω_{0 }= 2 × 10^{14 }s^{1}, in which the quantity L is less than , or ħω_{0 }is greater than , differing in both cases by one order of magnitude as can be seen in Figure 1. Because of this difference, several authors[5,6,8] used the strongconfinement limit as a satisfactory approximation in the case of small quantum dots. Accordingly with Equation 9, the frequencies ω_{0 }= 1 × 10^{13 }s^{1 }and ω_{0 }= 2 × 10^{14 }s^{1 }define a quantum dot size of L = 17.4 nm and L = 3.9 nm, respectively. This means that our results are suitable for the current quantum dot sizes that range from few nanometers to tens of nanometers.
We find numerically eigenenergies and eigenstates of Hamiltonian (8) by a finite elements method for the two frequencies mentioned above. We have used the software COMSOL Multiphysics, which offers the possibility of defining a geometry, in this case the upper half of a sphere, and to solve the timeindependent Schrödinger equation with appropriate boundary conditions.
The terms involving quantum states and energies in Equations 13 and 19 are found using the eigenstates and eigenenergies previously calculated. The OR and SHG coefficients are shown in Figures 2 and 3, respectively. Figures 2a and 3a correspond to ω_{0 }= 1 × 10^{13 }s^{1}, and Figures 2b and 3b to ω_{0 }= 2 × 10^{14 }s^{1}. In each figure, we present the corresponding nonlinear optical coefficient considering excitonic effects with and without Coulomb interaction. For comparative purposes, we also present the case without excitonic effects, i.e., when only one electron exists in the quantum dot.
Figure 2. The OR coefficient as a functions of the incident photon energy ħω for (a) ω_{0 }= 1 × 10^{13 }s^{1 }and (b) ω_{0 }= 2 × 10^{14 }s^{1}, considering excitonic effects with (red line) and without Coulomb (black line)interaction. The blue line corresponds to the case without excitonic effects.
Figure 3. The SHG coefficient as a functions of the incident photon energy ħω for (a) ω_{0 }= 1 × 10^{13 }s^{1 }and (b) ω_{0 }= 2 × 10^{14 }s^{1}, considering excitonic effects with (red line) and without Coulomb (black line)interaction. The blue line corresponds to the case without excitonic effects.
Figures 2 and 3 reproduce the reported results quite well [5,8] regarding the enhancement of the nonlinear optical coefficients due to the quantum confinement of the exciton. However, the same figures show additionally a significant blueshift of the OR and SHG peaks when both the harmonic and Coulomb potentials are taken into account. In Tables 1 and 2, we present the eigenenergies of the exciton and peak energies of the OR and SHG coefficients with and without Coulomb interaction for the two frequencies under study. The peak energies are estimated by Equations 17 and 22 for the OR and SHG coefficients, respectively.
Table 1. Eigenenergies of the exciton and peak energies of the OR and SHG coefficients with and without Coulomb interaction for ω_{0 }= 1 × 10^{13 }s^{1}
Table 2. Eigenenergies of the exciton and peak energies of the OR and SHG coefficients with and without Coulomb interaction for ω_{0 }= 2 × 10^{14 }s^{1}
We can see from Tables 1 and 2 that the eigenenergies obtained with Coulomb interaction are smaller than those obtained without that interaction. The explanation to this fact is that there is an attractive Coulomb potential between the electronhole pair that implies a reduction of the eigenenergies for the exciton. However, the eigenenergies are affected in different ways depending on the quantum state. For example, for the ground state ω_{0 }= 1 × 10^{13 }s^{1}, Table 1, we have an energy difference of 5.237 meV, while for the first and second excited states the differences are of 3.472 and 2.784 meV, respectively. We have a similar situation for ω_{0 }= 2 × 10^{14 }s^{1}, Table 2. This is because the mean spatial separation between the electron and the hole increases, and therefore the Coulomb interaction decreases, as the energy of the quantum state increases. The final result is a blueshift of the OR and SHG peaks of the order of meV for both ω_{0 }= 1 × 10^{13 }s^{1 }and ω_{0 }= 2 × 10^{14 }s^{1}.
In addition, the OR and SHG coefficients exhibit different peak intensities depending on the consideration of the Coulomb interaction, as it can be seen in Figures 2 and 3. This fact originates from the modification of the dipole matrix elements defined in Equations14, 15, and 20 when the Coulomb interaction is considered. According to Equation 18, the peak intensity of OR coefficient depends essentially on the product , while for SHG coefficient, Equation 21, the peak intensity depends on μ_{01 }μ_{12 }μ_{20}. Tables 3 and 4 show the values of these dipole matrix element products with and without Coulomb interaction for the two frequencies considered. As one can see, the product is greater with Coulomb interaction than without it for both confinement frequencies. Therefore, in Figure 2a, b, the OR intensity is higher in the former case than in the later one. In the case of SHG coefficient, the product μ_{01 }μ_{12 }μ_{20 }is smaller with Coulomb than without that interaction. This fact makes the SHG intensity smaller in the former case as can be seen in Figure 3a, b.
Conclusions
Contrary to the assumption that Coulomb interaction can be neglected when the quantum dot dimensions are smaller than the effective Bohr radius, we show that this interaction affects the excitonic effects of the secondorder nonlinear optical properties of semispherical quantum dots. We find that Coulomb interaction manifests itself in a blueshift of the energy peaks of the order of several meV in the studied spectra. These results were found for two quantum dot sizes, in the first one the characteristic quantities of the harmonic and Coulomb potentials are equals, and in the second one they differ by one order of magnitude. This means that the Coulomb interaction plays an important role even when the quantum dot sizes are smaller than the effective Bohr radius.
Therefore, we encourage experimentalists to carry out measurements of secondorder optical nonlinearities in asymmetrical quantum dots with the aim of to detect the magnitude of this effect.
Abbreviations
OR: optical rectification; SHG: second harmonic generation.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
JF: carried out the numerical calculations and drafted the manuscript. AC: performed ananalys and interpretation of results, and gave final approval of the version to be published. All authors read and approved the final manuscript.
Acknowledgements
This work was partially supported by Facultad de Ciencias of Universidad de los Andes.
References

Fu Y, Hellström S, Ågren H: Nonlinear optical properties of quantum dots: excitons in nanostructures.
J Nonlinear Opt Phys Mater 2009, 18:195. Publisher Full Text

Gotoh H, Kamada H: Coherent nonlinear optical properties in quantum dots.

Bimberg D, Meuer C, Lämmlin M, Liebich S, Kim J, Kovsh A, Krestnikov I, Eisenstein G: Nonlinear properties of quantum dot semiconductor optical amplifiers at 1.3 μm.
Chin Opt Lett 2008, 6:724. Publisher Full Text

Rosencher E, Vinter B: Optoelectronics. Cambridge: Cambridge University Press; 2003.

Yu YB, Zhu SN, Guo KX: Exciton effects on the nonlinear optical rectification in onedimesional quantum dots.
Phys Lett A 2005, 335:175. Publisher Full Text

Baskoutas S, Paspalakis E, Terzis AF: Effects of excitons in nonlinear optical rectification in semiparabolic quantum dots.
Phys Rev B 2006, 74:153306. Publisher Full Text

Zhang CJ, Guo KX, Lu ZE: Exciton effects on the optical absorptions in onedimensional quantum dots.
Phys E 2007, 36:92. Publisher Full Text

Karabulut İ, Şafak H, Tomak M: Excitonic effects on the nonlinear optical properties of small quantum dots.
J Phys D: Appl Phys 2008, 41:155104. Publisher Full Text

Brunhes T, Boucaud P, Sauvage S, Lemaître A, Gérard JM, Glotin F, Prazeres R, Ortega JM: Infrared secondorder optical susceptibility in InAs/GaAs selfassembled quantum dots.
Phys Rev B 2000, 61:5662. Publisher Full Text

Sauvage S, Boucaud P, Brunhes T, Glotin F, Prazeres R, Ortega JM, Gérard JM: Secondharmonic generation resonant with sp transition in InAs/GaAs selfassembled quantum dots.
Phys Rev B 2001, 63:113312. Publisher Full Text

Que W: Excitons in quantum dots with parabolic confinement.
Phys Rev B 1992, 45:11036. Publisher Full Text

Rosencher E, Bois P: Model System for optical nonlinearities: Asymmetric quantum wells.
Phys Rev B 1991, 44:11315. Publisher Full Text

Boyd R: Nonlinear Optics. 3rd edition. New York: Elsevier; 2008.