**Preface to "Semiconductor Laser Dynamics"**

It is my great pleasure to publish this book. All contents were peer-reviewed by multiple referees and published as papers in the Special Issue "Semiconductor Laser Dynamics: Fundamentals and Applications" in the journal *Photonics*.

These studies provide new and interesting results in different branches of semiconductor laser dynamics, dealing with the dynamics and stability of semiconductor lasers in a broad sense. This book offers a small window with a view of the present interests and developments in this lively field, which forms a fertile ground for innovative ideas.

> **Daan Lenstra** *Editor*

### *Editorial* **Special Issue "Semiconductor Laser Dynamics: Fundamentals and Applications"**

#### **Daan Lenstra**

Institute of Photonic Integration, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands; dlenstra@tue.nl; Tel.: +31-648-875-241

Received: 28 May 2020; Accepted: 10 June 2020; Published: 11 June 2020

**Abstract:** With the advent of integrated photonics, a crucial role is played by semiconductor diode lasers (SDLs) as coherent light sources. Old paradigms of semiconductor laser dynamics, like optical injection, external feedback and the coupling of lasers, regained relevance when SDLs were integrated on photonic chips. This Special Issue presents a collection of seven invited feature papers and 11 contributed papers reporting on recent advances in semiconductor laser dynamics.

**Keywords:** semiconductor laser; dynamics and stability; laser coupling; integrated lasers

#### **1. Introduction**

As one of the most widely used coherent light sources today, the semiconductor laser is an essential component of many optical systems, notably for communication, storage, sensing and metrological applications but nowadays mainly as parts of photonic integrated systems. They can be linear Fabry–Pérot or ring-type lasers, operating in narrow linewidth, single frequency or pulsed. Their numerous applications are ever increasing due to the unprecedented fabrication accuracy and reproducibility offered by photonic integration technology, allowing total control of the phase and intensity of the generated laser light. Many of these applications involve the nonlinear dynamics of the coupled photon inversion system in one way or another. We mention lasers for the generation of micro-waves or short mode-locked pulses and lasers for the generation of chaotic light in encrypted communication, as well as linewidth narrowing and frequency stabilization by external optical feedback and increased modulation bandwidth by optical injection.

In the well-defined embedded setting of integrated lasers, the issues of reproducibility and long-term dynamical stability are becoming ever more important and should be considered in the design and fabrication of such laser systems. Since precise control of quantities like optical distance, group velocity, wave-guide loss, gain and many other relevant parameters is very feasible, knowledge of the dynamical behaviour of semiconductor lasers in their dependence on parameter values can be successfully incorporated into the optimal design of these lasers and laser systems.

This Special Issue presents a collection of original state-of-the-art research articles dealing with the dynamics and stability of semiconductor lasers in a broad sense, sometimes with special emphasis on their operation in a photonic chip. Specifically, this issue comprises 18 papers dealing with semiconductor lasers coupled to various kinds of optical perturbations, such as delayed feedback, delayed coupling and optical injection, etc. Among these papers, seven are invited "feature" papers on the highly topical subjects of coupled lasers, reservoir computing, injection locking, external optical feedback and very narrow linewidth lasers. The feature papers are reviewed in Section 2 and the contributed papers in Section 3.

#### **2. Feature Papers**

A long-standing and central problem in semiconductor laser dynamics (SLD) is the influence of external delayed optical feedback [1]. This is the situation in which part of the output laser light is reflected from an external reflector and coupled back into the laser. The paper by A. Locquet [2] reviews various aspects of the routes to chaos that can occur under these circumstances. One important application of delayed optical feedback is found in reservoir computing [3], and the task-independent computational abilities are the subject of the paper by Harkhoe and Van der Sande [4]. The review paper by Boller et al. [5] presents an overview of their record-breaking results on linewidth narrowing in hybrid-integrated diode lasers with feedback from low-loss silicon nitride circuits.

