Vectorial polariton solitons in semiconductor ... - OSA Publishing

Vectorial polariton solitons in semiconductor microcavities W. L. Zhang and S. F. Yu* School of Electrical and Electronic Engineering, Nanyang Technological University, 639798 Singapore * [email protected]

Abstract: This paper presents numerical studies of vectorial polariton solitons in semiconductor microcavities. In the simulation, polarization degree of freedom of the polariton fields is taken into consideration. In the bistable regime, bright and/or dark solitons are found to bifurcate from the homogonous solutions of the two circular polarization modes. The combinations of solitons in the two polarization directions can be brightdark, dark-bright, bright-bright, and dark-dark. ©2010 Optical Society of America OCIS codes: (140.3945) microcavities; (240.5240) polaritons; (260.5430) polarization

References and links 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19.

C. Weisbuch, M. Nishioka, A. Ishikawa, and Y. Arakawa, “Observation of the coupled exciton-photon mode splitting in a semiconductor quantum microcavity,” Phys. Rev. Lett. 69(23), 3314–3317 (1992). A. Amo, D. Sanvitto, F. P. Laussy, D. Ballarini, E. del Valle, M. D. Martin, A. Lemaître, J. Bloch, D. N. Krizhanovskii, M. S. Skolnick, C. Tejedor, and L. Viña, “Collective fluid dynamics of a polariton condensate in a semiconductor microcavity,” Nature 457(7227), 291–295 (2009). A. Baas, J. Ph. Karr, H. Eleuch, and E. Giacobino, “Optical bistability in semiconductor microcavities,” Phys. Rev. A 69(2), 023809 (2004). D. Bajoni, E. Semenova, A. Lemaître, S. Bouchoule, E. Wertz, P. Senellart, S. Barbay, R. Kuszelewicz, and J. Bloch, “Optical bistability in a GaAs-based polariton diode,” Phys. Rev. Lett. 101(26), 266402 (2008). C. Ciuti, and I. Carusotto, “Quantum fluid effects and parametric instabilities in microcavities,” Phys. Status Solidi 242(11), 2224–2245 (2005). D. Sarchi, M. Wouter, and V. Savona, “Polariton parametric photoluminescence in spatially inhomogeneous systems,” Phys. Rev. B 79(16), 165315 (2009). K. G. Lagoudakis, M. Wouters, M. Richard, A. Baas, I. Carusotto, R. André, L. S. Dang, and B. DeveaudPlédran, “Quantized vortices in an exciton-polariton condensate,” Nat. Phys. 4(9), 706–710 (2008). R. Balili, V. Hartwell, D. Snoke, L. Pfeiffer, and K. West, “Bose-Einstein condensation of microcavity polaritons in a trap,” Science 316(5827), 1007–1010 (2007). Y. Larionova, W. Stolz, and C. O. Weiss, “Optical bistability and spatial resonator solitons based on excitonpolariton nonlinearity,” Opt. Lett. 33(4), 321–323 (2008). A. V. Yulin, O. A. Egorov, F. Lederer, and D. V. Skryabin, “Dark polariton solitons in semiconductor microcavities,” Phys. Rev. A 78(6), 061801 (2008). A. R. Davoyan, I. V. Shadrivov, and Y. S. Kivshar, “Self-focusing and spatial plasmon-polariton solitons,” Opt. Express 17(24), 21732–21737 (2009), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-24-21732. G. Tissoni, L. Spinelli, L. Lugiato, M. Brambilla, I. Perrini, and T. Maggipinto, “Spatio-temporal dynamics in semiconductor microresonators with thermal effects,” Opt. Express 10(19), 1009–1017 (2002), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-10-19-1009. Y. Tanguy, T. Ackemann, W. J. Firth, and R. Jäger, “Realization of a semiconductor-based cavity soliton laser,” Phys. Rev. Lett. 100(1), 013907 (2008). X. Hachair, G. Tissoni, H. Thienpont, and K. Panajotov, “Linearly polarized bistable localized structure in medium-size vertical-cavity surface-emitting lasers,” Phys. Rev. B 79, 011801 (2009). O. A. Egorov, D. V. Skryabin, A. V. Yulin, and F. Lederer, “Bright cavity polariton solitons,” Phys. Rev. Lett. 102(15), 153904 (2009). I. A. Shelykh, Y. G. Rubo, G. Malpuech, D. D. Solnyshkov, and A. Kavokin, “Polarization and propagation of polariton condensates,” Phys. Rev. Lett. 97(6), 066402 (2006). C. Leyder, T. C. H. Liew, A. V. Kavokin, I. A. Shelykh, M. Romanelli, J. Ph. Karr, E. Giacobino, and A. Bramati, “Interference of coherent polariton beams in microcavities: polarization-controlled optical gates,” Phys. Rev. Lett. 99(19), 196402 (2007). N. A. Gippius, I. A. Shelykh, D. D. Solnyshkov, S. S. Gavrilov, Y. G. Rubo, A. V. Kavokin, S. G. Tikhodeev, and G. Malpuech, “Polarization multistability of cavity polaritons,” Phys. Rev. Lett. 98(23), 236401 (2007). D. V. Skryabin, O. A. Egorov, A. V. Gorbach, and F. Lederer, “One-dimensional polariton solitons and soliton waveguiding in microcavities,” Superlattices Microstruct. 47(1), 5–9 (2010).

