Acta Polytechnica doi:10.14311/AP.2017.57.0399 Acta Polytechnica 57(6):399–403, 2017 © Czech Technical University in Prague, 2017 available online at http://ojs.cvut.cz/ojs/index.php/ap LAMBERT FUNCTION METHODS FOR LASER DYNAMICS WITH TIME-DELAYED FEEDBACK Yogesh N. Joglekar∗, Andrew Wilkey, Gautam Vemuri Department of Physics, Indiana University Purdue University Indianapolis (IUPUI), Indianapolis 46202 Indiana, USA ∗ corresponding author: yojoglek@iupui.edu Abstract. Time-delayed differential equations arise frequently in the study of nonlinear dynamics of lasers with optical feedback. Traditionally, one has resorted to numerical methods because the analytical solution of such equations are intractable. In this manuscript, we show that under some conditions, the rate equations model that is used to model semiconductor lasers with feedback can be analytically solved by using the Lambert W-function. In particular, we discuss the conditions under which the coupled rate equations for the intra-cavity electric field and excess carrier inversion can be reduced to a single equation for the field, and how this single rate equation can be cast in a form that is amenable to the use of the Lambert W-function. We conclude the manuscript with a similar discussion for two lasers coupled via time-delayed feedbacks. Keywords: time-delayed differential equations; Lambert function; lasers. 1. Introduction Time-delayed differential equations arise naturally in a wide variety of physical phenomena where one or more system parameters are fed back into the system after a certain amount of time. Such time-delayed feedbacks are seen in population behaviors in biology and ecology [1–3], chemical reactions [4], interactions between time-delayed non-Markovian laser fields and resonant media [5] and a host of nonlinear dynamical systems [6]. One of the more prominent examples of a physical system with time-delayed feedback is a laser wherein the light emitted by a laser is injected back into the laser by reflection from a distant mirror outside the laser cavity [7]. The mathematical model for a time-delayed feedback system often reduces to a first order differential equation with a time-delayed term, and the analytical solution of such differential equations can be difficult because one has to deal with an infinite-dimensional equation. In this article, we demonstrate that the Lambert W-function [8] can be invoked in some situations to obtain analytical solutions to time-delayed equations of physical interest, and explore some of the consequences of using this method. For the sake of concreteness, we focus on the prob- lem of a semiconductor laser that is subject to time- delayed feedback of light into the laser [9]. Lasers with time-delayed feedback are a paradigm for the study of time-delayed systems in part because the delay can be easily controlled, which allows one to study the behavior of the system for delays that are shorter than the intrinsic time-scales of the laser as well as for delays that are longer than the natural time scales of the isolated laser system. Such lasers are of fundamen- tal interest due to the variety of nonlinear dynamical behaviors that arise as a function of the time-delay and the strength of the feedback [10]. In particular, there are combinations of delay and feedback strengths that produce single-tone oscillations in the optical fre- quency of the laser, period doubling routes to chaos, and coherence collapse and line-narrowing. Each of these dynamical responses have been studied for a vari- ety of applications such as the development of stable, all-optical microwave frequency oscillators, chaotic synchronization for all-optical encryption, and stable, narrow line-width lasers [11]. Another system that has been of immense interest to the semiconductor laser community is the coupling of two lasers by mutual injection of light from each laser into the other [13]. These systems have a natural time-delay built into them due to the finite amount of time it takes for the light from one laser to reach the other laser due to the physical separation between the lasers. 2. Lang–Kobayashi equations Semiconductor lasers are usually modeled by the Lang– Kobayashi equations [14] which are known to describe the experimentally observed behavior of these lasers very well. For a single-mode laser, these equations describe the coupled time evolution of the electric field and the excess carrier inversion inside a laser cavity. In the slowly-varying envelope approximation, these equations in the non-dimensional form are dE dt = (1 + iα)ζN(t)E(t) + κE(t− τ), (1) T dN dt = P −N(t) − (1 + 2N(t))|E(t)|2. (2) Here E(t) is the complex, time-dependent, intra-cavity electric field, α is the line-width enhancement factor 399 http://dx.doi.org/10.14311/AP.2017.57.0399 http://ojs.cvut.cz/ojs/index.php/ap Yogesh, Andrew, Gautam Acta Polytechnica 0 20 40 60 80 100 Time (ns) 0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 D im e n si o n le ss e le ct ri c- fie ld in te n si ty P=5x10-4 P=10-3 P=10-2 Laser Power above threshold (P=0) 0 20 40 60 80 100 Time (ns) -10 -8 -6 -4 -2 0 2 4 D im e n si o n le ss e xc e ss c a rr ie r d e n si ty 10-3 Figure 1. Full temporal dynamics of Electric field intensity I(t) = |E(t)|2 and excess carrier inversion N(t) in a typical semiconductor laser with line-width enhancement factor α = 3 and T = 100, as a function of the pump power P above the threshold. Note that both dynamics have the same time scale, and change from overdamped to underdamped in their approach to equilibrium as the pump power increases. for the gain medium, N(t) is the time-dependent ex- cess carrier inversion (above the carrier inversion at the lasing threshold), ζ is the differential-gain coef- ficient, κ is the feedback coupling strength, τ is the time-delay, T is the ratio of excess carrier-inversion lifetime to the cavity photon lifetime, and P is the external pumping to the laser. Note that the rate equa- tion for the macroscopic polarization within the gain medium does not enter this model because it decays very rapidly, relative to the time scale at which E(t) and N(t) evolve in semiconductor lasers, and hence can be adiabatically eliminated. Figure 1 shows the results for the envelope field intensity |E(t)|2 and the excess carrier inversion N(t) as a function of the pump power above the lasing threshold. When the pump power is small, both intensity and excess inversion approach their equilibrium values rapidly in an over- damped manner. As the pump power increases, the approach to the equilibrium values changes to an un- derdamped manner. In dissipative systems, this tran- sition from overdamped to underdamped approach can be mapped onto a parity-time (PT ) symmetry breaking transition [15] where the overdamped region is associated with PT -broken phase and the under- damped region is associated with the PT -symmetric phase. For realistic parameters of standard semicon- ductor lasers, used in Figure 1, we note that the time scales for variation in the electric field and the excess inversion are the same. Thus, it is not possible to elimi- nate the excess carrier density dynamics in the current set up, and the Lambert function formalism cannot be directly applied to above set of equations (1)–(2). The Lambert W-function is defined by solutions of the equation wew = z. For a general complex number z, this equation has countably infinite number of solutions denoted by Wk(z) for integers k, out of which, by convention, only branches k = 0 and k = −1 are real-valued for any z [8]. Suppose we are able to reduce the Lang–Kobayashi equations to a single equation of the form dx dt = ax + bx(t− τ) (3) where x is complex and a,b are constants that could be real or complex. We note that such equations commonly arise in time-delayed population dynamics models in biology and ecology [1–3], but in those cases, the coefficients and the solution of the equation are both constrained to be purely real. Since (3 is a linear equation, its general solution will be given by the linear superposition of exponential-in-time terms. Assuming x(t) ∼ eλt leads to a transcendental characteristic equation (λ−a)τe(λ−a)τ = bτeaτ. (4) Equation 4 shows that the eigenvalues λ can be ex- pressed in terms of the Lambert W-function. Specif- ically, when if bτeτa > 0, the solution for is given by λ = a + W0(bτeaτ )/τ. If bτeaτ < 0, there are two possible solutions when −1/e < bτeaτ < 0, one smaller than 1 and the other larger than 1. Thus, in general, the eigenvalues λ that characterize the exponential-in-time behavior of x(t) are given by λk = a + 1 τ Wk(bτeaτ ) (5) where for bτeaτ < −1/e or complex, the general solu- tion is obtained by an analytical continuation of the function to the complex plane [8]. 3. Domain of validity of Lambert formalism It is generally the case that the modeling of semi- conductor lasers requires the coupled rate equations for the intracavity electric field, E(t) and the excess carrier inversion, N(t). However, the use of the Lam- bert equation formalism, as discussed above, requires 400 vol. 57 no. 6/2017 Lambert Function Methods 0 500 1000 1500 2000 Time (ns) -0.1 -0.05 0 0.05 0.1 0.15 T im e -d e p e n d e n t I( t) a n d N (t ) I(t) N(t) P=0.03, =1 Figure 2. Time evolution of the dimensionless electric field intensity I(t) = |E(t)|2 and the excess carrier inversion N(t) in typical, solitary semiconductor laser shows that the electric field rises exponentially in a very short window before nonlinearities in the gain medium become effective and saturate the intensity to a steady state value. that the two coupled rate equations be reduced to a single equation with time-delay. In this case, it means that the model must be reduced to a rate equation for the electric field only. This is, strictly speaking, not possible for semiconductor lasers because, as seen in Figure 1, the characteristic time scales on which the intracavity field and the excess carrier inversion evolve are of the same order of magnitude. However, with advancements in technology, it may well be possible to fabricate lasers in which the excess carrier inversion evolves much faster than the electric field, in which case the rate equation for the excess carrier inversion can be eliminated. This entails set- ting the time evolution of the inversion to zero, i.e., dN(t)/ dt = 0, solving for the steady state value of N(t), and substituting this expression for N(t) in the rate equation for the electric field (1). The model is thereby reduced to a single, time-delayed rate equation that is amenable to the use of the Lambert function. Another possibility is to fabricate a laser in which the carrier inversion relaxation time is much larger than the electric field decay time. In this case, the inversion does not evolve during the time that the electric field reaches a steady state, and one can, once again, focus on the single, time-delayed rate equation for the electric field. In addition to the elimination of the inversion equa- tion, the Lambert formalism is strictly applicable to the case of a semiconductor laser with time-delayed feedback only when the nonlinearities that arise from the cubic term in the electric field are also ignored. Physically, this means that the formalism is valid when two conditions are simultaneously satisfied - (i) the ex- cess carrier inversion N(t) has reached its steady state value, and (ii) the laser has reaches threshold and the intracavity electric field E(t) is starting to grow expo- nentially, but the intensity saturation that typically sets in due to the nonlinearities in the gain is yet to occur. For typical semiconductor laser parameters (with no time-delayed feedback), this time window between the inversion settling to a steady state value and the intensity getting saturated is negligibly small. Figure 2 shows the typical time evolution of a solitary semiconductor laser (α = 3, pumping P = 0.03, and differential gain coefficient ζ = 1) with initial excess carrier inversion, i.e., N(t = 0) = −0.1, and we see that the exponential electric field rise occurs in a very short time interval. This essentially means that the predictions of the analytic solutions (5), obtained via the Lambert function approach will be visible in a very tiny time window. What is desired however is a wider time window between the inversion reaching a steady state and the intensity of the laser not yet being saturated, i.e., the intensity is still in the exponential amplification regime. Since the intensity at very short times grows exponentially with the product of the differential gain coefficient ζ and N(t), and because the steady-state value of N(t) is a function of the pumping P, one can manipulate these two quantities to slow down the exponential growth of laser intensity and thereby enhance the time window in which both conditions are simultaneously met. Figure 3 shows the time-dependent laser intensities I(t) = |E(t)|2 when the inversion is adiabatically eliminated and only the electric field equation is solved, as a function of differential gain ζ and pump power P . We note that the time axis only contains range after which the carrier inversion has settled down. It is evident from the left-hand panel that for a fixed ζ, as the pump power is reduced from P = 0.09 to P = 0.07, the exponential growth of intensity is slowed, thereby providing a wider window of time in which to make any desired measurements of the laser intensity dynamics. The right-hand panel shows that for a fixed pumping, as ζ is reduced, the intensity growth is slowed. While the external pump is an easily varied parameter in experiments, the differential gain coefficient ζ is a material parameter and hence cannot be tuned in a given laser. On the other hand, reducing the pump power close to zero means the laser is operating very close to the threshold, and that leads to enhanced quantum fluctuations whose effects are not included in the present analysis. Thus, our results provide some guidance on the material parameters that are necessary for a laser to meet such that the Lambert formalism will be applicable to analytically study the dynamics of the laser. Under these conditions, the rate equation for the electric field can be written as dE dt = (1 + iα)ζN0E(t) + κE(t− τ) (6) where N0 ∼ P is the steady-state value of the carrier inversion. This equation is identical to (3) with a 401 Yogesh, Andrew, Gautam Acta Polytechnica 50 60 70 80 90 100 Time (ns) 0 500 1000 1500 2000 2500 D im e n si o n le ss e le ct ri c- fie ld I n te n si ty =0.06 =0.07 =0.08 Laser differential gain 50 60 70 80 90 100 Time (ns) 0 1 2 3 4 5 6 7 D im e n si o n le ss e le ct ri c- fie ld I n te n si ty 104 P=0.08 P=0.07 P=0.09 Laser pump above threshold Figure 3. Rise of the intracavity electric field E(t) when the carrier inversion is adiabatically eliminated. The left-hand panel shows that for a fixed ζ = 0.1, as the pump current is reduced, the exponential growth of intensity slows down. The right-hand panel shows that for a fixed pumping P = 0.08, as ζ is reduced, the time window during which the gain-nonlinearity can be ignored is widened. In both cases, the calculations are carried out using typical semiconductor laser parameters, α = 3. Note that the time axis only contains time window after the carrier inversion has reached steady-state value. manifestly complex a = (1 +iα)ζN0 and a purely real, positive b = κ. Thus, within the appropriate time- window, the electric field exponents are determined by the properties of Lambert W-function. 4. Bidirectionally coupled lasers Another interesting situation in which the Lambert formalism can be invoked is when two identical semi- conductor lasers, at optical frequencies ω1 and ω2 are mutually coupled to each other. Such systems have been extensively studied in the context of their nonlin- ear dynamics [10]. The four rate equations for such a system, two for the intracavity electric fields and two for the corresponding excess carrier inversions, are given by a modified form of the Lang-Kobayashi equa- tions wherein the bidirectional coupling is accounted for, i.e., dE1(t) dt = (1 + iα)ζN1(t)E1(t) + i∆ωE1(t) + κe−iΘτE2(t− τ), (7) dE2(t) dt = (1 + iα)ζN2(t)E2(t) − i∆ωE2(t) + κe−iΘτE1(t− τ), (8) T dN1(t) dt = J1 −N1(t) − (1 + 2N1(t))|E1(t)|2, (9) T dN2(t) dt = J2 −N2(t) − (1 + 2N2(t))|E2(t)|2. (10) Here the subscripts 1, 2 denote laser index,and these equations are written in a frame rotating at a fre- quency that is the average of the two laser frequencies, Θ = (ω1 + ω2)/2, so that each laser is detuned by an equal amount ±∆ω = ±(ω1 − ω2)/2. κ is the coupling coefficient between the two lasers, τ is the time-delay in the coupling that depends on the phys- ical separation between the two lasers, e−iΘτ is the phase accumulation due to light propagating from one laser to the other, and J1,2 are the injection currents above threshold for each laser. If these four coupled equations can be reduced to two by eliminating the carrier inversion equations as discussed for the sin- gle laser with feedback, one is left with two coupled, time-delayed rate equations for the intracavity fields E(t) = [E1(t),E2(t)]T , d dt E(t) = M ( d dt ) E(t) (11) where the 2 × 2 non-Hermitian, time-delay operator M is given by M = [ (1 + iα)ζN10 + i∆ω κe−iΘτe−τ∂t κe−iΘτe−τ∂t (1 + iα)ζN20 − i∆ω ] (12) and N10 and N20 are the steady-state carrier inver- sions. Equations (11) and (12) are also amenable to analytic solution via the Lambert W equation for- malism. An experimental study of the bidirectionally coupled lasers, and a detailed comparison between the predictions of the Lambert function solution and nu- merical solutions of the full Lang–Kobayashi equations will be reported elsewhere. In summary, we have shown that the Lambert W- function provides an hitherto unexplored, analytic method for studying the intensity dynamics of a semi- conductor laser with time-delayed optical feedback. The formalism is valid in a regime where the two rate equations given by the Lang–Kobayashi model are reduced to a single, time-delayed rate equation for the intracavity electric field. The Lambert function can be invoked when the nonlinearities that arise from gain saturation are neglected, which implies that the analytic results are valid at short times after the laser 402 vol. 57 no. 6/2017 Lambert Function Methods intensity crosses its threshold value, i.e., when the intensity is still in the exponentially amplifying stage. Furthermore, the analytic technique assumes that the carrier inversion has reached its steady state value while the intensity is till growing. To overcome the problem of a very narrow observation time window for the laser intensity dynamics, we have suggested some remedies that could be implemented at the laser fabrication step to modify the material parameters of the laser. In particular, reducing the differential gain coefficient, and modifying other properties of the laser so that a very weak pump will induce the desired population inversion, will enable a much wider time window between the time at which the inver- sion reaches a steady state and the time at which the laser intensity saturates. Once a wider time window is attained, the predictions of the Lambert function results can be tested. References [1] See e.g. H. Smith, An introduction to delay differental equations with applications to the life sciences, ISBN 978-1-4419-7646-8, Springer (2011). [2] S.P. Blythe, R. M. Nisbet, and W. S. C. Gurney, Instability and complex dynamic behavior in population models with long time delays. Theoretical Population Biology 22, 147âĂŞ176 (1982). [3] Y. Kuang, Delay differential equation with applications in population dynamics. Academic Press, Inc. (1993). [4] Marc R. Roussel, The Use of Delay Differential Equations in Chemical Kinetics, J. Phys. Chem. 100, 8323-8330 (1996). [5] M.H. Anderson, G. Vemuri, J. Cooper, P. Zoller and S.J. Smith, Experimental Study of Absorption and Gain from Two-Level Atoms in Non-Markovian, Phase Diffusing Optical Fields, Phys. Rev. A 47, 3202-3209 (1993). [6] B. Krauskopf and D. Lenstra, Fundamental issues of nonlinear laser dynamics, AIP Conference Proc. 548 (2000). [7] See e.g. D. Kane and K.A. Shore, Unlocking dynamical diversity: Optical feedback effects on semiconductor lasers, ISBN 0-470-85619-X. [8] R.M. Corless, G.H. Gonnet, D.E.G. Hare, D.J. Jeffery, and D.E. Knuth, On the Lambert W function, Adv Comput Math 5, 329 (1996). [9] G.H.M. van Tartwijk and D. Lenstra, Quant. and Semiclassical Optics 7, 87 (1995). [10] J. Mork, J. Mark and B. Tromborg, Route to chaos and competition between relaxation oscillations for a semiconductor laser with optical feedback, Phys. Rev. A 65, 1999-2002 (1990). [11] M. Soriano, J. Garcia-Ojalvo, C.R. Mirasso, and I. Fischer, Rev. Mod. Phys. 85, 421-470 (2013). [12] M. Sciamanna and K.A. Shore, Nature Photonics 9, 151-162 (2015). [13] J. Mulet, C. Masoller and C.R. Mirasso, Modeling bidirectionally coupled single-mode semiconductor lasers, Phys. Re. A 65, 063815 (2002). [14] R. Lang and K. Kobayashi, IEEE J. Quant. Elec. QE-16, 347-355 (1980). [15] J. Li, A.K. Harter, L. de Melo, Y.N. Joglekar, and L. Luo, arXiv:1608.05061. 403 Acta Polytechnica 57(6):399–403, 2017 1 Introduction 2 Lang–Kobayashi equations 3 Domain of validity of Lambert formalism 4 Bidirectionally coupled lasers References