Next Article in Journal
Two-Agent Pareto-Scheduling of Minimizing Total Weighted Completion Time and Total Weighted Late Work
Next Article in Special Issue
Oscillatory Solutions to Neutral Delay Differential Equations
Previous Article in Journal
Multiple Solutions for a Class of New p(x)-Kirchhoff Problem without the Ambrosetti-Rabinowitz Conditions
Previous Article in Special Issue
Analysis of the Sign of the Solution for Certain Second-Order Periodic Boundary Value Problems with Piecewise Constant Arguments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Spiking Laser Controlled by a Delayed Feedback

by
Anton V. Kovalev
1,†,
Evgeny A. Viktorov
1,† and
Thomas Erneux
2,*,†
1
Birzhevaya Liniya 14, ITMO University, 199034 Saint Petersburg, Russia
2
Optique Nonlinéaire Théorique, Université Libre de Bruxelles, Campus Plaine C.P. 231, 1050 Bruxelles, Belgium
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2020, 8(11), 2069; https://doi.org/10.3390/math8112069
Submission received: 13 September 2020 / Revised: 10 November 2020 / Accepted: 16 November 2020 / Published: 20 November 2020
(This article belongs to the Special Issue Recent Advances in Delay Differential and Difference Equations)

Abstract

:
In 1965, Statz et al. (J. Appl. Phys. 30, 1510 (1965)) investigated theoretically and experimentally the conditions under which spiking in the laser output can be completely suppressed by using a delayed optical feedback. In order to explore its effects, they formulate a delay differential equation model within the framework of laser rate equations. From their numerical simulations, they concluded that the feedback is effective in controlling the intensity laser pulses provided the delay is short enough. Ten years later, Krivoshchekov et al. (Sov. J. Quant. Electron. 5394 (1975)) reconsidered the Statz et al. delay differential equation and analyzed the limit of small delays. The stability conditions for arbitrary delays, however, were not determined. In this paper, we revisit Statz et al.’s delay differential equation model by using modern mathematical tools. We determine an asymptotic approximation of both the domains of stable steady states as well as a sub-domain of purely exponential transients.

1. Introduction

One of the unexpected features of early laser systems was the emergence of output pulsations. Even before the first laser has been built, researchers had undertaken studies of instabilities in masers (a maser is a device that produces coherent electromagnetic waves at microwave frequencies through amplification by stimulated emission. The laser works by the same principle as the maser but produces higher frequency coherent radiation at visible wavelengths) and pulsations were actually observed in ruby masers in 1958 [1,2]. The first working laser was a ruby laser made by Maiman in 1960 (the first laser was a solid-state laser: Ruby emitting at 694.3 nm. Ruby consists of the naturally formed crystal of aluminum oxide (Al 2 O 3 )) [3]. The complex spiking behavior of ruby lasers [4] was not anticipated and attracted a great deal of theoretical interest (see Figure 1). Similar pulsations have been reported for solid-state, liquid, and gas laser systems.
First derived by Statz and de Mars [6], two rate equations provide a fair description of the laser dynamical response. Computer simulations, however, predicted a train of regular and damped spikes at the output of the laser [7,8,9]. They are called “relaxation oscillations” (ROs). This discrepancy between theory and experiment has to do with the fact that the spiking behavior dies out very slowly. Mechanical and thermal shocks and disturbances present in many lasers, and especially the ruby laser, act to continually re-excite the spiking behavior [10].
It is desirable for some applications to have a smooth output with no or very small amplitude modulation. Non-spiking operation has been achieved in ruby lasers [11,12,13,14] by sampling a portion of the output beam with a photodiode and using the detector signal to control the voltage applied to an electro-optic shutter located inside the resonator. In order to obtain spike suppression, the delay time between the detector and the shutter must be much shorter than the spike width of approximately 0.1 to 1 μ s [10].
In 1965, Statz et al. [8] investigated theoretically and experimentally the conditions under which spiking in the laser output can be completely suppressed using a delayed optical feedback. They consider the laser rate equations supplemented by a delayed feedback term. The equations in their original form, as well as their formulation in modern notations, are documented in Appendix A. The same stabilization problem due to a delayed feedback has been later considered by Krivoshchekov et al. [15]. Their rate equations are documented in Appendix B. The delay differential equation (DDE) of Statz et al. [8] has been explored numerically, while the DDE problem of Krivoshchekov et al. [15] has been analyzed for small delays and small feedback gains.
To the best of our knowledge, the publication by Statz et al. [8] is the first one in laser physics, where a DDE problem is clearly formulated. The effect of a delayed feedback on a laser is an important topic in laser physics. In particular, semiconductor lasers (SLs) used in our everyday applications are extremely sensitive to optical feedback. Dynamical instabilities were identified in the early 1970s [16,17] and have been studied since then with various motivations [18,19,20].
The plan of the paper is as follows. Section 2 introduces Statz et al. delay equations and analyzes the stability properties of the non-zero intensity steady state. Section 3 reviews the analysis in [15] based on the short delay limit. Lastly, we discuss the role of the damping rate on the stability conditions.

