Abstract

In this paper, we have performed a theoretical study on nonlinear optical rectification (OR) and second harmonic generation (SHG) for three-level dome-shaped InAs/GaAs quantum dots (QDs) in the presence of wetting layer (WL). We used the compact density matrix framework and effective mass approximation to investigate the second order nonlinear phenomena on InAs/GaAs QD. It is demonstrated that second harmonic generation (SHG), optical rectification (OR), and their mutual absorption and refractive index changes are quite sensitive to the size of QDs. The size variations have profound irregular behavior owing to distribution of envelope function on WL and QD simultaneously. Moreover it is found that  nm is a critical radius where the regular variation takes place. It is shown that size variation causes blue shift until Critical radius ( nm) and after that, increasing the QD size lead to redshift in second order phenomena.

1. Introduction

In the last few years the nonlinear optical properties of intersubband transitions in semiconducting materials have attracted remarkable attention due to the potential application in electronics and optoelectronics devices [113]. Among the semiconductor materials, zero-dimensional quantum dots (QDs) are very interesting; tunable features of QDs are one of the main reasons for this increasing performance of semiconductor QDs. The rapid advances in nanotechnology techniques such as molecular beam epitaxy (MBE) [14, 15], metal-organic chemical vapor deposition [16], and Stranski-Krastanov (S-K) [17] methods have made it possible to prepare different shape and geometry of QDs [1418]. In QDs the carriers are confined by the confinement potential, named as quantum confinement effect in all spatial directions. Confinement provides the quantization of electronic energy levels based on the size of the dots; confinement leads to significant optical nonlinearities [11, 18]. Photons with proper energy can cause the intersublevel transitions involving large electric dipole moments [19, 20]. Among the nonlinear optical properties, the second-order nonlinear optical property plays a crucial role because it is the simplest and the lowest-order nonlinear effect; moreover, its magnitude is usually stronger than the other optical nonlinearities [21]. The second order nonlinear optical interaction of two incident fields with optical media leads to some phenomena’s like Second harmonic generation (SHG) and optical rectification (OR). To the best of our knowledge the early work on OR is dated back to 1962 by Xie and Bass et al. [21, 22] as well as the very first experimental observation of SHG reported by Franken et al. [23]. It is worthwhile to mention that the beginning of the field of nonlinear optics is taken by discovery of SHG [24]. SHG can be so efficient since nearly all incident beams at frequencytransformed to radiation at second harmonic frequency [24]. The OR inference on semiconductor appeared in literature of 1992 under the rotating-wave and effective-mass approximation by Chuang et al. [25]. After this demonstration this sequence keeps going forward for different types of semiconductor nanostructures like quantum wells and quantum dots [26, 27]. The OR has been studied in several QDs systems including the cubical, disc-like, and spherical shapes [2830]. The effect of exciton in OR for semiparabolic QDs for confined electron and hole in semiparabolic potential was studied by Baskoutas et al. [31]. Other studies show that the OR is a function of the applied electric field and confinement strength on QDs [30, 32]. In QDs the spin-orbit interaction decreases the magnitude of OR as well as the peaks which happen at lower energies but the donor impurity causes a big shift in the magnitude of the OR peaks position as reported by Chen et al. and Vaseghi et al. respectively [33, 34]. Recently it is found that the OR in lense-shaped QDs is quite sensitive to the applied hydrostatic pressure and temperature [35]. The Polaron effect in OR and SHG in cylindrical QDs was studied by Liu et al. [36]; indeed they made a focus on the dependence of both the OR and the SHG on confinement frequency, external magnetic field, and the linear potential coefficient .

The early experimental research on mid-infrared and near-field SHG in semiconductor QDs was reported by Brunhes et al. [37], Liu and Bryant [38], respectively. Indeed Brunhesa et al. [39, 40] showed that the resonant SHG in association with the intraband transitions of the P-type QDs grown by MBE can be achieved. Also the enhanced SHG with intersubband transition (specifically S-P transition) with conduction states of InAs/GaAs was demonstrated experimentally for the first time by Brunhesa et al. [41]. The size of the QDs influences the confinement potential and thus optical properties of low dimensional QDs face some alteration. The effect of size under different values of electric and magnetic fields on SHG for parabolic studied by Li et al. [42]; also the length-dependent SHG on cubical QDs as function of incident photon energy was presented by Shao et al. [43]. The spin-orbit interaction effects on the OR for different size of QDs were reported by Vaseghi et al. [34]. Sabeian and Khaledi-Nasab highlighted the size-dependent intersubband optical properties of dome-shaped InAs/GaAs quantum dots with wetting layer [44].

