Nonlinear topological symmetry protection in a dissipative system

We investigate experimentally and theoretically a system ruled by an intricate interplay between topology, nonlinearity, and spontaneous symmetry breaking. The experiment is based on a two-mode coherently-driven optical resonator where photons interact through the Kerr nonlinearity. In presence of a phase defect, the modal structure acquires a synthetic Möbius topology enabling the realization of spontaneous symmetry breaking in inherently bias-free conditions without fine tuning of parameters. Rigorous statistical tests confirm the robustness of the underlying symmetry protection, which manifests itself by a periodic alternation of the modes reminiscent of period-doubling. This dynamic also confers long term stability to various localized structures including domain walls, solitons, and breathers. Our findings are supported by an effective Hamiltonian model and have relevance to other systems of interacting bosons and to the Floquet engineering of quantum matter. They could also be beneficial to the implementation of coherent Ising machines.


I. INTRODUCTION
Topology, nonlinearity, and spontaneous symmetry breaking are key concepts of modern physics.Topology is concerned with the study of properties that are invariant under continuous deformations, hence inherently robust against fluctuations and disorder [1,2].Nonlinearity is the source of complexity, underpinning attractors, chaos, self-organization, and self-localization (solitons) [3,4].Spontaneous symmetry breaking (SSB) underlies much of the diversity observed in Nature [5][6][7], from ferromagnetism [8] to embryo development [9].Here we report on an experiment in which these three concepts are fundamentally intertwined.Nonlinearity and SSB combine to endow our system with a topological invariant; in return that invariance protects the underlying symmetry, enabling a robust realization of SSB that is immune to perturbations.
The discovery of topological states of matter in solidstate materials has revealed the central role of topology in the classification of quantum matter [1,10].Besides, recent progress in designing and controlling synthetic lattice systems has allowed for the exploration of topological properties in a broad class of physical contexts [11,12], including ultracold gases [11], photonic crystals [12,13], mechanics [14], and systems with synthetic dimensions [13,15].Topological robustness has potential wideranging applications, from fault-tolerant quantum computers [16,17] to lasers capable of single-mode operation at high power [18].While topological band structures concern the properties of single-particle Bloch states, interesting behaviors were observed by combining these with classical nonlinearities [19].This includes, nonexhaustively, edge solitons, nonlinearity-induced topological phase transitions [20,21], novel symmetry-protected phases [22], or quantized Thouless pumping of solitons [23].
The results presented in this Article demonstrate a new facet of the potent bridging of topology and nonlinearity, and how it can confer robustness to the realization of SSB.Upon SSB, a system bifurcates to states with a broken symmetry even though the equations of motion retain that symmetry [5][6][7].Intense interest has been paid to this process in recent years, especially in optical microresonators where it has been studied, e.g., for optical memories [24], nonreciprocal propagation in integrated photonic circuits [25,26], dynamical control of laser directionality and chirality [27], or sensors with divergent sensitivities [28,29].Experiments on SSB have also been used for fundamental studies of domain walls [30], of the In that regime (P2 SSB), the phase of the y-mode changes by π at each resonator roundtrip (top right), corresponding to a periodic roundtrip-to-roundtrip swapping of the hybrid mode occupations (bottom right).tR is the roundtrip time of the resonator.
universal unfolding of the pitchfork bifurcation [31,32], or of the Bose-Hubbard dimer model [33].These experiments most often rely on an exchange symmetry (Z 2 ) between two states.Invariably, however, this symmetry is affected by manufacturing imperfections, experimental biases, or other non-idealities.As a consequence, instead of a spontaneous, random selection between the two states, one state is statistically favored over the other [32,34], hindering the experiments.Here we show how the combination of nonlinearity and topology can solve that problem.Specifically, we report on the experimental realization of an optical resonator with a synthetic modal structure characterized by a Möbius topology protected by the Kerr nonlinearity.That protection leads to an exact exchange symmetry, which enables robust, bias-free SSB that does not require any fine tuning of the system parameters.The generation of truly random binary sequences based on this platform confirms the robustness of the exchange symmetry.Additionally, we reveal a new class of topological symmetry-broken localized structures characterized by a rapid periodic exchange of their symmetry.Our observations are supported by theory, which reveals in particular the key role played by degenerate parametric processes.This connects our work with some recent implementations of so-called coherent Ising machines [35,36] that could potentially benefit from our symmetry-protection scheme.These benefits could extend to experiments on Bose-Einstein condensates.Indeed, the Kerr nonlinearity stems from interaction processes among photons, in direct analogy with the nonlinearity inherent to Bose-Einstein condensates and described by the Gross-Pitaevskii equation [37].We also present an effective Hamiltonian approach, similar to that used for the Floquet engineering of quantum mat-ter [38][39][40], to explain aspects of the emergence of the symmetry in our system.Finally, our work could have relevance to the symmetry restoration techniques used in mean-field approaches of quantum many-body systems in nuclear, atomic, and molecular physics [41].

