Next Article in Journal
Comparison of Novel Raman Techniques for Fiber Forensics
Previous Article in Journal
Thermal Effects on Optical Chirality, Mechanics, and Associated Symmetry Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Temporal Cavity Soliton Interaction in Passively Mode-Locked Semiconductor Lasers

by
Andrei G. Vladimirov
Weierstrass Institute, Mohrenstr. 39, 10117 Berlin, Germany
Optics 2023, 4(3), 433-446; https://doi.org/10.3390/opt4030031
Submission received: 14 June 2023 / Revised: 2 July 2023 / Accepted: 13 July 2023 / Published: 26 July 2023
(This article belongs to the Section Nonlinear Optics)

Abstract

:
Weak interactions of temporal cavity solitons resulting from gain saturation and recovery in a delay differential model of a long cavity semiconductor laser were studied numerically and analytically using an asymptotic approach. This paper shows that in addition to the usual soliton repulsion leading to a harmonic mode-locking regime, soliton attraction is also possible in a laser with a nonzero linewidth enhancement factor. It is shown numerically that this attraction can lead either to pulse merging or to pulse bound-state formation.

1. Introduction

Temporal cavity solitons (TCSs) are short nonlinear optical pulses generated by mode-locked lasers and optical microresonators that preserve their shape in the course of propagation [1,2,3]. In lasers, unlike the usual self-starting mode-locked pulses generated above the linear laser threshold, TCSs coexist with a stable laser off regime and require a finite perturbation for their excitation. For example, when the cavity length of a laser with a semiconductor gain medium is sufficiently large, usual mode-locked pulses can be transformed into TCSs [4], corresponding to a non-self-starting mode-locking regime. In many practical situations, when more than one TCS is exited in an optical cavity, weak interactions among the TCSs may take place via their exponentially decaying tails. Spatial and temporal dissipative soliton interaction in lasers with saturable absorbers have been studied in many publications in cases in which the gain and absorption variables are eliminated adiabatically and interactions take place only via the overlapping electric fields of the pulses [5,6,7,8,9,10]. The interaction of mode-locked pulses in the presence of finite relaxation times of the gain and/or absorber media has been less investigated. In this case, the electromagnetic field saturates gain and absorption behind the pulse, and their slow relaxation can affect the position of the next pulse traveling in the cavity. This type of interaction was studied in refs. [11,12,13,14,15]. In particular, it was demonstrated theoretically and verified experimentally with solid state and fiber lasers [11] that the interaction due to gain depletion and very slow recovery can produce a repulsive force between adjacent pulses, leading to the formation of a harmonic mode-locking regime. A similar conclusion was made in ref. [13] using the delay differential equation (DDE) model [16,17,18] of a mode-locked monolithic semiconductor laser, where similarly to ref. [11], the gain recovery time was much longer than the cavity round-trip time. Here, using the same DDE model, I consider mode-locked pulse interaction in the TCS regime, where the cavity round-trip time is sufficiently long—much longer than the gain recovery time. Based on an asymptotic approach, the equations governing the slow evolution of the time separation and phase difference of the interacting TCSs are derived and analyzed. Asymptotic study of weak TCS interactions in DDE models of optical systems has already been previously carried out in [19,20,21]. However, in ref. [21], devoted to TCS interaction in nonlinear mirror mode-locked lasers, and in this paper, a closed analytical form of the interaction equations is derived. Furthermore, among the above-cited works, only refs. [13,15] are devoted to pulse interactions in semiconductor lasers. However, in ref. [13], a short cavity laser operating far away from the TCS regime is considered, while the analysis in ref. [15] is more or less empirical. Therefore, a more rigorous analysis of TCS interactions in semiconductor lasers that takes into account their characteristic features, such as the linewidth enhancement factor, is required. This work aimed to fill in this gap. Using interaction equations, I show that the TCS interaction scenarios depend on the “interaction coefficients” introduced below and can be more rich than those described in [11,13,15]. In addition to the pulse repulsion resulting in a harmonic mode-locking regime, TCS attraction leading to either pulse merging or bound-state formation can take place in a laser with a nonzero linewidth enhancement factor. Note that soliton attraction leading to bound-state formation was previously observed in ref. [12] in a complex Ginzburg–Landau equation-type mode-locked laser model with second-order dispersion and in the DDE model of a nonlinear mirror mode-locked laser [21]. However, unlike the present work, where the mode-locking mechanism is due to the interplay of gain and absorption saturation and recovery [16], in both those papers, Kerr nonlinearity played a decisive role in the process of mode-locked pulse formation. Furthermore, since ref. [12] considered the limit of infinitely large gain recovery time, the mechanism of pulse interaction in this paper is different and can be attributed to the saturation and slow recovery of absorption rather than gain. The results of this work might be useful for pulse spacing manipulation in multipulse mode-locked lasers.

2. Model Equations