The QDs in the practical growth method always are grown on wetting layer (WL) [45, 46]. In our study the InAs is sandwiched between two GaAs matrixes. In almost all the previous studies the second order susceptibility equations were simplified and some of the none resonant terms were removed based on some approximations. In this work we are presenting a deep understanding of the OR and SHG in a three-level dome-shaped InAs/GaAs QDs with WL using finite element method (FEM) and compact density matrix framework. The WL zone which is usually ignored in some of the preceding articles [4763] is fully adopted as well as its effect on second order nonlinear phenomena is investigated. By taking all resonant and nonresonant terms of second order nonlinear susceptibility under consideration, the effect of all expressions is being highlighted.

2. Theoretical Framework

This work is concerned with the investigation of dome-shaped InAs QDs. First of all we calculated the transition frequencies between the levels and envelop function in single InAs QDs by making use of finite element method (FEM) numerical solution for Schrödinger equation within the effective mass approximation. Throughout this work the InAs QDs are grown on InAs WL which is embedded in a large GaAs matrix and the GaAs cap cover the grown QDs as Figure 1(a) depicts.

By this definition we need to consider a dome-shaped InAs/GaAs QD with cylindrical symmetry whose Schrödinger wave equation in one-band envelop-function formalism is given by whereis the electron effective mass having two different values in QD region and in surrounding medium, is the electronic envelop function, and is the eigenenergy. Using the cylindrical symmetry of the presented geometry in Figure 1(a) gave the permission to separate the envelop function as where, , andare cylindrical coordinates. Substituting (2) in (1) and dividing both sides bygive whereis a separation constant. As the envelop function must be single-valued underrotation, the-part can be expressed asleading to integer values of. Equation (4) can be rearranged as As it is obvious the three-dimensional Schrödinger equation is being simplified to two variable equation. To solve the Schrödinger equation we must define adequate boundary conditions. According to Figure 1(b), for the boundaries 1 (rooftop) and 2 (bottom), the condition ofis considered, while for all other exterior boundaries the condition of is employed where is outward unit vector. For interface boundaries of InAs QD-WL and GaAs matrixes, the condition ofis adopted because of the finiteness of the potential barrier. In our numerical simulations, the effective electron mass was set as 0.023for InAs and 0.069for GaAs [64], whereis the free electron mass. The finite potential was taken as 0 and 0.697 eV for InAs and GaAs, respectively.

Assume that our system is excited by linear-polarized monochromatic electric field propagating along-direction as According to Boyd notation [24], . Based on nonlinear optics formalism [24] the refractive index changes and the absorption coefficient are given by: where,is the vacuum permeability, is the real part of permittivity, and is the vacuum permittivity. Consider is effective susceptibility expressed in term of linear and nonlinear parts. In this study we are looking only for second order nonlinear susceptibility, so the absorption coefficient and refractive index changes equations will rearranged be as: Based on the density matrix approach [24], the second order OR and SHG relationship could be calculated. According to this formalism, the OR and SHG for 3 levels of InAs/GaAs QD system are given by whereis the density of the carriers,is the component of electric dipole moment matrix,is the transition angular frequency, is the relaxation rate from levelto level , and andare the dummy indices, indicating the number of the levels which is 3 here.

In this case OR is purely real, so it has no contribution in absorption coefficient but SHG collaborates in both refractive index changes and absorption coefficient. It is noteworthy to indicate that we avoid opening the summation because both OR and SHG consist of almost 76 and 144 successive expression, respectively. We admit removing any nonresonant expressions since the effect of all expression will leads to more practical results.

3. Result and Discussion

Figure 1(b) depicts the geometry of models set as 30 nm and 100 nm for the width and height of the 2D modeling area, respectively. The height of WL has been set as 3 nm based on what presented in [44, 65, 66]. In the next step the radius of QDs was varied from 1 nm to 30 nm. To carry out our simulations, we set the carrier density toand the damping rate as all the relaxation rates[6769]. Also, the refractive index of GaAs was set as constant value .

