Dissipative Kerr solitons in a photonic dimer on both sides of exceptional point

Exceptional points are a ubiquitous concept widely present in driven-dissipative coupled systems described by a non-Hermitian Hamiltonian. It is characterized by the degeneracy of the Hamiltonian's eigenvalues and coalescence of corresponding eigenvectors. Recent developments demonstrated that exceptional points can play an important role in photonics. However, to date, exceptional points have been extensively examined in the systems supporting only a few optical modes, thereby leaving the observation of collective (multimode) effects outside of the scope of study. In the present paper, we analyze the role of exceptional points in nonlinear multimode photonics. Specifically, we provide insights into complex nonlinear dynamics arising in a continuous wave-driven pair of strongly coupled nonlinear micro-resonators (i.e. a nonlinear photonic dimer) operating in the multimode regime. Investigating this system, which is known to possess exceptional points, we find two fundamentally different nonlinear regimes of operation corresponding to effective parity-time symmetric and broken parity-time symmetry states. We demonstrate that the photonic dimer can be critically coupled to a bus waveguide, thereby, providing an efficient generation of the dissipative Kerr solitons on both sides of the exceptional point. The parity-time symmetric case, which corresponds to a pair of symmetrically split resonances, has been recently shown to exhibit a variety of emergent phenomena including gear soliton generation, symmetry breaking, and soliton hopping. Dissipative solitons generation in the parity-time symmetry broken case - leading to the dissipation splitting - up to now remains unexplored.


