Aug 20, 2014 - inj < Pon inj, i.e. a hysteresis region exists with an additional LS either turned on or off. The two ..... mobile cavity solitons, Phy...

0 downloads 10 Views 682KB Size

arXiv:1408.4667v1 [nlin.PS] 20 Aug 2014

1

We investigate a control of the motion of localized structures of light by means of delay feedback in the transverse section of a broad area nonlinear optical system. The delayed feedback is found to induce a spontaneous motion of a solitary localized structure that is stationary and stable in the absence of feedback. We focus our analysis on an experimentally relevant system namely the Vertical-Cavity Surface-Emitting Laser (VCSEL). In the absence of the delay feedback we present experimental evidence of stationary localized structures in a 80 µm aperture VCSEL. The spontaneous formation of localized structures takes place above the lasing threshold and under optical injection. Then, we consider the effect of the time-delayed optical feedback and investigate analytically the role of the phase of the feedback and the carrier lifetime on the self-mobility properties of the localized structures. We show that these two parameters affect strongly the space time dynamics of two-dimensional localized structures. We derive an analytical formula for the threshold associated with drift instability of localized structures and a normal form equation describing the slow time evolution of the speed of the moving structure. PACS numbers:

I.

INTRODUCTION

Transverse localized structures often called cavity solitons were observed experimentally in Vertical-Cavity SurfaceEmitting Laser (VCSEL) [1, 2]. Due to the maturity of the semiconductor technology and the possible applications of localized structures of light in all-optical delay lines [3] and logic gates [4], these structures in VCSELs have been a subject of an active research in the field of nonlinear optics. Moreover, the fast response time of VCSELs makes them attractive devices for potential applications in all-optical control of light. Localized structures appear as solitary peaks or dips on homogeneous background of the field emitted by a nonlinear microresonator with a high Fresnel number. Here, we report experimental evidence of spontaneous formation of localized structures in a 80 µ m diameter VCSEL biased above the lasing threshold and under optical injection. Such localized structures are bistable with the injected beam power and the VCSEL current. Various mechanisms have proven to be responsible for the generation of localized structures (LS) in VCSELs: coherent optical injection (holding beam, i.e. the part of the optical injection that is used to ensure bistability of the system) in combination with a narrow writing beam (the part of optical injection which is used to locally make the system evolve from the lower branch of the bistability curve to the higher branch and vice versa)[2], frequency selective feedback with [5] or without [6] a writing beam, saturable absorption [7], amplifiers [8], optically injected lasers [9, 10], “spatial translational coupling” introduced in [11, 12], and others (see [13] for a review). LS can start to move spontaneously in the laser transverse section in the presence of saturable absorption [17]. In particular, in the case when the pump beam has a circular shape LS move along the boundary on a circular trajectory [14]. It has been shown that they can undergo a spontaneous motion under the thermal effects [15, 16]. Delayed feedback control is a well documented technique that has been applied to various spatially extended systems in optics, hydrodynamics, chemistry, and biology. It has been demonstrated recently that a simple feedback loop provides a robust and a controllable mechanism for the motion of localized structures (LS) and localized patterns [18–24]. These works demonstrated that when the product of the delay time and the feedback rate exceeds some threshold value, localized structures start to move in an arbitrary direction in the transverse section of the device. In these studies, the analysis was restricted to the specific case of nascent optical bistability described by the real Swift-Hohenberg equation with a real feedback term. More recently, analytical study supported by numerical simulations revealed the role of the phase of the delayed feedback and the carrier lifetime on the motion of cavity solitons in a broad-area VCSEL structure, driven by a coherent externally injected beam. We show that for certain values of the feedback phase LS can be destabilized via a drift bifurcation leading to a spontaneous motion of a solitary two-dimensional LS [18]. We

2 demonstrate also that the slower is the carrier decay rate in the semiconductor medium, the higher is the threshold associated with the motion of cavity solitons. The paper is organized as follows: In Sec. 2, we discuss the experimental observation of stationary LS in medium size VCELS. In Sec. 3, we introduce model equations. In Sec. 4, we investigate the drift instability induced by delay feedback. We conclude in Sec. 5. II.

EXPERIMENTAL OBSERVATION OF STATIONARY LS IN MEDIUM SIZE VCELS