#131354 - $15.00 USD

(C) 2010 OSA

Received 13 Jul 2010; revised 16 Aug 2010; accepted 16 Aug 2010; published 22 Sep 2010

27 September 2010 / Vol. 18, No. 20 / OPTICS EXPRESS 21219

1. Introduction Cavity polaritons, which exhibit neither exciton nor photon characteristics, are half-matter half-light quasi-particles arising from the strong coupling between excitons and photons [1,2]. The nonlinear effect of polaritons in quantum well semiconductor microcavities is an interesting area of recent research. Related topics including polariton bistability [3, 4], parametric wave-mixing [5, 6], vortices and Bose-Einstein condensation [7, 8] have been investigated extensively. Furthermore, studies of spatial solitons in a strong-coupling regime have initialized the investigation of polariton solitons in microcavities [9–11]. This is because, when compared with the weak-coupling regime [12–14], the nonlinear effect and optical response speed of semiconductor microcavities can be significantly improved by 2 to 3 orders of magnitude in the strong-coupling regime. Hence, the response time and the required pumping intensity to operate polariton solitons can be efficiently reduced. In addition, it is expected that both lifetime and relaxation kinetics of polaritons can be easily modified by changing the photon-exciton detuning, which can be used to design a specific dispersion relation for soliton formation. Recently, dark and bright polariton solitons of semiconductor microcavities have been reported [10, 15] and verified experimentally [2, 9]. However, these studies only considered the scalar polariton solitons, and the polarization of polariton fields was ignored. It is known that the polariton fields in microcavities are vectorial polarization dependence which governed by the underlying spin-dependent polariton-polariton scattering processes. Two possible spins, where are + 1 and –1 correspond to the right circular polarization (RCP) and left circular polarization (LCP) respectively, are supported inside microcavities [16–18]. This polarization dependence of polaritons adds the complexity of their nonlinearities and provides more flexible ways to form solitons. Hence, it is of great importance to study the polariton solitons with polarization dependence. In this paper, vectorial polariton solitons is studied numerically in a semiconductor microcavity. A modified 4th-order Runge-Kutta method is used to find the soliton solutions. Based on this method, various combinations of polariton solitons in the two polarization directions are detected and analyzed for both vertical and oblique pumping. 2. Theoretical model The polariton fields are usually described by two coupled equations for their excitonic and photonic components. We start our discussions with the polarization-resolved dimensionless model for the intracavity photon field, E, and for the exciton field, ψ [10, 15–19]

∂ t E1,2 = i ( ∂ 2x + 2ikin ∂ x − kin2 ) E1,2 − ( γ p − iδ ) E1,2 + iψ 1,2 + Ein1,in 2

(

∂ tψ 1,2 = iE1,2 − ( γ e − iδ )ψ 1,2 − i ψ 1,2 + η ψ 2,1 2

2

)ψ

1,2

(1) (2)

where subscripts 1 and 2 correspond to the RCP and LCP modes, respectively. γp (γc) is the photon (exciton) decay rate normalized to the Rabi frequency, Ω. δ is the normalized frequency detuning, i.e., δ = (ωin – ω0)/Ω wherein ωin is the pump frequency and ω0 denotes the identical resonance frequency of exciton and cavity. The time, t, transverse coordinate, x, and wave number, kin, are normalized to 1/Ω, x0, and 1/x0 respectively. x0 = c 2 / 2n 2ωin Ω , where c is the vacuum light velocity and n is the refractive index. The nonlinear terms in (2) denote the contributions of exciton-exciton interctions in the same (|ψ1,2|2ψ1,2) and opposite (η|ψ2,1|2ψ1,2) polarization directions. For zero-pumping condition, if we neglect the nonlinear effect and assume E1,2 = E1s,2sexp(ikx), ψ1,2 = ψ1s,2sexp(ikx) and ω0 = 0, the dispersion relationship of polaritons can be obtained through the eigenvalues of coupled Eqs. (1) and (2). This paper only refer to the lower polariton branch, ELP, and it can be expressed as [10, 15]