INTRODUCTION
Symmetry is one of the main fundamental concepts in physics which underpins conservation laws, micro-and macro-properties of matter, rising of degeneracies, etc. Breaking of symmetry in spontaneous or deterministic way is responsible for a variety of phenomena with a particularly illustrative example of the Higgs mechanism explaining the appearance of masses of certain bosons [1]. Systems obeying Parity-Time (PT ) symmetry have been shown to provide a possibility to treat open quantum systems described by a non-Hermitian Hamiltonian and nonetheless retrieve a real spectrum of eigenvalues [2,3]. Eigenvalues of PT symmetric Hamiltonian have two typical regions on the parameters space corresponding to preserved and broken symmetry. The transient point where the eigenvalues and eigenvectors coalesce is called an exceptional point (EP) [4,5].
The effect has been observed in various physical systems [3]. Due to the well-controlled laboratory conditions and a wide range of possible applications, guided optics serves as one of the primary platforms for investigating effects that emerge in PT -symmetric systems [5][6][7][8]. These effects cover observation of strong non-reciprocity in banded waveguides, enhanced lasing and (classical noise limited) sensing in coupled ring resonators with loss and gain (for more examples see review article [5]). Remarkably, after a gauge transformation, PT symmetry can be unraveled in completely passive resonators having different loss rates [5,9]. Here, it is important to mention that the majority of these studies have been considering linear effects in single-mode arrangements or nonlinear ones in non-resonant systems [10].
Almost in parallel, and without apparent intersection, another field of study has been rapidly developing in photonics: design and fabrication of passive nonlinear coherent broadband light sources [11]. For this purpose, exactly the same optical platforms have been used: optical waveguides and resonators but in the opposite, strongly nonlinear regime. Nonlinear waveguides served as a sources of supercontinuum signals [12], while microand macro-resonators has been used for generation of stable and coherent frequency combs [13]. The latter has been achieved due to the observation of localized coherent structures in passive optical micro-resonators [14]. These structures generated in media with χ (3) nonlinearity are called dissipative Kerr solitons (DKSs). They have been widely investigated in nonlinear photonics over the last decades [13,[15][16][17]. The stability of DKSs relies on the balance between chromatic dispersion, Kerr-type nonlinearity, parametric gain and the intrinsic cavity losses [18]. DKSs are exact solutions of the damped-driven nonlinear Schrödinger equation known as Lugiato-Lefever equation (LLE) [19,20]. Their observation in passive macroscopic (fiber) cavities [17] and (integrated) microcavity arXiv:2101.09237v1 [physics.optics] 22 Jan 2021 systems [14,21] has spurred a vivid research effort unravelling a rich inherent dynamical behaviour [22][23][24][25][26][27][28]. The discovery of DKSs in microresonators revolutionized the field, bringing coherent frequency combs to outside-oflaboratory applications [13]. Remarkably, a pioneering theoretical proposal considering the DKS generation in the PT -symmetric system has been presented in [29].
Recently, DKSs have been discovered in a high-Q multimode photonic dimer (pair of strongly-coupled, almost identical nonlinear resonators) [30]. The photonic dimer has revealed a pleiad of emergent phenomena including soliton hopping, periodic appearance of commensurate and incommensurate dispersive waves (DWs), and symmetry breaking related to the discreteness of the system. Solitons have been generated in both resonators simultaneously and due to the underlying field symmetry were called gear solitons (GSs).
In this work, we investigate all-passive photonic dimer with a hidden PT -symmetry (further referred to as PTsymmetry for simplicity) in linear and nonlinear multimode regimes. In the linear regime, we analyse conditions for the critical coupling and demonstrate that the line of EPs is a demarcation of these conditions. The exceptional point line splits the parameter space into two parts which we refer to as split resonance (PT symmetric) and split dissipation (PT -symmetry broken). In the split resonance regime, we show novel insights into the effects previously reported in [30], using the supermode basis representation. Passing through an EP, we observe the divergence of the nonlinear interaction efficiency which hints at the enhanced sensitivity. Further, we investigate the nonlinear dynamics in the split dissipation regime which includes single-resonator DKS, dark-bright DKSs pairs, and highly-efficient perfect soliton crystals. Finally, we demonstrate switching of the soliton-generating cavity caused by the nonlinear alteration of the PT symmetry.

II. EXCEPTIONAL POINT AS A DEMARCATION OF THE CRITICAL COUPLING CONDITIONS
We consider the system of two multimode resonators [ Fig. 1(a)], with identical intrinsic loss rate κ 0 , mode spacing D 1 , and geometry, such that the dispersion and Kerr nonlinearity coefficient are also identical. A global offset between the resonant frequencies ω µ of their respective modes µ is introduced with the inter-resonator detuning δ. The two resonators are coupled to each other by the evanescent field with the rate J µ , which generally depends on the mode number. Each resonator is coupled to a waveguide (through and drop ports) with the rates κ ex,i , i = 1, 2. Resonator 1 is pumped by a continuous wave (cw) laser at frequency ω p . Nonlinear dynamics in the photonic dimer can be described by two coupled Lugiato-Lefever equations (LLEs), which in Fourier space is expressed as follows [30,31]: where g K = ω 2 0 cn2 n 2 0 V eff is the Kerr coefficient, c stands for the speed of light in vacuum, -the Planck constant, ω 0 -the frequency of the pumped mode, V eff -the effective mode volume, n 0 and n 2 are linear and nonlinear refractive indexes, respectively, δ µ,0 is the Kronecker delta, s in = Pin ω0 -the input pump field amplitude, A µ , B µ are the field amplitudes of the modes with index µ in the first and second resonator, respectively. The variables A, B are the slowly varying intra-resonator field envelops, and F[...] µ denotes the µ th -component of the discrete Fourier transform, which are defined in Appendix A.
Modes of each resonator with identical angular momentum are linearly coupled with the inter-resonator coupling rate J µ . In contrast, the Kerr nonlinearity couples all the modes within each resonator via four-wave mixing (FWM) processes. The interplay between the linear coupling in the spatial dimension and the nonlinear coupling in the frequency dimension is the source of the rich dynamics of the system.
In the present section, we restrict ourselves to a linear and single mode analysis by considering only the central mode µ = 0 with g K = 0. The critical coupling conditions are of particular interest for maximizing the pump transfer to the resonators and the eigenvalue analysis of the coupled mode matrix for finding the EP conditions, which separate two conceptually different dimer states.

A. Critical coupling conditions
Let us analyse the condition for critical coupling. In the linear single-mode representation, Eq. 1 is simplified as follows: In Eq. 2, we defined the external coupling mismatch ∆κ ex = κ ex,1 − κ ex,2 and the average external coupling κ ex = 1 2 (κ ex,1 + κ ex,2 ). The identity matrix is denoted as I. Critical coupling is achieved when the transmission via the through port [ Fig. 1(a)] vanishes, i.e. s out,1 = 0. In the case of a single resonator, critical coupling is achieved when the external coupling rate matches the loss, i.e. κ ex = κ 0 [32,33]. For two resonators, the conditions are easily found in case of δ = 0. There are two possibilities Eq. 3 is a natural generalization of the critical coupling conditions for a single resonator that can be achieved by setting J to zero. Eq. 4 satisfies the critical coupling condition at ω p = ω 0 ± 4J 2 − (κ 0 + κ ex,2 ) 2 . This condition requires strong coupling, i.e. J > 1 2 (κ 0 + κ ex,2 ). The critical coupling conditions are shown in Fig. 1(d) for κ ex,2 = 0. Typical cavity field intensities for both cases are plotted in Fig. 1(b,c) as a function laser detuning.
The first critical coupling condition given by Eq. 3 has a quadratic dependence on the inter-resonator coupling rate [ Fig. 1(d)]. It leads to a broad resonance with a dip in the first resonator (blue) and a narrow resonance in the second resonator (red) at the same resonance frequency [ Fig. 1(c)]. The second critical coupling condition (Eq. 4) branches off the first one at J/κ 0 = 0.5 and does not depend on the inter-resonator coupling rate. It features split resonances with identical linewidths [ Fig. 1(b)].
Experimental implementation of the multimode photonic dimer demonstrated the presence of the nonvanishing inter-resonator detuning δ caused by the fabrication imperfectness. Nonetheless, the possibility to control and manipulate δ, and thereby establish control over the output solitonic spectrum, has been demonstrated and efficiently implemented by imprinting heaters directly on the photonic device [30,34]. Critical coupling at non-vanishing δ is possible as well. The inter-resonator detuning introduces asymmetry in the distribution of the supermodes (eigenvectors) in each resonator. Thus, the supermode confined in the first (second) resonator requires smaller (larger) κ ex,1 to be critically coupled. It follows that in general when δ = 0 only one supermode can be critically coupled for a given value of κ ex,1 . For more details, see Appendix B.
The qualitative behavior of the photonic dimer can be anticipated by examining the eigenvalues of the system Eq. 2. Operating with a naturally Non-Hermitian system, we can exploit the concept of EP [4] to shed light on the nature of each critical coupling conditions.

B. Eigenvalues and exceptional points
The eigenvalues of the matrix M defined in Eq. 2 in case of δ = 0, κ ex,2 = 0, and ω 0 = ω p are given by where the real (imaginary) part corresponds to resonance frequency (loss rate). The eigenvalues are shown in Fig. 1(e) as a function of κ ex,1 for an inter-resonator coupling J = 1.5κ 0 . Two different regions of split resonance and split dissipation are identified and shaded in Fig. 1 with orange and purple, respectively. For κ ex,1 < 6κ 0 , the eigenvalues have degenerate imaginary part and split real parts, associated with the split resonances as depicted in Fig. 1(b). In contrast, κ ex,1 > 6κ 0 leads to degenerate real parts and split imaginary parts, i.e. to identical resonance frequencies but different loss rates, as can be seen in Fig. 1(c). The two regions correspond to the PT -symmetric and PT -symmetry broken states, respectively. An EP is found between the two regions at κ ex,1 = 6κ 0 , where the system eigenvalues become degenerate and the two eigenvectors coalesce because of the vanishing square root in Eq. 5. EPs lie along the line defined by κ ex,1 = 4J, which separates the two critical coupling conditions in the (J, κ ex,1 ) plane. Remarkably, the two critical coupling conditions and the EP line fork at J = 1 2 κ 0 , κ ex,1 = 2κ 0 . This particular point is highlighted by the blue star in Fig. 1(d). It is the only EP that satisfies a critical coupling condition. This point also marks the entry into the strong coupling regime (J > 1 2 κ 0 ). Above reasoning is valid when κ ex,2 = 0. In the general case (κ ex,2 = 0) the line of EPs can cross the line corresponding to critical coupling conditions.
Concluding, there are two types of critical coupling conditions in the photonic dimer. These conditions are found on both sides of the EPs, such that critical coupling can be achieved in the PT -symmetric as well as PT -symmetry broken states. In the next sections, we examine the versatile nonlinear dynamics and dissipative Kerr soliton generation in theses states.

III. CRITICALLY COUPLED RESONATORS: SPLIT RESONANCE (PT -SYMMETRIC)
In the present section we discuss the case of the split resonance (PT -symmetric). We revisit ideas presented earlier in [30] by looking at the nonlinear dynamics from the supermode perspective. We demonstrate a separability of the GS dynamics from DWs living in the S supermodes. Finally, we show how this representation explains the origin of soliton hopping effect. An essential part of the investigation of the dynamics inherent to the photonic dimer and described by Eq. 1 relies on numerical simulations.

A. Four-wave mixing pathways between supermodes
The linear part of Eq. 1 can be diagonalized by a linear transformation on each pair of modes with index µ. We define the complex inter-resonator detuning δ c = δ − i 1 2 ∆κ ex . If the inter-resonator coupling is independent of the wavelength, the complex frequency splitting ∆ω c = 4J 2 + δ 2 c , and the (complex) normalized inter-resonator detuning d c ≡ δ c /∆ω c are independent of the mode index. Therefore, the nonunitary transformation diagonalizing the linear part of Eq. 1 is given by A s,µ = αA µ +βB µ , A as,µ = βA µ −αB µ , where Here S and AS stand for the symmetric (S) and antisymmetric (AS) mode, as they are completely symmetric (antisymmetric) at d c = 0 in the split resonance regime. Then, by defining the spatial envelope of the field in the S and AS modes A s (θ) = µ A s,µ e iµθ , A as (θ) = µ A as,µ e iµθ , we can express Eq. 1 in the supermode basis (see details in Appendix C and alternative Hamiltonian formulation in [30]): As one can see, the linear anti-diagonal terms are eliminated in the supermode basis, while the nonlinear terms (diagonal in the resonator basis) induce nonlinear coupling between the supermodes. In particular, we identify FWM processes between the supermodes. The efficiencies of FWM pathways are associated with nonlinear coupling coefficients t i , i = 1, .., 4, defined in Eq. C14 of Appendix C. For example, the term 1 2 t 3 A 2 as A * s in Eq. 7a signifies the annihilation of two photons in the AS supermode and the creation of two photons in the S supermode. The rate of the process is proportional to g K t 3 . Each FWM process can be represented by a quantum Hamiltonian term, for examplê whereâ † i,µ (â i,µ ) is the annihilation (creation) operator for supermode i = s, as with longitudinal mode index µ, and µ 1 + µ 2 = µ 3 + µ 4 . The 9 different non-linear processes from Eq. 7 are depicted in Fig. 2

Intra-band Even
Inter-band Even Inter-band Odd are arranged in categories corresponding to intra-band even processes, inter-band even processes and inter-band odd processes.
We refer to a nonlinear process as intra-band when two annihilated and two created photons are from the same supermode family, while inter-band processes imply nonlinear mixing of photons belonging to different supermodes, inspired by the concept of Bloch bands in condensed matter Physics. The number parity of the process (even or odd) refers to the number of photons from each supermode family that is involved. We note that processes (2,4,7,9) are the counterparts of processes While the index µ has been omitted in Fig. 2(a) for readability, both the mode number and the energy have to be conserved in a FWM process. We employ the concept of integrated dispersion D int (µ) = ω µ − (ω 0 + D 1 µ) to depict the processes which satisfy the phase matching conditions.
The real and imaginary parts of the nonlinear coupling coefficients are shown in Fig. 2(b) as a function of the normalized inter-resonator detuning d = δ/ √ 4J 2 + δ 2 with solid and dashed dotted lines, respectively. The parameters are chosen in the split resonance regime with J = 20κ 0 . Vanishing and non-vanishing ∆κ ex are considered, emphasizing that the imaginary parts of all the nonlinear coupling coefficients originate from the external coupling mismatch. The imaginary part, however, generally constitutes only a small fraction of the absolute value of the nonlinear coupling coefficients.
Coefficient t 1 is responsible for the intra-band processes, that is, the usual FWM within the same mode family (S) or (AS). It has its lowest value equal to 0.5 at the maximum hybridization (d = 0). In contrast, coefficient t 3 is maximized at d = 0 and causes inter-band and even processes. The coefficients t 2 and t 4 are responsible for inter-band and odd processes. Their real parts are odd with respect to d. Therefore, there are no odd FWM processes at d = 0, unless an external coupling mismatch is present. In this case, coefficient t 4 has a non-vanishing absolute value.
B. Nonlinear dynamics and soliton generation in split resonance (PT -symmetric) regime The type of the critical coupling conditions corresponding to the split resonances allows for accessing dynamical states characterized by efficient generation of bright dissipative Kerr solitons in both cavities. Emergent dynamical effects described in [30] are found in this regime. In this section, we provide an additional (to the result already shown in [30]) and complementary description of these phenomena by representing the interresonator field in the hybridized supermodes basis.

Modulation Instability state
We restrict ourselves to the AS supermodes pumping scheme since we did not observe dynamics different from the single resonator case when exciting the S supermode family. Fig. 3(a,b) show the intracavity power evolution as a function of laser detuning in the resonator and supermode basis, respectively. It is numerically generated by exciting the system in a soft manner, i.e. adiabatically changing the laser detuning ξ = ω 0 − ω p from blue to red side of the AS hybridized resonance. Initial dynamics is found to be similar to the single resonator case. We observe the formation of primary combs in the AS supermode family followed by cnoidal waves (Turing rolls). The subsequent chaotic modulation instability stage [yellow area in Fig. 3(a,b)] already demonstrates a significant difference. Namely, the average intracavity power evolution in the second resonator, which is depicted by the red line, as a function of normalized detuning ξ/∆ω exhibits a local maximum inside the modulation instability area which corresponds to the efficient photon transfer to the S supermodes [violet curve in Fig. 3(b)]. At these values of detuning ( ξ/∆ω ≈ -1), we observe an enhancement of spectral components distinct from the modulation instability gain region. The mode number of the components correspond exactly to the distance from the pumped mode to the lower (S supermodes) parabola for a given value of the laser detuning, as described in [30]. This is a first signature of the interaction between the supermodes. Fig. 3(c,d) provide the underlying evolution of intracavity power (spatiotemporal diagram) in the supermode basis. The modulation instability region in the conventional basis does not differ for the single-particle dynamics. However, the supermode basis reveals that the transfer of photons to the S supermode family occurs after a certain detuning threshold. As follows from the spatiotemporal diagram of the AS state, it occurs in the developed AS supermodes modulation instability stage, where collision and annihilation of unstable coherent structures lead to the enhancement of wings in the optical spectrum [35] and thereby populates the modes in vicinity of the symmetric resonances.

Breathing state
As in conventional single resonator systems above a threshold pump power level, the modulation instability region is followed by the breathing solitons region [violet area in Fig. 3(a,b)]. Breathing originates from the Hopf bifurcation as demonstrated for the single resonator case [36]. It manifests itself as a periodic oscillation of localized coherent structures (similar to solitons on a finite background, such as Kuznetsov-Ma soliton [37,38] -a solution of the nonlinear Schrödinger equation), which radiates DWs at every cycle of oscillation. Fig. 4(a) shows the intracavity power evolution. Breathers in the photonic dimer exist in both resonators and oscillate in phase. However, the intracavity trace is found to be randomly deviating from the average because of the photon transfer to the S supermodes and, therefore, the generation of additional DWs.
The periodic oscillation of a coherent structure in slow time results in the appearance of a ladder of straight and equally-spaced lines on the nonlinear dispersion relation [39]. The presence of such ladder has been experimentally demonstrated (see [30] supplementary information) by reconstructing the comb spectrum with high resolution. Indeed, as follows from Fig. 4(c,e) the breathing frequency is given by the frequency offset between the lines. The same reasoning can be applied to the single resonator breathing states. There, it has been demonstrated experimentally that the breathing frequency linearly depends on the pump laser detuning [24]. Therefore, we can make a conjecture that the breathing occurs due to the photon transfer between the Kerr-shifted dispersion parabola and the first solitonic line given by the laser detuning, while the breathing frequency is the corresponding gap. In the points where the ladder crosses the AS supermodes parabola enhancement of the comb power is observed. Therefore, optical spectrum of a breather contains a set of sidebands [24].
In the supermode basis [see Fig. 4(b)] it becomes evident that the breathing occurs mostly in AS supermode families. Therefore, the breathing dynamics in the AS mode family does not show significant difference from the conventional breathing found in the single resonator case as follows from the nonlinear dispersion relation [ Fig. 4(d)], although it demonstrates significant differences in the resonator basis. Fig. 4(f)   ear dispersion relation for the S supermode family. The origin of the DWs which perturb the breathing state can be seen as an enhancement of the certain supermodes in the S family (µ ≈ ±70) in the places where the ladder from AS supermodes crosses the S parabola.

Soliton hopping state
The soliton hopping state recently predicted in the photonic dimer [30] is characterized by a periodic energy exchange between the coupled resonators in the presence of temporally-localized coherent structures. Interresonator oscillations have a frequency equal to the splitting between the supermode parabolas. The average power modulation is much stronger than in the breathing state, which leads to the enhanced sideband amplitudes in the optical spectrum [see Fig. 5(b)].
The spatiotemporal diagrams for the diagonalized sys-tem [see Fig. 3 Fig. 3(a,b)] coincides with the emergence of a localised coherent structure in the S supermodes family accompanied by a characteristic solitonic step in the average intracavity power evolution as follow from Fig. 3(b). This coherent structure is generated via the emerged FWM pathways depicted in Fig. 2(d). Soliton in AS supermode family acts in this case as a source of photons which triggers the parametric processes, thereby resonantly populating the S parabola in the vicinity of the 0 th mode (i.e. with the offset −∆ω) via the process #5 (â † as,µ1â † s,µ2âas,µ3âs,µ4 ). Energy conservation is ensured by populating supermodes offset by approximately +∆ω. Cascaded parametric process #2 (â † s,µ1â † s,µ2âs,µ3âs,µ4 ) populates the neighbouring S supermodes similarly to the CW-pumped single resonator. Therefore, we assume that the coherent structure generated in S supermodes is a GS (i.e. supermode dissipative Kerr soliton). Thus, the origin of the oscillatory behaviour can be seen as a time periodic interference of coherent structures living in different supermodes. Fig. 5(a) shows the dimer dynamics at fixed pump laser detuning, in the soliton hopping regime in AS and S supermodes representation. It can be obtained numerically by seeding the solitonic state in the AS supermode (see Appendix D) and further tune into the soliton hopping state. The average power exhibits small amplitude oscillations around a certain value. Periodic oscillations in slow time results into the series of sidebands (similar to Kelly-sidebands widely present in the modelocked lasers [40]) in the optical spectrum as has been shown in [30]. Corresponding nonlinear dispersion relation shows a ladder of lines similar to the breathing state discussed in Sec. III B 2, but the spacing between them is equal to the splitting between the DWs parabolas. The origin of the double maxima spectral sidebands is well seen in the supermode basis [see Fig. 5(b)]. They appear due to the different Kerr nonlinearity-induced shift of supermodes in the presence of inter-resonator detuning. White dots indicate the point where the ladder crosses dispersive parabolas and the continuing arrows indicate the corresponding spectral components enhancement. Both nonlinear dispersion relations depicted in Fig. 5(d,c) contain a corresponding DW parabola and the hopping ladder. Since the ladder crosses parabolas at slightly different mode numbers, sidebands have two maxima.

C. Experimental and numerical evidences of the deterministic single soliton generation
Deterministic generation of a single soliton state in optical microresonators is essential for a turn-key dissipative Kerr soliton-based broadband frequency combs generation. Indeed, passing the chaotic modulation instability stage, soliton arrangement inside the cavity can be arbitrary which leads to a non-homogeneous spectral profile due the interference of different solitonic components. One way to control and structure the soliton arrangement inside the cavity is to introduce a background modulation which leads to the generation of perfect soliton crystals [27]. However, the single soliton state is, nonetheless, difficult to achieve in this configuration.
Another way to naturally fall into the single soliton state has been described in [41]. It has been proposed to employ a strong avoided mode crossing with higher-order modes of the resonator, which leads to an extensive cavity Cerenkov radiation [21]. In this case, the soliton, being a line on the nonlinear dispersion relation [42], crosses the distorted cavity mode, which leads to effective photon transfer towards the higher-order modes [43]. In this way, every soliton acts as a "source" of DWs. Therefore, if the energy of the DWs is sufficient to perturb the solitonic states the number of solitons will decay towards unity, where the state will be stabilized.
Here, we present a deterministic version of this mechanism utilizing discovered properties of the photonic dimer [30]. Due to the more complex dispersion landscape, the single soliton generation process does not require any additional interaction with higher-order mode, even though it is shown to be enhanced for certain supermodes due to the underlying symmetry [44]. Indeed, the periodic intra-resonator field enhancement due to the crossing of the S supermode parabola is found to be sufficient to trigger the process discussed in [41]. In order to verify this claim, we investigate both numerically and experimentally the GSs generation. Fig. 6 shows the qualitative comparison of numerical and experimental phase diagrams. As follows from the numerical simulations of coupled LLEs 1, single soliton generation occurs when passing a threshold power of 0.9 W. A similar result follows from the experimental investigations. A schematic of the experimental setup is shown in Fig. 6(a). Single GSs are generated with an integrated Si 3 N 4 photonic dimer driven by external cavity diode laser. The CW pump is amplified by an erbium doped fiber amplifier to achieve the power level need for the investigation of the threshold of the process. After passing a fiber polarization controller needed to guarantee that the solitons are generated in a single polarization mode family, the light is injected in the photonic dimer. The generated light in filtered by a fiber Bragg grating and recorded with a fast oscilloscope. Figure 6(c) shows 50 traces of the generated combs power at different values of the pump power as a function of the laser detuning from the position of AS resonance. The central frequency of a widely tunable external cavity diode laser has been controlled by the piezo-tuning technique. Other details of the experimental measurements can be found in [30], Methods section.

IV. CRITICALLY COUPLED RESONATORS: SPLIT DISSIPATION (PT -BROKEN)
Passing through an EP, which exhibits a singularity of nonlinear interactions efficiency, we enter the domain of split dissipation (PT -symmetry broken phase), which exhibits drastically different dynamical features. We study soliton generation in this region and show that PTsymmetry breaking leads to soliton generation in either cavities. The soliton localization can be switched by increasing the pump power and thereby flipping the broken PT -symmetry. Four distinct dynamical states are identified, we observe among them on-demand perfect soliton crystals generation, which can be a promising alternative to the existing technology relying on the resonator's mode interaction [27].
A. Nonlinear coupling coefficients and their divergence at the exceptional point The system of Eq. 1 can be diagonalized to Eq. 7 in the PT -symmetry broken state in the same way as in the PT -symmetric state, by applying a transformation matrix T m such that T m M T −1 m is diagonal (see Eq. 2). The diagonalization enables the computation of the nonlinear coupling coefficients describing the nonlinear interaction between the supermodes. The non-vanishing values of the nonlinear coupling coefficients are displayed in Fig. 7(c) as a function of J, where the external coupling κ ex,1 is varied quadratically with J to satisfy the critical coupling condition (Eq. 3). The nonlinear coupling coefficients are normalized by a real factor equal to the square of the norm of T m (1, 0) T . Indeed, the system is non-Hermitian such that T m does not preserve the norm (the transformation is not Unitary). This means the variables |A s,µ | 2 , |A as,µ | 2 are not proportional to the number of photons. The normalization allows the interpretation of t i g K /N as a rate per photon.
The nonlinear coupling coefficients behave differently in the regime of split dissipation with d = 0 than in the regime of split resonance: t 3 , associated with even interband processes, vanishes completely, while t 2 , t 4 , (t 1 ) become purely imaginary (real). These are related to odd inter-band (even intra-band) processes. At the EP (J/κ 0 = 0.5) they diverge. The normalization corrects the divergence of t 1 and t 2 , but not of t 4 , which exhibits a singularity at the EP. The constant value t 1 /N = 1 signifies that the rate of the intra-band FWM does not vary with J in this critical coupling condition.
Although the linear concepts of PT symmetry are use-ful to understand the dynamics of the PT -symmetry broken state, we observe that the solitonic states presented in the next section are not distributed spatially according to the supermode basis like in the PT -symmetric regime. On the contrary, they are distributed in the resonator basis, which is diagonal with respect to nonlinearity.

B. Phase diagram: inter-resonator coupling vs pump power
We numerically explore the phase diagram under the condition of critical coupling in the nondegenerate dissipation regime [see Fig. 7(a)]. In Fig. 7(c), κ ex,1 is varied with J in the way that the critical coupling condition (Eq. 3) is satisfied across the phase diagram and the dimer is in a state of broken PT -symmetry (split dissipation). An EP is found at J = 0.5κ 0 [see Fig. 1(d)].
We selected 14 values of pump power distributed logarithmically from 0.01 W to 1.5 W and 13 values of interresonator detuning distributed linearly from 0κ 0 to 2.4κ 0 . For each set of parameters, we employ the conventional soliton generation scheme by scanning the resonance from blue to red-detuned side. The spatiotemporal and spectrum evolution diagrams in the resonator basis is used to identify the stable soliton state that is generated during the scan. Colors in Fig 7(a) correspond to different stationary states attainable in the split dissipation regime. Thus, for different points on the phase diagram, the value of detuning is not the same. If several stationary states are identified, we choose the first state in the soliton existence range. The phase diagram is averaged over 3 realizations and the pump laser frequency is swept at the speed 1 10 κ 2 0 /2π, corresponding to a change of frequency κ 0 every 10 photon lifetimes (2π/κ 0 ).

Dynamical states in split dissipation (PT -broken) regime
We differentiate four dynamical states in this regime: multi or single soliton in resonator 1 (S1, blue region), soliton crystal in resonator 1 (C1, cyan region), coexistence of periodic coherent structures in resonator 1 and soliton in resonator 2 (S2/C1, green region), and soliton in resonator 2 (S2, yellow region). The parameter regions enabling their generation are coloured on the phase diagram and their characteristic intracavity intensity profile are shown in the insets. White region refer to the absence of solitonic states.
At weak inter-resonator coupling (J < 1 2 κ 0 ), the system qualitatively follows the single resonator dynamics and features the S1 state, where DKSs exist in resonator 1 while resonator 2 only features their low-power projection. Increasing the inter-resonator coupling, dynamical regions corresponding to states C1, S2/C1, and eventually S2 are accessed.
The states are almost exactly partitioned in the resonator basis. For example, state S1 is confined in resonator 1 although a negligible amount ( 1%) is found in resonator 2. That is, the field amplitude distribution between the resonators does not follow the supermode distribution which is given by linear analysis. We assume that the nonlinearity changes the field distribution of the supermode, making them localized in the resonators. These supermodes are referred to as high-loss (confined in resonator 1) and low-loss (confined in resonator 2).

Parametric switching of the soliton-localisation
In the range of 1.2κ 0 J 1.8κ 0 , the four stable states can also be accessed by changing the pump power. Fig. 7(b) shows the intracavity power evolution as a function of the laser detuning for pump power levels corresponding to four dynamical states in this range of J. Therefore, gradually increasing the pump power, states S2, S2/C1, C1, and S1 can be sequentially accessed.
In this process, the soliton-localization switches from resonator 2 to resonator 1. Linear analysis predicts that the resonant soliton will be confined in the low-loss supermode while the red-detuned CW background will be in the high-loss supermode [see Fig. 1(c)] as is the case for state S2 [inset of Fig. 7(a)]. However, despite the field distribution predicted by the linear analysis, the presence of nonlinearity in the system introduces a mechanism allowing for the parametric switching between the cavities. Indeed, at low pump power, only the low-loss supermode has a quality factor sufficient for the soliton generation, resulting in state S2. At higher pump powers, both supermodes can sustain a coherent structure, leading to the coexistance of soliton and periodic coherent structure that has been observed in a limited intermediate range of parameters. Above a threshold, solitons are not generated in resonator 2. Moreover, in the C1 and S1 state, the parametric gain is able to compensate the difference of losses between the supermodes, and invert the PT symmetry: the parametric gain (via intra-band FWM) is larger in the supermode localized in resonator 1 than in the other supermode making the state of broken PT -symmetry flipped in comparison to the linear regime for longitudinal modes with µ = 0.
We note that no specific solitonic state was found at the EP. We suppose that the Kerr shift lifts the degeneracy between the two supermodes. However, an extensive investigation of the soliton generation in the close vicinity of the EP is beyond the scope of this study.

C. Deterministic soliton crystal and efficient comb generation
In this section, we study the generation of state C1, which corresponds to the deterministic generation of a soliton crystal [27]. Fig. 8 shows the numerical simulations of Eq. 1 setting the pump power to 0.2 W and J/2π = 75 MHz. Fig. 8(a) shows the intracavity power in both resonators as a function of laser detuning. The incident light couples into both high-and low-loss, supermodes of the system simultaneously, such that the low-loss supermode features a chaotic regime while the high-loss supermode remains in the cnoidal wave regime [ Fig. 8(c,e)]. After passing a critical detuning (∼ 7κ 0 ), resonator 2 leaves the chaotic regime without any coherent structures generated while cnoidal waves of resonator 1 transition into a soliton crystal state. Fig. 8(b) shows that the crystal state at detuning I [ Fig. 8(a)] is perfect [27] with more than 100 dB of extinction over almost the full existence range.
This state is known to exhibit a high conversion efficiency due to the high occupancy of the resonator 1 as shown in Fig. 8(b,d). Fig. 8(f) shows the output power in the pump mode (µ = 0) and comb modes (µ = 0). We observe that the perfect soliton crystal formation leads (a) Schematic phase diagram represented as a function of input power and J. Four states are identified: soliton in resonator 1 (S1, blue region), soliton crystal in resonator 1 (C1, cyan region), coexistence of soliton in resonator 2 and periodic coherent structure resonator 1 (S2/C1, green region), and soliton in resonator 2 (S2, yellow region). The insets show the intracavity field intensity in the resonators 1 (blue) and 2 (red), at the parameter indicated by the grey dot. White color refers to the absence of stable solitonic state. (b) Intracavity power traces at different input powers for J/2π = 70 MHz. (c) Nonvanishing values of the nonlinear coupling coefficients and normalizing factor N presented as a function of J along the critical coupling condition. The crossing of the EP at J/κ0 = 0.5 corresponds to a singularity in coefficients t1, t2, and t4. Parameters for the simulations are κ0/2π = 50 MHz, κex,2 = 0, δ = 0, D1/2π =180 GHz, D2/2π = 4 MHz. to a conversion efficiency higher than 75%. Also, we note that the pump is almost completely absorbed by the system, such that an effective nonlinear critical coupling is achieved. According to [27], the soliton crystal is generated deterministically when the pump power is below the threshold to avoid spatiotemporal chaos under the condition that modal crossings with higher-order modes trigger background modulation. Here, we observe deterministic soliton generation in the absence of modal crossings.

D. Bright-dark solitons coexistence and their interaction with periodic coherent structures
We perform and analyse a simulation with J/2π = 75 MHz, P in = 0.1 W, as shown in Fig. 9 in order to generate S2/C1 state depicted by green in the phase diagram [ Fig. 7(a)]. The power trace [ Fig. 9(a)] shows the presence of a step in each resonator. The spectrum and temporal intensity at detuning I are shown in Fig. 9(b). A soliton exists in resonator 2, while background modulation reminiscent of C1 state are present in both resonators. The comb modes in both resonators are excited in this state, hinting at nonlinearly-induced PT -transition that restores the PT symmetry in the comb modes [45,46]. Spatiotemporal diagrams shown in Fig. 9(c,e) as a function of the laser detuning indicate that the S2/C1 state decays into a S2 state after the end of the soliton existence range in resonator 1. After transitioning to the S2 state, the field in resonator 1 acts as a source for resonator 2, resonantly supplying additional energy to the solitonic state. This results into a coexistence of a bright and dark solitons synchronously rotating in the resonators. This situation is similar to a dual fiber loop arrangement presented in [47] but in the limit of equal cavities.
While the existence of a periodic coherent structure in resonator 1, bright-dark soliton pair can be generated as well. A simulation at a fixed detuning starting from the initial conditions I [ Fig. 9(b)] is shown in the spatiotemporal diagrams Fig. 9(d,f). We observe that the soliton pair is bounded by the effective nonlinear potential induced by the periodic structures in the neighbouring cavity. The spatiotemporal diagram depicts the possibility for the bright-dark soliton pair to tunnel from one potential unit cell to another interacting with their boundaries in an oscillatory manner. In addition to the fast oscillations, a random walk of the pair is observed at a slower timescale.

V. DISCUSSION AND CONCLUSION
In this article we investigate nonlinear dynamics in a driven-dissipative photonic dimer exhibiting an exceptional point. We analyse the generation of dissipative Kerr solitons on both sides of the exceptional point which acts as a demarcation of the dimer critical coupling conditions. These two regimes are the split resonance regime (with the preserved symmetry) and split dissipation regime (with broken symmetry), as found in conventional PT -symmetric systems with gain and loss. Each regime exhibits unique nonlinear dynamics not found in the single resonator.
In the split resonance regime, which has been substantially discussed in [30], we observe that dimer solitons can be generated in either supermodes, however, only the antisymmetric one exhibits non-conventional soliton dynamics related to the emerging efficient four-wave mixing pathways. The dynamics is conveniently expressed in the supermode basis, for which we developed the concepts of inter-band four wave mixing. Supermode representation reveal that, despite the complexity of the dimer system, we are able to separate conventional single-resonator soliton dynamics from the dispersive waves emerging in another supermode. Breathing state of the photonic dimer in the supermode basis appears remarkably similar to its single-resonator counterpart except for a small perturbation. In this case, the intra-resonator power of both cavities oscillates in phase. Rapid and counter-phase power oscillations (soliton hopping) are observed above a threshold pump power, originating from the generation of synchronized solitons in both supermodes. We highlight the fact that all the dynamics appearing in the split resonance regime can be well understood in the supermode representation. The same does not apply to the regime of split dissipation. The absence of the resonance splitting implies the pumping of both supermodes simultaneously. Therefore, the most convenient representation in this case is the conventional resonator basis which exhibits the broken PT -symmetry of the system. Satisfying the critical coupling condition, we impose different loss rates on the two resonators so resonator 1 becomes substantially overcoupled. We observe the generation of four different stable solitonic states localized in either or both resonators. In these states we observe: the synchronization of bright-dark soliton pairs (in resonators 2 and 1, respectively), the interaction of periodic coherent structures in resonator 1 with the bright-dark soliton pair, the deterministic generation of soliton crystal states with more than 75% pump conversion efficiency into the comb lines, and bright solitons in resonator 1. Thereby, the pump power enables the parametric switching of the soliton localization between resonator 1 and 2. The switching is induced by flipping the broken PT symmetry so the lossy resonator 1 has more gain in the comb modes than resonator 2. We suspect nonlinearly induced PT transition to take place in the intermediate regime [45,46]. Moreover, we compute the nonlinear coupling coefficients between the supermodes. We normalize the divergence of two of them. One however seems to exhibit a singu-larity at the exceptional point, which could demonstrate enhanced sensitivity in its vicinity [48].
Concluding, we would like highlight the abundance of nonlinear dynamics occurring in the simplest element of soliton resonator lattices -photonic dimer. Despite the formal similarities with the single-mode dimer systems extensively studied in the context of non-Hermitian photonics, massively multimode nonlinear dimer exploiting another degree freedom reveals unprecedented solitonic states and dynamics, which cannot be covered in depth in one study. The fundamental aspects of this systems -two ideal coupled resonators -can be of interest far beyond the photonics community. In this work, we define the Fourier Transform of a function f as

VI. ACKNOWLEDGEMENTS
Thereby, the amplitude A µ of the cavity mode with index µ is defined from the slowly-varying field amplitude corotating with angular frequency D 1 as This relation can be inverted to yield These definitions allows us to conveniently switch the nonlinear interaction from a temporal to a Fourier perspective [31]: Appendix B: Critical coupling with inter-resonator detuning A non-vanishing value of inter-resonator detuning does not allow exceptional point (unless the coupling is complex). However, it is still important to find critical coupling conditions in the split resonance case in the general case including inter-resonator detuning. Indeed, two resonators will generally have different resonance frequencies such that δ = 0. Moreover, the use of integrated heaters [30] or piezo-electric components [49] can make the inter-resonator detuning a control parameter for the system. Practically, a numerical method is easier to implement. Here, we would like to discuss qualitatively the dependence of the critical external coupling to the system parameters In single resonators, critical coupling is achieved when the coupling rate matches the loss. This principle is maintained when more resonators are added. Indeed, Eq. 4 originates in the fact that each supermode is equally distributed in each resonator. Thus, the effective coupling to the through port and drop port are halved, whereas the intrinsic loss remains constant, because the circulating photons see the each port only half of the time. Equalizing the loss and the effective coupling between the waveguide and the supermode results in 1 2 κ ex,1 = κ 0 + 1 2 κ ex,2 , as given by Eq. 4. Similarly, the critical external coupling in the case of split dissipation is relatively high because coupling is achieved in the state with a low loss, which is only sparsely confined in resonator 1, such that the effective coupling is equal to the loss.
When the inter-resonator detuning has a non-vanishing value, the supermodes are unequally distributed between the resonators. The effective coupling of each supermode to the through port depends on the hybridization, which strongly depends on δ. For δ > 0, the AS supermode, which has a higher resonance frequency, is rather confined in resonator 1. Thus, the critical coupling condition is achieved for κ ex,1 < 2κ 0 + κ ex,2 . Inversely, the S supermode requires higher κ ex,1 to achieve critical coupling. Fig. 10(a) shows example curves of the critical coupling condition to the AS supermode for different parameters in split resonance regime as a function of the normalized inter-resonator detuning d. The dashed line shows the critical coupling condition to the S mode for one set of parameter. Four critical coupling conditions are marked and the intracavity fields in resonator 1, resonator 2, and the transmission spectrum at these sets of parameters are shown in Fig. 10(b-d), respectively.
Appendix C: Projection into supermodes

Derivation with complex inter-resonator detuning
The same analysis can be adapted when drive and dissipation are included. However, the inter-resonator detuning takes a complex component accounting for different loss rates between the resonators. We therefore start with and define ∆κ ex to satisfy κ ex,1 = κ ex + 1 2 ∆κ ex κ ex,2 = κ ex − 1 2 ∆κ ex , such that we can define the complex inter-resonator detuning δ c = δ − i 1 2 ∆κ ex . We apply the following non-unitary transformation to Eq. C7 with rescaled effective parameters: Hence, the validity of the vanishing S (AS) supermodes power implies that strongly coupled mulimode dimer can sustain conventional single resonator dissipative Kerr solitons in the AS (S) supermode family. The rescaling of the parameters accounts for the larger effective cavity volume if we consider the parameter g K for example.