Anisotropy and nonlinearity in superlattices II

This paper uses the anisotropic medium approach to calculate nonlinear optical absorption in semiconductor superlattices. The study confirms that the nonlinearities increase with anisotropy suggesting that strongly anisotropic materials such as those used for solar cells may also play a relevant role for nonlinear optical applications.

new materials being currently investigated for solar cells and which are strongly anisotropic (Steinmann et al. 2015), may also be very useful for nonlinear optics applications, e.g. the controllable anisotropy in the superlattice case might be useful for applications such as power limiting (Poirier et al. 2002;Wu et al. 2003). Here the anisotropy-induced nonlinearity in the nanostructure is controlled per design in contrast with less controllable recent studies in which it depends on the nanoparticle shape (Hua et al. 2015). We show that the nonlinearities increase with the anisotropic character, but the results are only valid as long as sufficient tunnelling occurs. In other words, isolated quantum wells should not be addressed by this method. This paper extends and complements the study in Pereira (2016) with different well widths and a more detailed analysis of the population inversion and is organized as follows: first, the main equations used are summarized. Next, numerical applications for GaAs/AlGaAs superlattices are given and a short summary follows.

Mathematical approach and model equations
The carrier-induced optical nonlinearities of semiconductors can be obtained by the selfconsistent evaluation of Many Body Nonequilibrium Green's Functions (NEGF). Efficient numerical methods used here have been successfully applied to both inter-subband (Pereira and Faragai 2014;Pereira 2007Pereira , 2008Pereira , 2011Pereira et al. 2007;Pereira and Tomić 2011) and inter-band transitions (Grempel et al. 1996;Pereira et al. 1994;Pereira and Henneberger 1997;Chow et al. 1992) in quantum wells and superlattices. For completeness, we repeat here the main equations of Pereira (2016) for a guide to the reader, by starting with the interband polarization that has been used to describe superlattices as an effective anisotropic 3D material (Pereira 1995), where the meaning of the anisotropy dependent bandgap is explained below and n(h, c) = sin 2 h ? ccos 2 h. The anisotropy parameter c is given by the ratio between the inplane l k and perpendicular l ? reduced effective masses, c ¼ l k . l ? , with 1 More details are given in Pereira (1995). A phenomenological scattering rate C is used to simulate the average dephasing that stems from the electron-electron, electron-phonon and electronimpurity scattering (Chow et al. 1992;Schmielau and Pereira 2009a, b), in order to keep the approach as simple as possible without affecting the conclusions. We have neglected any k-dependence on the transition dipole moment } induced by the electric field E(x). In a superlattices there is a preferred direction determined by the growth (z-direction), leading to cylindrical symmetry. The next step is to perform an angle average, n h; c ð Þ h i¼ The resulting anisotropic mass M k ¼ 2l k = 1 þ c ð Þdetermines the exciton Bohr radius a 0 and corresponding 1S binding energy E 0 , where e 0 and e denote, respectively the background dielectric constant and the electron charge. The resulting exciton binding energies in superlattices are in very good agreement with experiments, as demonstrated in Pereira (1995). By approximating the full inversion factor by A we make sure that the cross-over from absorption to gain takes place exactly at the total chemical potential. Next, we replace the usual Yukawa potential describing screened carriers in 3D (which does not have analytical solutions for the corresponding Schrödinger equation) by the Hulthén potential W r ð Þ ¼ À2e 2 j 0À1 0 = expð2jrÞ À 1 ð Þ (Pereira 1995;Bányai and Koch 1986;Flügge 1974), which has successfully reproduced bulk nonlinear optical spectra (Schmielau and Pereira 2009c) and has well known analytical solutions.
The total chemical potential l = l e ? l h and screening wavelength j ¼ j e þ j h are given by where K 1 = 4.897, K 2 = 0.045 and K 3 = 0.133. The parameter m k is obtained from the particle density for electrons and holes, n = n e = n h by m k Pereira 1995).
It is beyond the scope of this paper to show the intermediate details that lead Eq. (6). The main steps are to combine the partly phenomenological approach of Bányai and Koch (1986) with the material parameters calculated with the anisotropic medium approach (Pereira 1995;Pereira et al. 1990). Furthermore, we use hypergeometric functions instead of the Jost function approach of Pereira (1995) and Flügge (1974). The resulting absorption spectrum then reads where a 0 c ð Þ Here, and c denote, respectively, the background refractive index and the speed of light in vacuum. In Eq. (6) above, the band gap renormalization stems from the Mott criterion. This expression for the bandgap renormalization is usually in good agreement with the full Green's function approach and the Single Plasmon Pole Approximation (SPPA) simplified under a quasi-static approximation (Pereira 1995).
where g c ð Þ ¼ 1=ja 0 c ð Þ and E gap ¼ E c 1 ð Þ þ E HH 1 ð Þ þ E gap;bulk is the bulk (temperature dependent) isotropic bandgap found in semiconductor material tables in the literature plus the confinement energies of the lowest conduction subband and the top heavy hole subband.
The broadened delta function representation is d C x ð Þ ¼ 1=pdgcoshðx=gÞ, where g ¼ C=E 0 c ð Þ which reproduces the Urbach tail very efficiently. The sum in the exciton part runs through the available states within the largest integer value of ffiffiffiffiffiffiffiffi ffi g c ð Þ p . In the low density limit, g(c) ? ? and we recover the Elliot formula for excitonic luminescence with the correct balance between bound and continuum states. Figure 1 depicts the nonlinear optical absorption of GaAs-Al 0.3 Ga 0.7 As superlattices at T = 300 K. The scaled photon energy (x-axis) highlights the Coulomb correction effects that would otherwise be mixed with the changes due to carrier confinement by increasing the well width, included in The increase in anisotropy (smaller c) by increasing the quantum well length, stems from reducing the wavefunction overlap through the barrier, here measured by a larger mass along the growth direction. The single quantum well case would be obtained for a very large barrier, with l ? ! 1 and correspondingly c ! 0. However, this method should not be used in the quasi-2D case and calculations dedicated to quantum wells (see e.g. Pereira et al. 1994) should be used. From (a) to (d) we do see an increase in nonlinearity, i.e. a larger reduction in absorption with increasing carrier density. Note that even though Pereira (1995) also has analytical expressions for the nonlinear spectra of superlattices, to the best of that author's knowledge numerical applications have never been given before as in the present paper. Figure 2 makes it more clear by showing the differential absorption Da Figure 3 allows a better understanding of the excitonic bleaching and gain development as a function of anisotropy and carrier density.