Another equally important and often encountered problem in SLD concerns the semiconductor laser with optical injection, usually from another laser. The slave laser may exhibit a large variety of dynamical features; for example, frequency locking to the injected signal, micro-wave oscillations, chaos and excitability [6]. The invited paper by Torre and Masoller [7] explores the combined effects of excitability and the emission of extreme pulses with promising applications to sensing. A problem which is intimately related to laser injection is laser coupling, that is, where each laser injects light into the other at the same time. The feature article by Perrott et al. [8] compares the cases of true injection and pure mutual coupling between semiconductor diode lasers in one photonic integrated circuit. The observed additional types of dynamics in the case of mutual coupling are general features of coupled lasers, which are studied in the invited paper by Erneux and Lenstra [9]. In the latter article, the synchronization of mutually delay-coupled quantum-cascade lasers with different pump strengths is theoretically analyzed. In all the above-mentioned cases of coupled lasers, the coupling was typically face-to-face. A different type of coupling is treated in the feature paper of Vaughan et al. [10], in which the dynamical behavior of two laterally coupled semiconductor lasers is theoretically analyzed.

#### **3. Contributed Papers**

The contributed papers reflect the importance of optical injection and feedback as the generic fundamental processes in semiconductor laser systems. The paper by Sortiss et al. [11] describes the use of injection locking for side-mode suppression with the application to optical communication in general and optical demultiplexing in particular. Jiang et al. [12] numerically investigate the dynamical properties of excited-state emitting quantum-dot lasers with optical injection.

In the numerical study by Ebisawa and Komatsu [13], an ingenious combination of three diode lasers with optical injection and feedback is investigated in order to quantify the orbital instability of the produced chaotic dynamics in terms of Lyapunov exponents. Jayaprasath et al. [14] numerically investigate the properties of the chaotic output light that is produced by a semiconductor laser with delayed external optical feedback, with consequences for the security of chaotic communication. The security theme is also addressed in the numerical study by Wang et al. [15], who consider the risk of the bias current as a key for secure communication.

Using the technology described in the invited paper by Ref. [5], the generation of tunable microwave oscillations by optical sideband injection is described in a paper by Khan and Hoque [16]. Microwave generation is also the theme of the paper by Qi et al. [17], in which a monolithically integrated laser-photodetector chip was designed and fabricated.

An interesting problem is external feedback in a ring laser since the feedback light from a clockwise mode will couple into the counterclockwise mode. The optical-feedback sensitivity of such a laser is studied, experimentally and numerically, by Verschaffelt et al. [18] by applying on-chip filtered optical feedback. The article by Zhang et al. [19] presents the design and performance of a compact, highly stable, external-cavity diode laser for use in an optical clock in space.

Vertical-cavity surface-emitting lasers (VCSELs) are well-suited for high-speed data communication. In the paper by Sanayeh et al. [20], an equivalent circuit model is presented that accurately describes the dynamic behavior of high-performance VCSELs and applies this to a simulation of their intrinsic modulation response. The article by Wilkey et al. [21] addresses the fundamental problem of whether a pair of coupled semiconductor lasers could possess Parity-Time (PT) symmetry. Based on a rate-equation model, they predict intensity dynamics like those in a PT-symmetric system.

#### **4. Outlook and Prospective Further Developments**

The collection of papers in this Special Issue on semiconductor dynamics offers only a small window with a view on the present interests and developments. The field is very much alive and forms a fertile ground for innovative ideas, of which we have seen a few examples only. Promising novel developments are to be expected for applications in the sensing of PT-symmetric photonic systems with exceptional points of operation [22], in photonic neural networks [23] and excitable laser systems [24], in the metrology of super-stable mode-locked pulse lasers and frequency combs [25] and in the search for feedback-resistant lasers [26] and integrated non-reciprocal devices [27].

**Funding:** This research received no external funding.

**Acknowledgments:** The author acknowledges the assistance from the editorial office of Photonics during the preparation of the special issue.

**Conflicts of Interest:** The author declares no conflict of interest.

#### **References**


© 2020 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

### *Review* **Routes to Chaos of a Semiconductor Laser Subjected to External Optical Feedback: A Review**

#### **Alexandre Locquet 1,2**


Received: 6 January 2020; Accepted: 27 February 2020; Published: 5 March 2020