The DDE model of a passively mode-locked semiconductor laser for the electric field amplitude ( A t ) at the entrance of the laser absorber section, the saturable gain ( G t ), and the saturable absorption ( Q t ), can be written in the form [16,17,18]:
γ 1 t A + 1 + i ω A = R ( t T ) A ( t T ) ,
t G = g 0 γ g G e Q e G 1 | A | 2 ,
t Q = q 0 γ q Q s 1 e Q | A | 2 ,
with
R ( t ) = κ e ( 1 i α g ) G ( t ) / 2 ( 1 i α q ) Q ( t ) / 2 + i ϕ i ω T .
where t is the time variable; κ is the attenuation factor describing linear non-resonant intensity losses per cavity round trip; and α g and α q are the linewidth enhancement factors in the gain and absorber sections, respectively. The time delay parameter (T) stands for the cold cavity round-trip time; γ is the spectral filtering bandwidth; γ g and γ q are the normalized carrier relaxation rates in the gain and absorber sections; and s is the ratio of the saturation intensities of these sections. The pump parameter ( g 0 ) depends on the injection current in the gain section, while q 0 is the unsaturated loss parameter, which depends on the inverse voltage applied to the absorber section. Parameter ϕ is the phase shift describing the detuning between the central frequency of the spectral filter and the closest cavity mode and ω is the reference frequency.
It is well known that in a certain parameter domain, Equations (1)–(3) demonstrate pulsed solutions corresponding to fundamental single pulse and harmonic multipulse mode-locking regimes [16,17,18]. Furthermore, it was shown in [4] that when the laser cavity is sufficiently long that the round-trip time is much larger than the gain relaxation time, these pulses can be transformed into TCSs sitting on the stable laser off solution. In this situation, two well-separated mode-locking pulses can interact only weakly via their exponentially decaying tails. Furthermore, when the pulses are sufficiently far away from one another, this interaction is mainly due to the gain component (G), which usually decays much slower than the electric field envelope (A) and saturable absorption (Q). Note, however, that when the distance between the TCSs becomes small enough, the interaction via absorption dynamics might also come into play and even lead to pulse bound-state formation; see ref. [12], where the case of infinitely large gain recovery time is considered.
In order to derive the TCS interaction equations, we rewrite the model equations in a more general real vector form
t U = F ω U + H ω U ( t T ) ,
where U = U 1 U 2 U 3 U 4 T is the real column vector with the components U 1 = A , U 2 = A , U 3 = G g 0 / γ g , and U 4 = Q q 0 / γ q . The two vectors in the right-hand side of Equation (4) are defined by
F ω U = γ U 1 ω U 2 γ U 2 + ω U 1 γ g U 3 e U 4 q 0 / γ q e U 3 + g 0 / γ g 1 U 1 2 + U 2 2 γ q U 4 s 1 e U 4 q 0 / γ q U 1 2 + U 2 2 ,
and
H ω U = R U 1 + i U 2 R U 1 + i U 2 0 0
with
R t = γ κ e ( 1 i α g ) U 3 + g 0 / γ g / 2 ( 1 i α q ) U 4 + q 0 / γ q / 2 i ω T .

3. Temporal Cavity Soliton