Numerical results and discussion
The combined figures show that as the anisotropy increases (smaller c, so does the inverse screening length j leading to a faster reduction of the screened Coulomb interaction and consequently larger optical nonlinearity, measured here directly by the differential absorption. This shows that an increase in nonlinearity is observed for all carrier densities, which are considered in the manuscript, and in the whole range of energies depicted in Fig. 1. Note however that it only makes sense to consider the range of linear to saturated absorption in the differential transmission calculations while two positive values are compared. The negative absorption (gain) confuses the comparison.
By increasing the anisotropy, the z-direction electron and hole masses become too different, with a much larger increase in hole mass. When the upper and lower (average) curvatures or equivalently, the average effective masses are too different the total chemical potential is relatively smaller (see Fig. 3b) and so is the inversion factor A (1) and (6), thus reducing the gain in Fig. 1. This influence of different electron and hole masses on the inversion/gain is fully consistent with the detailed analysis for isolated quantum wells seen e.g. in Pereira and Henneberger (1997). The results are complemented by a two-dimensional plot of the population inversion factor as a function of detuning and anisotropy parameter in Fig. 4. Note the larger inversion region and values (the factor gets more negative) for smaller anisotropy (larger c). The method used here can be used for a variety of other materials and superlattices as long as tunnelling between adjacent periods allows a 3D-like spread of carriers wavefunctions-either electrons or holes so that movement along the z-direction is possible and an effective 3D medium can be considered with corresponding effective masses. Thus the model is better suited for absorption in superlattices far from the 2D limit. In the quasi-2D limit of quantum wells, the full numerical solution (Grempel et al. 1996;Pereira et al. 1994;Pereira and Henneberger 1997) should be used. Thus, increase in nonlinearity is clearly demonstrated, which may be very important for truly 3D anisotropic (1) as a function of detuning and anisotropy parameter at 300 K. The left and right panels are respectively for population densities N = 0.02 and 2 9 10 18 cm -3 new materials, but the high gain in quantum wells cannot be described by the method presented here. In summary, the analytical expressions developed confirm the connection between an increase in optical nonlinearity with anisotropy by directly controlling the anisotropy and evaluating the resulting differential absorption and nonlinear spectra. This study illustrates the potential for three dimensional strongly anisotropic semiconductors materials for a role in nonlinear optics and a possible recipe to improve their efficiency by controlling the anisotropy.