Temporal Solitons in an Optically Injected Kerr Cavity with Two Spectral Filters

: We propose and analyze a theoretical scheme of an injected Kerr cavity, where the chromatic dispersion is induced by propagation of light through two Lorentzian spectral ﬁlters with different widths and central frequencies. We show that this setup can be modeled by a second order delay differential equation that can be considered as a generalization of the Ikeda map with included spectral ﬁltering, dispersion, and coherent injection terms. We demonstrate that this equation can exhibit modulational instability and bright localized structures formation in the anomalous dispersion regime.


Introduction
Optical frequency combs generated from a continuous wave laser output in microcavity Kerr resonators have revolutionized many fields of natural science and technology [1]. Of particular interest are the so-called soliton frequency combs [2][3][4] associated with the formation in the time domain of the temporal cavity solitons-nonlinear localized structures of light, which preserve their shape in the course of propagation [5]. Kerr cavity temporal dissipative solitons were reported experimentally in micro-cavity resonators [2,[6][7][8], and in driven fiber cavities [9]. Spectral filtering is commonly used to improve the output characteristics of multisection mode-locked lasers [10]. It was also demonstrated that a small spectral filtering can suppress the oscillatory instability of Kerr cavity solitons and stabilize their bound states by eliminating high frequency perturbations [11][12][13][14]. On the other hand, unlike conventional Kerr cavities, where spectral filtering is only a small perturbation, in Mamyshev oscillators a pair of spectral filters plays a crucial role in the process of the short pulse generation [15,16]. The nonlinear dynamics of such two filter systems is yet to be understood theoretically. In order to fill this gap we consider an externally driven Kerr cavity with two spectral filters and demonstrate that the effective dispersion created by the filters can lead to a temporal soliton generation even when the material dispersion is negligibly small.
Time-delay models of optical systems like Ikeda map [17], Lang-Kobayashi equations [18], delay differential equation (DDE) mode-locked laser models [19][20][21][22][23][24][25][26][27], frequency swept laser models [28,29], and others were successfully used to describe unidirectional propagation of light through linear and nonlinear optical elements in a ring cavity. Unlike simple complex Ginzburg-Landau (CGL)-type models of passive and active optical cavities based on the mean field approximation, these models are valid for arbitrary large gain and losses in the cavity. An important drawback of the time-delay models, however, is that the inclusion of an arbitrary second order chromatic dispersion of the intracavity media into these models is not a trivial task. Recently, it was shown that chromatic dispersion in photonic crystal mode-locked laser [30] and a SOA-fiber laser with fiber delay line [31] can be described using a distributed delay term, which arises from the transfer function of a detuned Lorentzian absorption line in frequency domain. Furthermore, under assumption of weak dispersion one can replace the distributed time delay model with an extended DDE model containing a single additional ordinary differential equation for the polarization variable [32]. Using this extended DDE model, the conventional combined effects of chromatic dispersion and nonlinearity such as modulational instability (MI) in the anomalous dispersion regime and bright localized structures formation were demonstrated [32]. Nevertheless, these assumptions and approximations limit our ability to describe accurately and characterize chromatic dispersion at all the frequencies that are important for dynamics of optical devices using DDE models. Another approach is to investigate rigorously derived DDEs where higher order dispersion arises naturally (e.g., coupled cavities [33]), however quantification of its magnitude for a given DDE is not a trivial task.
In this paper, we consider an externally injected ring Kerr cavity with two linear spectral filters introducing an effective chromatic dispersion and demonstrate the possibility to model arbitrary second-order dispersion near a chosen frequency within the DDE framework. We develop a second order DDE model of the system under consideration and demonstrate the appearance of MI and the formation of bright localized structures in this model in the anomalous dispersion regime. The system under consideration can be realized experimentally to generate temporal cavity solitons and the corresponding optical frequency combs [6,34]. We show that in a certain limit our model can be reduced to a generalized version of the well-known Lugiato-Lefever equation [35], which is widely used to describe optical microcomb generation [36][37][38], with an additional diffusion term. These results can be applied to qualitatively analyze any optical set-up that can be modeled using delay equations such as Fourier domain mode-locked [31], optically injected [32], and multisection mode-locked semiconductor lasers as well as any other system, where second-order chromatic dispersion is important (see references in [31]). We note also that two spectral filters with different central frequencies are used in Mamyshev oscillators, which employ active cavity to generate short optical pulses. Hence, this work not only presents the simplest dispersive second-order DDE model that can describe complicated phenomena like localized structures, but also could provide better understanding of the effect of two filters in the cavity and a basis for the theoretical investigation of Mamyshev oscillators [15,16].

Model Equations
We consider an optically injected passive nonlinear cavity with two linear Lorentzian filters inside it (see Figure 1). For this system using the lumped element method described in [19][20][21] we obtain the following set of delay differential equations where B and A represent electrical field envelopes after the first and the second filter, respectively, T is the cavity round-trip time, κ l is the intensity attenuation factor due to the linear cavity losses, ϕ is the phase shift, α is the Kerr coefficient which, without the loss of generality can be assumed to be positive, η 0 and ω 0 are the injection rate and frequency offset, ω 1 and ω 2 are the central frequencies of the two filters, while γ 1 and γ 2 are their bandwidths. The system (1) and (2) can be considered as an extension of the Ikeda map [17], which takes into consideration the spectral filtering introduced by two Lorentzian filters.

Transfer Function of the Filter
The transfer function of the two filters shown in Figure 1 can be written in the standard form [30,31] F where the complex function f (ω) can be expanded in power series near ω = 0 and D k represents the dispersion of the kth order. ( For example, for a single Lorentzian absorption line we have f =

−σL
Γ+i(ω+Ω) [31], and second-order dispersion coefficient takes the form D 2 = d 2 f dω 2 ≈ 2σL Ω 3 for 0 < Γ |Ω|, and the sign of the coefficient coincides with the sign of Ω.) From Equations (1) and (2) one obtains that the two filter transfer function is Since in the frequency domain the transfer functions of two filters are multiplied, see Equation (5), their ordering does not play any role. Therefore, without loss of generality we can assume that γ 1 ≥ γ 2 and the reference frequency is chosen in such a way that the maximum of the function |F(ω)| is at ω = 0, so that d|F(ω)| dω ω=0 = 0. The latter condition is equivalent to which has two solutions One can see that both ω 2± can be obtained for any values of ω 1 = 0 and 0 < γ 2 < γ 1 due to D F > 0, and, moreover, |ω 2 1 − γ 2 1 | ≤ D F ≤ γ 2 1 + ω 2 1 . Below we will assume that the condition ω 2 = ω 2+ is satisfied that corresponds to a situation when the largest of the two maximums of |F(ω)| is located at zero frequency, see Figure 2. In the left panel of this figure we fix γ 1 = γ 2 = 1 and demonstrate that for ω 1 < 1 we have one global maximum, whereas for ω 1 > 1 we can get two local maxima in our combined filter of the same magnitude. In the second case, or for values of γ 2 close to γ 1 , this kind of filter can lead to complex dynamical behaviour in the system, similarly to the case of the Mamyshev oscillator [15,16]. However, in this paper we are focused on the impact of the filter on the frequency of injection, which is assumed to be near or at the global maximum at zero frequency. In the right panel we can observe how the second maximum is suppressed with γ 2 decreasing from 1 to 1 8 . For ω 2 = ω 2+ we have and one can see that i.e., the difference between the two frequencies grows with |ω 1 |, and since for any of ω 1 = ±ω * 1 we have ω 2 = ∓ω * 2 , respectively, for any fixed 0 < γ 2 < γ 1 we can choose |ω 2 − ω 1 | and find corresponding value of ω 1 (negative or positive) and ω 2 = ω 2+ , which means that we can have any possible combinations of two filters (narrow filter to the left or to the right of the broad filter), disregarding frequency shift of the combined filter. Moreover, these inequalities immediately imply that the largest maximum of |F(ω)| is at ω = 0 for ω 2 = ω 2+ as claimed earlier, since more narrow Lorentzian filter with the width γ 2 has the central frequency ω 2 closer to zero frequency. The other root (6) ω 2 = ω 2− can correspond to a global maximum, a local maximum, or a local minimum, which do not provide any additional useful alternatives, if we are interested to fix the maximum of the combined filter |F(ω)| at zero frequency ω = 0.  Since f (ω) = ln F(ω), in the power series expansion (4) we obtain the second-order coefficient as One can check that D 2 can be represented in the following way with It follows from the conditions (6), where ω 2 = ω 2+ , and the relation γ 1 ≥ γ 2 that the second term in the right hand side (RHS) of (7) vanishes, ρ > 0, and signβ = signω 1 . Hence, the parameters β and σ 2 −2 2 represent the normalized second-order dispersion and the diffusion near zero frequency, respectively, whereas the parameter √ ρ can be considered as a scaling coefficient.

Normalized DDE Model
By substituting Equation (2) into Equation (1), using Equation (6), and rescaling time as t → t/ √ ρ, we obtain the following normalized dimensionless second-order DDE where losses and forcing take the form Here, similarly to the previous section one can see that dr d|ω 1 | < 0 for γ 1 > γ 2 and ω 2 = ω 2+ , hence for fixed κ l the parameter κ decreases with increasing detuning between the two frequencies |ω 2 − ω 1 |. With the condition (11) on dimensionless parameters σ and β the equation is equivalent to the system with two filters (1) and (2) for any κ < 1 with global maximum of the combined filter transfer function F(ω) (5) located at zero frequency ω = 0, |F(ω)| ≤ |F(0)| ≤ 1. We note that for √ 2 ≤ σ ≤ 2(1 + β 2 ) the global maximum of F(ω) is at another frequency and additional constraints on parameter κ are necessary, whereas for σ < √ 2 the diffusion coefficient in (7) becomes negative, equivalence with (1) and (2) is lost and zero solution A = 0 for η = 0 is unstable.
For β = 0 relation (7) gives zero dispersion coefficient and diffusion coefficient equal to (σ 2 − 2)/2, which agrees with previous analysis of a similar DDE [39] in case of σ > √ 2. We note that the direct application of truncated expansion (4) in the RHS of DDE models similar to (1) and (12) would lead to the appearance of the second derivative of the delayed variable A (t − T) in and spurious instability [31]. In contrast, non-delayed second order derivative A (t) in (12) appears without any expansions, and the term iβA (t) is responsible for the second-order chromatic dispersion similarly to the CGL equations. On the other hand, unlike the cubic CGL model, where the real coefficient by the second derivative is responsible for the diffusion (or parabolic spectral filtering), the filtering in (12) is performed by two Lorentzian filters, which are introduced by the presence of of both the first and the second order derivatives. Therefore, the generality of this kind of dispersion operator is still not directly comparable to the operators in CGL-type equations or systems with distributed delay [31]. However, in contrast to the approximate DDEs [32], the physical meaning of this filter is clear for any feasible parameters. For simplicity below we consider the case of non-detuned injection, ω 0 = 0.

Limit of Lugiato-Lefever Equation
Assuming large delay limit T = r/ε, where ε 1 and r = 1 + a 1 ε + a 2 ε 2 + . . . , taking ω 0 = 0, and rescaling the time variable t → t/ε we can rewrite (12) in the form Let the injection rate, the linear cavity losses and the phase shift be small, Then looking the solution in the form A(t) = εu(t, τ) + ε 2 v(t, τ) + . . . with τ = ε 2 t, applying multiscale analysis [40,41], collecting the first order in ε, and using solvability condition of the resulting equation we get the periodic boundary condition: Next, in the second order in ε we obtain which implies the periodicity of v, v(t − 1, τ) = v(t, τ) and the relation a 1 = σ with the periodicity condition (15). Finally collecting the third order terms in ε, using the relation a 2 = σ 2 , and applying solvability condition [40,41] we get the generalized Lugiato-Lefever equation (LLE) with the diffusion coefficient (σ 2 − 2)/2 and the boundary condition (15). It is well known that for β < 0 in Equation (16) corresponds to anomalous dispersion regime, and for θ < 0 this equation can demonstrate the formation of bright dissipative solitons [35]. Dimensionless dispersion and diffusion coefficients in Equation (16) coincide with those defined by Equation (7).

Continuous Wave (CW) State
The CW state of the Equation (12) with ω 0 = 0 takes the form A = A 0 e iψ , where the real quantities A 0 and ψ satisfy the system of the transcendental equations which leads to a single transcendental equation for Assuming which coincides with the equation for the uniform stationary solutions of the LLE (16). Therefore, in this limit Equation (12) can have up to three coexisting CW states similarly to LLE (see Figure 3, left), however out of this limit for strong injection there can be more coexisting CW states, see right panel of Figure 3. Here, the upper CW state looses stability via a modulational instability, and unstable CW states are shown by dotted red line.  (17) and (18) obtained by varying parameter η (left). Here, , which in the LLE limit (13) corresponds to ε = 0.01. In the right panel we choose κ = exp(−2), ϕ = −2, which corresponds to ε = 1, where solid lines correspond to stable CWs and dashed -to unstable.

Stability of CW Solution and MI
Here, we demonstrate how MI of an initially stable CW state can appear in the anomalous dispersion regime in the limit of large delay T 1. For that, we linearize the Equation (12) near the CW state A(t) = (A 0 + δAe λt )e iψ and calculate the determinant of the Jacobian of the linearised system to obtain the following characteristic equation for the eigenvalues λ describing the stability of the CW solution: where Y(λ) = exp(−λT). In the limit of large delay time T → ∞ the eigenvalues belonging to the pseudo-continuous spectrum can be represented in the form λ = iµ + Λ with real µ [42]. Thus, in this limit the characteristic equation is a quadratic equation for Y with the coefficients depending only on the imaginary part µ of the eigenvalue, and we can obtain from (21) two branches of pseudo-continuous spectrum given by For a stable CW solution we have Λ ± (µ) ≤ 0 and, in particular, Moreover, the second derivative of Λ ± (µ) at µ = 0 takes the form One can see from (25) that for small β ≈ 0 and σ > √ 2 we have Λ ± (0) < 0 for the CW solutions on the upper branch of S-shaped bifurcation curve depicted in Figure 3 (see top-left panel of Figure 4). It is known that strong anomalous dispersion can lead to the change of the curvature of one of the two branches Λ ± (µ) of the pseudo-continuous spectrum at µ = 0 (see Figure 4, bottom-left) [31]. Moreover, further increase of anomalous dispersion can lead to a MI of a CW state (see Figure 4, bottom-right) [32]. On the other hand, in the case of strong normal dispersion regime the sign change of the curvature of the curve with smaller Λ ± (0) does not lead to a MI, as it is seen from the top-right panel of Figure 4.  The stability condition of the CW state can be written in the form |Y ± (iµ)| > 1 for all the wavenumbers µ. In particular, at zero wavenumber µ = 0 we get A CW solution satisfying this condition is stable with respect to perturbations at zero wavenumber, but may be unstable with respect to MI at nonzero wavenumbers. A possible (but not unique) way how such a MI can develop is related to the change of the sign from negative to positive of one of the two quantities Λ ± (0) corresponding to the greater of the two values Λ ± (0). Note, that both the quantities Λ ± (0) are always negative at β = 0 due to the inequality (11). It follows from Equation (26) that such a sign change can take place only in the case when D > 0 and, hence, Y ± (0) must be real. Otherwise, B is purely imaginary and the last term in the RHS of Equation (25) vanishes. It is shown in Appendix A that in the LLE limit (14) we have Λ + (0) < Λ − (0) and for the CW state with the highest intensity we can have Λ − (0) > 0 only in the anomalous dispersion regime β < 0 . Out of LLE limit we can have positive Λ ± (0) in case of β > 0 as well (see Appendix B for numerical treatment). Note, that the development of MI on the stable upper part of the CW branch can be correlated with the appearance of stable localized structures. In the next section we study numerically how the dispersion coefficient β affects the existence range of these structures. Remarkably, there are also scenarios out of the LLE limit where localized structures can be observed for zero and small positive β.
We note that the condition Λ (0) > 0 precedes the appearance of MI of the CW for stronger dispersion if Y ± (0) is real-valued. If it is complex-valued, MI is still possible for sufficiently strong dispersion (see Figure 5). In the LLE limit, one can derive an approximate condition for MI development (see Appendix A), and demonstrate that for typical parameters in this limit strong anomalous dispersion can lead to the instability.  Figure 5. Curves of pseudocontinuous spectrum in the limit of large delay (21) of (12) of the CW with the largest intensity (|A 0 | ≈ 0.0178 out of another two |A 0 | ≈ 0.0156 and |A 0 | ≈ 0.0065) near the LLE limit with β = 1/2 (left) and β = 1/2 (right), and η ≈ 1.797 × 10 −6 . Other parameters are as in Figure 3

Numerical Results
In this section, we perform numerical bifurcation analysis of Equation (12) using DDE-BIFTOOL [43]. Let us start with the parameter set close to the LLE limit. We take σ = √ 2.5, β = −0.5, which corresponds to two filters (5) with γ 1 = γ 2 , ω 1 = − 1+ are chosen according to analysis of DDE (13) in the limit ε → 0. These values directly correspond to the parameters of LLE (16) where localized structures are known to exist (θ = ϕ/ε 2 < − √ 3). We also pick a reference injection strength where three CW states are guaranteed to coexist for ε 1, and vary η near η ε . The main numerical difficulty in this limit is that the delay time T and the width of the localized structures are proportional to the large quantity ε −1 . Figure 6 was obtained with ε = 0.05, which was the smallest value of ε used in our simulations and corresponded to κ ≈ 0.995. Similarly to the LLE [44,45]   Thus, close to the limit ε → 0 the normalized second-order DDE (12) demonstrates the bifurcation structure similar to that of the LLE. Note, however, that the magnitude of the parameters σ, β and 1 − κ in the normalized Equation (12) depends on the frequency detuning of the linear filters with respect to each other in the original equivalent system (1) and (2), and stronger detuning results for fixed σ in larger |β| as discussed in Section 2.1, and higher losses 1 − κ at the same time as discussed in Section 2.2. Therefore, it is necessary to investigate how larger values of ε corresponding to smaller values of the attenuation factor κ (larger losses) affect the properties of the localized solutions. For example, for the considered parameters σ = √ 2.5 and β = −0.5 from Equations (1) and (2) and the condition κ l ≤ 1 one obtains κ ≤ 0.4 in (12), which is satisfied for ε ≥ 0.7.
For larger ε = 0.2, which corresponds to κ ≈ 0.92, the stable part of the periodic solution branch is split into two parts separated by two fold bifurcation points, see left panel of Figure 7. Both these parts correspond to very asymmetric localized pulses, but the second part contains wider pulses then the first one, see right panel of Figure 7. Figures 8 and 9 show the branches of CW and periodic solutions with scaled intensity (In these figures η is shifted, more precisely,η = η − η ε − η 1,j + 10 −5 j for the left panel of Figure 8,η = η − η ε − η 2,j + 5jη 2,1 for the right panel,η = η − η ε − η 1,j + 10 −2 (j − 4) for the left panel of Figure 9,η = η − η ε − η 2,j + 2 ×  Other parameters are as in Figure 8.
Finally, let us study how the variation of the effective dispersion coefficient β influences the existence range of the temporal localized structures. One can see in the top left panel of Figure 11 that near the LLE limit (ε = 0.05) the localized solutions can be observed only for negative β < −0.05 corresponding to the anomalous dispersion regime. It is seen that when β increases the interval of the injection rates between two fold bifurcations, where stable localized solutions exist, shrinks so that for β ≥ −0.06 one can hardly see any stable localized structure in the region of S-shaped CW curve. However, by increasing ε out of the LLE limit first to ε = 0.3 (top right panel of Figure 11) and then to ε = 0.5 (bottom left panel), one can see that around ε ≈ 0.5 localized structures can be observed for β = 0 as well. Furthermore, for 0.5 ≤ ε ≤ 1 (bottom right panel) stable bright temporal dissipative solitons can be observed even with small positive β, 0 ≤ β < 0.1. Therefore, we conclude that out of the LLE limit temporal localized structures in the DDE model (12) could be observed not only in anomalous dispersion regime but also for positive β, although in a smaller range of the injection rates η. We have considered here only one possible way to exit the LLE limit in a continuous manner, however our preliminary numerical simulations suggest that there are many possible combinations of values of κ, η, ϕ, and β away from the LLE limit, where stable or unstable temporally localized structures can be observed. We leave this question to further studies, and, in particular, for experimentally justified parameters. Other parameters are as in Figure 8.
In this section we have demonstrated the existence of bright localized structures of the second oder DDE model (12). It follows from this result that such kind of structures should also exist in the original set of two DDEs (1) and (2). We can use this theory then to find parameters for (1) and (2) where localized structures exist (see Figure 12).

Conclusions
We have considered a DDE model of an optically injected ring Kerr cavity with two spectral filters having different widths and central frequencies. We have derived a normalized complex-valued second order DDE (12), which is similar to real-valued DDE reported in [39]. This equation can be considered as a generalization of the Ikeda map, which explicitly contains second-order dispersion coefficient at zero frequency as a parameter. We have derived an admissibility relation for these parameters, analyzed stability of the CW solutions of this model in the limit of large delay, and demonstrated the effect of strong dispersion on the development of modulational instability. We have shown that in the limit of small losses and weak injection the DDE model can be reduced to a generalized version of the well known LLE, which is known to have dissipative soliton solutions. We have performed numerical bifurcation analysis of CW solutions and temporal localized structures using DDE-BIFTOOL package [43] and demonstrated qualitative similarity of the solutions of the DDE model with those of the LLE in the regime of anomalous dispersion (β < 0) in the corresponding limit. Moreover, out of LLE limit one can observe numerically stable asymmetric bright localized structures not only in the anomalous dispersion regime, but also at zero and small positive values of β, albeit in a shrinkingly smaller interval of existence. This is in agreement with the results of Ref. [46] obtained using an extended LLE with the third order dispersion term included. Experimental observations of asymmetric temporal Kerr cavity solitons near zero dispersion point, where higher order dispersion comes into play, were reported in [47].
To summarize, we have proposed theoretically a Kerr cavity scheme, where the dispersion introduced by two spectral filters with different central frequencies can lead to the development of modulational instability and appearance of stable temporal cavity solitons even when the material dispersion is negligible. This scheme is different from the traditional Kerr cavity setups for frequency comb generation, where chromatic dispersion of the intracavity material plays a crucial role in soliton formation. The modeling approach we have developed is also suitable for the analysis of more complex systems with two spectral filters, such as Mamyshev oscillators. Investigation of short pulse generation in Mamyshev oscillators using the theoretical basis developed here might be a subject for future research. Acknowledgments: Not applicable.

Conflicts of Interest:
The authors have no conflicts of interest to declare.

Appendix A. MI in the LLE Limit
Linear stability of the CW solutions of the classical Lugiato-Lefever equations was studied alalytically and numerically in a number of works, see e.g., Refs. [48][49][50] for both the normal and anomalous dispersion regimes. In this Appendix we consider the DDE model (12) in the LLE limit, where this model can be reduced to the generalized LLE (16) with the additional diffusion (spectral filtering) term (σ 2 − 2)/2, which exerts a stabilizing effect on the CW solutions. A detailed linear stability analysis of the generaized LLE model (12) is beyond the scope of this paper. We present here only some analytical results concerning the sign change of one of the two quantities Λ ± (0) (see Equation (25)), which can be a precusor of the MI.
Without the loss of generality we can assume that α = 1 and k = 1. Substituting (14) and (25) and using Λ ± (0) = − ln Y ± (0) we obtain: where Λ ± (0) can be positive only whenD > 0, i.e., for − θ 3 < u 2 0 < −θ. We see from (A1) and (A2) that forD > 0 in the LLE limit we have Λ + (0) < Λ − (0) and the curvature Λ − (0) can be positive only for Let us consider the case when Equation (20) has three real roots, which similarly to the standard LLE takes place for θ < − √ 3. Taking a real root u 2 0 = a 2 of Equation (20) we can express the remaining two roots as where three distinct roots are real for θ < − √ 3 and From this condition it follows that 2a 2 + θ > 0 for all − 4 the condition that u 2 0± < a 2 for both the roots defined by (A4) is also possible only for 2a 2 + θ > 0. Therefore it follows from (A3) that for any θ < − √ 3 the CW state with the highest intensity can have Λ − (0) > 0 only in the anomalous dispersion regime β < 0.
Since the condition Λ − (0) > 0 is not sufficient for the development of MI of a CW state, we assume further in (21) with k = α = 1 that λ = εw. Then up to order ε 2 we obtain: where D(w) = − u 2 0 + βw 2 − θ 3u 2 0 + βw 2 − θ . We can find the frequencies w at which the MI can observed using the conditions In particular, in the normal dispersion regime (β > 0) solving (A6) with respect to w 2 and θ we obtain and two solitions for θ Since according to (A7)-(A9) both w 2 and −θ increase with u 0 , the value of −θ corresponding to w = 0 gives the lower bound of the detuning parameter, for which the development of MI is possible. We note that for w = 0 and θ < 0 we have D(0) < 0, and Λ − (0) = 0, hence MI is not possible, and the lower bound obtained using w = 0 is not tight and provides a necessary condition for the MI. It follows from Equation (A7) that for w = 0 we get u 2 0 = √ 4β 2 +(σ 2 −2) 2 β > 2. Substituting this expression into (A8) and (A9) yields the following expression for the lower bound of the MI: where the sign "+" ("−") corresponds to Equation (A8) [Equation (A9)]. The expression in the right hand side of Equation (A10) achieves its minimal value −θ ≈ 2 Therefore, for −θ < 2 √ 3 MI is not possible in the normal dispersion regime. From the previous paragraph we see that in normal dispersion regime the MI is possible only for sufficiently large −θ > 2 √ 3. In the anomalous dispersion case (β < 0) by solving Equation (A6) we get the following expression for the MI frequencies −β > 2 and two solutions for the detuning parameter θ: Here, the coefficients by u 2 0 can change sign depending on the parameters, hence there is no clear lower boundary for −θ in contrast to normal dispersion regime. Moreover, we could numerically observe MI for σ = √ 2.5, see Figures 4 and 5 Finally, we conclude that while it is possible to observe MI in the case of anomalous dispersion, in normal dispersion case that can be done only for sufficiently large −θ > 2 √ 3.

Appendix B. Examples of MI Out of LLE Limit
Finally, we consider some examples of specific CWs, where relation (26) can be simplified, so that conditions for the change of curvature of Λ ± (ξ) at ξ = 0 can be obtained explicitly, and then MI can be demonstrated numerically for close parameters.
Let us consider a simple case where ϕ = −αA 2 0 + π/2 in Equation (12) and, therefore, cos αA 2 0 + ϕ = 0 in Equation (19). Hence, we get A 0 = η √ 1+κ < η, which corresponds to the CW with the lowest intensity in case of bistable S-shaped CW branch. Furthermore, from (23) we obtain D = α 2 η 4 (κ+1) 2 − 1 and Y + (0) = 1 √ κ αη 2 κ+1 + √ D . Therefore, from (25) and (A12) with sufficiently large β we can have This condition is reminiscent of the MI condition of a CW in the anomalous dispersion regime in the case of CW solutions of a semiconductor laser model [31], though secondorder dispersion coefficient β is multiplied here not by α but by a function of α, η, κ. Similarly to the case of a laser under optical injection [32], this condition manifests the change of curvature of Λ(ξ) at ξ = 0 for some β > 0, and it precedes appearance of a MI for larger β, which can be observed numerically (see Figure A1, top). The most usual way to find a localized structure in form of a bright dissipative soliton, is to look for three coexisting CWs, where CW with the largest field intensity A 2 0 is modulationally unstable. Indeed, for ϕ = −αA 2 0 + π 6 , where ζ = √ 3 2 , the CW takes the form and the condition Λ (0) ≥ 0 takes the form Similarly, the change of curvature of Λ(ξ) at ξ = 0 can occur only for β > 0, where also a MI can be observed (see Figure A1, bottom). We can use the CW (A13) by choosing corresponding parameter ϕ to find parameters where localized structures exist, and obtain parameters for the original system (1) and (2) to obtain these structures in numerical simulations (see Figure 12).