Let us assume that the relation
γ 1 < γ q 1 < γ g 1 T ,
for the relaxation rates in the model Equations (1)–(3) is satisfied. This means that the round trip time in a multimode semiconductor laser cavity is sufficiently long, much longer than the gain relaxation time. In this case, the DDE model can have TCS solutions [4]. Let U = u = u 1 u 2 u 3 u 4 T and ω = ω 0 be a TCS solution of Equation (4) corresponding to a narrow mode-locked pulse with the duration τ p γ 1 . Here, u t = u t + T 0 is periodic in time with the period T 0 close to the delay time T. In terms of the original model Equations (1)–(3) we have
u = A 0 t A 0 t G 0 t g 0 / γ g Q 0 t q 0 / γ q T ,
where A 0 t , G 0 t , and Q 0 t is a T 0 -periodic TCS solution of these equations. A numerically calculated intensity time trace of the TCS solution is shown in Figure 1a.
The decay rates of the TCS tails are determined by the following linearization [22,23,24] of Equation (4) on the solution U = 0 :
γ 1 t a + 1 + i ω 0 a = R 0 a t + δ ,
t v 3 = γ g v 3 .
t v 4 = γ q v 4 ,
where a = v 1 + i v 2 , v = v 1 v 2 v 2 v 4 T is a small perturbation vector,
R 0 = κ e ( 1 i α g ) g 0 / ( 2 γ g ) ( 1 i α q ) q 0 / ( 2 γ q ) + i ϕ + i ω 0 δ ,
and the time advance parameter is δ = T 0 T . It follows from Equations (7) and (8) that the decay rates of the TCS gain and absorption components at large positive times t are determined by the two eigenvalues λ g , q = γ g , q , while Equation (6) has an infinite number of eigenvalues defined by
λ k = γ 1 + i ω 0 δ 1 W k γ δ e 1 + i ω 0 γ δ R 0 .
here W k is the Lambert function with the index k = 0 , ± 1 , ± 2 . For the parameter values of Figure 1 we get λ 0 = 3.741 and λ 1 = 16.892 , while all the other eigenvalues are complex and have positive real parts greater than λ 1 .
Assuming that the origin of the time coordinate, t = 0 , is located at the TCS power peak, we determine based on (7)–(9) that at sufficiently large t > 0 the TCS trailing edge can be represented as
u 1 , 2 b 1 , 2 e λ 0 t , u 3 b 3 e γ g t , u 4 b 4 e γ q t ,
where b 1 , 2 , 3 , 4 are real constants, which can be calculated numerically.
Further, let us consider the leading tail of the TCS at negative times, t < 0 . Since Equations (7) and (8) have no eigenvalues with positive real parts, gain and absorption components of the TCS leading edge, u 3 and u 4 , decay faster than exponentially in negative time [22]. The field component u 1 + i u 2 of the leading tail decays exponentially at t < 0 , with the decay rate determined by the eigenvalue λ 1 having the smallest positive real part. Since the inequality γ g , q < | λ 0 | < λ 1 is satisfied, the field component of the TCS decays faster in both time directions than the gain and absorption components in positive time. This means that the interaction via the electromagnetic field dynamics can be neglected when considering the interaction of two well-separated TCSs. Such a type of interaction is typical of lasers with slow gain and absorption and can be viewed as the long-range interaction [21], in contrast with the short-range interaction due to the overlap of the electric fields considered in [5,6,7,8,9,25]. Furthermore, since gain and absorption components of a TCS decay faster than exponentially in negative time, the leading tails of the TCSs can be neglected in the derivation of the interaction equations. Figure 1b shows the time evolution of the absolute values of field envelope U 01 2 + U 02 2 , gain U 03 , and absorption U 04 components of this solution in logarithmic scale. It is seen from this figure that the gain component dominates over the field and absorption ones during almost the whole time interval between the two consequent pulses.
The stability of the TCS depends on the spectrum of the linear operator L obtained via linearization of Equation (4) at the solution U = u and ω = ω 0 . Since the model Equations (1)–(3) are invariant under time translations, U t U t t 0 , and phase shifts, U 1 + i U 2 U 1 + i U 2 e i ϕ 0 , with arbitrary constants t 0 and ϕ 0 , the operator L has two zero eigenvalues corresponding to the neutral (Goldstone) modes given by θ = t u and φ = u 2 u 1 0 0 T , respectively, L θ = t θ + B θ + C T θ t T = 0 and L φ = 0 . Here, B = B u and C T = C u t T are the linearization matrices of F ω 0 U and H ω 0 U ( t T ) at U = u . Below, it is assumed that the TCS is stable and the rest of the spectrum of the operator L is located in the left half of the complex plane. The adjoint linear operator L also has two zero eigenvalues associated with the so-called adjoint neutral modes θ and φ , L θ = t θ + θ B + θ t + T C = 0 and L φ = 0 . Let the adjoint neutral modes be biorthogonal to the neutral modes, θ · φ = φ · θ = 0 and θ · θ = φ · φ = 1 , where x · y = 0 τ 0 x · y d t . Asymptotic behavior of the row vector adjoint neutral modes θ = θ 1 θ 2 θ 3 θ 4 and φ = φ 1 φ 2 φ 3 φ 4 at sufficiently large negative times t < 0 is given by
θ 1 , 2 c 1 , 2 e λ 0 t , θ 3 c 3 e γ g t , θ 4 c 4 e γ q t ,
φ 1 , 2 d 1 , 2 e λ 0 t , φ 3 d 3 e γ g t , φ 4 d 4 e γ q t ,
where c 1 , 2 , 3 , 4 and d 1 , 2 , 3 , 4 are real coefficients, which can be calculated numerically. Similar to the leading tail of the TCS solution, the gain and absorption components of the trailing tail of the adjoint neutral modes decay faster than exponentially at t > 0 . Therefore, trailing tails of the adjoint neutral modes can be neglected in the derivation of the interaction equations. The temporal evolution of the field, gain, and loss components of the translational adjoint neutral mode θ = θ 1 θ 2 θ 3 θ 4 are shown in Figure 2. One can see that, similar to Figure 1, the gain component θ 3 of the adjoint neutral mode dominates almost everywhere between the consequent mode-locked pulses. Therefore, one can conclude that the pulse interaction via the field and absorption dynamics can be neglected for the parameter values of these figures.

4. Interaction Equations

To derive the equations describing slow evolution of the time coordinates and phases of weakly interacting TCSs, let us look for the solution of Equation (4) in the form of a sum of two unperturbed TCS solutions plus a small correction, w ( t ) = O ϵ , due to the interaction:
U = k = 1 2 u k + w ,
where u k = u 1 k u 2 k u 3 k u 4 k T with u 1 k + i u 2 k = u 1 t τ k + i u 2 t τ k e i ϕ k , u 3 k = u 3 t τ k , and u 4 k = u 4 t τ k . Coordinates τ k and phases ϕ k of the interacting TCSs are assumed to be slow functions of time, t τ k , t ϕ k = O ϵ , k = 1 , 2 . The small parameter ϵ characterizes weak overlap of the TCSs. Similar to the case of dissipative soliton interaction in partial differential equation laser models [9,25,26,27], the right-hand side of the interaction equations obtained from our DDE model can be expressed in terms of the TCS solutions and their adjoint neutral modes evaluated at the point between the two TCSs [21]. The details of the calculations are given in the Appendix A, where it is shown that the interaction equations for the time separation Δ τ = τ 2 τ 1 and phase difference Δ ϕ = ϕ 2 ϕ 1 of a pair of interacting T 0 -periodic TCS take the form
t Δ τ θ 1 T 0 / 2 u 2 T 0 / 2 θ 2 0 u 1 0 ,
t Δ ϕ φ 1 T 0 / 2 u 2 T 0 / 2 φ 2 0 u 1 0 ,
where θ k = θ ( t t k ) and φ k = φ ( t t k ) are the adjoint neutral modes evaluated at the kth TCS ( k = 1 , 2 ) and without loss of generality, one can assume that t = 0 and t = T 0 / 2 correspond, respectively, to the middle point between the two interacting TCSs and the opposite point on a circle with the circumference T 0 .
Substituting asymptotic Expressions (10)–(12) into the interaction Equations (14) and (15) and neglecting the field components we get
t Δ τ = K τ g e γ g T 0 Δ τ e γ g Δ τ + K τ q e γ q T 0 Δ τ e γ q Δ τ ,
t Δ ϕ = K ϕ g e γ g T 0 Δ τ e γ g Δ τ + K ϕ q e γ q T 0 Δ τ e γ q Δ τ
with K τ g = b 3 c 3 , K τ q = b 4 c 4 , K ϕ g = b 3 d 3 , and K ϕ q = b 4 d 4 .
Interaction Equations (16) and (17) describe the long-range interaction of two well-separated TCS via gain and absorption dynamics and do not take into account the short-range interaction due to weak overlap of the electric field envelopes of the TCSs. They reflect the fact that in a ring cavity, the interaction of two TCSs is twofold. The trailing tail of the first (second) TCS overlaps with the leading tail of the adjoint neutral mode of the second (first) TCS, which is located at distance Δ τ ( T 0 Δ τ ) behind it. This is reflected by the presence of the two exponential terms in the square brackets of Equations (16) and (17). In the case of TCS repulsion ( K τ g < 0 ), this type of twofold interaction leads to a regime with two equally spaced pulses per cavity round trip corresponding to a harmonic mode-locking regime. As I have already noted, due to the inequality γ q > γ g typical of semiconductor lasers, the interaction force related to absorption dynamics decays much faster than that due to gain dynamics and the terms proportional to K τ q and K ϕ q can be neglected in the interaction equations.