In recent years, a considerable amount of experimental work has been realized on stationary localized structures in Vertical-Cavity Surface-Emitting Lasers (VCSELs). They were observed in very broad (aperture d > 100µm) [2], broad (d ∼ 80µm) [25], and medium (d ∼ 40µm) [26] size VCSELs. Here we present experimental results obtained with bottom-emitting InGaAs multiple quantum well VCSEL with d = 80µm and threshold current of 42.5 mA at 20◦ C. The holding (injection) beam is provided by a commercial tunable semiconductor laser (Sacher Lasertechnik TEC100-0960-60 External Cavity Diode Laser), isolated from the rest of the setup by an optical isolator (OFR IO5TiS2-HP). The long-term electrical and temperature stability of this laser are less than 20 mA RMS and 0.05 ◦ C, respectively. A half-wave plate is used to adjust the linear polarization of the holding beam to be the same as the one of the VCSEL. The injection beam power is tuned using a variable optical density filter. The detuning between the master laser and the VCSEL is defined as θ = νinj − νslave , where νinj is the frequency of the injection beam, and νslave the frequency of the strongest peak in the spectrum of the stand-alone VCSEL. It is experimentally tuned by changing the wavelength of the injection beam. The beam waist dinj is defined as the diameter of the smallest circle in the plane of propagation of the injection beam containing half of the beam power when it encounters the VCSEL. The power of the source is monitored by a Newport 818-SL photodiode connected to a Newport 2832-C powermeter. The near field is recorded by imaging it on a CCD camera. An example of stationary LSs is presented in the Fig. 1 illustrating the process of spontaneous creation and annihilation of two-dimensional localized structures. These experimental measurements have been performed when the VCSEL operated in an injection locked regime. The injection beam waist and the detuning are fixed to dinj = 50µm on and θ= -146 GHz. When increasing the injected beam power, a new LS appears at Pinj = Pinj as shown in the insets of Fig. 1. This results in a slight jump of the total output power as shown in light-versus-current characteristics. The process of switching-off is realized when decreasing the injection power: the recently created LS persists until of f of f on Pinj = Pinj with Pinj < Pinj , i.e. a hysteresis region exists with an additional LS either turned on or off. The two figures on the right show one dimensional scans along the vertical lines indicated in the near field images.

off

on

Figure 1: Total output power as a function of injection power displaying bistability when a new LS appears. The insets show the near field images of on the upper and the lower branch of the hysteresis curve. The two figures on the right show one dimensional scans along the vertical lines indicated in the near field images.

3 III.

MODEL EQUATIONS

The mean field model describing the space-time evolution of the electric field envelope E and the carrier density N in a VCSEL subjected to optical injection is given by the following set of dimensionless partial differential equations ∂E = − (µ + iθ) E + 2C(1 − iα)(N − 1)E ∂t + Ei − ηeiϕ E(t − τ ) + i∇2 E , h i ∂N = −γ N − I + (N − 1) |E|2 − d∇2 N . ∂t

(1)

(2)

Here the parameter α describes the linewidth enhancement factor, µ and θ are the cavity decay rate and the cavity detuning parameter, respectively. Below we will assume η to be small enough, so that we can neglect the dependence of the parameters µ and θ on ϕ. The parameter Ei is the amplitude of the injected field, C is the bistability parameter, γ is the carrier decay rate, I is the injection current, and d is the carrier diffusion coefficient. The diffraction of light and the diffusion of the carrier density are described by the terms i∇2 E and d∇2 N , respectively, where ∇2 is the Laplace operator acting in the transverse plane (x, y). Below we consider the case when the laser is subjected to coherent delayed feedback from an external mirror. To minimize the effect of diffraction on the feedback field we assume that the external cavity is self-imaging [18]. The feedback is characterized by the delay time τ = 2Lext /c, the feedback rate η ≥ 0, and phase ϕ, where Lext is the external cavity length, and c is the speed of light. The link between dimensionless and physical parameters is provided in [19]. Using the expression for the feedback rate 1/2