II. PRINCIPLE OF SYMMETRY PROTECTION
Our system consists of a passive nonlinear (Kerr) optical resonator coherently-driven by a single-frequency field and that presents two distinct polarization mode families.Hereafter ψ x and ψ y denote the modes' classical electric field (complex) amplitudes inside the resonator scaled such that |ψ x | 2 and |ψ y | 2 represent the number of photons in each mode.In this configuration, the nonlinearity is well known to induce a polarization SSB bifurcation contingent on the existence of an exchange symmetry between the two modes [42,43].That symmetry is typically realized by being as close as possible to conditions where the modes are equally driven and equally detuned from the driving field frequency [25,26,32].Here, we deviate from this paradigm.Specifically, we introduce a localized π-phase-shift defect between the two mode families of the resonator.As a result, their resonance frequencies are separated by half a free-spectral-range (FSR).In these conditions, the exchange symmetry cannot be realized as described above: the single-frequency driving field can only be resonant with one of the modes, assumed here to be the x mode, hence only that mode is driven and the roundtrip phase detunings differ by π.These different aspects are illustrated schematically in Fig. 1(a).
Consider now an alternative description of this res-onator in terms of the hybrid mode amplitudes ψ + and ψ − defined by the unitary transformation Because of the π-phase defect, the amplitude of the y mode flips its sign with respect to the driven mode at each roundtrip in the resonator, ψ y → ψ y e iπ = −ψ y .
From the above relations, it is clear that this is associated with a roundtrip-to-roundtrip periodic swapping of the + and − hybrid mode amplitudes, ψ + ψ − [right side in Fig. 1(c)].The field must circulate two roundtrips in the resonator to be restored to its original state, revealing that the hybrid modes embody the topology of a Möbius strip [44][45][46][47].In terms of those hybrid modes, the resonator shown in Fig. 1(a) is thus equivalent to that represented in Fig. 1(b).Note that the Möbius topology only exists in the synthetic modal dimension rather than in real physical space [48].A consequence of the single-sided, non-orientable nature of the Möbius topology is that the two hybrid mode amplitudes will necessarily accumulate the same changes over two resonator roundtrips.The resonator therefore exhibits an exchange symmetry, ψ + ψ − .That symmetry, coupled with the nonlinearity, enables SSB between the hybrid modes.
Here, it is important to realize that SSB in our system is fundamentally intertwined with the generation of the y field.As can be seen from Eqs. (1), a symmetric state for the hybrid modes, ψ + = ψ − , is necessarily associated with a y field of zero amplitude.By contrast, a symmetry-broken hybrid state, ψ + = ψ − , requires a nonzero y field.By symmetry, the mirror state, ψ + ψ − , must also exist, corresponding to a y field with a phase differing by π, ψ y −ψ y .Keeping in mind that the yfield is undriven, it must be generated internally through nonlinear interactions past a certain threshold of detuning and/or driving strength.What appears as an SSB bifurcation for the populations |ψ + | 2 , |ψ − | 2 of the + and − hybrid modes can therefore be seen as the threshold of parametric oscillation for the y mode [49] [left side in Fig. 1(c)].However, our system is not just a standard parametric oscillator.Above the bifurcation threshold, the underlying Möbius topology imparted by the π phase defect forces the system to alternate between the two symmetry-broken states one roundtrip to the next, and to display an overall periodicity of two roundtrips.As a reminder to these specific characteristics, we will refer to this regime as P2 SSB.One can also note that, in the P2 SSB regime, the amplitudes of the two hybrid modes over one resonator roundtrip become indistinguishable from the amplitude of one of those modes over two subsequent roundtrips.An alternative interpretation is thus that SSB now occurs between subsequent resonator roundtrips, in essence breaking the time translation symmetry, in a way reminiscent of period doubling instabilities [50][51][52].

