Nonlinear topological symmetry protection in a dissipative system

Stéphane Coen, Bruno Garbin, Gang Xu, Liam Quinn, Nathan Goldman, Gian-Luca Oppo, Miro Erkintalo, Stuart Murdoch, Julien Fatome

► To cite this version:

Stéphane Coen, Bruno Garbin, Gang Xu, Liam Quinn, Nathan Goldman, et al. Nonlinear topological symmetry protection in a dissipative system. Nature Communications, 2024, 15 (1), pp.1398. 10.1038/s41467-023-44640-x. hal-04076324

HAL Id: hal-04076324

https://hal.science/hal-04076324v1

Submitted on 7 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Nonlinear topological symmetry protection in a dissipative system

Stéphane Coen,<sup>1,2,*</sup> Bruno Garbin,<sup>1,2,†</sup> Gang Xu,<sup>1,2,‡</sup> Liam Quinn,<sup>1,2</sup> Nathan Goldman,<sup>3</sup>
Gian-Luca Oppo,<sup>4</sup> Miro Erkintalo,<sup>1,2</sup> Stuart G. Murdoch,<sup>1,2</sup> and Julien Fatome<sup>1,2,5</sup>

<sup>1</sup>*Physics Department, The University of Auckland,*
*Private Bag 92019, Auckland 1142, New Zealand*

<sup>2</sup>*The Dodd-Walls Centre for Photonic and Quantum Technologies, Dunedin, New Zealand*

<sup>3</sup>*Center for Nonlinear Phenomena and Complex Systems,*
*Université Libre de Bruxelles, CP 231, B-1050 Brussels, Belgium*

<sup>4</sup>*SUPA and Department of Physics, University of Strathclyde, Glasgow G4 0NG, Scotland, European Union*

<sup>5</sup>*Laboratoire Interdisciplinaire Carnot de Bourgogne (ICB), UMR 6303 CNRS,*
*Université de Bourgogne Franche-Comté, 9 Avenue Alain Savary, BP 47870, F-21078 Dijon, France*

We report an experimental and theoretical investigation of a system whose dynamics is dominated by an intricate interplay between three key concepts of modern physics: topology, nonlinearity, and spontaneous symmetry breaking. The experiment is based on a two-mode coherently-driven optical resonator in which photons interact through the Kerr nonlinearity. In presence of a phase defect between the modes, a nonlinear attractor develops, which confers a synthetic Möbius topology to the modal structure of the system. That topology is associated with an inherently protected exchange symmetry between the modes, enabling the realization of spontaneous symmetry breaking in ideal, bias-free, conditions without any fine tuning of parameters. The dynamic manifests itself by a periodic alternation of the modes from one resonator roundtrip to the next reminiscent of period-doubling. This extends to a range of localized structures in the form of domain walls, bright solitons, and breathers, which have all been observed with remarkable long term stability. A rigorous testing of the randomness of the symmetry-broken state selection statistics has also confirmed the robustness of the exchange symmetry in our experiment. Our results and conclusions are supported by an effective Hamiltonian model explaining the symmetry protection in our system. The model also shows that our work has relevance to other systems of interacting bosons and to the Floquet engineering of quantum matter. Our work 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–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 solid-state 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 wide-ranging 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, non-exhaustively, 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–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

<footnote>* Corresponding author, s.coen@auckland.ac.nz</footnote>

<footnote>† Present address: Université Paris-Saclay, CNRS, Centre de Nanosciences et de Nanotechnologies, 91120, Palaiseau, France</footnote>

<footnote>‡ Present address: School of Optical and Electronic Information, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan, China</footnote>
<figure><img src="image_2.png" /><figcaption>FIG. 1. (a) Schematic representation of a two-mode coherently driven Kerr resonator with a π-phase defect. Inset on the right shows the linear resonances of the two modes and how only one mode is driven (here the x mode). (b) Equivalent Möbius-like resonator described in terms of the + and − hybrid modes defined by Eqs. (1). (c) Parametric oscillation of the y-mode, $|\psi_y|^2 > 0$ (top left) corresponds to an SSB bifurcation for the occupations of the hybrid modes, $|\psi_+|^2 \neq |\psi_-|^2$ (bottom left). In that regime (P2 SSB), the phase of the y-mode changes by $\pi$ at each resonator roundtrip (top right), corresponding to a periodic roundtrip-to-roundtrip swapping of the hybrid mode occupations (bottom right). $t_R$ is the roundtrip time of the resonator.</figcaption></figure>

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–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 $\psi_x$ and $\psi_y$ denote the modes’ classical electric field (complex) amplitudes inside the resonator scaled such that $|\psi_x|^2$ and $|\psi_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 $\pi$-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 $\pi$. 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 $\psi_+$ and $\psi_-$ defined by the unitary transformation

$$ \psi_+ = \frac{1}{\sqrt{2}} (\psi_x + i\psi_y), \quad \psi_- = \frac{1}{\sqrt{2}} (\psi_x - i\psi_y). \quad (1) $$