5. Results of Numerical Simulations

For the parameter values of Figure 1 and corresponding to zero linewidth enhancement factors, α g = α q = 0 , numerically, we get K τ g = 1.120 and K τ q = 2.145 in Equation (16). Due to the relation K ϕ g = K ϕ q = d 3 = d 4 = 0 , which is the consequence of ω 0 = A 0 = 0 , the second interaction Equation (17) transforms into t Δ ϕ = 0 . For negative values of K τ g the TCS interaction is repulsive, while positive K τ q corresponds to TCS attraction via absorption dynamics. However, as discussed above, for the parameter values of these figures, the interaction via gain dynamics dominates for almost all sufficiently large soliton separations, and the soliton attraction due to absorption dynamics is hardly possible to observe. The repulsive interaction is illustrated in Figure 3 obtained via numerical integration of Equations (1)–(3) using the RADAR5 code [28]. The initial condition was taken as a sum of two or more well-separated unperturbed TCSs. Figure 3a shows the standard mechanism of the harmonic mode-locking regime formation as a result of the repulsion of a pair of TCSs due to the interaction via gain dynamics. Figure 3b was obtained using the same parameter values as in Figure 3a, but with smaller initial separation of the two TCSs. It is seen that during the first stage of the interaction, there is still repulsion between the TCSs, but later, the second TCS becomes smaller and finally disappears. The equation t Δ ϕ = 0 implies that the TCS phase difference remains almost constant in the course of the interaction. This difference is affected only by a very weak overlap of the field components which were neglected in the derivation of the interaction Equations (16) and (17). Repulsive interaction of three and four TCSs leading to the development of harmonic mode-locking regimes with three and four pulses per cavity round trip is illustrated in Figure 3c and Figure 3d, respectively. Here, the transformation of a train of initially non-equidistant pulses into an equidistant pulse train clearly indicates the repulsive nature of the interaction. Note that in order to produce Figure 3 with minimized common TCSs drift, which is close to the single TCS drift T 0 T per cavity round trip, the following procedure was used. The intensity time trace obtained through numerical solution of Equations (1)–(3) was split into consecutive intervals of equal duration T 0 . Then, these intervals were associated with consecutive round trip numbers.
The dependence of the interaction coefficient K τ g on the linewidth enhancement factor in the gain section α g is shown in Figure 4a. It is seen that this dependence is non-monotonic and has a pronounced resonant character. The interaction coefficient is negative (TCS repulsion) when the linewidth enhancement factor is sufficiently small, and it becomes positive (TCS attraction) with the increase in α g , showing a sharp peak around α g 0.94 . A further increase in the α g leads to a non-monotonic gradual decrease in the interaction coefficient, which becomes negative again at α g 2.37 . The results of numerical simulations of Equations (1)–(3) with α g = 2.0 corresponding to a small positive values of the interaction coefficient are illustrated in Figure 5. It is seen that the interaction is very asymmetric; see refs. [15,21,26] and Appendix A. Figure 5b corresponding to q 0 = 4.0 and positive K τ g 0.854 × 10 2 shows the TCS attraction leading to the merging of two pulses when one of them disappears after collision. In Figure 5a obtained for q 0 = 5.0 and K τ g 1.076 × 10 2 the soliton attraction leads to a formation of a pulse bound state. Since for α g 0 , the relation d 3 = d 4 = 0 does not hold any more, the TCS phases are evolving with time in the course of interaction. Therefore, the bound state shown in Figure 5a is similar to the “incoherent” bound state described in [21] with the phase difference Δ ϕ between two pulses growing monotonically with time (see Figure 6 illustrating the intensity time trace and the evolution of the TCS phase difference of the incoherent bound state). It was demonstrated in [21] that due to the electric field overlap of the interacting TCS, such types of bound states are characterized by slightly oscillating time separation Δ τ . However, since the interaction via electric fields is extremely weak, for the bound state shown in Figure 5a, such oscillation is hardy possible to detect. Figure 4b shows the evolution of the inter-soliton time separation Δ τ as a function of the round trip number obtained via direct numerical simulation of the laser model (1)–(3). The parameter values are the same as in Figure 4a. It is seen that for α g = 0.5 , when the interaction coefficient K τ g is negative, the TCS interaction is repulsive and leads to a harmonic mode-locking regime. On the contrary, for α g = 1.0 , 1.5 , 2.0 , which correspond to K τ g > 0 , the interaction results in the formation of a pulse bound state. Furthermore, comparing Figure 4b with Figure 4a, we can see that the smaller the interaction coefficient, the weaker the interaction force and, hence, the longer the transient time before the equilibrium inter-pulse time separation is achieved. The final inter-pulse distance in the bound state is, however, almost independent of α g and K τ g . Note that the time separation of the pulses in the incoherent bound state shown in Figure 5a is of the same order of magnitude as the gain relaxation time. This is why the pulses in this bound state have significantly different peak powers (see Figure 6a) and cannot be considered as individual TCSs any more. Therefore, the interaction Equations (16) and (17) are not valid when the pulses are so close to one another. Indeed, in order for a bound state to be formed, the attraction predicted by the interaction equations should be compensated by a repulsion at sufficiently small inter-pulse distances. This repulsion dominating at small pulse separations might be related to the pulse repulsion in a laser with a cavity round trip time shorter or much shorter than the gain relaxation time discussed in [11,13].

6. Conclusions

The interactions of two well-separated TCSs in a long cavity semiconductor laser were studied numerically and analytically using the DDE mode-locking model. Interaction equations governing the slow evolution of the time separation and phase difference of the TCSs were derived and analyzed in the parameter range typical of semiconductor lasers, where the interaction via the gain saturation and recovery dominates over the interaction via absorption and field dynamics. Analytical results were compared to direct numerical simulations of the DDE mode-locking model. It was demonstrated that in addition to usual pulse repulsion predicted in [11,13], an attractive TCS interaction is also possible in a laser with a nonzero linewidth enhancement factor. This attractive interaction can result either in pulse merging or in a formation of an incoherent pulse bound state. In the latter case, the repulsion force counteracting the soliton attraction force might be attributed to the standard mechanism of the mode-locking pulse repulsion described in [11,13], which acts beyond the TCS limit. The incoherent bound pulse state discussed here is similar to that observed experimentally [29] and described theoretically [21] in a nonlinear mirror mode-locked laser. It is also similar to the “type A” pulse bound states reported in [12]. The mechanism of the latter bound states formation is, however, different from that described here and can be related to the TCS attraction due to the interaction via the absorption component of the pulsed solution in a laser with an infinitely large gain relaxation time.

Funding

This research was funded by the Deutsche Forschungsgemeinschaft (DFG projects No. 445430311 and No. 491234846).

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Derivation of the Interaction Equations