III. THEORETICAL ANALYSIS
The conclusions of the previous Section are based on the premise that the phase defect seen by the y mode at each resonator roundtrip is exactly π.To understand how the Möbius topology and the associated symmetry can persist in presence of deviations from this condition, we must examine some aspects in more detail.To that end, we introduce the Hamiltonian ĤNL describing the Kerrmediated interactions taking place along the resonator roundtrip, Here g is the single photon-induced Kerr frequency shift of the resonator (g > 0) [53] while â † σ (resp.âσ ) represent creation (annihilation) operators for the mode σ = x, y.These operators satisfy bosonic commutation relations, [â σ , â † σ ] = δ σ,σ .The three terms on the right-hand side of Eq. ( 2) describe, respectively, on-site (Bose-Hubbard) interactions (self-phase modulation), inter-site interactions (cross-phase modulation), and pair hopping (parametric four-photon mixing).The latter corresponds to two photons in the x mode being converted to two photons in the y mode (or vice-versa); since the y-mode is undriven, this term is required as the source of y photons.Without loss of generality, we will assume A = 1, B = 2/3 and C = 1/3.These are the values found in silica (which is the Kerr material used in our experimental demonstration presented below) for modes that are linearly polarized [54,55].Generalization to other polarization states will be discussed in Appendix.
To understand specifically the evolution of the phase of the y mode, we examine the equation of motion for ây , i.e., dâ y /dt = i[ ĤNL , ây ]/ .In the classical limit, â † y → ψ y , we get where the right-hand side has the form of the interaction terms found in the Gross-Pitaevskii/nonlinear Schrödinger equation ( * denotes complex conjugation).Separating amplitude and phase, ψ σ = |ψ σ |e iφσ , we find that the phase φ y of the y mode obeys (4) The three terms in the above equation represent phase shifts induced by the Kerr nonlinearity.If we restrict ourselves to high-Q resonators, these can all be assumed 1.Indeed, nonlinear effects will only have an important role when close to resonance.In these conditions, we can use a single Euler step to approximate the integration of Eq. ( 4) over one resonator roundtrip time t R .To obtain the total change in the phase of the y mode from one roundtrip to the next, we still need to consider the role of the boundary conditions, which add a linear contribution due to the resonator detuning (see Appendix for more detail), leading to Here m is the roundtrip index, δ 0 = δ x is the roundtrip phase detuning of the driving field (modulo 2π), also assumed 1 (δ 0 > 0 means red-detuned), and the corresponding detuning for the y mode is defined through δ x − δ y = π + δ π , with δ π representing deviations from an exact value of π.We can observe that all terms above are in fact π, except π itself (a similar argument for the x mode would show that φ . Moreover, the pair hopping (parametric) contribution (∝ C), which is the only phase sensitive term, is π periodic in φ y .This is a well known characteristics of degenerate parametric oscillators [56].It follows from Eq. ( 5) that the only stable configuration is for the phase of the y mode to step by exactly π, one roundtrip to the next.That solution acts as an attractor for the system, irrespective of the value of δ π .Effectively, the intracavity fields will self-adjust through the dissipative nonlinear dynamics to guarantee a roundtrip-to-roundtrip y mode phase step of π.In this way, the Möbius topology of the resonator is intrinsically robust and protected by the nonlinearity.
As a second step, we examine how the Möbius topology induces an exact exchange symmetry ψ + ψ − .To that end, we introduce an explicit asymmetry and show how it cancels out.Specifically, we consider a difference in detunings between the + and − hybrid modes.The Hamiltonian Eq. ( 2) can be expressed in terms of these modes with creation (annihilation) operators â † ± (â ± ); these are related to the corresponding x, y operators by the same unitary transformations as Eqs.(1).We find where A = 2/3 and B = 4/3.For convenience, let us define a set of angular momentum (Schwinger) operators as ) These operators satisfy commutation relations [ Ĵµ , Ĵν ] = i µνλ Ĵλ similar to spins ( µνλ is the Levi-Civita symbol), while N counts the total number of photons.In this part of our analysis, we assume a lossless resonator (conservative limit) so N is constant.With these notations, a difference in angular frequency detuning δΩ between the two hybrid modes enter the Hamiltonian as δΩ Ĵz , so that the total Hamiltonian of the system reads where a = A − B , b = (A + B )/4, and c = A /2.We now proceed by seeking an effective (Floquet) Hamiltonian Ĥeff [38,57] which captures stroboscopically the dynamics of a full period of the Möbius configuration depicted in Fig. 1(b), i.e., every two actual roundtrips of the resonator.The swapping of the two hybrid modes at each roundtrip is introduced through the discrete application of the operator Ûswap = σx = e i(π/2)(σx−1) (with σx a Pauli matrix), leading to The first three factors on the right hand side of this expression define Ĥ1 which can be obtained exactly as [58] Ĥ1 = e i π 2 σx Ĥasym e −i π 2 σx .
Noting that the swapping operation is a single-particle process, we can substitute Ĵx / for σx /2 in the exponentials above and then use the Baker-Campbell-Hausdorff formula.Specifically we have e iπ Ĵx/ Ĵz e −iπ Ĵx/ = − Ĵz while all other terms of Ĥasym are left unchanged, leading to Comparing the expressions for Ĥasym , Eq. ( 8), and Ĥ1 , above, we observe that the sign of the asymmetric term is reversed.It is also clear that Ĥasym commutes with Ĥ1 , which from Eq. ( 9) leads to Ĥeff = ( Ĥ1 + Ĥasym )/2.It follows that the asymmetric term cancels out exactly over two resonator roundtrips, and we have Ĥeff = ĤNL .Our analysis therefore shows that the effective Hamiltonian, which describes propagation every two resonator roundtrips, including two swaps of the hybrid modes, only contains the original -symmetric -nonlinear interactions; it is immune to fluctuations and asymmetries in detunings as foreseen in our qualitative description of the previous Section.This confirms the existence of an exact exchange symmetry for the hybrid modes, ψ + ψ − , and the topological robustness of the P2 SSB regime.
We are now in a position to derive mean-field equations of motion for the classical fields ψ + , ψ − based on the effective Hamiltonian Ĥeff .We have already established that the dissipation guarantees an exact overall roundtrip phase jump of π for the y mode, so we do not need to consider the influence of δ π , the deviation from an exact π phase defect, on the swapping operator Ûswap .Still, δ π contributes to the linear detunings through the relation δ 0 − δ y = π + δ π .This effect can be taken into account by the following roundtrip phase detuning operator, Note that, for the y mode, we have discarded the π phaseshift that is already incorporated in Ĥeff .The contribution of the detunings can be distributed along the resonator roundtrip and combined with the effective Floquet Hamiltonian, leading to, in the classical limit, To these equations, we have added phenomenologically losses (dissipation), associated with the linewidth ∆ω of the resonances, and driving, with strength S and external coupling fraction θ.In the conservative limit applied above, these terms are not included directly.An alternative classical derivation of the mean-field equations ( 13)-( 14), taking into account these dissipative effects, is presented in Appendix for completeness.That derivation confirms in particular that the Möbius topology leads, as for the detunings, to a robust symmetrization of the driving.Note that the symmetry of the driving can also be interpreted as stemming from the absence of driving of the y mode; under these conditions, Eqs.(1) indeed entails S + = S − = S = S x / √ 2, where S x is the driving strength of the x mode only.Additionally, we can observe that a deviation δ π from an exact π phase defect simply causes a common shift in detuning as well as some linear coupling between the modes, none of which alters the exchange symmetry of the equations.
Apart from the absence of kinetic terms, the role of which will be discussed in Section VI, Eqs. ( 13)-( 14) have the form of coupled, driven, damped, Gross-Pitaevskii/nonlinear Schrödinger equations.These equations are well known to exhibit SSB, provided the detunings and the driving strengths of the two modes are identical [25,26,29,32,42,43].Contrary to previous implementations, in our system, and thanks to the presence of the π phase defect, this condition is automatically and robustly satisfied.Finally, while the equations above describe the stroboscopic mean-field evolution of ψ + and ψ − two resonator roundtrips at a time, it is important to keep in mind that the amplitudes of the two fields swap at each resonator roundtrip due to the underlying Möbius topology.