Because of the $\pi$-phase defect, the amplitude of the $y$ mode flips its sign with respect to the driven mode at each roundtrip in the resonator, $\psi_y \to \psi_y e^{i\pi} = -\psi_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, $\psi_+ \Leftrightarrow \psi_-$ [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–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, $\psi_+ \Leftrightarrow \psi_-$. 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, $\psi_+ = \psi_-$, is necessarily associated with a $y$ field of zero amplitude. By contrast, a symmetry-broken hybrid state, $\psi_+ \neq \psi_-$, requires a non-zero $y$ field. By symmetry, the mirror state, $\psi_+ \Leftrightarrow \psi_-$, must also exist, corresponding to a $y$ field with a phase differing by $\pi$, $\psi_y \Leftrightarrow -\psi_y$. Keeping in mind that the $y$-field 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 $|\psi_+|^2, |\psi_-|^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 $\pi$ 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–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 $\pi$. 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 $\hat{H}_{\text{NL}}$ describing the Kerr-mediated interactions taking place along the resonator roundtrip,

$$ \frac{1}{\hbar g} \hat{H}_{\text{NL}} = \frac{A}{2} (\hat{a}_x^\dagger \hat{a}_x^\dagger \hat{a}_x \hat{a}_x + \hat{a}_y^\dagger \hat{a}_y^\dagger \hat{a}_y \hat{a}_y) + B (\hat{a}_x^\dagger \hat{a}_y^\dagger \hat{a}_x \hat{a}_y) + \frac{C}{2} (\hat{a}_x^\dagger \hat{a}_x^\dagger \hat{a}_y \hat{a}_y + \hat{a}_y^\dagger \hat{a}_y^\dagger \hat{a}_x \hat{a}_x). \quad (2) $$

Here $g$ is the single photon-induced Kerr frequency shift of the resonator ($g > 0$) [53] while $\hat{a}_\sigma^\dagger$ (resp. $\hat{a}_\sigma$) represent creation (annihilation) operators for the mode $\sigma = x, y$. These operators satisfy bosonic commutation relations, $[\hat{a}_\sigma, \hat{a}_{\sigma'}^\dagger] = \delta_{\sigma,\sigma'}$. 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 the Methods.

To understand specifically the evolution of the phase of the $y$ mode, we examine the equation of motion for $\hat{a}_y$, i.e., $d\hat{a}_y/dt = i[\hat{H}_{\text{NL}}, \hat{a}_y]/\hbar$. In the classical limit, $\hat{a}_y^\dagger \to \psi_y$, we get

$$ \frac{d\psi_y}{dt} = ig \left[ (A|\psi_y|^2 + B|\psi_x|^2) \psi_y + C \psi_y^* \psi_x^2 \right], \quad (3) $$

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, $\psi_\sigma = |\psi_\sigma| e^{i\phi_\sigma}$, we find that the phase $\phi_y$ of the $y$ mode obeys

$$ \frac{d\phi_y}{dt} = g [A|\psi_y|^2 + B|\psi_x|^2 + C|\psi_x|^2 \cos(2\phi_y - 2\phi_x)]. \quad (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 $\ll 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 Methods for more detail), leading to

$$ \phi_y^{(m+1)} - \phi_y^{(m)} = gt_R (A|\psi_y|^2 + B|\psi_x|^2) + gt_R C |\psi_x|^2 \cos (2\phi_y^{(m)} - 2\phi_x^{(m)}) \underbrace{-\delta_0 + \pi + \delta_\pi}_{-\delta_y} . \quad (5) $$

Here m is the roundtrip index, $\delta_0 = \delta_x$ is the roundtrip phase detuning of the driving field (modulo $2\pi$), also assumed $\ll 1$ ($\delta_0 > 0$ means red-detuned), and the corresponding detuning for the y mode is defined through $\delta_x - \delta_y = \pi + \delta_\pi$, with $\delta_\pi$ representing deviations from an exact value of $\pi$. We can observe that all terms above are in fact $\ll \pi$, except $\pi$ itself (a similar argument for the x mode would show that $\phi_x^{(m+1)} - \phi_x^{(m)} \ll 1$). Moreover, the pair hopping (parametric) contribution ($\propto C$), which is the only phase sensitive term, is $\pi$ periodic in $\phi_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 $\pi$, one roundtrip to the next. That solution acts as an attractor for the system, irrespective of the value of $\delta_\pi$. Effectively, the intracavity fields will self-adjust through the dissipative nonlinear dynamics to guarantee a roundtrip-to-roundtrip y mode phase step of $\pi$. 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 $\psi_+ \Leftrightarrow \psi_-$. 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 $\hat{a}_\pm^\dagger (\hat{a}_\pm)$; these are related to the corresponding x, y operators by the same unitary transformations as Eqs. (1). We find

$$ \frac{1}{\hbar g} \hat{H}_{\text{NL}} = \frac{A'}{2} (\hat{a}_+^\dagger \hat{a}_+^\dagger \hat{a}_+ \hat{a}_+ + \hat{a}_-^\dagger \hat{a}_-^\dagger \hat{a}_- \hat{a}_-) + B' (\hat{a}_+^\dagger \hat{a}_-^\dagger \hat{a}_+ \hat{a}_-) \quad (6) $$

where $A' = 2/3$ and $B' = 4/3$. For convenience, let us define a set of angular momentum (Schwinger) operators as

$$ \begin{aligned} \hat{J}_x &= \frac{\hbar}{2} (a_+^\dagger a_- + a_-^\dagger a_+), & \hat{J}_y &= \frac{\hbar}{2i} (a_-^\dagger a_+ - a_+^\dagger a_-), \\ \hat{J}_z &= \frac{\hbar}{2} (a_-^\dagger a_- - a_+^\dagger a_+), & \hat{N} &= a_+^\dagger a_+ + a_-^\dagger a_- \end{aligned} \quad (7) $$

These operators satisfy commutation relations $[\hat{J}_\mu, \hat{J}_\nu] = i\hbar\epsilon_{\mu\nu\lambda}\hat{J}_\lambda$ similar to spins ($\epsilon_{\mu\nu\lambda}$ is the Levi-Civita symbol), while $\hat{N}$ counts the total number of photons. In

this part of our analysis, we assume a lossless resonator (conservative limit) so $\hat{N}$ is constant. With these notations, a difference in angular frequency detuning $\delta\Omega$ between the two hybrid modes enter the Hamiltonian as $\delta\Omega \hat{J}_z$, so that the total Hamiltonian of the system reads

$$ \begin{aligned} \hat{H}_{\text{asym}} &= \hat{H}_{\text{NL}} + \delta\Omega \hat{J}_z \\ &= \hbar g \left[ a \hat{J}_z^2 / \hbar^2 + b \hat{N}^2 - c \hat{N} \right] + \delta\Omega \hat{J}_z \end{aligned} \quad (8) $$

where $a = A' - B'$, $b = (A' + B')/4$, and $c = A'/2$.

We now proceed by seeking an effective (Floquet) Hamiltonian $\hat{H}_{\text{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 $\hat{U}_{\text{swap}} = \hat{\sigma}_x = e^{i(\pi/2)(\hat{\sigma}_x - 1)}$ (with $\hat{\sigma}_x$ a Pauli matrix), leading to

$$ e^{-i \frac{1}{\hbar} 2t_{\text{R}} \hat{H}_{\text{eff}}} = \underbrace{\hat{U}_{\text{swap}} e^{-i \frac{1}{\hbar} t_{\text{R}} \hat{H}_{\text{asym}}} \hat{U}_{\text{swap}}}_{e^{-i \frac{1}{\hbar} t_{\text{R}} \hat{H}_1}} e^{-i \frac{1}{\hbar} t_{\text{R}} \hat{H}_{\text{asym}}}. \quad (9) $$

The first three factors on the right hand side of this expression define $\hat{H}_1$ which can be obtained exactly as [58]

$$ \hat{H}_1 = e^{i \frac{\pi}{2} \hat{\sigma}_x} \hat{H}_{\text{asym}} e^{-i \frac{\pi}{2} \hat{\sigma}_x}. \quad (10) $$

Noting that the swapping operation is a single-particle process, we can substitute $\hat{J}_x/\hbar$ for $\hat{\sigma}_x/2$ in the exponentials above and then use the Baker-Campbell-Hausdorff formula. Specifically we have $e^{i\pi\hat{J}_x/\hbar} \hat{J}_z e^{-i\pi\hat{J}_x/\hbar} = -\hat{J}_z$ while all other terms of $\hat{H}_{\text{asym}}$ are left unchanged, leading to

$$ \hat{H}_1 = \hat{H}_{\text{NL}} - \delta\Omega \hat{J}_z. \quad (11) $$

Comparing the expressions for $\hat{H}_{\text{asym}}$, Eq. (8), and $\hat{H}_1$, above, we observe that the sign of the asymmetric term is reversed. It is also clear that $\hat{H}_{\text{asym}}$ commutes with $\hat{H}_1$, which from Eq. (9) leads to $\hat{H}_{\text{eff}} = (\hat{H}_1 + \hat{H}_{\text{asym}})/2$. It follows that the asymmetric term cancels out exactly over two resonator roundtrips, and we have $\hat{H}_{\text{eff}} = \hat{H}_{\text{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, $\psi_+ \Leftrightarrow \psi_-$, 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 $\psi_+$, $\psi_-$ based on the effective Hamiltonian $\hat{H}_{\text{eff}}$. We have already established that the dissipation guarantees an exact overall roundtrip phase jump of $\pi$ for the y mode, so we do not
need to consider the influence of $\delta_\pi$, the deviation from an exact $\pi$ phase defect, on the swapping operator $\hat{U}_{swap}$. Still, $\delta_\pi$ contributes to the linear detunings through the relation $\delta_0 - \delta_y = \pi + \delta_\pi$. This effect can be taken into account by the following roundtrip phase detuning operator,

$$ \hbar [\delta_0(\hat{a}_x^\dagger \hat{a}_x) + (\delta_0 - \delta_\pi)(\hat{a}_y^\dagger \hat{a}_y)] = \hbar \left(\delta_0 - \frac{\delta_\pi}{2}\right) \hat{N} + \delta_\pi \hat{J}_x. \quad (12) $$

Note that, for the $y$ mode, we have discarded the $\pi$ phase-shift that is already incorporated in $\hat{H}_{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,

$$ t_R \frac{d\psi_+}{dt} = \left[ -\frac{t_R \Delta\omega}{2} + i g t_R \left( A' |\psi_+|^2 + B' |\psi_-|^2 \right) - i \left( \delta_0 - \frac{\delta_\pi}{2} \right) \right] \psi_+ - i \frac{\delta_\pi}{2} \psi_- + \sqrt{\theta} S, \quad (13) $$

$$ t_R \frac{d\psi_-}{dt} = \left[ -\frac{t_R \Delta\omega}{2} + i g t_R \left( A' |\psi_-|^2 + B' |\psi_+|^2 \right) - i \left( \delta_0 - \frac{\delta_\pi}{2} \right) \right] \psi_- - i \frac{\delta_\pi}{2} \psi_+ + \sqrt{\theta} S. \quad (14) $$

To these equations, we have added phenomenologically losses (dissipation), associated with the linewidth $\Delta\omega$ of the resonances, and driving, with strength $S$ and external coupling fraction $\theta$. 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 the Methods 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/\sqrt{2}$, where $S_x$ is the driving strength of the $x$ mode only. Additionally, we can observe that a deviation $\delta_\pi$ from an exact $\pi$ 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 $\pi$ phase defect, this condition is automatically and robustly satisfied. Finally, while the equations above describe the stroboscopic mean-field evolution of $\psi_+$ and $\psi_-$ 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 topologically-protected 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 $\pi$ 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 $\Delta = \delta_0/\alpha$ and the normalized driving power $X = gt_R\theta S_x^2/\alpha^3$, where $\alpha = t_R\Delta\omega/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 ($\theta = 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

<figure><img src="image_3.png" /><figcaption>FIG. 2. Schematic of the experimental setup. FPC: Optical fiber polarization controller. PBS: Polarizing beam-splitter.</figcaption></figure>
<figure><img src="image_4.png" /></figure>

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, \Delta = 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 $\psi_x, \psi_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 $\pi$ at every roundtrip. The spread of each density is magnified four times for clarity.

scalar parametric instabilities [55]. Overall, the resonator is 12 m-long, corresponding to a roundtrip time $t_R$ of 57 ns (FSR = 17.54 MHz), and has a resonance linewidth $\Delta\omega/(2\pi) = 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<sub>in</sub>) 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<sub>π</sub>), defines the $\pi$ phase-shift defect. We adjust it for $\delta_\pi$ 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<sub>out</sub>) 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 $\delta_0 = 1.16$ rad ($\Delta = 10$). Figures 3(a)–(b) reveal that under these conditions the hybrid mode intensity levels exhibit a broken symmetry ($|\psi_+|^2 \neq |\psi_-|^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 $\psi_+ \leftrightarrow \psi_-$ 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 $\Delta$ 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 $\Delta \approx 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 $\pi$ (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 *y*-mode should be entirely contained into a roundtrip-to-round phase step of $\pi$, 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 \propto \text{Re}(\psi_x \psi_y^*)$, and quadrature, $Q \propto \text{Im}(\psi_x \psi_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 $\pi$ at each roundtrip. The average phase step is found to be $(0.995 \pm 0.020)\pi$. Again this has been obtained without fine tuning of parameters and confirms that the system is nonlinearly attracted to a roundtrip phase jump of $\pi$ 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.6$ ns, and resonance linewidth $\Delta\omega/(2\pi) = 825$ kHz).

V. TEST OF THE EXCHANGE SYMMETRY

The results presented in the previous Section provide compelling evidence that a $\pi$-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-

<figure><img src="image_5.png" /><figcaption>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 \approx 12$.</figcaption></figure>
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 two-roundtrip 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 exciton-polariton 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 symmetry-broken 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 roundtrip-to-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) shows 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 symmetry-broken 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 \simeq +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-
<figure><img src="image_6.png" /></figure>

FIG. 5. P2 SSB localized structures. (a)–(c), P2 domain walls observed for repulsive interactions ($\phi_2 > 0$) with $X = 50$ and $\Delta = 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 ($\phi_2 < 0$) with $X = 30$ and $\Delta = 14.5$. (g)–(k) Breathing dynamics of the P2 bright symmetry-broken cavity solitons observed with $X = 10$ and $\Delta = 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).

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 $\Delta = 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 $\Delta = 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 symmetry-broken states of our resonator alternate roundtrip-to-roundtrip, 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 phase-defect-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\pi/3$ and $4\pi/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 $\alpha, \beta, \gamma$, such that the hybrid mode amplitudes would be converted at each roundtrip $\alpha \to \beta \to \gamma \to \alpha$, and display a generalized Möbius topology over a cycle of three resonator roundtrips. This will realize a cyclic $\mathbb{Z}_3$ symmetry. The nonlinear four-wave interactions enabled by energy conservation are of the form $x + x \to y + z$ (and permu-

tations). Using the same approach as in Section III, it is easy to show that these interactions make the $2\pi/3$ and $4\pi/3$ roundtrip phase steps robust attractors of the system, making the $\mathbb{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].

METHODS

**Classical derivation:** We present here a classical derivation of the mean-field equations (13)–(14). This derivation complements the Hamiltonian description of Section III, which was derived in the conservative limit, by including dissipation and driving explicitly. For completeness, we also consider modes with arbitrary polarization states, i.e., arbitrary values of the A, B, and C coefficients [see Eq. (2)]. We start from the classical equations of motion for the two modes, i.e., Eq. (3) and a similar equation for $\psi_x$. These equations are complemented by boundary conditions,

$$ \begin{pmatrix} \psi_x \\ \psi_y \end{pmatrix}^{(m+1)} = e^{-\alpha} M \begin{pmatrix} \psi_x \\ \psi_y \end{pmatrix}^{(m)} + \sqrt{\theta} S_{\text{in}} \begin{pmatrix} \cos \chi \\ \sin \chi \end{pmatrix}. \quad (A.1) $$

Here, $\alpha = t_R \Delta\omega/2$ is the dissipation rate of the resonator, as introduced before, $S_{\text{in}}$ is the total driving amplitude, with the angle $\chi$ describing a split of the driving between the two modes ($\chi = 0$ corresponds to only driving the x mode as previously considered, and $S_x = S_{\text{in}} \cos \chi$), and the matrix M represents the action of the detunings,

$$ M = \begin{pmatrix} e^{-i\delta_0} & 0 \\ 0 & e^{i\pi} e^{-i(\delta_0 - \delta_\pi)} \end{pmatrix}. \quad (A.2) $$
We proceed by restricting ourselves to high-Q resonators operated close to resonance and by following an approach similar to that of Refs. [85, 86]. In these conditions, $\alpha$, $\theta$, $\delta$, and $\delta_{\pi}$ can be assumed $\ll 1$. Similarly, all the nonlinear terms in the equations of motion are small

and these equations can be integrated over one roundtrip time $t_R$ by a single Euler step, as we have done with Eq. (4). Combining the result with the boundary conditions, Eq. (A.1), and expanding all small terms at first order, we can express the fields at roundtrip $m + 1$ in terms of those at roundtrip $m$ as

$$ \psi_x^{(m+1)} = [1 - \alpha - i\delta_0 + igt_R (A|\psi_x^{(m)}|^2 + B|\psi_y^{(m)}|^2)] \psi_x^{(m)} + igt_R C \psi_x^{(m)*} \psi_y^{(m)} + \sqrt{\theta} S_{\text{in}} \cos \chi, \quad (A.3) $$

$$ \psi_y^{(m+1)} = -[1 - \alpha - i(\delta_0 - \delta_\pi) + i gt_R (A|\psi_y^{(m)}|^2 + B|\psi_x^{(m)}|^2)] \psi_y^{(m)} - i gt_R C \psi_y^{(m)*} \psi_x^{(m)} + \sqrt{\theta} S_{\text{in}} \sin \chi. \quad (A.4) $$

Above, we can observe that the leading contribution for the $y$ mode is a change of sign, as expected from the $\pi$ phase defect.

The next step is to iterate Eqs. (A.3)–(A.4) and seek expressions for $\psi_{x,y}^{(m+2)}$ in terms of the $m$th roundtrip fields. At first order, we can assume that, in the nonlinear terms that result, $\psi_x^{(m+1)} \approx \psi_x^{(m)}$ and $\psi_y^{(m+1)} \approx -\psi_y^{(m)}$. We can also observe that the $y$ component of the driving will cancel out. Because of the double action of the $\pi$ phase defect over two roundtrips, we then find that the changes in the fields over two roundtrips, $\psi_{x,y}^{(m+2)} - \psi_{x,y}^{(m)}$, are first order quantities. This enables us to introduce a slow-time derivative over two resonator roundtrips, defined as $d./dt = [(m+2) - (m)]/(2t_R)$, which transforms the difference equations we get into the following differential equations (this latter step is analogous to the stochastic Floquet approach used in Section III),

$$ t_R \frac{d\psi_x}{dt} = [-\alpha - i\delta_0 + i gt_R (A|\psi_x|^2 + B|\psi_y|^2)] \psi_x + i gt_R C \psi_x^{*} \psi_y^2 + \sqrt{\theta} S_{\text{in}} \cos \chi, \quad (A.5) $$

$$ t_R \frac{d\psi_y}{dt} = [-\alpha - i(\delta_0 - \delta_\pi) + i gt_R (A|\psi_y|^2 + B|\psi_x|^2)] \psi_y + i gt_R C \psi_y^{*} \psi_x^2. \quad (A.6) $$

where we have dropped the roundtrip index $m$.

The last step is to express the equations above in terms of the + and – hybrid mode amplitudes defined by Eqs. (1). This leads to Eqs. (13)–(14), with

$$ A' = \frac{1 + B - C}{2} A, \quad B' = (1 + C)A, \quad (A.7) $$

but with additional parametric terms, $i gt_R C' \psi_+^* \psi_-^2$ for Eq. (13) and $i gt_R C' \psi_-^* \psi_+^2$ for Eq. (14), where

$$ C' = \frac{1 - B - C}{2} A. \quad (A.8) $$

We can observe that these terms do not break the $\psi_+ \leftrightarrow \psi_-$ exchange symmetry of the equations, and overall do not change the phenomenology. When the $x$ and $y$ modes

correspond to linear polarization states, $C' = 0$ as in Eq. (6).

Note that the above analysis assumes a $\pi$ phase defect that is purely localized. Additional considerations show that distributing the phase defect over a certain fraction of the resonator length amounts to a phase-mismatch for the parametric interaction (non conservation of momentum). This leads to less efficient generation of the $y$ mode but P2 SSB remains possible. When parametric generation cannot overcome the losses of the resonator anymore, the P2 SSB regime eventually disappears.

Kinetic terms (chromatic dispersion) can also be included in our derivation in a straightforward manner.

## DATA AVAILABILITY

The data that support the plots within this paper and other findings of this study are available from the corresponding author upon reasonable request.

## ACKNOWLEDGMENTS

We acknowledge financial support from The Royal Society of New Zealand, in the form of Marsden Funding (18-UOA-310) as well as M.E.’s Rutherford Discovery Fellowship (RDF-15-UOA-015). J.F. thanks the Conseil régional de Bourgogne Franche-Comté, mobility (2019-7-10614), financial support from the CNRS, the IRP WallIN project, and the FEDER, Solstice project. N.G. is supported by the FRS-FNRS (Belgium), the EOS program (CHEQS project), and the ERC (TopoCold and LATIS projects).

## AUTHOR CONTRIBUTIONS

BG and JF set up the experiment and made the initial observations of P2 SSB. JF carried on subsequent numerical modeling and made the measurements shown in Fig. 3 as well as the P2 domain walls [Fig. 5(a)–(c)].
GX performed the experiment related to the P2 bright symmetry-broken cavity solitons and associated breather [Fig. 5(d)–(k)]. LQ experimentally implemented and carried on the randomness statistical test of Fig. 4. NG derived the Hamiltonian model discussed in Section III. GLO, ME and SM provided extra theoretical and/or experimental support. SC suggested the topological interpretation of P2 SSB and derived the classical mean-field

description. All authors helped in interpreting the results. SC wrote the paper with feedback from all the authors.

## COMPETING INTERESTS

The authors declare no competing interests.

[1] M. Z. Hasan and C. L. Kane, Colloquium: Topological insulators, Rev. Mod. Phys. **82**, 3045 (2010).

[2] L. Lu, J. D. Joannopoulos, and M. Soljačić, Topological photonics, Nat. Photon. **8**, 821 (2014).

[3] G. Nicolis and I. Prigogine, Exploring Complexity: An Introduction (St. Martin's Press, New York, 1989).

[4] N. N. Akhmediev and A. Ankiewicz, Solitons — Nonlinear Pulses and Beams, 1st ed. (Chapman & Hall, London, 1997).

[5] P. W. Anderson, More is different, Science **177**, 393 (1972).

[6] F. Strocchi, Symmetry Breaking (Springer, Berlin New York, 2008).

[7] Y. Nambu, Nobel Lecture: Spontaneous symmetry breaking in particle physics: A case of cross fertilization, Rev. Mod. Phys. **81**, 1015 (2009).

[8] L. D. Landau and E. M. Lifshitz, Statistical Physics (Elsevier, 1980).

[9] E. Korotkevich, R. Niwayama, A. Courtois, S. Friese, N. Berger, F. Buchholz, and T. Hiiragi, The apical domain is required and sufficient for the first lineage segregation in the mouse embryo, Dev. Cell **40**, 235 (2017).

[10] X.-L. Qi and S.-C. Zhang, Topological insulators and superconductors, Rev. Mod. Phys. **83**, 1057 (2011).

[11] N. R. Cooper, J. Dalibard, and I. B. Spielman, Topological bands for ultracold atoms, Rev. Mod. Phys. **91**, 015005 (2019).

[12] T. Ozawa, H. M. Price, A. Amo, N. Goldman, M. Hafezi, L. Lu, M. C. Rechtsman, D. Schuster, J. Simon, O. Zilberberg, and I. Carusotto, Topological photonics, Rev. Mod. Phys. **91**, 015006 (2019).

[13] M. Segev and M. A. Bandres, Topological photonics: Where do we go from here?, Nanophot. **10**, 425 (2021).

[14] S. D. Huber, Topological mechanics, Nat. Phys. **12**, 621 (2016).

[15] T. Ozawa and H. M. Price, Topological quantum matter in synthetic dimensions, Nat. Rev. Phys. **1**, 349 (2019).

[16] A. Y. Kitaev, Fault-tolerant quantum computation by anyons, Ann. Phys. **303**, 2 (2003).

[17] R. Raussendorf and J. Harrington, Fault-tolerant quantum computation with high threshold in two dimensions, Phys. Rev. Lett. **98**, 190504 (2007).

[18] M. A. Bandres, S. Wittek, G. Harari, M. Parto, J. Ren, M. Segev, D. N. Christodoulides, and M. Khajavikhan, Topological insulator laser: Experiments, Science **359**, eaar4005 (2018).

[19] D. Smirnova, D. Leykam, Y. Chong, and Y. Kivshar, Nonlinear topological photonics, Appl. Phys. Rev. **7**, 021306 (2020).

[20] L. J. Maczewsky, M. Heinrich, M. Kremer, S. K. Ivanov, M. Ehrhardt, F. Martinez, Y. V. Kartashov, V. V. Kono-

top, L. Torner, D. Bauer, and A. Szameit, Nonlinearity-induced photonic topological insulator, Science **370**, 701 (2020).

[21] S. Xia, D. Kaltsas, D. Song, I. Komis, J. Xu, A. Szameit, H. Buljan, K. G. Makris, and Z. Chen, Nonlinear tuning of PT symmetry and non-Hermitian topological states, Science **372**, 72 (2021).

[22] S. de Léséleuc, V. Lienhard, P. Scholl, D. Barredo, S. Weber, N. Lang, H. P. Büchler, T. Lahaye, and A. Browaeys, Observation of a symmetry-protected topological phase of interacting bosons with Rydberg atoms, Science **365**, 775 (2019).

[23] M. Jürgensen, S. Mukherjee, and M. C. Rechtsman, Quantized nonlinear Thouless pumping, Nature (London) **596**, 63 (2021).

[24] P. Hamel, S. Haddadi, F. Raineri, P. Monnier, G. Beaudoin, I. Sagnes, A. Levenson, and A. M. Yacomotti, Spontaneous mirror-symmetry breaking in coupled photonic-crystal nanolasers, Nat. Photon. **9**, 311 (2015).

[25] Q.-T. Cao, H. Wang, C.-H. Dong, H. Jing, R.-S. Liu, X. Chen, L. Ge, Q. Gong, and Y.-F. Xiao, Experimental demonstration of spontaneous chirality in a nonlinear microresonator, Phys. Rev. Lett. **118**, 033901 (2017).

[26] L. Del Bino, J. M. Silver, S. L. Stebbings, and P. Del'Haye, Symmetry breaking of counter-propagating light in a nonlinear resonator, Sci. Rep. **7**, 43142 (2017).

[27] Q.-T. Cao, R. Liu, H. Wang, Y.-K. Lu, C.-W. Qiu, S. Rotter, Q. Gong, and Y.-F. Xiao, Reconfigurable symmetry-broken laser in a symmetric microcavity, Nat. Commun. **11**, 1136 (2020).

[28] A. E. Kaplan and P. Meystre, Enhancement of the Sagnac effect due to nonlinearly induced nonreciprocity, Opt. Lett. **6**, 590 (1981).

[29] M. T. M. Woodley, J. M. Silver, L. Hill, F. Copie, L. Del Bino, S. Zhang, G.-L. Oppo, and P. Del'Haye, Universal symmetry-breaking dynamics for the Kerr interaction of counterpropagating light in dielectric ring resonators, Phys. Rev. A **98**, 053863 (2018).

[30] B. Garbin, J. Fatome, G.-L. Oppo, M. Erkintalo, S. G. Murdoch, and S. Coen, Dissipative polarization domain walls in a passive coherently driven Kerr resonator, Phys. Rev. Lett. **126**, 023904/1 (2021).

[31] T. B. Benjamin, Bifurcation phenomena in steady flows of a viscous fluid. II. Experiments, Proc. R. Soc. Lond. A **359**, 27 (1978).

[32] B. Garbin, J. Fatome, G.-L. Oppo, M. Erkintalo, S. G. Murdoch, and S. Coen, Asymmetric balance in symmetry breaking, Phys. Rev. Res. **2**, 023244/1 (2020).

[33] B. Garbin, A. Giraldo, K. J. H. Peters, N. G. R. Broderick, A. Spakman, F. Raineri, A. Levenson, S. R. K. Rodriguez, B. Krauskopf, and A. M. Yacomotti, Spon-
taneous symmetry breaking in a coherently driven
nanophotonic Bose-Hubbard dimer, Phys. Rev. Lett.
**128**, 053901 (2022).

[34] H. Crauel and F. Flandoli, Additive noise destroys a pitchfork bifurcation, J. Dyn. Differ. Equ. **10**, 259 (1998).

[35] T. Inagaki, K. Inaba, R. Hamerly, K. Inoue, Y. Yamamoto, and H. Takesue, Large-scale Ising spin network based on degenerate optical parametric oscillators, Nat. Photon. **10**, 415 (2016).

[36] T. Honjo, T. Sonobe, K. Inaba, T. Inagaki, T. Ikuta, Y. Yamada, T. Kazama, K. Enbutsu, T. Umeki, R. Kasahara, K.-i. Kawarabayashi, and H. Takesue, 100,000-spin coherent Ising machine, Sci. Adv. **7**, eabh0952 (2021).

[37] I. Carusotto and C. Ciuti, Quantum fluids of light, Rev. Mod. Phys. **85**, 299 (2013).

[38] N. Goldman and J. Dalibard, Periodically driven quantum systems: Effective Hamiltonians and engineered gauge fields, Phys. Rev. X **4**, 031027 (2014).

[39] A. Eckardt, Colloquium: Atomic quantum gases in periodically driven optical lattices, Rev. Mod. Phys. **89**, 011004 (2017).

[40] N. Goldman, Floquet engineering of optical nonlinearities: A quantum many-body approach (2022), arXiv:2203.05554 [cond-mat, physics:physics, physics:quant-ph].

[41] J. A. Sheikh, J. Dobaczewski, P. Ring, L. M. Robledo, and C. Yannouleas, Symmetry restoration in mean-field approaches, J. Phys. G: Nucl. Part. Phys. **48**, 123001 (2021).

[42] A. E. Kaplan and P. Meystre, Directionally asymmetrical bistability in a symmetrically pumped nonlinear ring interferometer, Opt. Commun. **40**, 229 (1982).

[43] M. Haelterman, S. Trillo, and S. Wabnitz, Polarization multistability and instability in a nonlinear dispersive ring cavity, J. Opt. Soc. Am. B **11**, 446 (1994).

[44] J. Pond, Möbius dual-mode resonators and bandpass filters, IEEE Trans. Microw. Theory Tech. **48**, 2465 (2000).

[45] E. L. Starostin and G. H. M. van der Heijden, The shape of a Möbius strip, Nat. Mater. **6**, 563 (2007).

[46] D. J. Ballon and H. U. Voss, Classical Möbius-ring resonators exhibit fermion-boson rotational symmetry, Phys. Rev. Lett. **101**, 247701 (2008).

[47] C. Maitland, M. Conforti, A. Mussot, and F. Biancalana, Stationary states and instabilities of a Möbius fiber resonator, Phys. Rev. Research **2**, 043195 (2020).

[48] X.-B. Xu, L. Shi, G.-C. Guo, C.-H. Dong, and C.-L. Zou, "Möbius" microring resonator, Appl. Phys. Lett. **114**, 101106 (2019).

[49] G. K. L. Wong, S. G. Murdoch, R. Leonhardt, J. D. Harvey, and V. Marie, High-conversion-efficiency widely-tunable all-fiber optical parametric oscillator, Opt. Express **15**, 2947 (2007).

[50] H. J. Carmichael, Optical bistability and multimode instabilities, Phys. Rev. Lett. **52**, 1292 (1984).

[51] M. Haelterman and M. D. Tolley, Pure polarization period-doubling instability in a Kerr-type nonlinear ring cavity, Opt. Commun. **108**, 165 (1994).

[52] H. G. Schuster and W. Just, *Deterministic Chaos: An Introduction*, 4th ed. (Wiley-VCH, Weinheim, 2005).

[53] Y. K. Chembo and N. Yu, On the generation of octave-spanning optical frequency combs using monolithic whispering-gallery-mode microresonators, Opt. Lett. **35**, 2696 (2010).

[54] C. Menyuk, Pulse propagation in an elliptically birefringent Kerr medium, IEEE J. Quantum Elec. **25**, 2674 (1989).

[55] G. P. Agrawal, *Nonlinear Fiber Optics*, 5th ed. (Academic Press, 2013).

[56] J. Woo and R. Landauer, Fluctuations in a parametrically excited subharmonic oscillator, IEEE J. Quantum Elec. **7**, 435 (1971).

[57] C. Weitenberg and J. Simonet, Tailoring quantum gases by Floquet engineering, Nat. Phys. **17**, 1342 (2021).

[58] We use the relation $e^{i\bar{u}}e^{i\bar{v}}e^{-i\bar{u}} = \exp(e^{i\bar{u}}\bar{v}e^{-i\bar{u}})$, and the fact that $\hat{U}_{swap}$ is unitary.

[59] H. Lefevre, Single-mode fibre fractional wave devices and polarisation controllers, Electron. Lett. **16**, 778 (1980).

[60] A. J. Barlow, J. J. Ramskov-Hansen, and D. N. Payne, Birefringence and polarization mode-dispersion in spun single-mode fibers, Appl. Opt. **20**, 2962 (1981).

[61] A. U. Nielsen, B. Garbin, S. Coen, S. G. Murdoch, and M. Erkintalo, Coexistence and interactions between nonlinear states with different polarizations in a monochromatically driven passive Kerr resonator, Phys. Rev. Lett. **123**, 013902/1 (2019).

[62] A. U. Nielsen, B. Garbin, S. Coen, S. G. Murdoch, and M. Erkintalo, Invited Article: Emission of intense resonant radiation by dispersion-managed Kerr cavity solitons, APL Photonics **3**, 120804/1 (2018).

[63] T. Steinle, J. N. Greiner, J. Wrachtrup, H. Giessen, and I. Gerhardt, Unbiased All-Optical Random-Number Generator, Phys. Rev. X **7**, 041050 (2017).

[64] L. E. Bassham, A. L. Rukhin, J. Soto, J. R. Nechvatal, M. E. Smid, S. D. Leigh, M. Levenson, M. Vangel, N. A. Heckert, and D. L. Banks, *A Statistical Test Suite for Random and Pseudorandom Number Generators for Cryptographic Applications*, Special Publication (NIST SP) 800-22 Rev 1a (National Institute of Standards and Technology, Gaithersburg, MD, 2010).

[65] A. Amo, D. Sanvitto, F. P. Laussy, D. Ballarini, E. del Valle, M. D. Martin, A. Lemaitre, J. Bloch, D. N. Krizhanovskii, M. S. Skolnick, C. Tejedor, and L. Viña, Collective fluid dynamics of a polariton condensate in a semiconductor microcavity, Nature (London) **457**, 291 (2009).

[66] A. Amo, T. C. H. Liew, C. Adrados, R. Houdré, E. Giacobino, A. V. Kavokin, and A. Bramati, Exciton-polariton spin switches, Nat. Photon. **4**, 361 (2010).

[67] L. A. Lugiato and R. Lefever, Spatial dissipative structures in passive optical systems, Phys. Rev. Lett. **58**, 2209 (1987).

[68] A. J. Scroggie, W. J. Firth, G. S. McDonald, M. Tlidi, R. Lefever, and L. A. Lugiato, Pattern formation in a passive Kerr cavity, Chaos, Solitons & Fractals **4**, 1323 (1994).

[69] S. Barland, J. R. Tredicce, M. Brambilla, L. A. Lugiato, S. Balle, M. Giudici, T. Maggipinto, L. Spinelli, G. Tisoni, T. Knödl, M. Miller, and R. Jäger, Cavity solitons as pixels in semiconductor microcavities, Nature (London) **419**, 699 (2002).

[70] F. Leo, S. Coen, P. Kockaert, S.-P. Gorza, Ph. Emplit, and M. Haelterman, Temporal cavity solitons in one-dimensional Kerr media as bits in an all-optical buffer, Nat. Photon. **4**, 471 (2010).

[71] S. Coen, H. G. Randle, T. Sylvestre, and M. Erkintalo, Modeling of octave-spanning Kerr frequency combs us-
ing a generalized mean-field Lugiato-Lefever model, Opt. Lett. 38, 37 (2013).

[72] T. Herr, V. Brasch, J. D. Jost, C. Y. Wang, N. M. Kondratiev, M. L. Gorodetsky, and T. J. Kippenberg, Temporal solitons in optical microresonators, Nat. Photon. 8, 145 (2014).

[73] W. Weng, R. Bouchand, E. Lucas, E. Obrzud, T. Herr, and T. J. Kippenberg, Heteronuclear soliton molecules in optical microresonators, Nat. Commun. 11, 2402 (2020).

[74] G. Xu, A. U. Nielsen, B. Garbin, L. Hill, G.-L. Oppo, J. Fatome, S. G. Murdoch, S. Coen, and M. Erkintalo, Spontaneous symmetry breaking of dissipative optical solitons in a two-component Kerr resonator, Nat. Commun. 12, 4023 (2021).

[75] G. Xu, L. Hill, J. Fatome, G.-L. Oppo, M. Erkintalo, S. G. Murdoch, and S. Coen, Breathing dynamics of symmetry-broken temporal cavity solitons in Kerr ring resonators, Opt. Lett. 47, 1486 (2022).

[76] M. Haelterman and A. P. Sheppard, Polarization domain walls in diffractive or dispersive Kerr media, Opt. Lett. 19, 96 (1994).

[77] M. Gilles, P.-Y. Bony, J. Garnier, A. Picozzi, M. Guasoni, and J. Fatome, Polarization domain walls in optical fibres as topological bits for data transmission, Nat. Photon. 11, 102 (2017).

[78] J. K. Jang, M. Erkintalo, S. Coen, and S. G. Murdoch, Temporal tweezing of light through the trapping and manipulation of temporal cavity solitons, Nat. Commun. 6, 7370/1 (2015).

[79] M. Wang, F. Zhou, X. Lu, A. McClung, M. Davanco, V. A. Aksyuk, and K. Srinivasan, Fractional optical angular momentum and multi-defect-mediated mode renormalization and orientation control in photonic crystal microring resonators, Phys. Rev. Lett. 129, 186101 (2022).

[80] D. B. Kitchen, H. Decornez, J. R. Furr, and J. Bajorath, Docking and scoring in virtual screening for drug discovery: Methods and applications, Nat. Rev. Drug. Discov. 3, 935 (2004).

[81] P. J. Carrington, J. Scott, and S. Wasserman, eds., *Models and Methods in Social Network Analysis*, Structural Analysis in the Social Sciences No. 28 (Cambridge University Press, 2005).

[82] F. Barahona, On the computational complexity of Ising spin glass models, J. Phys. A: Math. Gen. 15, 3241 (1982).

[83] A. Lucas, Ising formulations of many NP problems, Front. Phys. 2, 5 (2014).

[84] M. Honari-Latifpour and M.-A. Miri, Optical Potts machine through networks of three-photon down-conversion oscillators, Nanophot. 9, 4199 (2020).

[85] M. Haelterman, S. Trillo, and S. Wabnitz, Dissipative modulation instability in a nonlinear dispersive ring cavity, Opt. Commun. 91, 401 (1992).

[86] M. Haelterman, Simple model for the study of period-doubling instabilities in the nonlinear ring cavity, Appl. Phys. Lett. 61, 2756 (1992).