2. Rate Equations

The basic idea of Statz et al. [8] was to provide additional losses in the laser cavity through a negative delayed feedback. Their objective was to determine the feedback conditions for which damped oscillations are replaced by pure exponential decays. In the dimensionless form, the laser rate equations for the intensity output of the laser I and the population inversion D are given by (see Appendix A for details)
d I d t = I 1 + D δ I ( t τ ) ,
d D d t = γ A D ( 1 + I )
where A > 1 is the pump parameter and is one control parameter. γ = O ( 10 6 ) is the ratio of the intensity and population time scales. δ > 0 and τ denote the gain and the delay of the feedback, respectively. The values of the parameters used in [8] are documented in Appendix A. We note that γ is small. The limit γ 0 is singular because γ multiplies the right hand side of Equation (2). We eliminate this singularity by a change of variables described in Appendix A. The new dependent variables x and y represent deviations of I and D with respect to their steady state values I s and D s [21] . The new time s = ω t takes into account the basic time scale of the laser damped ROs when the steady state is slightly perturbed. These equations are
d y d s = ( 1 + y ) x η y ( s θ ) ,
d x d s = y ε x ( 1 + I s ( 1 + y ) ) .
The new parameters η , θ , and ε are defined in Appendix A. The term multiplying ε in Equation (4) describes the slow damping of the free-running laser intensity oscillations. The delayed term in Equation (1) is a second source of dissipation if the delay θ is small. Assuming η > > ε , we neglect the term multiplying ε in (4). We obtain
d y d s = ( 1 + y ) x η y ( s θ ) ,
d x d s = y .
The non-zero intensity steady state is ( x , y ) = ( 0 , 0 ) . From the linearized equations, we determine the characteristic equation for the growth rate λ . It is given by
λ 2 + λ η exp ( λ θ ) + 1 = 0 .
Introducing λ = u + i v into Equation (7) and taking the real and imaginary parts provide two equations for u and v . We find
u 2 v 2 + η exp ( u θ ) u cos ( v θ ) + v sin ( v θ ) + 1 = 0 ,
2 u v + η exp ( u θ ) u sin ( v θ ) + v cos ( v θ ) = 0 .
The case of a purely real λ corresponds to v = 0 . From Equation (8), we determine u = u ( η ) in implicit form as
η = u 2 + 1 u exp ( u θ )
and since η > 0 , u < 0 . The case of complex λ has been investigated by solving numerically Equations (8) and (9). Figure 2 represents the real part u as a function of η . The full lines correspond to real λ for θ = 0 (black) and θ = 0.1 (red). The broken lines correspond to complex λ for θ = 0 (black) and θ = 0.1 (red). The circles indicate limit points (LPs) where there is a transition between real and complex λ . If θ 0 , the domain of pure exponential decays lies between the two red LPs. They are determined by using (10) and by applying the condition d η / d u = 0 . We find that u = u ( θ ) , in implicit form, is
θ = 1 u 2 u ( 1 + u 2 ) .
Together with (10), we have a parametric solution for u ( θ ) and η ( θ ) . We proceed as follows. We progressively decrease u from a low negative value and determine θ from (11). We then use (10) with θ = θ ( u ) and determine η . The solution η = η ( θ ) is shown in Figure 3. There is a critical point θ = θ c above which no real roots are possible. To find its value, we use (11) and determine u = u c from the condition d θ / d u = 0 . We find u c = 2 + 5 . We substitute this value into (11) and obtain
θ c 0.3 .
If θ 0 , we find from (10) and (11), the following approximations for the LPs
η 1 2 2 θ ,
η 2 θ 1 exp ( 1 ) .
The first LP moves to the fixed value η 1 = 2 , while the second LP moves to infinity in the limit of small delays. The later case will be further analyzed in Section 3.
A Hopf bifurcation instability is possible and corresponds to u = 0 . From (8) and (9), we determine the conditions
v 2 + η v sin ( v θ ) + 1 = 0 ,
η v cos ( v θ ) = 0 .
They admit the solutions
( 1 ) : η = 0 ,
( 2 ) : η = θ n 2 θ 2 θ n θ > 0 , ( n = 0 , 2 , . . ) ,
( 3 ) : η = θ 2 θ n 2 θ n θ > 0 ( n = 1 , 3 , . . )
where
θ n π 2 ( 1 + 2 n )
and n is detailed in (18) and (19). The dots in Figure 4 are the values considered by Statz et al. [8] in their numerical simulations (listed in Table A1). As we can see, five dots are in the domain of unstable steady states. The simulations in [8] were limited in integration time but clearly showed growing oscillations for those values (except for the middle point at η = 4 ) .
The predictions of the linear stability analysis are verified numerically by simulating Equations (3) and (4). Figure 5a corresponds to η = 0 (no feedback) and exhibits the slowly damped ROs. Figure 5b is for η = 1 , which is still part of the domain of decaying oscillations although they are quickly damped. Figure 5c is for η = 2.2 , where the decay to the steady state x = 0 is purely exponential. We note that already for η = 1 , oscillations are strongly damped. This can be explained by considering Equations (3) and (4), and by analyzing the linear stability of ( x , y ) = ( 0 , 0 ) . The characteristic equation for the growth rate λ is
λ 2 + λ η exp ( λ θ ) + ε ( 1 + I s ) + 1 + η ε ( 1 + I s ) exp ( λ θ ) = 0 .
If η = 0 , the damping rate of the RO oscillations is given by
Γ = ε 2 ( 1 + I s )
as ε 0 . Assuming η = O ( ε ) , the damping rate is
Γ = 1 2 ε ( 1 + I s ) + η cos ( θ )
in first approximation as ε 0 . If cos ( θ ) > 0 , as in Figure 5, Γ increases linearly with η .