IV. EXPERIMENTAL OBSERVATION OF P2 SSB
To demonstrate experimentally the topologicallyprotected P2 SSB regime, we have used, as modes x and y, the two polarization modes of a ring resonator constructed from single-transverse mode silica optical fiber.This platform enables straightforward implementation of both the π phase-shift defect and the projection of the x and y modes onto the + and − hybrid modes through manipulations of the polarization state of light with optical fiber polarization controllers (FPCs) [59].Note that, for x and y modes that are orthogonally linearly polarized, the + and − hybrid modes defined by Eqs.(1) correspond to circular polarization states of opposite handedness.
A schematic of the experimental setup is depicted in Fig. 2. As we used slightly different resonators with different specifications for different measurements (in part, to demonstrate the universality of the methodology), we list here the parameters corresponding to the results presented in this Section.Variations between setups will be highlighted as appropriate.To facilitate comparisons, experimental results will be presented in terms of normalized parameters.Specifically, we will be referring to the normalized roundtrip phase detuning ∆ = δ 0 /α and the normalized driving power X = gt R θS 2 x /α 3 , where α = t R ∆ω/2 is the dissipation rate of Eqs. ( 13)- (14).
Our first resonator is built around a fiber coupler (beam-splitter) that recirculates 90 % of the intracavity light.This coupler is also used for the injection of the coherent driving field (θ = 0.1).Another 1 % tap coupler extracts a small fraction of the intracavity field for analysis.Both couplers are made up of standard SMF28 (Corning) optical fiber.The rest of the resonator is made up of a highly nonlinear, low birefringence spun fiber (iXblue Photonics) [60].At the 1550 nm driving wavelength, that fiber exhibits normal group-velocity dispersion (corresponding to repulsive interactions, or a positive effective mass), which has been selected to avoid  scalar parametric instabilities [55].Overall, the resonator is 12 m-long, corresponding to a roundtrip time t R of 57 ns (FSR = 17.54MHz), and has a resonance linewidth ∆ω/(2π) = 650 kHz.The resonator is driven by a 1550 nm-wavelength continuous-wave laser (NKT Photonics) with an ultra-narrow linewidth < 1 kHz, ensuring coherent driving.The laser frequency can be varied with a piezoelectric transducer, which provides the ability to scan the detuning.Alternatively, we can use the transducer to lock the detuning at a set value via a feedback loop using the technique discussed in Ref. [61].The driving laser is intensity-modulated with a Mach-Zehnder amplitude modulator driven by a pulse pattern generator to generate flat-top 1.1 ns-long pulses with a period matching the resonator roundtrip time t R [62].Before injection into the cavity, the driving pulses are amplified with an erbium-doped fiber amplifier followed by an optical bandpass filter, which reduces amplified spontaneous emission noise.
The setup includes three FPCs.The first one is placed before the input coupler (FPC in ) to align the polarization state of the driving field with one of the two principal polarization modes of the resonator, henceforth defining the x mode.Note that, because of residual birefringence, these principal modes correspond in practice to polarization states which evolve along the fiber but map onto themselves at each roundtrip.The second FPC, placed inside the resonator (FPC π ), defines the π phase-shift defect.We adjust it for δ π to be as close as possible to zero, which is controlled by monitoring the linear cavity resonances while scanning the driving laser frequency.Finally, the third FPC, placed on the output port of the resonator (FPC out ) and followed by a polarization beam splitter (PBS), can be set so as to split the output field either in terms of the x and y modes, or in terms of the + and − hybrid modes.The intensity of the two components are then recorded with individual 12.5 GHz photodiodes.The total intensity is recorded separately.
To proceed, we record the output intensity levels across our driving pulses over subsequent resonator roundtrips using a real-time oscilloscope.A typical measurement is shown as color plots in Figs.3(a)-(c) in the form of a vertical concatenation (bottom-to-top) of a 20-roundtrip sequence of real-time oscilloscope traces plotted against time.Panels (a) and (b) correspond to the intensities of the + and − hybrid modes while (c) is the total intensity.The measurements have been obtained with 11 W of peak driving power (X = 50) and a roundtrip phase detuning locked to δ 0 = 1.16 rad (∆ = 10).Figures 3(a)-(b) reveal that under these conditions the hybrid mode intensity levels exhibit a broken symmetry (|ψ + | 2 = |ψ − | 2 ) as well as an anti-phase alternating dynamics: the hybrid modes appear to exchange their intensities one roundtrip to the next.This alternation is a key signature of the modal Möbius topology.At the same time, the total intensity [Fig.3(c)] appears constant and presents no sign of the periodic alternation of the hybrid modes.This points to a high level of symmetry between the two hybrid modes and confirms the realization of the ψ + ψ − exchange symmetry.We must note that this remarkable level of symmetry is achieved without any fine tuning of parameters and can persist for hours.
Further insights can be gained by measuring the high and low intensity levels of the hybrid modes for a range of resonator detunings.The results of such measurements are plotted in Fig. 3(d) for ∆ ranging from −4 to 10 and for the same driving power level as Figs.3(a)-(c).Blue and orange markers correspond to the + and − hybrid modes, respectively.For low values of the detuning, the two hybrid modes are found to have the same intensity and the intracavity field appears symmetric in terms of these modes.In contrast, beyond ∆ 2, we observe instead two mirror-like, asymmetric states, whose relative contrast increases with the detuning.The emergence of these asymmetric states is found to correspond to the start of the roundtrip-to-roundtrip alternating dynamics described above.These measurements agree with numerically-calculated steady-state solutions of Eqs. ( 13)-( 14) and which are shown as black curves in Fig. 3(d) (dashed lines correspond to unstable solutions).Overall, the data in Fig. 3(d) exhibit the behavior of a pitchfork bifurcation characteristics of SSB.This result is remarkable in that it is obtained without having to make any effort to realize an exchange symmetry.In fact, we observe this behavior as soon as the y-mode phase defect is sufficiently close to π (to within about the resonance linewidth).
To complete these observations, we have also made measurements in terms of the x and y modes.The emergence of the P2 SSB regime is characterized by the y-mode intensity rising above zero through parametric generation, but the intensities of the x and y modes do not otherwise reveal the alternating dynamics characteristics of the underlying Möbius topology.As described above, the associated dynamics for the ymode should be entirely contained into a roundtrip-toround phase step of π, which calls for a phase sensitive measurement.To that end, we have used homodyne detection by mixing the x and y output fields into a single-polarization 90 • optical hybrid (Kylia).A pair of 40 GHz-bandwidth balanced detectors then provide in phase, I ∝ Re(ψ x ψ * y ), and quadrature, Q ∝ Im(ψ x ψ * y ), components.Measurements of I and Q over 1000 resonator roundtrips in the P2 SSB regime are reported in Fig. 3(e), with cyan (respectively, purple) points corresponding to odd (even) roundtrips.The arrangement of these two groups of points in diametrically opposite positions in the complex plane reveal unequivocally that the relative phase between the x and y fields varies in step of π at each roundtrip.The average phase step is found to be (0.995 ± 0.020)π.Again this has been obtained without fine tuning of parameters and confirms that the system is nonlinearly attracted to a roundtrip phase jump of π as predicted in Section III.Note that, in comparison with the rest of the observations reported in Fig. 3, the data in Fig. 3(e) was obtained with a slightly shorter but otherwise identical resonator (resonator length of 10.5 m, roundtrip time t R = 50.6ns, and resonance linewidth ∆ω/(2π) = 825 kHz).