**Abstract:** This paper reviews experimental investigations of the route to chaos of a semiconductor laser subjected to optical feedback from a distant reflector. When the laser is biased close to threshold, as the feedback strength is increased, an alternation between stable continuous wave (CW) behavior and irregular, chaotic fluctuations, involving numerous external-cavity modes, is observed. CW operation occurs on an external-cavity mode whose optical frequency is significantly lower than that of the solitary laser. The scenario is significantly different for larger currents as the feedback level is increased. At low feedback, the laser displays periodic or quasiperiodic behavior, mostly around external-cavity modes whose frequency is slightly larger than that of the solitary laser. As the feedback level increases, the RF and optical frequencies involved progressively lock until complete locking is achieved in a mixed external-cavity mode state. In this regime, the optical intensity and voltage oscillate at a frequency that is also equal to the optical frequency spacing between the modes participating in the dynamics. For even higher feedback, the locking cannot be maintained and the laser displays fully developed coherence collapse.

**Keywords:** semiconductor laser; optical feedback; nonlinear dynamics; bifurcations; chaos

#### **1. Introduction**

In this article, the dynamical behavior of semiconductor lasers subjected to optical feedback from an external mirror, in the long cavity case [1], based on the experimental observations of the research group I belong to are reviewed. External optical feedback is known to lead to a wealth of dynamical regimes [1,2], some of which have been exploited in diverse applications such as laser feedback interferometry [3], reservoir computing [4], physical-layer secure communications [5], and random-number generation [6]. A classification of the different dynamical regimes of a laser diode with optical feedback has been proposed as early as 1986 by Tkach and Chraplyvy [7], and is still being referred to. The classification features five regimes, four of which involve CW dynamics, and only one, regime IV, corresponds to all other possible dynamics. It has been shown since then that regime IV actually contains a great variety of dynamical regimes. The sequence of regimes experimentally observed within regime IV and leading to chaotic behavior as the feedback level is increased will be focused on, and, when possible, agreement or disagreement with the Lang and Kobayashi rate equation model will be indicated.

The paper is organized as follows: Section 2 reviews previous experimental studies of routes to chaos, Section 3 presents the experimental setup, Section 4 discusses modeling considerations, and Sections 5 and 6 present our observations when the laser is biased close to and far from threshold, respectively; finally, Section 7 summarizes and discusses the main conclusions.

#### **2. State of the Art**

Laser diodes subjected to external optical feedback have been the subject of a large number of publications in the last three decades, focusing either on dynamical behavior or on their use in a variety of applications. We refer the reader to a book [1] and a review paper [2] for extensive information. We focus here on experimental investigations of the sequence of dynamical regimes experienced by the laser as the feedback strength is increased, from CW to chaotic behavior. These routes reveal the way in which intrinsic time scales of a laser with optical feedback interplay and lead to a variety of sustained periodic or quasiperiodic oscillations and eventually chaos. Quasiperiodic [8–10], period-doubling [11], and subharmonic [12] routes to chaos have been reported. Contrary to the quasiperiodic route, which is reported to occur for a wide range of operating conditions, the period-doubling and subharmonic routes have been observed for specific, restricted conditions. Of note, the routes have typically been studied based on observations of a discrete set of feedback levels, and not for continuous tuning. Hohl and Gavrielides have also observed [13], both experimentally and numerically, an alternating sequence of CW and chaotic behavior, referred to as a bifurcation cascade, for a laser biased close to threshold. In their experiment, the optical spectrum was monitored while the feedback level was continuously tuned.

Previous work from our group has revisited the various routes to chaos observed in the literature, confirming and complementing, in the case of a laser being biased close to threshold, the bifurcation cascade route but also providing a different interpretation of the route observed for larger bias currents. In particular, we show that the route that has been named "quasiperiodic" does not contain the sequence of regimes expected in such a case as it involves a number of different attractors and their interplay.

#### **3. Experimental Setup**