3. Small Delay Limit

Statz et al. [8]’s strategy for eliminating undesired laser oscillations is based on the observation that the roots of the characteristic equation are real and negative if η > 2 in the case of zero delay. In other words, the decay towards the stable steady state is purely exponential if the feedback gain is sufficiently large. Provided that the delay is small, the authors reasonably expected the same transient evolution to the steady state. Krivoshchekov et al. [15] expanded the delayed variable for small τ as
I ( t τ ) = I ( t ) τ I ( t ) + O ( τ 2 ) .
They then analyzed the stability of the steady state. Instead of (7), the characteristic equation is
( 1 η θ ) λ 2 + λ η + 1 = 0
The stability condition is
1 η θ > 0 .
Using Table A2 in Appendix B, (26) can be rewritten in terms of the original parameters as
1 2 α 2 I s 2 σ 0 T d > 0
which is the inequality (5) in [15]. From (25), we determine the conditions for real roots given by
η > η 1 = 2 2 θ + O ( θ 2 ) .
The expansion (24) however does not lead to the second LP because η 2 = O ( θ 1 ) is large . The limit where η θ 1 as θ 0 is however interesting because it suggests that a small delay may lead to sustained oscillations provided the negative feedback rate is sufficiently large. A careful analysis of the linearized equations evaluated at the Hopf bifurcation point suggests the following scalings as θ 0
η ¯ = θ η , s ¯ = θ 1 s , and x ¯ = θ 1 x .
Inserting (28) into Equations (5) and (6) leads to the following equations for x ¯ and y
d y d s ¯ = ( 1 + y ) ( θ 2 x ¯ η ¯ y ( s ¯ 1 ) ) ,
d x ¯ d s ¯ = y .
In the limit θ 0 , Equation (29) reduces to
d y d s ¯ = η ¯ y ( s ¯ 1 ) ( 1 + y )
which we recognize as Wright equation [22] (or equivalently, the delayed logistic equation, after the change of variable y ¯ = 1 + y ). This specific limit is documented in [23] (p. 117). Oscillations emerge at η ¯ = π / 2 with a period P = 4 and quickly becomes pulsating as η ¯ is further increased.