The variation of the first three transition frequency amplitudes in the sublevels of conduction band as a function of QDs size is shown in Figure 2. The first three subbands are named s, p, and d and thus these intersubband transitions belong to these specific subbands. The diagonal and off-diagonal elements dipole momentum matrixes are illustrated in Figures 3(a) and 3(b), respectively. According to Figures 2 and 3 approximately until nm the parameters take a single value with slight ignorable fluctuation. At  nm the amplitude rises up. After nm the situation becomes different in all the specific levels. The envelop function of ground state (GS, ) consists of single part but the envelop function on first excited state (FES, ) and second excited state (SES, ) consists of 2 and 3 individual units (lobes), respectively. When the radii of the QDs are below 5 nm, the envelop function for all states totally holds on WL, in other words the envelop functions are localized inside the WL. At nm GS envelop function and first lobes of FES and SES begin to leave the WL and enter the QD. Figure 2 reveals the dropdown for all frequencies from 5 to 8 nm. For upper dimensions, more complicated behavior is being observed.and() are related to GS and FES.also has a constant value untilnm. When the entry of second lobe of FES begins atnm, simultaneouslyslope changes from negative to positive, in the same mannershows a jump on its amplitude but it still takes close values. The second lobe gets inside the dot atnm. After that dimension, thesuddenly takes negative slope and for upper radii reduces with semiexponential manner. For radii larger than 15 nm, the GS and FES become more localized inside the dot; also their its slope vanishes gradually. Moreover thetakes linear reduction but thelinear behavior increase which makes one think that thecomponent partially controls these elements.and() are multifaceted parameters since the SES has three lobes. The SES has 3 lobes which leave the WL successively.generally has higher absolute value thanfor every specific dimension. Foruntil nm everything is just likebut after that when the third lobe of SES begins to enter the dot and Elevation occurs between nm and nm. Afterward the expected regular sequence take place.afteruntilnm increases but a ripple happens because entrance of first lobe in both GS and SES ended and the Second lobe initiate its journey toward the dot. The slope is positive until nm. However the afterward the slop changes to negative. Since  nm has slight changes on the negative sequence this leads the entire behavior to have the same slope for upper radii.

and() are calculated based on FES and SES. Thus we anticipate more complicated behavior. Indeed has smaller value than other transitions. Until nm everything is like previous transitions. At nm a little elevation occurs but soon after that the sequence keep moving down. The second lobes of FES and SES collaborate with each other cause this reduction until nm; this is due to the simultaneous movement of these lobes to the dot which intensify the decline and its slope. Afterward a promotion follows,  nm and in just one nanometer greater than the regular sequence shows up.has slight fluctuation until nm which is not so weird because the lobes are not so synchronized. A dropdown happens at nm; just one nanometer beyond a significant elevation appears and keeps being like that until nm. For higher radii the envelop functions for all levels only become more localized and the sequence regularly just rises.

Figure 3(a) shows the intersubband dipole momentums (ISDM). These elements are discussed only for single level for transition between s,p, and d sublevels so mainly they are so-called intersubband dipole momentums. Due to Simultaneous, Confinement effect and size promotion the ISDM have greater values than those of off-diagonal elements.

As Figure 3(a) shows the absolute value of(); when only lobe of the envelop function come inside the dot, the absolute value just keeps moving forward which looks like that the component take the lead.

The() has a ripple on nm because the first unit comes inside and the journey of the next one just begins. In the dimensions which the envelop function is part in/part out the dipole momentums take close value with almost zero slope. When the FES envelop function becomes localized inside the dot completely at nm it looks like the component takes the control and the absolute amplitude just keeps rising in the quasi-linear sequence.

For(), the situation is same with somehow more complexities since the SES consists of 3 lobes. Again like the other discussed parameters, the linear increasing behavior shows up on nm.