V. TEST OF THE EXCHANGE SYMMETRY
The results presented in the previous Section provide compelling evidence that a π-phase defect combined with the Kerr nonlinearity does confer a robust Möbius topology and exchange symmetry to a driven resonator, enabling SSB.In particular, the reproducibility of this behavior and our ability to observe it over long periods of time without any fine tuning of parameters strongly hint that the exchange symmetry is exact and protected by the nonlinearity as predicted by theory.In order to confirm this hypothesis further, we now present a more stringent and quantitative test of the quality of the exchange symmetry realized in our system.
Our test is based on the properties of SSB.Specifically, upon crossing the SSB bifurcation threshold in our resonator, an initially symmetric state will lose its stability in favor of one of two states with broken symmetry.Under perfectly symmetrical conditions, this selection should be truly random.Conversely, a deviation from symmetry will lead to a bias in the state selection process [32,34,63].The degree of randomness of the state selection process is therefore a measure of the quality of the exchange symmetry.To test for such randomness in practice, we have sinusoidally modulated the resonator detuning at 5.2 kHz around a mean (locked) value, so as to cross the SSB threshold repeatedly, and measured the state selection outcome.This measurement is based on the same resonator as that used to obtain the results of Fig. 3(e) but with a driving beam modulated so as to have 19 pulses simultaneously present in the resonator.Each pulse undergoes SSB independently, providing multiple simultaneous realizations of the experiment, thus enabling acquisition of a larger amount of data for bet-  ter statistical significance.The outcomes of the state selection process are determined by analyzing the output pulses captured with our real-time oscilloscope.
Over multiple runs, we have accumulated a sequence of 2.4 million individual events and expressed them as zero and one binary values.The randomness of that sequence was then rigorously assessed using the NIST Statistical Test Suite for random number generation (NIST STS-2.1.2) which consists of a series of statistical tests used to detect non-randomness within a given data set [64].In each test, a p-value is generated which gives the probability of this result occurring under the null hypothesis.We consider a p-value of 0.01 or lower to be evidence of non-randomness, or a failure.NIST recommends that the minimum pass rate for each statistical test should be 96 % (the whole sequence is partitioned into 100 individual sub-sequences for testing).The results are presented in Fig. 4 which shows that all tests are passed with a proportion well over 96 %, highlighting that we have no evidence of any bias in the SSB state selection process.This confirms the robustness of the exchange symmetry provided by the Möbius topology in our experiment.