#131354 - $15.00 USD

(C) 2010 OSA



2 1 2 1  k − i (γ p + γ e )  −  k 2 + i (γ e − γ p )  + 4 (3)  2 2 In order to find the bright soliton solutions for a value of Ein, we first guess the initial values of excitons, ψ1,2(0, x), and photons, E1,2(0, x). This can be done by assuming that the analytical expressions of ψ1,2(0, x), and E1,2(0, x) to be ψ1,2(0, x) = ψ1s,2s[1 + A1,2exp(−B1,2x2)] and E1,2(0, x) = E1s,2s[1 + C1,2exp(−D1,2x2)] respectively, where ψ1s,2s and E1s,2s are the homogenous solutions in the lower bistable branch, and A1,2, B1,2, C1,2, and D1,2 are real constants. The differential terms ∂xψ1,2(0, x) and ∂ 2xψ 1,2 (0, x ) can be deduced from central difference method. Hence, ψ1,2(∆t, x) and E1,2(∆t, x) can be obtained for the next time step by 4th-order Runge-Kutta method (∆t denotes the time step). Repeating the above steps, we can obtain the values of ψ1,2(m∆t, x) and E1,2(m∆t, x) at any time step where m is an integer. As unstable solitons are sensitive to the initial conditions, ψ1,2(0, x) and E1,2(0, x) should be carefully selected, i.e., different values of A1,2 are tested until the profiles of ψ1,2 and E1,2 remain unchanged with the increase of t. On the other hand, the soliton solutions are stable if they are insensitive to the initial conditions. Using similar method, dark solitons can also be excited. In the following simulations, we concentrate our investigation on the bistable regime and the excitonic component is taken as an example to discuss polariton solitons. Typical parameters of a GaAs/AlAs microcavity with InGaAs/GaAs quantum wells as the gain region are used to model cavity polaritons [10, 15–17]. Hence, the values of the mentioned parameters are assumed to be: η = –0.1, γp = 0.1, γe = 0.1, Ein1 = Ein, and Ein2 = 0.8Ein. δ, kin and Ein are variables to be specified latter.

ELP (k ) =

3. Pumping at normal incidence

In this section, we consider the case of normal incidence, i.e., kin = 0. Figure 1 shows the homogonous solutions of the excitonic components of polaritons in the two polarization directions. The solid (dotted) curves mark the stable (instable) states. We can see that the two polarization modes show a distinct ‘S’ shape as a function of pump strength, Ein. Because the pumping into the two polarization modes is asymmetric, the hysteresises of the two modes are not overlapped. A vertical-broken line separates Fig. 1 into two regions. In region I (II), the RCP (LCP) mode exhibits a dominate hysteresis. Its counterpart, the LCP (RCP) mode, also shows a reverse (similar) hysteresis due to nonlinear interaction between the two modes.

Fig. 1. The homogenous solutions of the excitonic components of the polaritons as a function of pumping strength, Ein, for kin = 0 and δ = −0.7.

Next, the solutions of vectorial solitons are searched along the bistable curves of the homogenous solutions, and the results are given in Fig. 2. The squares (dots) denote the peak (dip) values of bright (dark) solitons, and the stable solitons are also marked by solid lines. Figures 2(a) and (c) correspond to region I. In Fig. 2(a), a branch of bright and a branch of dark solitons are found to bifurcate from the homogonous solutions of the RCP mode. The

#131354 - $15.00 USD

(C) 2010 OSA



entire branch of bright solitons is unstable, and the branch of dark solitons is separated into a stable and an unstable part by a fold bifurcation point.

Fig. 2. Maxima (minima) of the excitonic components of bright (dark) polariton solitons as a function of pumping strength, Ein, for kin = 0 and δ = −0.7. The solid-dotted curves represent the background amplitude, and the squares (dots) denote the peak (dip) values of the bright (dark) solitons.

Fig. 3. The excitonic components of polaritons in time and real space domains. (a) and (b) correspond to, respectively, points P1 and P1’ marked in Fig. 2. For clarity, the dark soliton is plot using negative value of |ψ2|.

Fig. 4. Profiles of the excitonic components of polaritons. (a) – (f) correspond to, respectively, points P2, P3, P4, P2’, P3′, and P4’ marked in Fig. 2.

Furthermore, due to nonlinear interaction, the bright (dark) solitons of RCP mode also arouse small-amplitude dark (bright) solitons in the LCP mode, see Fig. 2(c). Hence, vectorial bright-dark and dark-bright solitons are found in region I. In region II, bright-bright and darkdark combinations of solitons are found, see Fig. 2(b) and (d). Figure 3 gives an example for the excitation of polariton solitons through the proposed method. Figure 3(a) and (b) correspond to points P1 and P1’ marked in Fig. 2 respectively