(1−R) η = r R1/2 given in [27], where r (R) is the power reflectivity of the feedback (VCSEL top) mirror and τin is the τin VCSEL cavity round trip time, we see that the necessary condition for the appearance of the soliton drift instability Rτ 2 ητ > 1 [18] can be rewritten in the form r > (1−R)in2 τ 2 . In particular, for R = 0.3 and τ = 20τin the latter inequality becomes r > 1.5 · 10−3 .

IV.

DRIFT INSTABILITY THRESHOLD

When the delayed feedback is absent, η = 0, Eqs. (1) and (2) are transformed into the well-known mean field model [28], which supports stable stationary patterns and LSs [2, 26, 29, 30]. It was demonstrated recently that when the feedback rate η exceeds a certain threshold value, which is inversely proportional to the delay time τ , LS starts to move in the transverse direction [18]. Example of moving two-dimensional LS are shown in Fig. 2. The single and the three moving peaks are obtained from numerical simulations of Eqs. (1) and (2). The boundary condition are periodic in both transverse dimensions.

a)

t= 12 000

20 40

b)

20 40

60

60

t= 10 400

80

80

100

100

120

120 20

40

60

80

100

120

20

40

60

80

100

120

Figure 2: Field intensity illustrating a moving a single (a) and three (b) peak LS. Parameter values are C=0.45, θ = −2., α = 5., µ = 1., Feedback parameters are η = 0.135, τ = 100, ϕ = 0.5. Maxima are plain white.

In the case when the system is transversely isotropic, the velocity of the LS motion has an arbitrary direction. The self-induced motion of the LS is associated with a pitchfork bifurcation where the stationary LS looses stability and a branch of stable LSs uniformly moving with the velocity v = |v| bifurcates from the stationary LS branch. The bifurcation point can be obtained from the first order expansion of the uniformly moving LS in power series of the

4 small velocity v. Close to the pitchfork bifurcation point this expansion reads: E(x − vt, y) = E0 (x − vt, y) + vE1 (x − vt, y) + ... N (x − vt, y) = N0 (x − vt, y) + vN1 (x − vt, y) + ...,

(3) (4)

where without the loss of generality we assume that the LS moves along the x-axis on the (x, y)-plane. Here E0 (x, y) = X0 (x, y) + iY0 (x, y) and N0 (x, y) describes the stationary axially symmetric LS profile, which corresponds to the time-independent solution of Eqs. (1) and (2) with τ = 0. Although formally this solution depends on the feedback parameters η and ϕ we neglect this dependence assuming that the feedback rate is sufficiently small, η ≪ 1. Substituting this expansion into Eqs. (1) and (2) and collecting the first order terms in small parameter v we obtain: ℜ[∂x E0 (1 − ητ eiϕ )] ℜE1 L ℑE1 = ℑ[∂x E0 (1 − ητ eiϕ )] N1 γ −1 ∂x N0

(5)

where the linear operator L is given by µ − 2C(N0 − 1) ∇2 − θ − 2C α (N0 − 1) −2C(X0 + αY0 ) µ − 2C(N0 − 1) −2C(Y0 − αX0 ) . L = −∇2 + θ + 2C α (N0 − 1) 2(N0 − 1)X0 2(N0 − 1)Y0 −d∇2 + 1 + |E0 |2

By applying the solvability condition to the right hand side of Eq. (3), we obtain the drift instability threshold 1 + γ −1 (b/c) ητ = p 1 + (a/c)2 cos[ϕ + arctan (a/c)]

(6)

with a = hψ1† , ψ2 i − hψ2† , ψ1 i,

b = hψ3† , ψ3 i,

c = hψ1† , ψ1 i + hψ2† , ψ2 i.

(7)

Here T

T

ψ = (ψ1 , ψ2 , ψ3 ) = ∂x (X0 , Y0 , N0 )

(8)