VI. TOPOLOGICAL LOCALIZED STRUCTURES
In the experimental results discussed so far, we have only considered cases where the intracavity field is homogeneous across the nanosecond driving pulses.The spatial extension of the driving provides however the opportunity to support more complex, inhomogeneous field configurations.The P2 SSB regime naturally extends to such structures because the mechanism giving rise to the Möbius topology of our resonator and to the tworoundtrip alternating dynamics is purely local.We now discuss this aspect in more detail.
Inhomogeneous field structures can be accounted for by introducing a kinetic term into the Hamiltonian.This corresponds to adding terms of the form −i(φ 2 /2)∂ 2 ψ ± /∂τ 2 to, respectively, the mean-field equation for the + [Eq.( 13)] and − [Eq.( 14)] hybrid modes.Here τ is a fast time variable describing variations of the fields along the resonator.In the optical context relevant to our experiments, φ 2 represents chromatic dispersion integrated along the roundtrip of the resonator.With these kinetic terms, Eqs. ( 13)-( 14) take the form of coupled driven-dissipative Gross-Pitaevskii equations.Such equations have been used, e.g., to describe excitonpolariton condensates [37,65,66].In optics, these equations are known as Lugiato-Lefever equations (LLEs).The LLE was initially introduced as a paradigmatic model of pattern formation in dissipative optical systems [67].It is also known to support localized structures referred to as cavity solitons [68,69].In the temporal domain [70], these solitons underlie the new key technology of microresonator broadband optical frequency combs [71,72].In their coupled form, Eqs. ( 13)- (14) exhibit an even richer dynamics [61,66,73].Of particular relevance here are the recent experimental observations of two-component localized structures emerging through SSB, namely domain walls between homogeneous symmetry-broken states [30] and bright symmetrybroken vector cavity solitons [74,75].Given the formal equivalence of the equations describing these structures and our system, it should be clear that these structures have a P2 counterpart; in fact the only difference with what was previously reported should be the roundtripto-roundtrip alternation of the two components.
We start by considering P2 domain walls, which exist in presence of repulsive interactions (φ 2 > 0).Experimental evidence are presented in Figs.5(a)-(c).Panel (a) show intensity profiles of the + (blue) and − (orange) hybrid modes at a selected odd (left) and even (right) roundtrip, while panels (b) and (c) present the roundtrip-to-roundtrip evolution of those intensities over 20 roundtrips as pseudo-color plots similar to those used in Fig. 3.These plots reveal a field structure whose components are perfectly anti-correlated and, at the same time, alternate roundtrip-to-roundtrip.The driving pulses are subdivided into domains realizing, at any one time, one or the other of the two homogenous symmetrybroken states.By symmetry, the other component has the opposite domain structure.The domains are segregated by sharp kink-like temporal transitions, which correspond to dissipative (polarization) domain walls or PDWs [30,76,77].A close up of a domain wall obtained with a 70 GHz-bandwidth photodiode and sampling oscilloscope is shown in inset (this measurement is bandwidth limited with simulations predicting a rise/fall time of 4 ps).We must also point out that the total intensity (not shown) does not reveal any structure, neither the domain walls, nor the roundtrip-to-roundtrip alternation.These observations were performed with the same resonator and in the same conditions (X = 50, ∆ = 10) as that used for the results of Figs.3(a)-(c) (for which we have φ 2 +0.5 ps/THz).We have only added a shallow 10 GHz sinusoidal phase modulation imprinted onto the driving so as to align and trap the PDWs onto a controlled temporal reference grid [78].The domain walls were excited spontaneously by ramping up the detuning through the SSB bifurcation threshold.
Next we consider P2 bright symmetry-broken temporal cavity solitons.Because these require attractive interactions, i.e., a negative effective mass or, equivalently, anomalous chromatic dispersion (φ 2 < 0), we have used here a different resonator than for all the other observations reported so far in this Article.The resonator was 86 m-long and entirely made up of standard SMF28 (Corning) optical fiber with an input coupler that recirculates 95 % of the light (θ = 0.05).The roundtrip time t R = 420 ns, the resonance linewidth ∆ω/(2π) = 57 kHz, and φ 2 = −1.7 ps/THz.The driving and detection schemes were the same as before.Experimental results are presented in Figs.5(d)-(f) using the same layout as for the P2 domain walls.Here we can observe a se- quence of five bright solitons, which are aligned on a 2.87 GHz driving phase modulation grid [78].Spectral measurements (not shown) indicate a pulse duration of 1.6 ps (not resolved by the oscilloscope).These solitons are symmetry-broken with one hybrid mode component dominating over the other.Due to the high driving level (X = 30, corresponding to 2.5 W peak driving power, and with ∆ = 14.5), the contrast between the components is almost 100 %.Again, the symmetry implies the existence of two mirror-like states, i.e., solitons exhibiting a different dominant component, and our data confirms that these simultaneously coexist in the resonator.Finally, we observe a regular alternation of the components, with + mode-dominant solitons becoming − dominant and vice versa, each roundtrip to the next, as expected from the Möbius topology.Bright symmetry-broken cavity solitons are also known to go through a Hopf bifurcation for lower detuning values, leading to a breathing dynamics.This dynamics is characterized by slow in-phase oscillations of the intensities of the two components of the soliton [75].In Figs.5(g)-(k), we provide evidence of this dynamic in the P2 SSB regime.As for the other structures discussed above, panels (i) and (j) show color plots of the evolution of the + and − hybrid mode temporal intensity profiles.Additionally, panel (k) displays the cor-responding evolution for the total intensity.The latter clearly reveals the breathing dynamics, occurring here with a period of about 9 resonator roundtrips, while the roundtrip-to-roundtrip alternation characteristics of the P2 SSB regime is only detected in the individual components.The superposition of breathing and alternation is better seen in panel (h), with curves highlighting the evolution at the center of the soliton in terms of + (blue), − (orange), and total (green) intensities; these correspond to cross-sections through the center of panels (i)-(k), respectively.Note that we have selected here a lower driving level (X = 10) than in Figs.5(d)-(f), with the detuning set to ∆ = 5.2.This choice leads to less contrast between the soliton components [compare Figs.5(d) and (g)], and provides a clearer illustration of the overall complexity of the P2 symmetry-broken breather dynamics.
Finally, we must point out that all the structures reported in Fig. 5 could consistently be maintained and observed for long periods of time, from 30 to 60 minutes.This confirms that the robustness of the Möbius topology extends to complex inhomogeneous structures.

VII. DISCUSSION
We have described theoretically and experimentally a system characterized by a complex interplay between topology, nonlinearity, and SSB.Specifically, we have shown that a π phase-defect can force a driven Kerr resonator with two modes x and y towards an attractor arising through nonlinear degenerate parametric interactions, so that the system displays a modal Möbius topology.That topology is associated with an exchange symmetry, which together with the nonlinearity, enables SSB.Because of the Möbius topology, the two symmetrybroken states of our resonator alternate roundtrip-toroundtrip, which we refer to as P2 SSB.In the experiment, this alternation is directly seen in the intensity of the + and − hybrid modes.Homodyne measurements have also enabled us to detect the corresponding phase steps in the evolution of the x and y modes of the resonator.Because the Möbius topology arises through a purely local mechanism, it extends to inhomogeneous field distributions, giving rise to a range of P2 localized structures which we have also been able to observe.Thanks to the robustness associated with the topology, SSB is realized without any bias or fine tuning of parameters.This fact has been theoretically confirmed by an effective Floquet Hamiltonian model and, experimentally, through rigorous testing of the randomness of the symmetry-broken state selection statistics.The inherent robustness is also seen in the fact that all the different behaviors we have reported can be obtained without any sensitive adjustments of parameters, and are consistently observed over long period of times, easily exceeding 30 minutes (typically only limited by our detuning locking feedback loop).As the realization of SSB is often hampered by biases and non-idealities, we believe that our findings can have interest for a number of applications, in photonics, but more broadly in other bosonic systems.
As a final remark, we would like to mention that phasedefect-induced topological protection can be generalized to systems with more than two modes, potentially leading to P3, P4, . . .regimes of SSB.In what follows, we present brief arguments to illustrate the general idea, restricting our attention to three modes for simplicity.Consider a resonator with modes x, y, and z, where the y and z modes present phase defects of about, respectively, 2π/3 and 4π/3 with respect to the x mode.This could be implemented, e.g., using etched waveguide-based adiabatic mode converters or photonic crystals [48,79].The resonance frequencies of the three mode families would be shifted in steps of a third of the FSR.With these modes, it is possible to define hybrid modes α, β, γ, such that the hybrid mode amplitudes would be converted at each roundtrip α → β → γ → α, and display a generalized Möbius topology over a cycle of three resonator roundtrips.This will realize a cyclic Z 3 symmetry.The nonlinear four-wave interactions enabled by energy conservation are of the form x + x → y + z (and permu-tations).Using the same approach as in Section III, it is easy to show that these interactions make the 2π/3 and 4π/3 roundtrip phase steps robust attractors of the system, making the Z 3 symmetry exact, and protecting the topology.We have further confirmed that parameters exist where this system exhibits three symmetry-broken states, leading to P3 SSB, where these states cycle over three roundtrips.
As a particular application, we must note that the symmetry-broken states described in our Article could potentially be used to realize artificial spins.In the P2 SSB regime, these would correspond to spin-1/2 and could be exploited to implement an Ising machine.Such a machine aims to find the ground state of interacting spins described by the Ising Hamiltonian with an appropriate physical system.The interest rests on the fact that many complex combinatorial optimization problems relevant to modern society, such as drug discovery [80] or the analysis of social networks [81], can be mapped onto an Ising Hamiltonian [82,83].We must note that coherent optical Ising machines based on degenerate parametric oscillation have been particularly successful in recent years [35,36].A P2 SSB Ising machine would work along the same principle, but the inherent topological robustness of the Möbius topology could provide additional benefits.In this context, using more than two modes would offer even more tantalizing possibilities.For example, P3 symmetry-broken states would correspond to spin-1 particles, and open the way to the realization of a Potts machine [84].