4. Discussion

The interest of purely negative real eigenvalues of a DDE is rarely emphasized in the DDE literature essentially devoted to the conditions of oscillatory instabilities. They do, however, play a role in applications, and they strongly depend on the delay.
It is quite remarkable that in 1965, a DDE was analyzed in order to understand how to eliminate undesirable oscillations by using a delayed feedback. Statz et al. [8] were exploring the conditions for purely real eigenvalues as means to obtain exponentially decaying transients. From their experiments and their numerical simulations, they conclude that the delay should be small, and that the feedback gain needs to be sufficiently large. In this paper, we substantiate their findings analytically by determining the conditions
η 1 ( θ ) < η < η 2 ( θ ) and θ < θ c
where η 1 and η 2 admit the approximations (13) and (14) for small delays, respectively, and θ c is given by (12).
The damping rate ε of the natural laser oscillations is small like 10 3 for the ruby laser, but also for other solid state and semiconductor lasers. In our analysis, we neglect its contribution because the delay term dominates the stability conditions as the feedback rate η is increased. However, there are other problems where the damping rate could be moderate to large. Consider, for example, the delayed Duffing equation, which has benefited from recent mathematical interests [24,25,26]. Duffing equation models an elastic pendulum whose spring’s stiffness is nonlinear. Supplemented by a delayed velocity control term, it is given by
x + ε x + x + x 3 + c x ( t τ ) = 0 .
The linearized problem leads to the characteristic equation
λ 2 + ε λ + 1 + c λ exp ( λ τ ) = 0
which is Equation (7) supplemented by the ε λ term. The analysis of the real and purely imaginary roots of Equation (34) is similar to our previous analysis with ε = 0 . In Figure 6, we contrast the stability diagrams of ε = 0 and ε = 1 . The domain of purely exponential decay increases with ε , while the change of stability through a Hopf bifurcation moves to higher delays. This is expected because the damping term contributes to stabilizing the zero solution.

Author Contributions

All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

A.V.K. and E.A.V. were supported by the Government of Russian Federation (grant 08-08).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The DDE Problem of Statz et al.

Equations (9) and (10) of Statz et al. [8] for the intensity of the laser field u and the inverted population v are
d u d T = u Q v 1 b u ( T T d ) ,
d v d T = 1 T 0 ( v α ) + ( 1 α ) u v .
u and v are normalized such that u = v = 1 is the non-zero intensity steady state in the absence of feedback ( b = 0 ). The dimensionless time T = ω 0 T is the original time variable T multiplied by the angular optical frequency ω 0 = 2.7 × 10 15 s 1 . Q = 3.6 × 10 7 is the quality factor of the cavity and T 0 = ω 0 t 0 = 6 × 10 12 where t 0 is approximately the fluorescent life time. T d = ω 0 T d is the dimensionless delay and α = 2 is the estimated value of the pump parameter. The values of the feedback amplitude b and delay T d used in the simulations are [8]
b = 10 3 ( Figure 1 ) , 5 × 10 3 ( Figure 2 ) , 10 2 ( Figure 3 ) , T d = 0 , 0.24 × 10 9 , 1.2 × 10 9 , 25 × 10 9 .
In the experiments, a feedback circuit is used with a minimum delay of 4 × 10 7 s . In terms of the original delay T d , the delays considered in the simulations are 0 , 0.9 × 10 7 s , 4.4 × 10 7 s , 92 × 10 7 s and are of the same range as the delay of the feedback control circuit. The experiments considered a feedback amplitude of approximately 5 × 10 3 and motivated the values of b detailed in (A3).
One parameter can be eliminated by reformulating Equations (A1) and (A2). Specifically, we introduce the new variable t , D, I and parameter A defined as
t T / Q , D v , I ( 1 α ) u , A α 1 ,
Inserting (A4) into Equations (A1) and (A2), we obtain
d I d t = I 1 + D δ I ( t τ ) ,
d D d t = γ A D ( 1 + I )
where
γ Q / T 0 = 6 × 10 6 , τ T d / Q , and δ b / ( α 1 ) .
The values of the feedback amplitude δ are identical to the values of b shown in (A3). The values of the delay τ corresponding to the values of T d listed in (A3) are
τ = 0 , 6.7 , 33 , and 694 .
Equations (A5) and (A6) are the classical rate equations for a Class B laser exhibiting a small γ [21] . The non-zero intensity steady state ( I , D ) = ( I s , D s ) is given by
I s = ( 1 + δ ) + ( 1 + δ ) 2 + 4 δ ( A 1 ) 2 δ = A 1 + O ( δ ) ,
D s = 1 + δ I s = 1 + O ( δ ) .
We next introduce the basic time scale of the laser equations given by
s ω t
where ω is the laser relaxation oscillation frequency defined by
ω γ D s I s
After introducing (A11) into Equations (A5) and (A6), we introduce the new variables x and y defined as the deviations from D and I from their steady state values. Specifically, x and y are related to D D s and I I s as
D D s = ω x and I = I s ( 1 + y )
In terms of s, x and y, Equations (A5) and (A6) become
d y d s = ( 1 + y ) x η y ( s θ ) ,
d x d s = y ε x ( 1 + I s ( 1 + y ) )
where the values of the parameters are listed in Table A1. Equations (A14) and (A15) are now in an ideal form for analysis.
Table A1. Dimensionless parameters and their values.
Table A1. Dimensionless parameters and their values.
ω = γ D s I s 2.45 × 10 3
η = δ I s ω 0.41 , 2.04 , 4.1
ε = ω D s I s 2.45 × 10 3
θ = ω τ 0 , 1.64 × 10 2 , 8.2 × 10 1 , 1.7