The experimental setup is represented in Figure 1. The laser diodes (LD) considered in this manuscript are a range of 1550 nm DFB lasers: packaged (different Mitsubishi ML925B11F diodes) and unpackaged quantum well and quantum dash-based diodes have been used. The temperature of the laser is stabilized +/− 0.01K and its current +/− 0.01A. The LD is subjected to optical self-feedback coming from an external mirror (M) placed at distance *L* from the LD. A variable attenuator, composed of a linear polarizer (LP) and a quarter-wave plate (QWP), is placed in the external cavity. Fine-grained rotation of the QWP allows for a quasi-continuous adjustment of the feedback level η. The optical intensity *I* is monitored with a fast photodetector, and a multimeter is used to determine the DC component, *VDC*, of the laser voltage. In the case of unpackaged lasers, the AC voltage across the laser diode, *VAC*, is measured with a real-time oscilloscope (OSC) and enables the monitoring of the charge carrier density [14,15]. The optical spectrum is tracked with a high-resolution optical spectrum analyzer. Finally, a heterodyne technique, exploiting the beating of the LD with a stable reference laser, is used to measure the optical phase. A description of the principles and implementation of the heterodyne technique can be found in Refs. [14,16].

**Figure 1.** Experimental Setup. LD: laser diode, M: mirror, QWP: quarter-wave plate, P: polarizer; BS: beam splitter, OI: optical isolator, PD: photodetector, BT: bias tee, Amp: amplifier. MM: multimeter, OSA: high-resolution optical spectrum analyzer, OSC: real-time oscilloscope. The model numbers are given in Refs. [17,18]. Not represented: the heterodyne scheme used to measure the optical phase (please refer to Refs. [14,16]).

#### **4. Modeling Considerations**

Even though experimental results are our focus, I will also refer to the Lang and Kobayashi (LK) model [19], which is widely used to interpret the nonlinear dynamics of single-mode laser diodes subjected to optical feedback. It is based on standard semi-classical rate equation modeling, and no spatial effects within the laser cavity are taken into account explicitly. The dynamics involve the total carrier population *N*(*t*), an intra-cavity electric field that is only time-dependent and represented as *E*(*t*)exp[*i*ω0*t* + *i*φ(*t*)], where *E* is the amplitude, φ the slowly-varying phase, and ω<sup>0</sup> the angular frequency of the solitary laser. The terms of the rate equations take into account sources of carrier and photon gains and losses, as well as a coupling between the amplitude and the phase represented by the linewidth enhancement factor α. Lang and Kobayashi have added, in the field equation, a term proportional to the delayed optical feedback. The LK model has proven to be useful in interpreting numerous experimentally observed dynamical behaviors of a LD, and has also been used for prediction (e.g., Refs. [20,21]). In particular, the model shows that, as feedback level is increased, potentially stable CW solutions, named external-cavity modes, and unstable CW solutions, referred to as antimodes, appear in pairs. The equilibria (ECMs) are spaced in frequency by ~*f*τ. They are located on an ellipse in the (*N*(t), φ(*t*) − φ(*t* − τ)) plane, where τ is the round-trip time in the external cavity. ECMs are located on the lower part of the ellipse and antimodes on the upper part of it, as represented in Figure 2. The mode that is closest in frequency to that of the solitary laser is called the minimum linewidth mode (MLM), and denoted ECM 0. Positively shifted ECMs with respect to ECM 0 use positive numbering

(1, 2, 3 ... ), while negatively shifted ECMs use negative numbering. The mode with the lowest optical frequency is the maximum gain mode (MGM).

**Figure 2.** Locations of the equilibria (ECMs) (circles) and antimodes (crosses) in the (*N*(*t*),φ(*t*) − φ(*t* − τ)) plane according to the Lang and Kobayashi (LK) model.

Finally, two time scales are of crucial importance. The first is the relaxation oscillation period, τ*RO*, which is intrinsic to the laser and represents the period of transient oscillations appearing in a LD as a result of the interaction between the carrier and photon populations. The second is the delay introduced by the optical feedback. The frequency of the relaxation oscillations is denoted *fRO* = 1/τ*RO,* and the inverse of the delay is called here the delay frequency *f*<sup>τ</sup> = 1/τ.

#### **5. Route to Chaos When the Laser Is Biased Close to Threshold**

In this section, I present a review of our observations in the case of a laser biased relatively close to threshold [22,23]. In this case, the sequence of bifurcations displays regular or irregular alternation between different regimes; this type of sequence will be referred to as a cascade of bifurcations [13].

Hohl and Gavrielides have reported in Ref. [13], for a current of *J* = 0.99*Jth*, where *Jth* is the solitary laser threshold current, an alternating sequence of CW and chaotic behaviors as the feedback level is increased. Figure 3 represents three experimental bifurcation diagrams for different currents and cavity lengths. The probability density function of the extrema of the optical intensity *I* is represented, using a color map, as a function of the feedback strength η. In panel (b), we observe a regular alternation between two distinct regimes: one is characterized by small-intensity fluctuations, while in the other fluctuations are much larger. This regular alternation is consistent with the optical spectra that have been observed in Ref. [13]. Hohl and Gavrielides also provide an interpretation, based on LK, in which slips toward newly created stable maximum gain modes (MGMs) occur regularly as the feedback level increases and the ellipse grows in size. These slips correspond to abrupt switches to a CW regime, which itself leads, as η is increased, to more complex behavior, including low frequency fluctuations (LFF) and fully developed coherence collapse (CC), involving a number of ECMs. The experimental bifurcation diagrams we have obtained confirm this interpretation and show the robustness of the alternation between regimes for a range of currents and cavity lengths. Specifically, we have found that regular or irregular alternations are consistently observed for currents *J* -1.6*Jth* [23].

**Figure 3.** Experimental bifurcation diagrams of a Mitsubishi ML925B11F diode with (**a**) *J* = 1.58 *Jth* and *L* = 30 cm, (**b**) *J* = 1.21 *Jth* and *L* = 15 cm, and (**c**) *J* = 1.21 *Jth* and *L* = 65 cm. From Ref. [22].

As the current is increased above threshold, we find that the bifurcation cascade progressively disappears. Regions of CW and of large fluctuations are still observed, but not in regular alternation, as illustrated in panel (a). Above 1.6*Jth* approximately, no alternation can be observed [22,23], and the bifurcation structure progressively becomes the one described in the next section.

An increase of the cavity length also leads to a degradation of the regularity of the alternation [23], as illustrated in panel (c). A possible explanation is that, as the cavity length increases, ECMs become more closely spaced in frequency and attractor merging is facilitated. This makes it more difficult for independent attractors to develop, with a significant basin, around a single ECM, and no slip toward a stable CW regime occurs.

Finally, I would like to point out that numerical simulations based on the Lang and Kobayashi model lead to bifurcation cascades for a significantly narrower range of parameter values than experiments do, indicating a possible limit of the model. Comparisons between numerical and experimental bifurcation diagrams can be found in Ref. [23].

#### **6. Route to Chaos at Larger Bias Current**

In this section, I focus on a review of the experimentally observed bifurcation scenario when the laser is biased significantly above threshold. Of note, numerous simulated bifurcation diagrams, based on the LK model, can be found in the literature ([1] and references therein). The first bifurcations observed in simulations usually correspond, in the long cavity case, to an undamping of the relaxation oscillations followed by quasiperiodic behavior in which a second timescale, close to the round-trip time in the external cavity, comes into play. The sequence of bifurcations observed after that is strongly dependent on the choice of the model parameter values. Experimentally, however, we have found that a consistent and robust scenario occurs when the current *J* 2*Jth*. Specifically, we have investigated systematically the route from CW behavior to fully developed coherence collapse, for a range of laser diodes, both packaged and unpackaged, quantum well- or quantum dash-based. Figure 4 represents the bifurcation diagram of a LD biased at *J* = 2.28 *Jth*, for a cavity length *L* =30 cm (*f*<sup>τ</sup> = 500 MHz, *fRO ~* 7.8 GHz). The bifurcation diagram is significantly different from the ones reported in the previous section, for lower current, as no alternation between CW and more complex regimes takes place. We observe a sequence of different regimes leading from CW (region α) to fully developed coherence collapse (region θ), going through quasiperiodic-like (QP) behavior (region β), limit cycle (LC) periodic behavior (region γ), a region of intermittency (region δ) involving a subharmonic (SH) regime, a period-doubled (PD) regime (region ε), and an intermittency region (ζ) between PD and fully developed CC [17,24].