T is the corresponding solution is a translational neutral mode of the operator L, Lψ = 0, while ψ † = ψ1† , ψ2† , ψ3† R +∞ † † of the homogeneous adjoint problem L ψ = 0. The scalar product h·i is defined as hψj† , ψk i = −∞ ψj† ψk dxdy. To estimate the coefficients a and b we have calculated the function ψ † numerically using the relaxation method in two transverse dimensions, (x, y). The results of these calculations are shown in Fig. 3 together with the axially symmetric profile E0 of the stationary LS. It is seen from this figure that similarly to the neutral mode ψ defined by (8) the neutral mode ψ † of the adjoint operator L† is an even function of the coordinate y and an odd function of the coordinate x, which is parallel to the LS direction of motion. The dependence of the critical feedback rate η corresponding to the drift instability threshold defined by Eq. (6) on the feedback phase ϕ and carrier relaxation rate γ is illustrated by Fig. 4. In this figure the curves labeled by different numbers correspond to different values of γ. Considering the fact that the feedback in Eq. (1) is introduced with the minus sign, we see that the drift instability takes place only for those feedback phases when the interference between the cavity field and the feedback field is destructive, i.e. when cos function in the denominator of the right hand side of Eq. (6) is positive. On the contrary, when this interference is constructive the feedback has a stabilizing effect on the LS. Furthermore, the slower is the carrier relaxation rate, the higher is the drift instability threshold. Since the stationary LS solution does not depend on the carrier relaxation rate γ, the coefficients a and b in the threshold condition (6) are also independent of γ. Therefore, (6) gives an explicit dependence of the threshold feedback rate on the carrier relaxation rate. In particular, in the limit of very fast carrier response, γ ≫ 1, and zero feedback phase, ϕ = 0, we recover from (6) the threshold condition ητ = 1 which was obtained earlier for the LS drift instability induced by a delayed feedback in the real Swift-Hohenberg equation [18]. Note that at γ → ∞, a 6= 0, and ϕ = − arctan a the critical feedback rate appears to be smaller than that obtained for the real Swift-Hohenberg −1/2 equation, ητ = (1 + a2 ) < 1.

5

Figure 3: Left panels: real and imaginary parts of the stationary soliton profile, X0 = ℜE0 (a) Y0 = ℑE0 (c). Right panels: real and imaginary parts of the neutral mode of the adjoint operator L† , ψ1† = ℜψ † (b) and ψ2† = ℑψ † (d). Parameters values: µ = 1.0, θ = −2.0, C = 0.45, α = 5.0, γ = 0.05, τ = 100, d = 0.052, Ei = 0.8, I = 2.

Figure 4: Critical value of the feedback rate η corresponding to the drift bifurcation vs feedback phase ϕ calculated for different values of the carrier relaxation rate γ. The values of the parameter γ are shown in the figure. Other parameters are the same as in Fig. 3.

As it was demonstrated above, the bifurcation threshold responsible for self-induced drift of LS in the VCSEL transverse section is obtained by expanding the slowly moving localized solution in the small velocity v, substituting this expansion into the model equations (1),(2), and equating the first order terms in v. In order to describe the slow evolution of the LS velocity a similar procedure with P3 slightly above the bifurcation threshold one needs to perform P3 E = E0 (x − x0 (t), y) + k=1 ǫk Ek (x − x0 (t), y, t) + ... and N = N0 (x − x0 (t), y) + k=1 ǫk Nk (x − x0 (t), y, t) + ..., where dx/dt = v(t) = O(ǫ), dv/dt = O(ǫ3 ) and ǫ is a small parameter characterizing the distance from the bifurcation point. Then, omitting detailed calculations, in the third order in ǫ, we obtain the normal form equation for the LS velocity: p dv = v(δηq − ητ 2 rv 2 ), 2 dt

(9)

6 where δη is the deviation of the feedback rate from the bifurcation point. The coefficients q, p, and r are given by q = a sin ϕ + c cos ϕ, p = q + b, and r = f sin ϕ + g cos ϕ + O(τ −1 ), respectively. Here a, b, and c are defined by Eq. (7) and f = hψ1† , ∂xxx Y0 i − hψ2† , ∂xxx X0 i, h = hψ3† , ∂xxx N0 i, g = hψ1† , ∂xxx X0 i + hψ2† , ∂xxx Y0 i. The stationary √ LS velocity above the drift instabilitypthreshold is obtained by calculating the nontrivial steady state of Eq. (9), v = δηQ, where the coefficient Q = (1/τ ) q/(rη) determines how fast the LS velocity increases with the square root of the deviation from the critical feedback rate. The dependence of this coefficient on the feedback phase is illustrated by Fig. 5.