Appendix B. The DDE Problem of Krivoshchekov et al.

Equations (1) and (2) of Krivoshchekov et al. [15] are given by
d I d T = σ 0 β N 1 α 2 I 2 ( T T d ) ,
d N d T = Γ 1 ( 1 + I ) N
where
σ 0 = 1.5 × 10 8 s 1 , Γ = 2.5 × 10 2 s 1 , β = 1.25 or 2 , and T d ( max ) = 10 7 s .
Introducing
t = σ 0 T , D = β N
into Equations (A16) and (A17), we obtain
d I d t = D 1 α 2 I 2 ( t τ ) I
d D d t = γ A ( 1 + I ) D
where the new parameters are
γ = Γ σ 0 = 1.67 × 10 6 , A = β , τ ( max ) = σ 0 T d = 15 .
The steady state satisfies
D 1 α 2 I 2 = 0 ,
A ( 1 + I ) D = 0 .
Introducing the new variables x, y and s defined by
I = I s ( 1 + y ) ,
D = D s + ω x ,
s = ω t
we obtain
d y d s = ( 1 + y ) x η ( y ( s θ ) + 1 2 y 2 ( s θ ) )
d x d x = y ε x ( 1 + I s ( 1 + y ) )
where the definitions and the values of the parameters are listed in Table A2. Except for the feedback amplitude η , the definitions of ω , ε , and θ are the same as for Statz et al. equations (see Table A1).
Table A2. Dimensionless parameters and their values.
Table A2. Dimensionless parameters and their values.
ω = γ D s I s 6.5 × 10 4
η = 2 α 2 I s 2 ω
ε = ω D s I s 2.6 × 10 3
θ = ω τ
θ ( max ) = ω τ ( max ) 9.75 × 10 3