Substituting Equation (13) into Equation (4), collecting the first order terms in small parameter ϵ , and applying the solvability conditions [30] to the resulting equation yield
t τ k = θ k · P , t ϕ k = φ k · P ,
P = t u Σ + F ω 0 u Σ + H ω 0 u Σ t T ,
where u Σ = u 1 + u 2 , · = 0 T 0 · d t and θ k ( φ k ) is the adjoint translational (phase) neutral mode evaluated on u k , k = 1 , 2 .
Since u k is the solution of Equation (4), the equality
k = 1 2 t u k + F ω 0 u k + H ω 0 u k t T = 0
is satisfied. Subtracting this equality from (A2) we get
P = F ω 0 u Σ k = 1 2 F ω 0 u k + H ω 0 u Σ t τ k = 1 2 H ω 0 u k t T .
Therefore, the equation for τ 2 in (A1) is
t τ 2 = θ 2 · P = θ 2 · F ω 0 u Σ k = 1 2 F ω 0 u k + H ω 0 u Σ t T k = 1 2 H ω 0 u k t T .
Using T 0 periodicity of θ 2 and u 1 , 2 , Equation (A3) can be rewritten as
t τ 2 = θ 2 · F ω 0 u Σ k = 1 2 F ω 0 u k θ 2 t + T · H ω 0 u Σ k = 1 2 H ω 0 u k .
Next, we split the integral · = 0 T 0 · d t into two parts · = · 1 + · 2 , where · 1 = T 0 / 2 0 · d t and · 2 = 0 T 0 / 2 · d t are the integrals over the intervals T 0 / 2 , 0 and 0 , T 0 / 2 , respectively:
t τ 2 = j = 1 2 θ 2 · F ω 0 u Σ k = 1 2 F ω 0 u k j j = 1 2 θ 2 t + T · H ω 0 u Σ k = 1 2 H ω 0 u k j .
On the first interval T 0 / 2 , 0 , where u 2 is small, one obtains
F ω 0 u Σ F ω 0 u 1 B 1 u 2 , H ω 0 u Σ H ω 0 u 1 C 1 u 2 ,
and
F ω 0 u 2 B 0 u 2 , H ω 0 u 2 C 0 u 2
where B 1 = B u 1 and C 1 = C u 1 [ B 0 = B 0 and C 0 = C 0 ] are the linearization matrices of F ω 0 U and H ω 0 U at U = u 1 ( U = 0 ). Similarly, on the second interval 0 , T 0 / 2 , where u 1 is small, one gets
F ω 0 u Σ F ω 0 u 2 B 2 u 1 , H ω 0 u Σ H ω 0 u 1 C 2 u 1 ,
where B 2 = B u 2 and C 2 = C u 2 and
F ω 0 u 1 B 0 u 1 , H ω 0 u 2 C 0 u 1 .
Consequently, one obtains
t τ 2 θ 2 · B 1 B 0 u 2 1 θ 2 t + T · C 1 C 0 u 2 1 θ 2 · B 2 B 0 u 1 2 θ 2 t + T · C 2 C 0 u 1 2 ,
where the first two terms in the right hand side containing the product of two small quantities θ 2 and u 2 on the first interval T 0 / 2 , 0 can be neglected. This yields
t τ 2 θ 2 · B 2 B 0 u 1 2 θ 2 t + T · C 2 C 0 u 1 2 .
Since u 1 is the solution of Equation (4), it satisfies the equation t u 1 + F ω 0 u 1 + H ω 0 u 1 ( t T ) = 0 . Using the relations (A5) valid on the second interval 0 , T 0 / 2 , one can rewrite it in the form
t u 1 + B 0 u 1 + C 0 u 1 t T 0 .
The adjoint neutral mode θ 2 satisfies the equation
t θ 2 + θ 2 B 2 + θ 2 t + T C 2 = 0 .
Multiplying Equation (A7) by u 1 , subtracting from the resulting equation Equation (A6) multiplied by θ 2 , and integrating over the second interval 0 , T 0 / 2 , one gets
θ 2 · B 2 B 0 u 1 2 t θ 2 · u 1 + θ 2 · t u 1 2 θ 2 t + T C 2 · u 1 θ 2 · C 0 u 1 t T 2 .
Substitution of this relation into (A6) gives
t τ 2 t θ 2 · u 1 + θ 2 · t u 1 2 + θ 2 t + T C 2 · u 1 θ 2 · C 0 u 1 t T 2 θ 2 t + T · C 2 C 0 u 1 2 .
Finally, integrating the full derivative t θ 2 · u 1 over the interval 0 , T 0 / 2 leads to
t τ 2 θ 2 T 0 / 2 u 1 T 0 / 2 θ 2 0 u 1 0 + θ 2 t + T · C 0 u 1 θ 2 C 0 · u 1 t T 2 .
Note, that all the elements of the 4 × 4 matrix C 0 are equal to zero, except those in the upper 2 × 2 diagonal block. Hence, the last term in Equation (A8),
θ 2 t + T · C 0 u 1 θ 2 C 0 · u 1 t T 2 = 0 δ + T / 2 T / 2 + δ θ 2 t + T C 0 u 1 d t ,
contains only the asymptotic expressions for the field components, which are assumed to be small and are neglected in this study. Therefore, we can drop the last term in Equation (A8).
The equation for slow evolution of τ 1 is derived in a similar way to Equation (A8):
t τ 1 θ 1 0 u 2 0 θ 1 T 0 / 2 u 2 T 0 / 2 + θ 1 t + T · C 0 u 2 θ 1 C 0 · u 2 t T 1 .
Note, that the terms θ 2 T 0 / 2 u 1 T 0 / 2 and θ 1 0 u 2 0 in Equations (A8) and (A9), respectively, can be neglected due to the fast decay of the leading tail of the TCS solution and trailing edge of the adjoint neutral mode. The remaining terms, θ 2 0 u 1 0 in Equation (A8) and θ 1 T 0 / 2 u 2 T 0 / 2 in Equation (A9), have very different magnitudes except for the case, where the TCSs are nearly equidistant in the cavity, Δ τ = τ 2 τ 1 τ 0 / 2 . This means that except for this case the TCS interaction is strongly asymmetric and does not satisfy Newton’s third law [15,21,26]. Thus, keeping only the second terms in the right hand sides of Equations (A8) and (A9), one gets the following equation for the TCS time separation Δ τ = τ 2 τ 1 :
t Δ τ θ 1 T 0 / 2 u 2 T 0 / 2 θ 2 0 u 1 0 .
The equation for the slow evolution of the phase difference Δ ϕ = ϕ 2 ϕ 1 can be derived in a similar way. This equation reads:
t Δ ϕ φ 1 T 0 / 2 u 2 T 0 / 2 φ 2 0 u 1 0 .