**Figure 4.** Experimental bifurcation diagrams of a Mitsubishi ML925B11F diode with *J* = 2.28 *Jth* and *L* = 30cm. Greek letters indicate regions of existence of various dynamical regimes. α: CW; β: quasiperiodic-like (QP); γ*:* limit cycles (LC); δ: multistate intermittency including subharmonic (SH) behavior, ε*:* period-doubled (PD), ζ: intermittency between PD and coherence collapse (CC), θ: CC.

Reproduced from Ref. [24], with the permission of AIP Publishing.

The traditional analysis of the optical intensity *I* alone is insufficient to unravel this complex sequence of dynamics. For this reason, we have analyzed simultaneously the laser voltage *V*, the optical spectrum, and the optical phase φ. Detailed experimental reports and interpretation of the intensity, phase, and optical spectrum in the various regimes can be found in Refs. [16,17,24,25]; the main points are focused on here.

The first regime that can be identified experimentally looks quasiperiodic in the time domain (region β). A study of the optical spectra and optical phase [16] reveals that it actually involves an alternation in time between a periodic oscillation located around ECM 1 and another ECM, which depends on feedback level (e.g., ECM −3 or ECM −4). The quasiperiodic appearance in time therefore does not result from a torus that would have developed from two successive Hopf bifurcations of a given ECM, as would be expected in a traditional quasiperiodic (Ruelle-Takens or Curry-Yorke) route [26], but is rather the result of the interaction between two equilibria (ECMs). In addition, it is interesting to note that the high frequency of the QP regime is equal to a multiple of the delay frequency *f*<sup>τ</sup> that is close to *fRO* [25] and it remains locked to that multiple if the current, and thus *fRO*, is slightly varied. As the feedback is further increased, and region γ is reached, periodic LC dynamics, located on a single ECM, are observed, again with frequencies that are (different) multiples of *f*τ. Specifically, the last limit cycle of region γ is located around ECM 2 and its RF frequency, measured both from *I*(*t*) and *V*(*t*) is equal to 7 GHz, corresponding to 14 times *f*τ. The previous observations show that there appears to be locking at the RF level between the intrinsic frequencies *f*<sup>τ</sup> and *fRO*, from the very early stages of the dynamics.

As the feedback level is increased, the optical frequencies involved in the dynamics also tend to lock. Indeed, when the PD regime is reached in region ε, the RF frequencies still display locking, as *I* and *V* show period-doubled oscillation at 3.5 GHz, corresponding to a halving of the frequency (7 GHz) of the last limit cycle of region γ. In addition, a locking occurs at the optical frequency level since the ECMs participating in the dynamics, as revealed by the optical spectrum and phase [16], are also separated by 7 GHz. Specifically, ECMS 3, −4, and −11 participate in the period-doubled dynamics: 3 and −11 are separated by 7 GHz, while 3 and −4, and −4 and −11, are separated by 3.5 GHz. Before this complete locking occurs, for lower feedback in region δ, a partial locking is observed to which 3 ECMS participate: two are separated by 3.5 GHz (e.g., ECMs 0 and −7), and thus exhibit locking, while the third, ECM 3, does not. The corresponding dynamical regime (region δ) shows a regular alternation in time between LC and PD oscillations. The duration of the LC and PD oscillations varies with the feedback level, but the sum of the main frequencies in the RF spectra always adds up to ~*fRO* (e.g., *fRO*/3 and 2*fRO*/3), corresponding to a subharmonic regime [12,16,24]. In summary, as the feedback level is increased, the locking between the dynamical frequencies involved in the laser dynamics progresses until full locking is attained. This locked regime is not maintained indefinitely, however, as it is lost to CC. In Ref. [16], the disappearance of the PD regime is interpreted as resulting from a crisis. To illustrate this point, the experimentally measured optical phase φ(*t*) and intensity *I(t)*, in the PD regime, are displayed in Figure 5. We see that the dynamical state moves, as a function of time, from ECM 3, to ECM −4, to ECM −11, then endures an abrupt repulsion toward ECM 3 again. A possible interpretation is that an antimode, located close ECM −11, provides the necessary repelling force, in the direction of its unstable manifold, and thus connects in phase space the distant ECMs −11 and +3. As the feedback level is raised, the ellipse grows in size, leading to an increase in the distance between ECM −11 and the closest antimode. When, for some feedback level, the ellipse becomes too wide, the connection breaks and a boundary crisis to CC behavior occurs. In the CC regime, as reported in numerous publications, a large number of ECMs are involved [1,2] and the intensity, voltage, and phase are observed to vary chaotically. Of note, the sequence of bifurcations described in this section has been observed consistently for a range of quantum well lasers, as well as with a quantum dash laser, illustrating the generality of the results for these types of quantum confinement. Finally, let me mention that the most common route to chaos of a semiconductor laser subjected to optical feedback, which is the one reported in this section, has often been described as a quasiperiodic route. Even though dynamical behaviors of quasiperiodic appearance are indeed observed, we have shown that the route actually differs significantly from a traditional quasiperiodic route [26,27] in which a single equilibrium point undergoes a series of Hopf bifurcations leading to periodic then quasiperiodic behavior, and finally chaos. In Ref. [16], we have proposed to name the sequence of bifurcations a crisis route to chaos.