References

  1. Makhov, G.; Kikuchi, C.; Lambe, J.; Terhune, R.W. Maser action in ruby. Phys. Rev. 1958, 109, 1399–1400. [Google Scholar] [CrossRef]
  2. Kikuchi, C.; Lambe, J.; Makhov, G.; Terhune, R.W. Ruby as a maser material. J. Appl. Phys. 1959, 30, 1061–1066. [Google Scholar] [CrossRef]
  3. Maiman, T.H. Stimulated optical radiation in ruby. Nature 1960, 187, 493–494. [Google Scholar] [CrossRef]
  4. Collins, R.J.; Nelson, D.F.; Schawlow, A.L.; Bond, W.; Garrett, C.G.B.; Kaiser, W. Coherence, narrowing, directionality, and relaxation oscillations in the light emission from ruby. Phys. Rev. Lett. 1960, 5, 303–305. [Google Scholar] [CrossRef]
  5. Nelson, D.F.; Boyle, W.S. A continuously operating ruby optical maser. Appl. Opt. 1962, 1, 181–183. [Google Scholar] [CrossRef]
  6. Statz, H.; DMars, G. Transients and Oscillation Pulses in Masers. In Quantum Electronics; Townes, C.H., Ed.; Columbia University Press: New York, NY, USA, 1960; pp. 530–537. [Google Scholar]
  7. Dunsmuir, R. Theory of relaxation oscillations in optical masers. J. Elerctron. Control 1981, 10, 453–458. [Google Scholar] [CrossRef]
  8. Staz, H.; DeMars, G.A.; Wilson, D.T. Problem of spike elimination in lasers. J. Appl. Phys. 1965, 30, 1510–1514. [Google Scholar] [CrossRef]
  9. Strozyk, J. Observations on a concentric spherical cavity laser oscillator. IEEE J. Quantum. Electr. 1967, 3, 343–348. [Google Scholar] [CrossRef]
  10. Koechner, W. Solid-State Laser Engineering, 6th ed.; Springer: New York, NY, USA, 2006. [Google Scholar]
  11. Marshall, F.R.; Roberts, D.L. Use of electro-optical shutters to stabilize ruby laser operation. Proc. IRE 1962, 50, 2108. [Google Scholar]
  12. Keller, D.V.; Davis, B.I. High power non-spiking operation of ruby laser. IEEE J. Quantum. Electr. 1966, QE-2, 179–181. [Google Scholar] [CrossRef]
  13. Thomas, C.H.; Price, E.V. Feedback control of a Q-switched ruby laser. IEEE J. Quantum. Electr. 1966, QE-2, 617–623. [Google Scholar] [CrossRef]
  14. Lovberg, R.V.; Wooding, E.R.; Yeoman, M.L. Pulse stretching and shape control by compound feedback in a Q-switched ruby laser. IEEE J. Quantum. Electr. 1975, QE-11, 17–21. [Google Scholar] [CrossRef]
  15. Krivoshchekov, G.V.; Makukha, V.K.; Tarasov, V.M. Stabilization of ruby laser radiation by an external negative feedback. Sov. J. Quantum Electron. 1975, 5, 394–396. [Google Scholar] [CrossRef]
  16. Broom, R.F. Self modulation at gigahertz frequencies of a diode laser coupled to an external cavity. Electron. Lett. 1969, 5, 571–572. [Google Scholar] [CrossRef]
  17. Broom, R.; Mohn, E.; Risch, C.; Salathe, R. Microwave self-modulation of a diode laser coupled to an external cavity. IEEE J. Quantum Electr. 1970, 6, 328–334. [Google Scholar] [CrossRef]
  18. Kane, D.M.; Shore, K.A. Unlocking Dynamical Diversity: Optical Feedback Effects on Semiconductor Lasers; Wiley: New York, NY, USA, 2005. [Google Scholar]
  19. Soriano, M.C.; Garcia-Ojalvo, J.; Mirasso, C.R.; Fischer, I. Complex photonics: dynamics and applications of delay-coupled semiconductor lasers. Rev. Mod. Phys. 2013, 85, 421–470. [Google Scholar] [CrossRef] [Green Version]
  20. Chembo, Y.K.; Brunner, D.; Jacquot, M.; Larger, L. Optoelectronic oscillators with time-delayed feedback. Rev. Mod. Phys. 2019, 91, 035006. [Google Scholar] [CrossRef]
  21. Erneux, T.; Glorieux, P. Laser Dynamics; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  22. Smith, H.L. An Introduction to Delay Differential Equations with Applications to the Life Sciences; Springer: New York, NY, USA, 2011. [Google Scholar]
  23. Grigorieva, E.V.; Kaschenko, S.A. Asymptotic Representation of Relaxation Oscillations in Lasers; Springer: Cham, Switzerland, 2017. [Google Scholar]
  24. Davidow, M.; Shayak, B.; Rand, R.H. Analysis of a remarkable singularity in a nonlinear DDE. Nonlear Dyn. 2017, 90, 317–323. [Google Scholar] [CrossRef] [Green Version]
  25. Sah, S.M.; Fiedler, B.; Shayak, B.; Rand, R.H. Unbounded sequences of stable limit cycles in the delayed Duffing equation: An exact analysis. arXiv 2019, arXiv:1908.06533. [Google Scholar]
  26. Fiedler, B.; Ĺopez Nieto, A.; Rand, R.H.; Sah, S.M.; Schneider, I.; de Wolff, B. Coexistence of infinitely many large, stable, rapidly oscillating periodic solutions in time-delayed Duffing oscillators. J. Differ. Equ. 2020, 268, 5969–5995. [Google Scholar] [CrossRef] [Green Version]