#131354 - $15.00 USD

(C) 2010 OSA



(i.e., vectorial bright-dark solitons). For clarity, the dark solitons of LCP mode are plot using negative values of |ψ2|. We can see that the profiles of ψ1 and ψ2 become unchanged after a short time interval of oscillations from t = 0. It is also noted that the solitons are stationary, i.e., the peak position of solitons in the x direction does not change with time. Figure 4 shows the profiles of other types of vectorial solitons, i.e., dark-bright, brightbright, and dark-dark solitons that correspond to points P2 – P2’, P3 – P3′ and P4 – P4’ respectively marked in Fig. 2. These entangled states of vectorial solitons induce additional flexibilities and are thought to provide a promising way for information processing. 4. Pumping at oblique incidence

This section considers the case at oblique incidence. We choose kin = 1.7, δ = –0.03, and the microcavity maintains within the bistable regime. In this case, the bistable curves and the maxima (minima) of bright (dark) solitons are given in Fig. 5. In Fig. 5(a), we find a branch of bright solitons for the RCP mode. Different from the case of normal incidence, a part of the bright solitons becomes stable. In the branch of solitons, the solutions are stable (instable) above (below) the fold bifurcation point, which is in accordance with the result of [15, 19]. Corresponding to the bright solitons of RCP mode, a branch of dark solitons of LCP mode is found in Fig. 5(c). Thus, vectorial bright-dark solitons are formed. In region II, vectorial bright-bright solitons are found, see Fig. 5(b) and (d).

Fig. 5. Maxima (minima) of the excitonic components of bright (dark) polariton soliton as a function of pump strength, Ein, for kin = 1.7 and δ = −0.03. Other descriptions are the same as those of Fig. 2.

Figure 6 repeats the plot of Fig. 3 except that the pumping is oblique, i.e., kin = 1.7. In this case, the profiles of the solitons become asymmetric and the solitons is nonstationary. From Fig. 6, the velocity of the solitons is deduced to be 0.31. It is close to the first order dispersion of the polariton, ∂kELP, which equals 0.3 for kin = 1.7. From the above results, we can see that the vectorial solitons are the combinations of a strong and a weak solitons in the two polarization directions. In region I (II), the strong and weak solitons are of the same (reverse) type. However, these results are obtained from a relative small value of |η|. In some cases (i.e., the active region has multiple quantum wells), the value of |η| may be even larger. Strong nonlinear interactions can greatly change the entangled solitons in the two polarization directions. In Fig. 7, η = –0.5 is considered and the results correspond to region II. Different from the case for η = –0.1 (i.e., vectorial brightbright solitons are found in region II), the combination of solitons is dark-bright, and the amplitudes of solitons are large for both polarizations. We also find that all the soliton solutions are stable.

#131354 - $15.00 USD

(C) 2010 OSA



Fig. 6. The excitonic components of polaritons in time and real space domains. (a) and (b) correspond to, respectively, points P1 and P1’ marked in Fig. 5. For clarity, the dark soliton is plot using negative value of |ψ2|.

Fig. 7. Maxima (minima) of the excitonic components of bright (dark) polariton solitons as a function of pump strength, Ein, for kin = 1.7, δ = −0.03, and η = −0.5. Other descriptions are the same as those of Fig. 2.

5. Conclusion

In conclusion, we have analyzed the formation of vectorial polariton solitons in semiconductor microcavities. In the simulation, RCP and LCP polaritons are in different bistable regimes due to an asymmetric pumping. For normal pumping, bright and dark solitons can bifurcate from homogonous solutions of one polarization mode which shows a dominant bistability. In the counterpart polarization mode, reverse (in region I) or same (in region II) types of solitons with small amplitudes are aroused due to nonlinear interaction. Thus, vectorial bright-dark, dark-bright, bright-bright, and dark-dark solitons are found. For oblique pumping, bright-dark and bright-bright solitons are found in regions I and II respectively. Moreover, there is a fold bifurcation point separating the soliton branchs into a stable and an instable part. The strength of nonlinear interaction between the two plarization modes also plays an important role in the formation of solitons, i.e., when the two polarization modes have strong interactions, they have reverse types of solitons in region II, and the amplitudes of solitons for both modes are large. Acknowledgments

The authors gratefully acknowledge the support of A*STAR SERC grant no. 082-101-0016 and LKY PDF 2/08 startup grant.

#131354 - $15.00 USD

(C) 2010 OSA