Figure 5: Coefficient Q describing the grows rate of the LS velocity with the square root of the deviation from the critical feedback rate. The values of the parameter γ are shown in the figure. Other parameters are the same as in Fig. 3.

V.

CONCLUSION

To conclude, we have studied experimentally the formation of transverse LS in a medium size bottom-emitting InGaAs multiple quantum well VCSEL operated in an injection locked regime. Creation and annihilation of a single localized structure have been demonstrated by changing the injection beam power. Theoretically we have analyzed the effect of delayed feedback from an external mirror on the stability of transverse LS in a broad area VCSEL. We have shown that depending on the phase of the feedback it can have either destabilizing or stabilizing effect in the LS. In particular, when the interference between the LS field and the feedback field is destructive, the LS can be destabilized via a pitchfork bifurcation, where a branch of uniformly moving LS bifurcates from this the stationary one. We have calculated analytically the threshold value of the feedback rate corresponding to this bifurcation and demonstrated that the faster the carrier relaxation rate in the semiconductor medium, the lower is the threshold of the spontaneous drift instability induced by the feedback. Finally, we have derived the normal for equation (9) governing the slow dynamics of the LS velocity. Acknowledgment

A.G.V. and A.P. acknowledge the support from SFB 787 of the DFG. A.G.V. acknowledges the support of the EU FP7 ITN PROPHET and E.T.S. Walton Visitors Award of the Science Foundation Ireland. M.T. received support from the Fonds National de la Recherche Scientifique (Belgium). This research was supported in part by the Interuniversity Attraction Poles program of the Belgian Science Policy Office under Grant No. IAP P7-35.

[1] V. B. Taranenko, I. Ganne, R. J. Kuszelewicz, and C. O. Weiss, Patterns and localized structures in bistable semiconductor resonators, Phys. Rev. A 61, 063818 (2000). [2] S. Barland, J. R. Tredicce, M. Brambilla, L. A. Lugiato, S. Balle, M. Giudici, T. Maggipinto, L. Spinelli, G. Tissoni, T. Knodl, Cavity solitons as pixels in semiconductor microcavities, Nature 419, 699 (2002).

