Dynamics on invariant tori emerging through forced symmetry breaking in phase oscillator networks

Abstract

We consider synchrony patterns in coupled phase oscillator networks that correspond to invariant tori. For specific nongeneric coupling, these tori are equilibria relative to a continuous symmetry action. We analyze how the invariant tori deform under forced symmetry breaking as more general network interaction terms are introduced. We first show in general that perturbed tori that are relative equilibria can be computed using a parametrization method; this yields an asymptotic expansion of an embedding of the perturbed torus, as well as the local dynamics on the torus. We then apply this result to a coupled oscillator network, and we numerically study the dynamics on the persisting tori in the network by looking for bifurcations of their periodic orbits in a boundary-value-problem setup. This way we find new bifurcating stable synchrony patterns that can be the building blocks of larger global structures such as heteroclinic cycles.

Figures (6)

• Figure 1: Dynamics of system \ref{['eq:nf_SDD']} for $\alpha=\beta=\frac{\pi}{2}$, $r=0.2$ and $\delta=0.01$. Panel (a) shows solutions to the reduced system \ref{['eq:SDD_red']} on the torus $\mathbb{T}^2$ described by coordinates $(\phi_2,\phi_3)$ representing the '$\mathrm{D}$' populations; these periodic orbits wind around $\mathbb{T}^2$. Shown are two stable periodic orbits (blue) at $\psi\in\{\frac{\pi}{2}, \frac{3\pi}{2}\}$ and two unstable periodic orbits (red) at $\psi\in\{0, \pi\}$. Panel (b) shows a periodic orbit in the full system that lies on $\mathrm{SDD}$ for $\phi_3-\phi_2 = 0$. Note that the two '$\mathrm{D}$' populations have approximately the same phase.
• Figure 2: Bifurcations can stabilize periodic orbits of the full system \ref{['eq:Dynamics3x2']} stemming from $\mathrm{SDD}$ for a perturbation determined by \ref{['eq:Zj']} and \ref{['eq:h2']} that is sufficiently large $\delta>0$. The rest of the parameters are as in \ref{['eq:par']}. Panel (a) shows the four branches of hyperbolic periodic orbits $\gamma_0, \gamma_{\frac{\pi}{2}}, \gamma_{\pi}$ and $\gamma_{\frac{3\pi}{2}}$ emanating from $\psi^* \in \{0, \frac{\pi}{2},\pi,\frac{3\pi}{2}\}$ at $\delta=0$, bifurcations that change their stability in the full system, and secondary solution branches. The vertical green line at $\delta=0$ represents the unperturbed case, in which the torus $\mathrm{SDD}$ is foliated by neutrally stable periodic orbits. Panel (b) highlights bifurcations at $\delta=\delta^*$ that stabilize the branches $\gamma_0$ and $\gamma_{\pi}$ in the full system and corresponding secondary branches of periodic orbits forming isolas.
• Figure 3: Stable synchrony patterns bifurcating from primary branches of periodic orbits. Panels (a)--(d) show the phase evolution of oscillators in '$\mathrm{D}$' populations along one period on coexisting stable periodic orbits lying on the secondary branches emanating from $\psi=0$ (bottom) and $\psi=\pi$ (top) for $\delta\approx 0.40372$; see the magnified parameter range in Figure \ref{['fig:pobranches_sdd']}(b). Shading indicates the deviation from $\pi$ of each phase oscillator, where black indicates that $\theta_{\sigma,k}=\pi$ and white $\theta_{\sigma,k}=0$ or $2\pi$.
• Figure 4: Some codimension-two bifurcations of system \ref{['eq:Dynamics3x2']} with the perturbation driven by \ref{['eq:Zj']} and \ref{['eq:h2']}. Panel (a) shows curves of saddle-node bifurcation arising on the primary (red) and secondary (magenta) branches that organize the parameter plane $(r,\delta)$. See the main text for an elaboration on the codimension-two points shown. Panel (b) shows the origin of the primary saddle-node bifurcation as one continues a primary branch of periodic orbits in parameter $r$, with $r$ decreasing. This bifurcation gives rise to the curve $\mathsf{SN}_p$ in Panel (a).
• Figure 5: Evolution of periodic orbits along the overlapping principal saddle-node curves $\mathsf{SN}_p$ in Figure \ref{['fig:R_delta_sdd']}. Panels (a)--(d) correspond to the points $r_a,r_b, r_c, r_d$, respectively. Each panel shows (1) a time series of the coordinate $\theta_{2,1}$ with respect to the rescaled Auto integration time for orbits on each principal saddle-node curve in green and red, (2) a relative position between the two '$\mathrm{D}$' populations on each principal saddle-node curve in green and red, and (3) how each '$\mathrm{D}$' populations deviate from $\pi$ as $(r,\delta)$ moves away from the origin.

Theorems & Definitions (12)

• Proposition 2.1
• proof : Proof of Proposition \ref{['prop:Embed']}
• Proposition 2.2
• proof : Proof of Proposition \ref{['lemmarelativeequilibrium']}
• Proposition 2.3
• proof : Proof of Proposition \ref{['prop:tangentialdynamics']}
• Proposition 2.4
• proof : Proof of Proposition \ref{['prop:splittingprop']}
• Proposition 2.5
• proof : Proof of Proposition \ref{['prop:Ansatz']}