1. Introduction
Optical frequency combs (OFCs) can be generated using mode-locked lasers, electro-optical modulation of continuous-wave light, or integrated microresonators [
1]. When these combs are produced in an integrated microresonator, they are called Kerr combs. Kerr combs have recently attracted a large amount of attention in the photonics research community; they have applications in areas such as communications, metrology, optical clocks, photonic radars, artificial intelligence, and quantum information [
2,
3,
4,
5,
6,
7]. Ensuring their reliability, robustness, and low cost is a key challenge.
Dissipative Kerr solitons (DKSs) are a special type of Kerr comb. DKS yield a series of discrete, equally spaced frequency lines in the output spectrum with a constant frequency offset. DKS generation can be performed using a single microring resonator (MRR) [
1,
2,
3,
8,
9,
10], defining paths in which the laser pump power and detuning are controlled [
8,
11]. To understand the dynamics of this structure, the different operating regions of the MRR in the parameter space
have been identified [
8,
12,
13]. These regions correspond to cnoidal waves, including Turing rolls, perfect soliton crystals, and single solitons or single DKSs, as well as chaotic solitons, soliton molecules, and imperfect soliton crystals. Specifically, Jose Jaramillo et al. defined a deterministic path for single DKS generation in a Si
3N
4 MRR, based on the above-mentioned regional classification [
8]. However, it is difficult to apply this technique in experiments because it requires the pump laser power and wavelength to be changed rapidly along a defined path.
Other approaches include thermal control, pulse-triggers, multiple laser pumps, and engineered spatial mode interactions, among others [
1,
2,
14,
15]. Although these approaches can generate solitons deterministrically, they involve difficult and expensive experimental setups. In [
16,
17], it was demonstrated that mode-crossings due to spatial mode-interactions can allow deterministic single soliton generation in a single MRR. However, its approach is conditioned by the structure’s geometry, limiting achievable dispersion profile and constraining generation to the paths proposed by the study, which may not be optimal. One of the limitations of a single MRR is the inability to control mode interactions.
Recently, the use of dual-coupled microring resonators (DCM) has been proposed to facilitate single DKS generation [
1,
18,
19,
20,
21]. This structure allows the formation of specific dispersion profiles that depend on the coupling of the auxiliary ring, thus facilitating the manufacturing and engineering of dispersion. DCM can be controlled to modify mode interactions and enlarge the region of the parameter space within which single solitons are stable, making them easier to obtain deterministrically. Approaches such as [
18,
19,
20,
22,
23,
24,
25] demonstrate mathematical models and the dynamics governing light propagation in this structure. While these approaches demonstrate the formation of a DKS, those structures do not control mode interaction. As a result, their impact is fixed and confined to the specific case studies.
Other work studies the stability and accessibility of DKSs in dual-coupled microresonators (DCMs) as the system parameters vary [
26,
27,
28,
29]. These approaches have led to mathematical frameworks that facilitate the study of frequency combs, in particular enabling the determination of parameter regimes (within which stationary waveforms such as Turing rolls and DKSs exist) and their stability. However, a drawback in these studies is that the models are exclusively designed for microresonators with normal dispersion, overlooking structures with anomalous dispersion. Additionally, the structures analyzed do not control the mode interaction, restricting the generalization of the obtained results. In [
22], the authors propose a DCM to generate Turing rolls and single DKS. This structure considers anomalous dispersion and establishes a configuration of mode interactions, simplifying the acquisition of DKS based on the characterization of the parameter space
. Additionally, they perform a stability analysis to identify regions in the parameter space in which DKSs are stable. However, the effect of mode interaction on DCM waveforms has not been thoroughly examined, nor their role in potentially limiting the accessibility of single DKS regions in the parameter space
.
In this study, we numerically demonstrate the generation of a single DKS using a straightforward constant pump power with an adiabatic frequency sweep in a dual-ring Si3N4 structure, which is sufficiently slow that the state of the optical waveform is not affected by the rate of the sweep. In the DCM, the mode interactions are electrically controlled by a microheater in the auxiliary microring. We chose strength and frequency values that facilitated access to the soliton and avoided the chaotic region, based on these characterizations. Finally, we verified the deterministic generation of the DKS. This document has the following structure: First, we characterize within the parameter space by fixing values of a and b. These values are experimentally chosen to generate DKSs in the lower region of the space, thereby circumventing the chaotic region. Next, we fixed values for four points of the space and computed the probability of DKS generation for each point. These points are manually selected, considering the position facilitating the DKS generation with a constant and a sweep from blue to red. Following this, we select four sets of values for a and b where the DKSs generation probability is high and compute the parameter space for these configurations. This process allows us to visualize potential regions where DKSs may generate and to ascertain which configurations exhibit a more pronounced presence of these waves. Finally, we analyze the stability of each zone to determine if reliable generation of DKSs is feasible with these values.
2. Characterization of Pump Power and Detuning Parameter Space
To understand the influence of mode interactions in the DCM, we explored the operating regions in the parameter space
. This study is necessary to identify the region where DKSs are generated. The characterization of the operating regions depends on applying the Lugiato–Lefever equation (LLE) [
8,
15]. The LLE serves as a mathematical framework to simulate and characterize the propagation of light within a microring resonator, modeled through the parameters of loss, coupling factor, and chromatic dispersion. This framework is decisive when responding to input pump and damping parameter variations, providing a comprehensive tool for our analysis.
The DCM was modeled by incorporating an avoided mode crossing (AMX) in the LLE [
17]. This AMX accounts for the optical effects caused by an auxiliary ring with a smaller radius within the main ring, as illustrated in
Figure 1. The AMX on the physical device can be controlled using a microheater positioned on top of the auxiliary ring. The observed two peaks in the DCM (see
Figure 1e) relative to the MRR result from mode interactions caused by AMX, where the resonance modes of the primary microring and the auxiliary ring couple strongly, manifesting as peaks in the spectral response. These peaks represent resonant states with altered mode density and enhanced light confinement, enabling specific dispersion profiles that support DKS generation. The LLE equation can be expressed as follows:
where
is the complex envelope of the total intracavity field,
t is the time variable,
is the fast time variable,
is the round-trip time,
is half of the total loss per round trip,
is the frequency detuning,
L is the cavity length,
is the second-order dispersion coefficient,
is the Kerr coefficient,
is the coupling coefficient between the waveguide and the microring resonator,
is the pump field,
a and
b are the strength and frequency of the mode interaction, respectively, and
is the Fourier transform, defined as
The LLE equation incorporates AMX effects with the term
The parameter values for
a and
b were selected based on previous studies that validate this approach [
16,
30,
31]. These studies highlight the role of AMX in facilitating DKS generation by lowering the threshold for soliton formation under certain dispersion conditions [
30,
31]. This incorporation of AMX into the LLE framework allows us to simulate and predict the influence of pump power and dissipation variations on modal dynamics and DKS formation. This approach is essential for understanding how device parameters impact DKS generation, providing insights for optimizing design in specific applications.
Equation (
1) can be written in terms of three normalized parameters. One standard way to normalize Equation (
1) is to let
in which case Equation (
1) becomes
where
To solve this equation numerically, we used the split-step Fourier method. Note that the normalized parameter definitions proposed in Equation (
8) facilitate our analyses in a more general framework, and allow for comparison with other methods presented in the literature [
8,
15]. To perform the simulation of the LLE model, we used the same parameters as in [
8] which correspond to a Si
3N
4 structure. The parameters used were a main radius of 100
m,
GHz,
,
,
,
, and
. The parameters chosen in the LLE model are founded on experimentally validated values from real physical devices [
8,
32,
33]. Additionally, we replicated real conditions in the simulation by initializing the intracavity field with Gaussian noise, set at a standard deviation
, which is added at each iteration. This approach ensures that our simulations are closely aligned with practical implementations, accurately reflecting the behavior of physical photonic systems. Due to the complexity and computational cost, thermal shift was not included in the model.
Additionally, we set the pump power for the single MRR at a detuning (
) and a pump power of
(57.1 mW), ensuring a well-defined behavior in the cnoidal wave region. Similarly, for the DCM, we tuned the pump power at a detuning (
) and
(47.6 mW). These values were chosen because this is the starting point for obtaining stable solitons [
8]. Then, a single step jump was made to a point in the parameter space
and maintained for 1.5
s
. We chose a configuration
that shifted the single soliton region to be accessible from the blue side of the comb generation region in the parameter space
. These values are
[GHz] and
, in which we observed the shift of DKS region taking into account that the should potentially be replicated in experiments. To emphasize, generating a mode interaction strength greater than
[GHz] could be physically unfeasible in chip design because achieving this interaction with the geometric designs proposed in the literature poses a challenge [
1,
34]. Hence, exploring values below this threshold is warranted to align with experimental limitations.
Figure 1c,f show the simulation results for a single MRR and DCM, respectively. In these simulations, we show the number of peaks in the time domain for each endpoint
after 3
s. In the first simulation with a single ring (see
Figure 1c), we emphasize the presence of the thin black DKS region on the right side of the comb generation zone. The location of this region in the parameter space poses challenges for accessing DKS. The difficulty occurs because the path must avoid the chaotic zone, requiring the proposal of a non-linear route. Next, we simulated the DCM by introducing the mode interaction with the parameters
we previously proposed. The shift of the DKS area to the left and bottom of the region where combs are generated is apparent. This new position of the single soliton region makes it accessible in a path at constant pump power and with detuning that sweeps directly from blue to red, as shown by the black arrow in
Figure 1f.
We next study the impact of parameters on the generation of single DKS. Instead of fixed values, we aim to identify effective combinations of a and b, increasing the probability of single DKS generation. This approach deepens system understanding and helps pattern recognition, aiding future research. To define , we analyzed the single DKS generation probability at four key points in the parameter space , producing an explicit function of a and b. The choice of these specific points was deliberate, guided by their strategic location below the chaotic region, avoiding crossing this region and making it accessible in a path with detuning that sweeps directly from blue to red. The probability calculation at each point was determined by systematically varying the values across the specified ranges, where a ranged from 0 to [GHz], and b ranged from 0 to 20. The number of iterations used in the probability calculation was 30 simulations, introducing variability in the initial field and in the integration path through Gaussian noise with .
Figure 2 presents the probability of single DKS generation as a function of different parameter values
, exploring these variations at the four specified points. Notably, various configurations of
exist, enabling reliable single DKS generation. However, many of these configurations are unattainable in a physics device, so we prefer those with a value of
a less than
[GHz]. To maximize the generation probability, we have deliberately chosen specific values corresponding to the regions of the graph marked with red, magenta, green, and black rectangles. These values, proven to be potential in terms of DKS generation, are as follows:
,
,
,
. This selection is based on the maximum generation probability associated with the chosen points, aiming to ensure reliable device performance under consideration.
Figure 3 displays the parameter space
for the previously proposed parameters
. It is important to note that all four figures show a single DKS area in the lower right part of the parameter space, highlighted in black. This signifies direct access to the soliton region without passing through the chaotic zone. In
Figure 3d, the
parameter space exhibits a larger DKS region than with other choices of
a and
b (see
Figure 3). However, the crossing strength is relatively large with
, which may pose challenges in a physical experiment. Conversely,
Figure 3a,c reveal a considerable area and offer potential values for a potential physical experiment. Although
Figure 3 demonstrates the presence of single DKS in these spaces, verifying if the generation in these areas is deterministic is required. There is a possibility that a single DKS may not always be obtained, potentially compromising the reliability of the path. To verify path accuracy, it is necessary to conduct multiple simulations, as proposed in [
17], and confirm the repeatability of the path under different noise conditions.
3. Deterministic Single Soliton Generation
To demonstrate the deterministic access of single DKS in the four previously proposed configurations of
a and
b, we conducted 200 repetitions of the LLE simulation, employing distinct initial random fields for each setting. In the case of configuration
, we maintained a constant normalized power
of 1.04. The normalized detuning
was swept from
to
over a simulation period of 3
s, as depicted in
Figure 4a. For parameters
, the magnitude
was set to 1.1538, and the detuning
ranged from
to
. Similarly, for
parameters,
equaled 1.08, with
spanning from
to
. Finally, in the case of
parameters,
was fixed at 1.09, and
varied from
to
. The values of
and
chosen for DKS generation along each path are determined by the trajectories indicated by the magenta arrows in
Figure 3.
We emphasize that these trajectories do not cross the chaotic region.
Figure 4a–d show that we obtain a single DKS state; see the temporal shape (center) and spectrum (right). It is noteworthy that when the value of
a is
, the contour of the spectrum deviates from the a hyperbolic-secant shape. However, increasing
a causes this shape to align more closely with the sech(x) shape. Note that configuration
does not conform to the sech(x) shape, unlike the other configurations.
Figure 4a–d show that
remains stable and noise-free along the selected paths, supporting controlled DKS generation. The absence of abrupt fluctuations confirms that the trajectory avoids chaotic regions, essential for stability. Additionally, the convergence of
to a constant value indicates that the DKS reaches a stable configuration, which is crucial for maintaining coherence in different applications.
Figure 5 shows the number of peaks (
) for the 200 simulations in the four proposed
a and
b configurations. We obtained in a single peak in all the realizations, indicating that a single DKS was always generated, demonstrating that the paths are deterministic. Although these configurations enable the generation of a single DKS, these simulations do not demonstrate that the DKS is stable. For this reason, we perform a stability analysis.
To validate DKS generation under variations in the DCM caused by changes in aspect ratio and manufacturing tolerances, which affect the radius, loss, coupling coefficient, and dispersion of the device, we conducted a Monte Carlo simulation. In this simulation, we varied the radius, loss, coupling coefficient, and dispersion randomly within the ranges specified in
Table 1. This approach allowed us to evaluate the robustness of DKS generation under realistic manufacturing fluctuations, ensuring reliable device performance within practical tolerance levels.
As shown in
Figure 6, we conducted 200 Monte Carlo iterations across the four paths defined in
Figure 4, achieving DKS generation in 175 cases (
Figure 6a), 130 cases (
Figure 6b), 160 cases (
Figure 6c), and 181 cases (
Figure 6d). These results indicate a DKS generation probability exceeding 80% for
Figure 6a,c,d and 65% for
Figure 6b. While variations in the aspect ratio influenced the results, the method still showed robustness and adaptability to specific structural changes. Although DKS generation is not entirely deterministic, our approach can reliably be applied post-fabrication using precise parameters obtained after manufacturing, once the aspect ratio and device adjustments are confirmed. This characterization provides key parameters such as quality factor (Q), FSR,
, and
, enabling precise DKS generation analysis and parameter mapping that ensures robustness and adaptability in fabricated devices.
4. Single DKS Stability Analysis
The selection of and values is essential for achieving a stable DKS. Conducting a stability analysis is crucial, as it identifies the specific ranges of and within which the DKS state can be maintained, thereby avoiding unintended frequency comb patterns or the absence of comb generation.
To analyze the stability region of the single DKS in the DCM, we propose a new method based on previous work [
8,
15]. Our proposed method first solves Equation (
10) while allowing
S or
to vary. To solve Equation (
10), we employ the split-step Fourier evolution method, which yields a highly stable and stationary solution. To obtain this stable solution, we use the paths proposed above, which sweeps the detuning from
to
while keeping
S constant. Subsequently, we jump to a final point in the range of
with a resolution of 0.025 in both parameters. At each point, we wait for 3
s to evaluate the stability of DKS.
In parallel with solving Equation (
10) to find the wave solutions, we determine their stability. We write
where
and
are perturbation of the solution
and its complex conjugate
. The perturbations
and
obey the linearized equation
where
and
The terms
and
correspond to the Fourier transform and its inverse function. Those functions enable mathematical modeling, as these transforms have a matrix representation, and the term
can be assigned to a diagonal matrix, resulting in a series of matrix operations of dimension
that operate on the solution vector
. Next, We convert Equation (
12) into an eigenvalue equation
, where
I is the identity operator, and we discretize this equation. We then proceed to computationally find the set of all eigenvalues
of this equation. If any eigenvalues are positive in its real component, the solution
is unstable for
and
S values.
is a drift term introduced in the model because the avoided crossing introduces a slight drift in the solution, causing it to become non-stationary. To find
, we calculate the soliton centroid in two consecutive round trips, determining the
difference between both to obtain this value. We repeat this process 30 times and calculate the average value. To evaluate the stability of single DKS at each point, initially we determine the number of peaks by computing the optical intensity [
8]. If the value of the number of peaks is one, it indicates that the wave may be stable. Otherwise, the wave will be unstable and will converge to a different type of wave. When the number of peaks is one, it indicates that the wave is a single DKS. Its stability is evaluated by computing the eigenvalue solution of the eigenvalue of the
matrix and validating that
. This process is necessary because a single soliton wave may not be inherently stable and can exhibit oscillations or breathing, which are considered unstable.
The red area in
Figure 7 shows the regions where the single DKS is still stable at the end of the simulation for the four parameter sets
that we investigated proposed. This result shows that the deterministic generation of a single DKS is correlated with a greatly enhanced stable region. We also observe that the parameter space
exhibits a small stable region, suggesting that while a substantial number of single DKSs are produced, they are unstable, rendering them unreliable for deterministic generation (see
Figure 7a). Furthermore, the other parameter spaces have a potential area for soliton generation (see
Figure 7b–d). All the stable zones lie below the chaotic zone, which facilitates the generation of single DKS, as proposed in this research. Note that configuration
presents a value greater than
. This situation poses a challenge in the design and manufacturing process of DCM. Therefore, it is more feasible to utilize the configurations proposed in
and
.