References

  1. Grelu, P.; Akhmediev, N. Dissipative solitons for mode-locked lasers. Nat. Photonics 2012, 6, 84–92. [Google Scholar] [CrossRef]
  2. Herr, T.; Brasch, V.; Jost, J.D.; Wang, C.Y.; Kondratiev, N.M.; Gorodetsky, M.L.; Kippenberg, T.J. Temporal solitons in optical microresonators. Nat. Photonics 2014, 8, 145–152. [Google Scholar] [CrossRef] [Green Version]
  3. Kippenberg, T.J.; Gaeta, A.L.; Lipson, M.; Gorodetsky, M.L. Dissipative Kerr solitons in optical microresonators. Science 2018, 361, 6402. [Google Scholar] [CrossRef] [Green Version]
  4. Marconi, M.; Javaloyes, J.; Balle, S.; Giudici, M. How lasing localized structures evolve out of passive mode locking. Phys. Rev. Lett. 2014, 112, 223901. [Google Scholar] [CrossRef] [Green Version]
  5. Rosanov, N.; Khodova, G. Diffractive autosolitons in nonlinear interferometers. J. Opt. Soc. Am. B 1990, 7, 1057–1065. [Google Scholar] [CrossRef]
  6. Malomed, B.A. Bound states of envelope solitons. Phys. Rev. E 1993, 47, 2874. [Google Scholar] [CrossRef] [PubMed]
  7. Akhmediev, N.; Ankiewicz, A.; Soto-Crespo, J. Multisoliton solutions of the complex Ginzburg-Landau equation. Phys. Rev. Lett. 1997, 79, 4047. [Google Scholar] [CrossRef] [Green Version]
  8. Akhmediev, N.; Rodrigues, A.S.; Town, G.E. Interaction of dual-frequency pulses in passively mode-locked lasers. Opt. Commun. 2001, 187, 419–426. [Google Scholar] [CrossRef]
  9. Vladimirov, A.G.; Khodova, G.V.; Rosanov, N.N. Stable bound states of one-dimensional autosolitons in a bistable laser. Phys. Rev. E 2001, 63, 056607. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Ablowitz, M.J.; Horikis, T.P.; Nixon, S.D. Soliton strings and interactions in mode-locked lasers. Opt. Commun. 2009, 282, 4127–4135. [Google Scholar] [CrossRef] [Green Version]
  11. Kutz, J.N.; Collings, B.; Bergman, K.; Knox, W. Stabilized pulse spacing in soliton lasers due to gain depletion and recovery. IEEE J. Quantum Electron. 1998, 34, 1749–1757. [Google Scholar] [CrossRef] [Green Version]
  12. Soto-Crespo, J.M.; Akhmediev, N.N. Multisoliton regime of pulse generation by lasers passively mode locked with a slow saturable absorber. J. Opt. Soc. Am. B 1999, 16, 674–677. [Google Scholar] [CrossRef] [Green Version]
  13. Nizette, M.; Rachinskii, D.; Vladimirov, A.; Wolfrum, M. Pulse interaction via gain and loss dynamics in passive mode locking. Phys. D 2006, 218, 95–104. [Google Scholar] [CrossRef]
  14. Zaviyalov, A.; Grelu, P.; Lederer, F. Impact of slow gain dynamics on soliton molecules in mode-locked fiber lasers. Opt. Lett. 2012, 37, 175–177. [Google Scholar] [CrossRef]
  15. Camelin, P.; Javaloyes, J.; Marconi, M.; Giudici, M. Electrical addressing and temporal tweezing of localized pulses in passively-mode-locked semiconductor lasers. Phys. Rev. A 2016, 94, 063854. [Google Scholar] [CrossRef] [Green Version]
  16. Vladimirov, A.G.; Turaev, D. Model for passive mode locking in semiconductor lasers. Phys. Rev. A 2005, 72, 033808. [Google Scholar] [CrossRef] [Green Version]
  17. Vladimirov, A.G.; Turaev, D.; Kozyreff, G. Delay differential equations for mode-locked semiconductor lasers. Opt. Lett. 2004, 29, 1221–1223. [Google Scholar] [CrossRef]
  18. Vladimirov, A.G.; Turaev, D. New model for mode-locking in semiconductor lasers. Radiophys. Quant. Electron. 2004, 47, 857–865. [Google Scholar] [CrossRef]
  19. Puzyrev, D.; Vladimirov, A.G.; Pimenov, A.; Gurevich, S.V.; Yanchuk, S. Bound pulse trains in arrays of coupled spatially extended dynamical systems. Phys. Rev. Lett. 2017, 119, 163901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Munsberg, L.; Javaloyes, J.; Gurevich, S.V. Topological localized states in the time delayed Adler model: Bifurcation analysis and interaction law. Chaos 2020, 30, 063137. [Google Scholar] [CrossRef]
  21. Vladimirov, A.G. Short-and long-range temporal cavity soliton interaction in delay models of mode-locked lasers. Phys. Rev. E 2022, 105, 044207. [Google Scholar] [CrossRef]
  22. Vladimirov, A.G.; Kovalev, A.V.; Viktorov, E.A.; Rebrova, N.; Huyet, G. Dynamics of a class-A nonlinear mirror mode-locked laser. Phys. Rev. E 2019, 100, 012216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Yanchuk, S.; Ruschel, S.; Sieber, J.; Wolfrum, M. Temporal dissipative solitons in time-delay feedback systems. Phys. Rev. Lett. 2019, 123, 053901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Stöhr, M.; Wolfrum, M. Temporal dissipative solitons in the Morris–Lecar model with time-delayed feedback. Chaos 2023, 33, 023117. [Google Scholar] [CrossRef]
  25. Turaev, D.; Vladimirov, A.G.; Zelik, S. Long-range interaction and synchronization of oscillating dissipative solitons. Phys. Rev. Lett. 2012, 108, 263906. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Vladimirov, A.G.; Gurevich, S.V.; Tlidi, M. Effect of Cherenkov radiation on localized-state interaction. Phys. Rev. A 2018, 97, 013816. [Google Scholar] [CrossRef] [Green Version]
  27. Vladimirov, A.G.; Tlidi, M.; Taki, M. Dissipative soliton interaction in Kerr resonators with high-order dispersion. Phys. Rev. A 2021, 103, 063505. [Google Scholar] [CrossRef]
  28. Guglielmi, N.; Hairer, E. Implementing Radau IIA methods for stiff delay differential equations. Computing 2001, 67, 1–12. [Google Scholar] [CrossRef]
  29. Kokhanovskiy, A.; Kuprikov, E.; Kobtsev, S. Single-and multi-soliton generation in figure-eight mode-locked fibre laser with two active media. Opt. Laser Technol. 2020, 131, 106422. [Google Scholar] [CrossRef]
  30. Halanay, A. Differential Equations: Stability, Oscillations, Time Lags; Academic Press: New York, NY, USA, 1966; Volume 23. [Google Scholar]