Figure 1. High-frequency pulsating oscillations in the maser output. Adapted with permission from Figure 2b of [5] © The Optical Society.
Figure 1. High-frequency pulsating oscillations in the maser output. Adapted with permission from Figure 2b of [5] © The Optical Society.
Mathematics 08 02069 g001
Figure 2. Real part u as a function of the feedback rate η for θ = 0 (black) and θ = 0.1 (red). The full and broken lines correspond to real roots ( v = 0 ) and complex roots ( v 0 ) of the characteristic Equation (7), respectively. Transitions between real and complex roots occur at limit points (encircled).
Figure 2. Real part u as a function of the feedback rate η for θ = 0 (black) and θ = 0.1 (red). The full and broken lines correspond to real roots ( v = 0 ) and complex roots ( v 0 ) of the characteristic Equation (7), respectively. Transitions between real and complex roots occur at limit points (encircled).
Mathematics 08 02069 g002
Figure 3. Domain in the ( θ , η ) parameter space for negative real roots of the characteristic Equation (7). The red dots at θ = 0.1 correspond to the limit points of the S-shaped curve in Figure 2.
Figure 3. Domain in the ( θ , η ) parameter space for negative real roots of the characteristic Equation (7). The red dots at θ = 0.1 correspond to the limit points of the S-shaped curve in Figure 2.
Mathematics 08 02069 g003
Figure 4. First Hopf bifurcation line delimiting the domain of stable steady states (black). It is given by (18) with n = 0 . It contains a smaller domain where the decay to the steady state is purely exponential (red). The dots are the values of the parameters used by Statz et al. [8] in their numerical simulations. They are documented in Table A1.
Figure 4. First Hopf bifurcation line delimiting the domain of stable steady states (black). It is given by (18) with n = 0 . It contains a smaller domain where the decay to the steady state is purely exponential (red). The dots are the values of the parameters used by Statz et al. [8] in their numerical simulations. They are documented in Table A1.
Mathematics 08 02069 g004
Figure 5. Stabilizing effect of the delayed feedback. Equations (3) and (4) have been simulated numerically with increasing values of η and with the initial conditions x ( 0 ) = 0.1 and y = 0 ( θ < s 0 ) . The fixed parameters are θ = 0.1 , ε = 2.45 × 10 3 , I s = 1 . (a) η = 0 , (b) η = 1 , (c) η = 2.2 .
Figure 5. Stabilizing effect of the delayed feedback. Equations (3) and (4) have been simulated numerically with increasing values of η and with the initial conditions x ( 0 ) = 0.1 and y = 0 ( θ < s 0 ) . The fixed parameters are θ = 0.1 , ε = 2.45 × 10 3 , I s = 1 . (a) η = 0 , (b) η = 1 , (c) η = 2.2 .
Mathematics 08 02069 g005
Figure 6. Stability diagram for the delayed Duffing equation. The first Hopf bifurcation line delimits the domain of stable steady states (black). It contains a smaller domain where the decay to the steady state is purely exponential (red).
Figure 6. Stability diagram for the delayed Duffing equation. The first Hopf bifurcation line delimits the domain of stable steady states (black). It contains a smaller domain where the decay to the steady state is purely exponential (red).
Mathematics 08 02069 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kovalev, A.V.; Viktorov, E.A.; Erneux, T. Non-Spiking Laser Controlled by a Delayed Feedback. Mathematics 2020, 8, 2069. https://doi.org/10.3390/math8112069

AMA Style

Kovalev AV, Viktorov EA, Erneux T. Non-Spiking Laser Controlled by a Delayed Feedback. Mathematics. 2020; 8(11):2069. https://doi.org/10.3390/math8112069

Chicago/Turabian Style

Kovalev, Anton V., Evgeny A. Viktorov, and Thomas Erneux. 2020. "Non-Spiking Laser Controlled by a Delayed Feedback" Mathematics 8, no. 11: 2069. https://doi.org/10.3390/math8112069

APA Style

Kovalev, A. V., Viktorov, E. A., & Erneux, T. (2020). Non-Spiking Laser Controlled by a Delayed Feedback. Mathematics, 8(11), 2069. https://doi.org/10.3390/math8112069

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop