A theoretical study of the propagation of light soliton produced by semiconductor quantum
dots through optical fibers
O. P. Swami1, V. Kumar2, B. Suthar3, A. K. Nagar2
department of Physics, Government College, Loonkaransar, Bikaner, Rajasthan, 334001, India
2Department of Physics, Government Dungar College, Bikaner, Rajasthan, 334001, India 3Department of Physics, MLB Government College, Nokha, Bikaner, Rajasthan, 334803, India
[email protected], [email protected], [email protected], [email protected] PACS 42.65.-k, 42.65.Tg, 78.67.Hc DOI 10.17586/2220-8054-2019-10-3-273-281
In this paper, the propagation of light soliton is studied in nonlinear optical fiber. We propose the external excitation of semiconductor quantum dot (SQD) waveguides through an optical source which allows the generation of solitonic pulses. These solitonic pluses are propagated through a nonlinear optical fiber. Soliton formation is investigated by the interaction between the external optical excitation and SQDs. Here, the SQDs are considered as a quantum system of three energy levels. In this study, the Fourier Split-Step (FSS) method is used to solve the numerically continuous nonlinear Schrdinger equation (NLSE) for evolution of the soliton pulse emitted by the SQDs inside an optical fiber with real physical parameters. The effect of a SQD's density and electric field on the pulse width is also studied. Phase plane portraits are drawn to observe the stability of soliton in fiber and SQDs.
Keywords: quantum dots, nonlinear optics, optical solitons, nonlinear guided waves. Received: 13 March 2019 Revised: 1 April 2019 Final revision: 10 May 2019
1. Introduction
Semiconductor quantum dots (SQDs) are model systems for investigating nonlinear light-matter interaction that have attracted much interest. SQDs are referred to zero-dimensional systems that restrict the movement of charge carriers in the three spatial dimensions, which results in atom-like discrete energy spectra and strongly enhanced carrier lifetimes [1]. These atoms like structures have been proposed for the use as qubits in quantum information processing [2] as well as for laser devices [3]. Some of the most recent investigations indicate that these types of heterostructures can undergo abrupt changes in the spectral response with minimal variations in their size and morphology, offering important applications for optics, among which are next-generation lasers, diodes, light emitters, optical multiplexers, biosensors, spectral tuners, quantum computing, logic gates, among others [4-12]. Due to their large dipole moments, reaching values on the order of 10-17 esu cm [13], the interaction between SQDs and optical light fields is strongly enhanced in comparison with atomic systems, making them good candidates for the study of nonlinear optical propagation effects. Recently, it has become possible to combine nanostructures with other polymeric materials such as optical fibers, giving rise to nanocomposites, which are generally composed of several phases
such as SiO2, where one or more of its dimensions are found at the nanoscale [14-22]. The optical fibers are important elements for the propagation of light wave signals with water and for the propagation of electromagnetic waves in the visible, ultraviolet and infrared range. In an optical fiber, it is necessary to consider a theory of wave propagation in dispersive media. The nonlinear Schrodinger equation (NLSE) gives a complete description of a variety of the localized non-linear effects that have been extensively studied in various contexts of the sciences and that the theory can directly link with the propagation of intense optical pulses in non-linear optical fibers that give rise to the optical solitons. Although there are many experimental advances in the propagation of solitons in optical fibers, in practice, the propagation of solitary pulses across large distances at a commercial level has not been possible due to the different technical functions [23-26]. The study allows us to propose a general model which permits the coupling of solitons from a set of quantum dots with specific characteristics in an optical fiber that has non-linear optical characteristics. As a result, a numerical simulation is developed to study the evolution of the soliton inside a non-linear optical fiber with the Fourier Split-Step numerical technique with the real parameters associated with the SQD and the optical fiber. The effect of soliton density and electric field on the full width at half maximum (FWHM) of re-emitted solitonic pulse is studied. Phase planes represent the stability/instability of solitons.
2. Theoretical analysis
SQDs can confine the movement of charge carriers in all dimensions due to their zero dimensional nanosized structures, which exhibits a discrete energy system. We consider a two-dimensional sheet of inhomogeneously broadened SQDs which forms a transition layer on one of the surfaces of the planar waveguide. The waveguide soliton satisfies the Maxwell-Bloch equations with nonlinear boundary conditions. The purpose of the present article is to theoretically investigate the formation processes of optical solitons under the condition of SIT in a semiconductor waveguide.
A SQD of three-level energy system is considered, whose ground state is and has energy = 0. The states
5x
|^2) and 1^3} have energies = hv0 = £x
+ — and
hv0 = 2ex — Sxx, respectively. Where the quantities
ex = + £^2 )/2 and describe the energies of the single-excitonic and biexcitonic states, respectively.
Fig. 1. Schematic diagram of energy structures of semiconductor quantum dots
The energy of exciton fine structure splitting is denoted by Sx = - and biexciton binding energy is denoted by 5xx. h is Planck constant. The Hamiltonian of system is described by:
H
Ho + H,
(1)
1
where H0 = hu01^2X^21 +2 hu'o l^3}(^3| is called the Hamiltonian of single excitonic state } and biexcitonic state The additional term H in the Hamiltonian is present due to the interaction of the light pulse with SQDs and H = ^ . During the excitation of the SQDs, an external light source is considered, composed of a linearly-polarized, high intensity optical wave from a laser and we will study the formation of non-linear optical waves that are re-emitted by the SQDs. Fig. 2 represents a schematic of a linearly-polarized plane wave incident on SQDs. In the process, due to the interaction of the light with the SQDs, the light is re-emitted in the form of a non-linear wave with special characteristics that will depend on the morphology of the quantum dot and the incident wave.
Fig. 2. Schematic arrangement of SQDs and optical fiber with an external light source
In the analytical treatment the electric field vector of incident light beam is taken as below:
e(z,t) eEivi,
i=±i
where:
i
+ i"1 ),
(2)
(3)
is the complex polarization vector. lX and V are unit vector along the x and y axes. <p l = e
il(kz-uit)
is a periodic
wave function and E is the slowly changing amplitude of the electric field. The incident light is considered linearly-polarized in TE mode and has a width T and angular frequency w > T — 1 propagating along z direction. When an external light beam interacts with SQDs, then a secondary dipolar field is produced in same direction as incident field and can be written as:
_c2 d2E(z,t) + d2E (z,t)
dz2
dt2
-An
d2P (z,t) dt2 ,
where the polarization component of ensemble of SQDs is the following:
11 (z,t) = l no
g (vol) e(p!2P2i + P23P32)(u, z,t)dv + c.c.,
(4)
(5)
which is determined by the interband transitions occurring in the quantum dots between the three energy levels. The quantities p12 and p23 are the dipole moments for the corresponding transitions which we assume to be parallel to each other. g (u01) is an inhomogeneous broadening function which arises due to the quantum dot size fluctuations. Because the energy levels depend on the size of dots, the exciton frequency u0 and detuning quantity u01 = u0 — u1 are also size dependent. n is the refractive index of semiconductor and n0 is the density of Quantum Dots. pj are the matrix elements of the density matrix and p can be determined from Liouville equation:
ihPmn = ^2 ((n \H| I) pim _ Pnl {I \H| m)) ,
where n, m, I = 1, 2, 3.
The density matrix of the pulse is considered of the form:
(6)
(7)
P21 = P21<f3, P32 = P32<fl, P31 = P31P2,
where pnm is the slowly varying complex amplitudes of the elements of density matrix. <p1t2,3 is taken in for envelope function. To guarantee that E is a real number we choose slow envelope approach, it is considered that amplitude El of pulse is very smooth in space and time and can be written as:
dEi
and
dt
dEi
dz
< ui
< k
Ei
Ei
(8)
(9)
then nonlinear carrier wave equation can be written as following form:
^ Vi
i=±i
dEi
dEi
I2 (c2k2 _ r]2v21) Ei _ 2iklc2-^— _ 2ilri2v1-^-
dt
= -2nl2-2
u1 no
P12P21 +P23P32
ei(kx-Ult) +c.c. (10)
The diagonal elements p12, p22 and p33 give rise to the populations in the states |^1), |^2) and |^3) respectively. Non-diagonal elements pmn(n = m) contain the relative phase between the states that describe the atomic coherence.
In the absence of phase modulation Ei amplitudes:
Ep-i = Ei = E. Therefore, the system of equations for the slowly varying
ihpii = (p*2i _ P21) P12E, ifrp22 = (P21 _ P21 + 5 P*32 _ à P 32 ) P12E,
ifrp33 = (p32 _ P*32) P12E, ihp2i = h (uo _ u) P21 _ [(pii _ P22) + àP31 ] P12E, ihp32 = h(y0 _ uo _ u) P32 + [à (p33 _ P22) + P3i] P12E, ihP3i = h(u'o _ 2u) P31 + [p32 _ àP21] P12E.
In the above equations, the rotating wave approximation has been applied and S = —. The solution of the
M12
system of equations (11), allows the calculation of the elements of the right side of equation (10). As the diagonal elements of matrix give rise to the populations in the states, for ground state p11 = 1, p11 = 0 and p33 = 0, the equations (11) take the form:
P21 = ^3 (sin 2Ad + 2S2 sin Ad) , iS
p32 = 2d3 (sin 2Ad - 2 sin Ad),
(12)
Ott I ft
where d = ^r+SS2) and A = ^ E ( z, t)dt, which is called the area of nonlinear optical pulse.
h J-oo
(2V0« + z^VOzt +-sin(2^) = 0, (13)
h
Applying the dispersion law for linear wave guide modes, wave equation (4) will take following form:
— \2w)zt +--r-2
n hrj2
where
^ = Ad. (14)
Equation (13) can be written in the form of electric field as below:
(d,E\ _E2 nV22£4
d£ T2 h2
(15)
z
with £ = t - — and constant phase velocity V. The solution of equation (15) for the envelope function has the form:
EE= —FJh-sec h(-), (16)
which is well-known solitonic solution, with width of the pulse remitted by SQD:
T y 4n2uno^2^ -V )' (17)
These soliton became stable when the speed of light (phase velocity) in SQD media c/n remains greater than the phase velocity V.
3. Result and discussion
The theoretical result of wave equation (4) is presented as a soliton. In other words, the excitation of SQDs from an intense nonlinear wave can produce optical solitons as a result of light-SQDs interaction. During the nonlinear interaction process, the SQDs are considered as a three-level energy quantum system, in which the optical transitions are given from the ground state to the biexcitonic and excitonic states. The allowed transitions between the ground state to the biexcitonic and excitonic states have very low dipole moment compared to the transition between the ground state and the background of the exciton band. On the other hand, the characteristics of the light re-emitted by the sheets of quantum dots in the form of optical solitons depend on the intensity of the incident light, whose minimum value to form the optical solitons can be specifically determined depending on the nature of the SQDs. In this way, the soliton remitted will depend on the refractive index of the semiconductor n, the dipole moment corresponding to the transitions between the ground state and the excitonic state or between the biexcitonic state and the excitonic state ^i2. In order to study the propagation of short optical pulses through nonlinear optical fibers, the non-linear Schrodinger equation (NLSE) is used, which takes into account the effects of the length of the fiber, the dispersion effect of group speed and non-linear optical effects as a consequence of the high intensity of light. The Fourier Split-Step method is a pseudo-spectral technique that is extremely useful due to its rapidity and high accuracy in calculations. In general, this method obtains an approximate solution of the propagation equation, assuming that dispersion effects and non-linear effects act independently along the fiber in very small steps. This technique was used to simulate the evolution of solitons that are re-emitted by the SQD.
In the results, it is observed that as the density of SQD increases, the peak intensity also increases. Because as the number of QDs per unit volume increases, absorption will also increases. Therefore the intensity of re-emitted beam will increase. Due to this the pulse width (FWHM) will also decrease. Same concept is applies for electric field verses FWHM. In the mathematical modelling of proposed work, the interactions of QDs in the Hamiltonian are omitted. For the simulation, we have proposed SQDs with pyramidal morphology of In as manufactured with a cylindrical symmetry with the parameters n = 2.11, ^12 = 1.92 x 10-28 C-m and V = 1.37 x 108 m/s.
Figures 3(a-d) represent the simulation of solitons created by SQDs arrays with different densities of SQDs according to following table. It is clear from equation (17) that the pulse width decreases as n0 increases. In many optical fiber communications systems, FWHM is required 1 ps or below, which can be achieved by increasing the density of SQDs. Simulated data are summarized in Table 1.
(a)
(b)
-4 -4
(c)
(d)
Fig. 3. Simulation of soliton generated by SQDs system with different densities of QDs.(a) nQ 2.59 x 1012 cm-3; (b) nQ = 8.00 x 1012 cm-3; (c) nQ = 23.35 x 1012 cm-3 and (d) nQ 99.12 x 1012 cm-3
Table 1. Relation between the density of SQDs and FWHM of solitons
S.No. Density of SQD(1012 cm-3) FWHM(ps)
1. 2.59 6.18
2. 8.00 3.52
3. 23.35 2.06
4. 99.12 1.00
Figures 4(a-d) represent that the soliton width as determined from equation (16) is decreased as electric field E is increased. According to Table 2, it is clear that when the amplitude of the electric field is increased, then the FWHM of TM mode of wave decreases, which is a basic requirements for many optical fiber communications systems.
Figures 5(a-d) represent the evolution of soliton profiles in SQDs for different values of refractive index of material of SQDs as taking V = 1.37 x 108 m/s.
Figures 6(a-d) represent the phase plane portraits for Figs. 5(a-d) respectively for stability analysis. Closed limit cycles of phase plane structures clearly indicate a stable soliton, while spiral-like phase plane structures indicate an unstable soliton. It is clear from below Table 3 that as c/n become less than V, the soliton become unstable, which is summarised in Table 3.
(a)
(b)
(c)
(d)
Fig. 4. Transverse view of soliton in SQD with electric field (a) E = 0.96 x 108 V/m; (b) E = 1.73 x 108 V/m; (c) E = 2.86 x 108 V/m and (d) E = 6.00 x 108 V/m
Table 2. Relation between the amplitude of Electric field and FWHM of TM mode of solitons
S.No. Amplitude of electric field(108 V/m) FWHM(ps)
1. 0.96 6.26
2. 1.73 3.46
3. 2.86 2.10
4. 6.00 1.00
Table 3. Relation between refractive index of material of SQDs and (in)stability of soliton
S.No. Refractive index(n) of SQDs c/n (m/s) Nature of soliton
1 2.05 1.46 x 108 Stable
2 2.11 1.A2 x 108 Stable
3 2.40 1.25 x 108 Unsable
4 3.50 0.86 x 108 Unstable
(a)
(b)
(c)
(d)
Fig. 5. Simulation of soliton generated by SQDs system with different refractive indices of QDs. (a) 2.05; (b) 2.11; (c) 2.40 and (d) 3.50 taking V = 1.37 x 108 m/s
Figure 7(a) represents the soliton profiles for input pulse in optical fiber and the Fig. 7(b) represents the pulse in optical fiber after travelling the distance of z = 1000 m. Simulations are made using split-step Fourier method with the parameters: nonlinear fiber parameter y = 0.32 W/m, fiber attenuation constant a = 0.2 dB/km and second order dispersion coefficient 3 = —2 x 10-26 s2/m.
4. Conclusions
In this paper, we have investigated the formation of optical solitonic pulses in three level semiconductor quantum dots embedded in nonlinear optical fibers for propagation of fields without losses over long distances. Light solitons are produced in semiconductor quantum dots by the interaction of laser field with SQDs. The FWHM of solitonic pulses can be modulated according to densities of SQDs and field intensity. If the density of SQD is increased, the peak intensity is also increased. Because as the number of QDs per unit volume is increased, absorption will also be increased. Therefore the intensity of re-emitted beam will be increased. Due to this, the pulse width (FWHM) will also be decreased. The same concept is applied for electric field verses FWHM. The stability/instability of generated solitons is dependent on the refractive indices of SQDs material. Likewise, it is proposed the manufacture of these types of nanostructures with different morphologies inserted in optical fibers that allow the propagation of solitons without losses over the long distances.
Fig. 6. Phase plane portraits for solitons for different refractive indices of QDs.(a) 2.05; (b) 2.11; (c) 2.40 and (d) 3.50 taking V = 1.37 x 108 m/s
(a) (b)
Fig. 7. Evolution of soliton through nonlinear optical fiber at (a) z = 0 m and (b) z = 1000 m References
[1] Pavesi L., Dal Negro L., et al. Optical gain in silicon nanocrystals. Nature, 2000, 408, P. 440-444.
[2] Hu Y.M., Yang W.L., Feng M., Du J.F. Distributed quantum-information processing with fullerene-caged electron spins in distant nanotubes. Phys. Rev. A, 2009, 80, P. 022322(1-13).
[3] Maximo L., Mendez V.H. Autoensamblado de puntos cunticos. Cinvestav, 2008, 27, P. 44-49.
[4] Hanewinkel B., Knorr A., Thomas P., Koch S.W. Optical near-field response of semiconductor quantum dots. Phys. Rev. B, 1997, 55, P. 13715(1-10).
[5] Aknmanov S.A., Visloukh V.A., Chirkin A.S. Optics of Femtosecond Laser Pulses. New York, American Institute of Physics, 1992.
[6] Klimov V.I., Mikhailovsky A.A., et al. Optical gain and stimulated emission in nanocrystal quantum dots. Science, 2000, 290, P. 314-317.
[7] Yang W.X., Chen A.X., Lee R.K., Wu Y. Matched slow optical soliton pairs via biexcition coherence in quantum dots. Phys. Rev. A, 2011, 84, P. 013835(1-11).
[8] Hu L., Wu H., et al. The effect of annealing and photoactivation on the optical transitions of band-band and surface trap states of colloidal quantum dots in PMMA. Nanotechnology, 2011, 22, P. 125202(1-8).
[9] Adamashvili G.T., Weber C., Knorr A. Optical nonlinear waves in semiconductor quantum dots: Solitons and breathers. Phys. Rev. A, 2007, 75, P. 1-9.
[10] Newell A.C. Solitons in Mathematics and Physics. Philadelphia, Society for Industrial and Applied Mathematics, 1985.
[11] Ajayan P.M., Schadler L.S., Braun P.V. Nanocomposite science and technology. Weinheim, Germany, Wiley, 2003.
[12] Rothenberg J.E., Grischkowsky D., Balant A.C. Observation of the Formation of the 0n Pulse. Phys. Rev. Lett., 1984, 53, P. 552-555.
[13] Chen M., Kaup D.J., Malomed B.A. Three-wave solitons and continuous waves in media with competing quadratic and cubic nonlinearities. Phys. Rev. E, 2004, 69, P. 056605(1-17).
[14] Bimberg D., Grundmann M., Ledentsov N.N. Quantum Dot Heterostructures. Chichester, Wiley, 1999.
[15] Schneider S., Borri P., et al. Self-induced transparency in InGaAs quantum-dot waveguides. Appl. Phys. Lett., 2003, 83, P. 3668-3670.
[16] Panzarini G., Hohenester U., Molinari E. Self induced transparency in semiconductor quantum dots. Phys. Rev. B, 2002, 65, P. 165322(1-6).
[17] Klimov V.I. Nanocrystal quantum dots. Los Alamos Science, 2003, 28, P. 214-220.
[18] Chen Y., Herrnsdorf J., et al. Colloidal quantum dot random laser. Opt. Express, 2011, 19, P. 2996-3003.
[19] Akhmediev N., Ankiewicz A. Solitons, Nonlinear pulses and Beams. London, Chapman and Hall, 1997.
[20] Kevrekidis P.G., Rasmussen K.O., Bishop A.R. Two-dimensional discrete breathers: Construction, stability, and bifurcations. Phys. Rev. E, 2000, 61, P. 2006(1-7).
[21] Blum K. Density Matrix Theory and Applications, 2nd ed. New York, Plenum Press, 1996.
[22] Adamashvili G.T., Weber C., Knorr A. Optical solitons in semiconductor Quantum dots. Eur. Phys. J. D, 2008, 47, P. 113-117.
[23] Malomed B.A., Kevrekidis P.G. Discrete vortex solitons. Phys. Rev. E, 2001, 64, P. 026601(1-6).
[24] Lamb Jr.G.L. Elements of soliton theory. New York, Wiley, 1980.
[25] Lamb Jr.G.L. Optical waves in crystals. New York, Wiley, 1980.
[26] Oster M., Johansson M. Stable stationary and quasiperiodic discrete vortex breathers with topological charge S =2. Phys. Rev. E, 2006, 73, P. 066608(1-6).