**Figure 5.** Period-evolution of the experimentally measured optical intensity and optical phase in the period-doubled regime (region ε of Figure 4). The top-right panel gives a visual representation of the ECMs involved. From Ref. [16].

#### **7. Discussion**

In this article, I have reviewed the experimentally observed routes to CC of a laser diode subjected to optical feedback as the feedback level is increased.

When the laser is biased close to threshold, stable dynamics around a single ECM are hardly observed, with the notable exception of CW behavior on the MGM, which displays a negative frequency shift with respect to solitary laser frequency. A typical bifurcation diagram consists of an alternation between complex regimes involving numerous ECMs and regions of stable behavior that occur when the MGM becomes accessible.

For currents significantly above threshold, the picture is different as some stable attractors develop around individual ECMs. In particular, for minor feedback, stable limit cycles develop around ECMs that have a positive frequency shift, which is consistent with the predictions by Masoller and Abraham [28]. Of note, the RF frequencies, measured from the intensity and voltage, always display a locking between the relaxation oscillation frequency and the delay frequency. Specifically, we consistently observe that the first dominant RF frequency that can be identified experimentally is not the relaxation oscillation frequency *fRO* but rather a multiple of the delay frequency *f*<sup>τ</sup> that is close to *fRO*. Of note, by dominant frequency, we either mean the fast frequency, when quasiperiodic-like behavior is the first observed in the bifurcation sequence, or the actual oscillation frequency, when LC dynamics are observed first (as reported in Ref. [18], different ECMs can be experimentally selected as starting states of bifurcation diagrams, resulting in different initial instabilities). Interestingly, similar locking has also been reported in quantum dot lasers subjected to optical feedback [25].

As the feedback level increases and the ellipse grows, allowing for the coexistence of ECMs that are distant in frequency, partial then complete locking of the optical frequencies of the ECMs also occurs. Partial locking takes place in the subharmonic regime. Complete locking occurs in the period-doubled regime, when the optical intensity and voltage oscillate at a frequency that is also equal to the optical frequency spacing between the ECMs involved. This type of regime appears to be a mixed ECM solution, as described by Pieroux et al. in Ref. [29]. Finally, for larger feedback, the locking is lost as a simple regime involving a limited number of ECMS becomes impossible and the laser dynamics become chaotic.

**Funding:** This research was co-funded by Région Grand Est.

**Acknowledgments:** The author thanks Byunchil Kim, Michael J. Wishon, Daeyoung Choi, and David S. Citrin. **Conflicts of Interest:** The author declares no conflict of interest.

#### **References**


© 2020 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

*Article*