7 [3] F. Pedaci, S. Barland, E. Caboche, P. Genevet, M. Giudici, J. R. Tredicce, T. Ackemann, A. J. Scroggie, W. J. Firth, G. L. Oppo, G. Tissoni, and R. Jager, All-optical delay line using semiconductor cavity solitons, Appl. Phys. Lett. 92, 011101 (2008). [4] A. Jacobo, D. Gomila, M. A. Mat´ aas, and P. Colet, Logical operations with localized structures, New J. Phys. 14 , 013040 (2012). [5] Y. Tanguy, T. Ackemann, W. J. Firth, and R. Jager, Realization of a semiconductor-based cavity soliton laser, Phys. Rev. Lett. 100, 013907 (2008). [6] N. Radwell and T. Ackemann, Characteristics of laser cavity solitons in a vertical-cavity surface-emitting laser with feedback from a volume bragg grating, IEEE J. Quantum Electron. 45, 1388 (2009). [7] P. Genevet, S. Barland, M. Giudici, and J. R. Tredicce, Cavity soliton laser based on mutually coupled semiconductor microresonators, Phys. Rev. Lett. 101 , 123905 (2008). [8] Y. Menesguen, S. Barbay, X. Hachair, L. Leroy, I. Sagnes, and R. Kuszelewicz, Optical self-organization and cavity solitons in optically pumped semiconductor microresonators, Phys. Rev. A 74, 023818 (2006). [9] Y. Larionova and C. O. Weiss, Spatial semiconductor resonator solitons with optical pumping, Opt. Express 13, 10711 (2005). [10] X. Hachair, F. Pedaci, E. Caboche, S. Barland, M. Giudici, J. R. Tredicce, F. Prati, G. Tissoni, R. Kheradmand,L. A. Lugiato, I. Protsenko, and M. Brambilla, Cavity solitons in a driven vcsel above threshold, IEEE J. Sel. Top. Quantum Electron. 12, 339 (2006). [11] A. F. Haudin, R. G. Rojas, U. Bortolozzo, M. G. Clerc, and S. Residori , Vortex emission accompanies the advection of optical localized structures , Phys. Rev. Lett. 106, 063901 (2011). [12] F. del Campo, F. Haudin, R. G. Rojas, U. Bortolozzo, M. G. Clerc, and S. Residori, Effects of translational coupling on dissipative localized states, Phys. Rev. E 86, 036201 (2012). [13] S. Barbay, R. Kuszelewicz, and J. R. Tredicce, Cavity solitons in vcsel devices, Adv. Opt. Technol. 2011, 628761 (2011) [14] F. Prati, G. Tissoni, L. A. Lugiato, K. M. Aghdami, and M. Brambilla, Spontaneously moving solitons in a cavity soliton laser with circular section, Eur. Phys. J. D 59, 73-79 (2010). [15] L. Spinelli, G. Tissoni, L. A. Lugiato, and M. Brambilla, Thermal effects and transverse structures in semiconductor microcavities with population inversion, Phys. Rev. A 66 , 023817 (2002). [16] A. J. Scroggie, J. M. McSloy, and W. J. Firth, Self-propelled cavity solitons in semiconductor microcavities, Phys. Rev. E 66, 036607 (2002). [17] S. V. Fedorov, A. G. Vladimirov, G. V. Khodova, and N. N. Rosanov, Effect of frequency detunings and finit relaxation rates on laser localized structures, Phys. Rev. E 61, 5814 (2000). [18] M. Tlidi, A. G. Vladimirov, D. Pieroux, and D. Turaev, Spontaneous motion of cavity solitons induced by a delayed feedback, Phys. Rev. Lett. 103, 103904 (2009). [19] K. Panajotov and M. Tlidi, Spontaneous motion of cavity solitons in vertical-cavity lasers subject to opticalinjection and to delayed feedback, Eur. Phys. J. D 59, 67 (2010). [20] M. Tlidi, E. Averlant, A. Vladimirov, and K. Panajotov, Delay feedback induces a spontaneous motion of two- dimensional cavity solitons in driven semiconductor microcavities, Phys. Rev. A 86, 033822 (2012). [21] S. V. Gurevich and R. Friedrich, Instabilities of localized structures in dissipative systems with delayed feed-back, Phys. Rev. Lett. 110, 014101 (2013). [22] M. Tlidi, A. Sonnino, and G. Sonnino, Delayed feedback induces motion of localized spots in reaction-diffusion systems, Phys. Rev. E 87, 042918 (2013). [23] S. V. Gurevich, Dynamics of localized structures in reaction-diffusion systems induced by delayed feedback, Phys. Rev. E 87, 052922 (2013). [24] A. Pimenov, A. G. Vladimirov, S. V. Gurevich, K. Panajotov, G. Huyet, and M. Tlidi, Delayed feedback control of selfmobile cavity solitons, Phys. Rev. A 88, 053830 (2013). [25] E. Averlant, M. Tlidi, H. Thienpont, T. Ackemann, and K. Panajotov, Experimental observation of localized structures in medium size VCSELs, Opt. Express 22, 762 (2014). [26] X. Hachair, G. Tissoni, H. Thienpont, and K. Panajotov, Linearly polarized bistable localized structure inmedium-size vertical-cavity surface-emitting lasers, Phys. Rev. A 79, 011801 (2009). [27] G. H. M. van Tartwijk and D. Lenstra, Semiconductor lasers with optical injection and feedback, J. Opt. B Quantum Semiclassical Opt. 7, 87-144 (1995) [28] L. Spinelli, G. Tissoni, M. Brambilla, F. Prati, and L. A. Lugiato, Spatial solitons in semiconductor microcavities, Phys. Rev. A 79, 2542 (1998). [29] X. Hachair, S. Barland, L. Furfaro, M. Giudici, S. Balle, J.R. Tredicce, M. Brambilla, T. Maggipinto, I.M. Perrini, G. Tissoni, and L. Lugiato, Cavity solitons in broad-area vertical-cavity surface-emitting lasers below threshold, Phys. Rev. A 69, 043817 (2004). [30] X. Hachair, L. Furfaro, J. Javaloyes, M. Giudici, S. Balle, J. Tredicce, G. Tissoni, L.A. Lugiato, M.Brambilla, and T. Maggipinto, Cavity-solitons switching in semiconductor microcavities, Phys. Rev. A 72, 013815 (2005).