FIG. 3 .
FIG. 3. Evidence of P2 SSB.Oscilloscope recordings of the intensity of the (a) + and (b) − hybrid modes and (c) total intensity over 20 subsequent resonator roundtrips (X = 50, ∆ = 10).The alternation of the hybrid mode intensities is a key signature of the underlying Möbius topology while the constant total intensity indicates a high level of symmetry between the hybrid modes.(d) High and low intensity levels of the + (blue) and − (orange) hybrid modes versus resonator detuning, revealing a pitchfork bifurcation (X = 50).The black lines are numerical predictions (dashed lines correspond to unstable states).(e) Complex plane showing the real (in-phase, I) and imaginary (quadrature, Q) components of ψxψ * y measured over 1000 resonator roundtrips in the P2 SSB regime.Cyan and purple density plots correspond, respectively, to odd and even roundtrips and their diametrically opposite placement signals that the relative phase between the x and y fields changes by π at every roundtrip.The spread of each density is magnified four times for clarity.

FIG. 4 .
FIG. 4. Test results for randomness of the P2 SSB state selection process using the NIST Statistical Test Suite.The bars indicate the percentage of sequences (total of 2.4 million events partitioned into 100 sub-sequences) that pass the tests with a significance level of 0.01.The recommended minimum pass rate of 96 % (black dotted line) is satisfied for all tests.The normalized driving power was set to X 12.

FastFIG. 5 .
FIG. 5. P2 SSB localized structures.(a)-(c), P2 domain walls observed for repulsive interactions (φ2 > 0) with X = 50 and ∆ = 10.(a) Temporal intensity profiles of the + (blue) and − (orange) hybrid modes at a selected odd (left) and even (right) roundtrip.The evolution of these profiles over 20 subsequent resonator roundtrips is illustrated in (b), (c), respectively, for the + and − hybrid modes.Inset: close up on one of the domain walls.(d)-(f), same as (a)-(c) but for P2 bright symmetry-broken cavity solitons observed for attractive interactions (φ2 < 0) with X = 30 and ∆ = 14.5.(g)-(k) Breathing dynamics of the P2 bright symmetry-broken cavity solitons observed with X = 10 and ∆ = 5.2.Same format as the other two examples with the addition of (k) an extra pseudo-color plot for the evolution of the total intensity and (h) a plot providing details of the evolution of the intensity at the center of the structure (blue/orange, +/− hybrid mode; green, total intensity).