1. Introduction
Optical frequency combs generated from a continuous wave laser output in micro-cavity 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].
2. 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,
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,
and
are the injection rate and frequency offset,
and
are the central frequencies of the two filters, while
and
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.
2.1. 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]
where the complex function
can be expanded in power series near
and
represents the dispersion of the
kth order. ( For example, for a single Lorentzian absorption line we have
[
31], and second-order dispersion coefficient takes the form
for
, 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
and the reference frequency is chosen in such a way that the maximum of the function
is at
, so that
. The latter condition is equivalent to
which has two solutions
,
. One can see that both
can be obtained for any values of
and
due to
, and, moreover,
Below we will assume that the condition
is satisfied that corresponds to a situation when the largest of the two maximums of
is located at zero frequency, see
Figure 2. In the left panel of this figure we fix
and demonstrate that for
we have one global maximum, whereas for
we can get two local maxima in our combined filter of the same magnitude. In the second case, or for values of
close to
, 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
decreasing from 1 to
. For
we have
and one can see that
i.e., the difference between the two frequencies grows with
, and since for any of
we have
, respectively, for any fixed
we can choose
and find corresponding value of
(negative or positive) and
, 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
is at
for
as claimed earlier, since more narrow Lorentzian filter with the width
has the central frequency
closer to zero frequency. The other root (
6)
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
at zero frequency
.
Since
, in the power series expansion (
4) we obtain the second-order coefficient as
One can check that
can be represented in the following way
with
It follows from the conditions (
6), where
, and the relation
that the second term in the right hand side (RHS) of (
7) vanishes,
, and
. Hence, the parameters
and
represent the normalized second-order dispersion and the diffusion near zero frequency, respectively, whereas the parameter
can be considered as a scaling coefficient.
One can see that under these conditions
and since
grows with
, second-order dispersion grows with the difference between the frequencies. Moreover, using (
6) to express
one obtains
hence due to
we see that
, and
decreases with
increasing from 0 to
. Therefore, the lower bound on
can be estimated assuming
, where
and
2.2. Normalized DDE Model
By substituting Equation (2) into Equation (
1), using Equation (
6), and rescaling time as
, we obtain the following normalized dimensionless second-order DDE
where losses and forcing take the form
,
with
, and the phase shift is
. Here, similarly to the previous section one can see that
for
and
, hence for fixed
the parameter
decreases with increasing detuning between the two frequencies
. With the condition (
11) on dimensionless parameters
and
the equation is equivalent to the system with two filters (
1) and (2) for any
with global maximum of the combined filter transfer function
(
5) located at zero frequency
,
. We note that for
the global maximum of
is at another frequency and additional constraints on parameter
are necessary, whereas for
the diffusion coefficient in (
7) becomes negative, equivalence with (
1) and (2) is lost and zero solution
for
is unstable.
For
relation (
7) gives zero dispersion coefficient and diffusion coefficient equal to
, which agrees with previous analysis of a similar DDE [
39] in case of
.
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
in and spurious instability [
31]. In contrast, non-delayed second order derivative
in (
12) appears without any expansions, and the term
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,
.
2.3. Limit of Lugiato-Lefever Equation
Assuming large delay limit
, where
and
, taking
, and rescaling the time variable
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
with
, 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,
and the relation
with the periodicity condition (
15). Finally collecting the third order terms in
, using the relation
, and applying solvability condition [
40,
41] we get the generalized Lugiato-Lefever equation (LLE)
with the diffusion coefficient
and the boundary condition (
15). It is well known that for
in Equation (
16) corresponds to anomalous dispersion regime, and for
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).
3. Continuous Wave (CW) State
The CW state of the Equation (
12) with
takes the form
, where the real quantities
and
satisfy the system of the transcendental equations
which leads to a single transcendental equation for
Assuming
,
,
,
∼
, and
, this equation can be approximated by a cubic equation for
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.
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
. For that, we linearize the Equation (
12) near the CW state
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
. In the limit of large delay time
the eigenvalues belonging to the pseudo-continuous spectrum can be represented in the form
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
and, in particular,
, where
Moreover, the second derivative of
at
takes the form
One can see from (
25) that for small
and
we have
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
(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
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
for all the wavenumbers
. In particular, at zero wavenumber
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
corresponding to the greater of the two values
. Note, that both the quantities
are always negative at
due to the inequality (
11). It follows from Equation (26) that such a sign change can take place only in the case when
and, hence,
must be real. Otherwise,
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
and for the CW state with the highest intensity we can have
only in the anomalous dispersion regime
. Out of LLE limit we can have positive
in case of
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
precedes the appearance of MI of the CW for stronger dispersion if
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.
4. 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
,
, which corresponds to two filters (
5) with
,
,
. Let
,
and
,
, where
are chosen according to analysis of DDE (
13) in the limit
. These values directly correspond to the parameters of LLE (
16) where localized structures are known to exist (
). We also pick a reference injection strength
where three CW states are guaranteed to coexist for
, 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
.
Figure 6 was obtained with
, which was the smallest value of
used in our simulations and corresponded to
. Similarly to the LLE [
44,
45] the bifurcation diagram in the left panel of this figure shows a typical S-shaped CW branch with stable lower part and modulationally unstable upper part. The branch of unstable periodic solutions bifurcates from the unstable middle part of the CW branch and it becomes stable after a fold bifurcation at
. The stable periodic solution has only a slight asymmetry in its time profile and resembles the temporal dissipative solitons of the LLE, see right panel of
Figure 6.
Thus, close to the limit
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
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
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
and
from Equations (
1) and (2) and the condition
one obtains
in (
12), which is satisfied for
.
For larger
, which corresponds to
, 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.
Figure 8 and
Figure 9 show the branches of CW and periodic solutions with scaled intensity (In these figures
is shifted, more precisely,
for the left panel of
Figure 8,
for the right panel,
for the left panel of
Figure 9,
for the right panel, and
,
,
,
,
) obtained by increasing gradually
up to 1, which corresponds to
. The width of the S-shaped area of the CW branches as well as the width of the branch of the localized solutions increases with
up to
and then decreases. For
the upper branch of CW states stabilizes for the chosen parameter values (
Figure 9, left). Asymmetric localized structures corresponding to the stable parts of the periodic solution branches with
in
Figure 10 look similar to the those shown in the right panel of
Figure 7 obtained for
.
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 (
) the localized solutions can be observed only for negative
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
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
(top right panel of
Figure 11) and then to
(bottom left panel), one can see that around
localized structures can be observed for
as well. Furthermore, for
(bottom right panel) stable bright temporal dissipative solitons can be observed even with small positive
,
. 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.
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).
5. 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 (
) 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.