Optical rectification is purely real so according to (9) it only contributes to refractive index changes. Figure 4 shows the mutual refractive changes versus the field frequency for different sizes. Figure 4(a) shows that the significant changes on the OR coefficient begin to appear at nm in which right peaks take smaller values. At nm changes become more significant in which the positive right peaks take negative values and also the left peaks take smaller value. Figure 4(b) shows that in nm the negative take approximately 5 times larger value with redshift and the left peaks become negative. In nm suddenly the right peak faces new exchange and both peaks turn to become positive. The sudden change is mutual to the entrance of GS to dot completely. The same profile is repeated for upper radii with blue shift as well as for the height of the peaks untilnm as Figure 4(c) depicts. According to Figure 4(d) at nm suddenly a left peak rises and takes the higher value than the previous radius. At nm two peaks take close values and get closer to each other. At nm the right peaks turn to become negative but the left peaks remain positive with smaller absolute value. For upper radii the same profile repeats itself but the absolute amplitude of the peaks rises; for instance, in nm as shown in Figure 4(e) it get the value almost equal to 10 times larger than the absolute peaks ofnm.

Figure 5 represents the absorption due to SHG. According to Figure 5(a) changes appear on nm and become remarkable on nm where the absolute value of the peaks vanished to relatively small value. At  nm the left peaks as well as the right peaks take higher negative amplitudes; mainly only amplification shows up as Figure 5(b) illustrates. At  nm the left peaks take smaller absolute values but in positive region also the right peaks turn to become negative. nm and 8 nm reveal a vice versa profile. At  nm everything changes; double positive peaks appear on the right side alongside with other negative peaks on left side as Figure 5(c) depicts. According to Figure 5(d) in the couple of next steps forand 11 nm the situation is almost the same for the two positive peaks on right side and two negative in left side. On the other hand the peaks have relatively close absolute values. nm shows a reduction on the values as well as the entire change on the profile, in which one peak and valley occur on the right side and the left peaks become positive with so slight relative value. nm according to Figure 5(e) looks so weird since suddenly the right peak vanished and a big valleys takes place with two other right peaks. At nm the two successive valley on the left side of the graph take place while only one positive peak shows up on right side. The sequence keeps repeating itself with intensifying the values of the peaks versus the upper radii according to Figure 5(e).

Figure 6 shows the refractive index changes due to SHG. Just likethe remarkable impact of size promotion appears on nm as Figure 6(a) illustrates. For, 8, and 9 nm in Figures 6(b), 6(c), and 6(d) the refractive index changes (RFC) take a chaotic successive peaks and valleys and the values for these three mentioned radii take extraordinary variations. Atand 11 nm (Figure 6(e)) the RFC graphs take a regular behavior for these two steps with blue shift. Figure 6(d) shows the RFC forand 13 nm. The irregular result takes place for 12 nm without agreement with the preceding dimensions. nm is the beginning of new behavior for overall shape of RFC curves. In all upper researched radii (Figures 6(e) and 6(f)) 3 successive peaks and valleys are available. The second and third peaks-valley come closer to each other in the upper radii. On the other hand the absolute values of the biggest peaks face incredible intensification.

Figure 7 shows the biggest peak value for refractive index changes because of optical rectification (Figure 7(a)), for the absorption coefficient because of SHG (Figure 7(b)), and the refractive index changes owing to SHG (Figure 7(c)), respectively. The amplitude of the peaks is almost constant comparatively until nm. For upper radii the changes in the absolute value of the peaks somehow take hundred time larger values than those of smaller radii.

4. Conclusion

In this work the effect of size on transition frequencies, all elements of dipole momentum matrix, ORs and SHG are investigated deeply. All resonant and none-resonant expressions of OR and SHG are taken into consideration, which makes the result more useful for practical applications. The WL caused weird behavior in all discussed parameters. The transition frequencies between subbands satisfy the general expectation, which quotes.

At nm in does not face same circumstance on difference transition and in every subband it is the beginning for new sequence. For the absorption and refractive index changes due to optical rectification and SHG,  nm is always turning point. As inference it could be concluded that in QDs system with WL there is always critical radius after which the result takes regular behavior sequence. The order of the SHG both in real and imaginary part is relatively smaller than OR by almost the order of. The presented result quotes that, in InAs/GaAs QDs, between the second order nonlinear phenomena, the OR plays a dominant role, moreover it is possible to engineer the amplitudes by choosing proper and optimum size for QD.