Figure 1. Intensity time trace of the periodic TCS solution of Equations (1)–(3) (a). Temporal evolution of the logarithm of the absolute values of the field envelope (blue), gain (red) and loss (green) components of the TCS solution (b). The parameters are: α g = α q = 0 , g 0 = 0.5 , q 0 = 4.0 , κ = 0.8 , s = 10.0 , γ = 5.0 , γ g = 0.2 , γ q = 1.0 , and T = 50.0 . The solution period is T 0 = 50.138425 and ω 0 = 0 .
Figure 1. Intensity time trace of the periodic TCS solution of Equations (1)–(3) (a). Temporal evolution of the logarithm of the absolute values of the field envelope (blue), gain (red) and loss (green) components of the TCS solution (b). The parameters are: α g = α q = 0 , g 0 = 0.5 , q 0 = 4.0 , κ = 0.8 , s = 10.0 , γ = 5.0 , γ g = 0.2 , γ q = 1.0 , and T = 50.0 . The solution period is T 0 = 50.138425 and ω 0 = 0 .
Optics 04 00031 g001
Figure 2. Field intensity of the adjoint neutral mode θ as a function of time (a). Temporal evolution of the logarithm of the absolute values of the field (blue), gain (red), and loss (green) components of the translational adjoint neutral mode θ (b). Parameters are the same as in Figure 1.
Figure 2. Field intensity of the adjoint neutral mode θ as a function of time (a). Temporal evolution of the logarithm of the absolute values of the field (blue), gain (red), and loss (green) components of the translational adjoint neutral mode θ (b). Parameters are the same as in Figure 1.
Optics 04 00031 g002aOptics 04 00031 g002b
Figure 3. TSC repulsion due to the interaction via gain dynamics leading to a harmonic mode-locking regimes with two (a), three (c), and four (d) pulses per cavity round trip. Colorful lines indicate the trajectories of the interacting TCSs. Panel (b) illustrates repulsive interaction resulting in the annihilation of the second pulse. (ac)— g 0 = 0.5 . (d)— g 0 = 0.8 . Cold cavity round trip time is T = 50 . Other parameters are the same as in Figure 1.
Figure 3. TSC repulsion due to the interaction via gain dynamics leading to a harmonic mode-locking regimes with two (a), three (c), and four (d) pulses per cavity round trip. Colorful lines indicate the trajectories of the interacting TCSs. Panel (b) illustrates repulsive interaction resulting in the annihilation of the second pulse. (ac)— g 0 = 0.5 . (d)— g 0 = 0.8 . Cold cavity round trip time is T = 50 . Other parameters are the same as in Figure 1.
Optics 04 00031 g003
Figure 4. Interaction coefficient K τ g as a function of α g (a) and pulse time separation as a function of the round trip number (b). Different curves in panel (b) correspond to different linewidth enhancement factor values, α g = 0.5 , 1.0 , 1.5 , 2.0 , as indicated in the panel. g 0 = 0.8 . Other parameters are the same as in Figure 1.
Figure 4. Interaction coefficient K τ g as a function of α g (a) and pulse time separation as a function of the round trip number (b). Different curves in panel (b) correspond to different linewidth enhancement factor values, α g = 0.5 , 1.0 , 1.5 , 2.0 , as indicated in the panel. g 0 = 0.8 . Other parameters are the same as in Figure 1.
Optics 04 00031 g004
Figure 5. TCS interaction resulting in pulse bound state formation at q 0 = 5.0 (a) and pulse merging at q 0 = 4.0 (b). Colorful lines indicate the trajectories of the interacting TCSs. g 0 = 0.8 , α g = 2.0 , α q = 0 . The cold cavity round trip time is T = 50 . Other parameters are as shown in Figure 3.
Figure 5. TCS interaction resulting in pulse bound state formation at q 0 = 5.0 (a) and pulse merging at q 0 = 4.0 (b). Colorful lines indicate the trajectories of the interacting TCSs. g 0 = 0.8 , α g = 2.0 , α q = 0 . The cold cavity round trip time is T = 50 . Other parameters are as shown in Figure 3.
Optics 04 00031 g005
Figure 6. Intensity time-trace (a) and pulse phase difference (b) of the TCS bound state. Parameters are the same as in Figure 5a.
Figure 6. Intensity time-trace (a) and pulse phase difference (b) of the TCS bound state. Parameters are the same as in Figure 5a.
Optics 04 00031 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vladimirov, A.G. Temporal Cavity Soliton Interaction in Passively Mode-Locked Semiconductor Lasers. Optics 2023, 4, 433-446. https://doi.org/10.3390/opt4030031

AMA Style

Vladimirov AG. Temporal Cavity Soliton Interaction in Passively Mode-Locked Semiconductor Lasers. Optics. 2023; 4(3):433-446. https://doi.org/10.3390/opt4030031

Chicago/Turabian Style

Vladimirov, Andrei G. 2023. "Temporal Cavity Soliton Interaction in Passively Mode-Locked Semiconductor Lasers" Optics 4, no. 3: 433-446. https://doi.org/10.3390/opt4030031

APA Style

Vladimirov, A. G. (2023). Temporal Cavity Soliton Interaction in Passively Mode-Locked Semiconductor Lasers. Optics, 4(3), 433-446. https://doi.org/10.3390/opt4030031

Article Metrics

Back to TopTop