--- title: Secondary Periodic Orbits url: https://www.emergentmind.com/topics/secondary-periodic-orbits type: topic --- # Secondary Periodic Orbits Secondary periodic orbits are periodic motions that arise beyond the simplest closed trajectories of a dynamical system and are attached to a primary structure by resonance, bifurcation, continuation, symmetry, or hierarchical modulation. Across recent and classical arXiv literature, the expression does not denote a single universally fixed object. In nonlinear vibration it can denote the periodic orbit of a reduced-order model whose lift is a two-frequency torus created in a secondary Hopf bifurcation; in fluid mechanics it can denote periodic states born from secondary bifurcations of coherent equilibria; in celestial mechanics it can denote continued or bifurcated orbit families; in billiards it aligns with additional non-Birkhoff branches beyond the ordered Birkhoff polygons; and in relativistic timing problems it can refer to a secondary orbital clock superposed on a primary recurrence [2412.02100] [1205.0347] [1811.05897] [2504.03574] [2604.09788]. ## 1. Terminological range and dynamical role The literature shows that “secondary periodic orbit” is field-dependent rather than universal. In the forced \(1\!:\!2\) internally resonant mechanical systems of "Bifurcation analysis of quasi-periodic orbits of mechanical systems with 1:2 internal resonance via spectral submanifolds" [2412.02100], a **secondary Hopf bifurcation** is a bifurcation of a periodic orbit of the forced system, not of an equilibrium, and the resulting object in the full system is a **2D invariant torus**. In that setting, the paper’s “secondary periodic orbit” language refers to the periodic orbit of the \(4\)D reduced-order model whose lift is quasi-periodic motion in the full system. In "Periodic orbits near onset of chaos in plane Couette flow" [1205.0347], secondary periodic orbits are periodic states born from secondary bifurcations of the upper-branch coherent equilibrium. In "A Family of Periodic Orbits in the Three-Dimensional Lunar Problem" [1811.05897], the term is used in the broad celestial-mechanics sense of a family continued from a simple limit family rather than a primary Keplerian one. By contrast, the galactic-potential and irregular-body papers do not formally define the term, but they study resonant or bifurcated periodic-orbit families that play the same structural role [1202.4084] [1610.09440]. | Context | Meaning of “secondary” | Representative paper | |---|---|---| | Forced internal resonance | Torus-producing orbit born from a periodic orbit via secondary Hopf | [2412.02100] | | Transitional shear flow | Periodic states from secondary bifurcations of coherent equilibria | [1205.0347] | | Celestial mechanics | Continued family beyond a primary/simple limit family | [1811.05897] | | Irregular-body dynamics | Additional branches/families generated during continuation and bifurcation | [1610.09440] | | Symmetric billiards | Non-Birkhoff periodic branches beyond ordered Birkhoff orbits | [2504.03574] | | Strong-field timing | Secondary orbital clock or super-period modulating a primary recurrence | [2604.09788] | A recurrent feature is that these objects are identified not merely by closure in time, but by their relation to another invariant structure: a resonance, an equilibrium branch, a parent periodic orbit, a symmetry class, or an underlying orbital clock. This suggests that “secondary” is best understood relationally rather than absolutely. ## 2. Resonance-supported families in smooth Hamiltonian potentials In "A semi-numerical method for periodic orbits in a bisymmetrical potential" [1202.4084], the underlying model is the two-dimensional harmonic oscillator plus a Plummer term, \[ V(x,y)=\frac{1}{2}\left(\omega_1^2 x^2+\omega_2^2 y^2\right)-\frac{\mu}{\sqrt{x^2+y^2+\alpha^2}}, \] with \[ \mu=0.001,\qquad \alpha=0.25, \] and resonances \[ \omega_1:\omega_2 = 1:1,\quad 2:1,\quad 2:3,\quad 4:3. \] The paper studies periodic orbits intersecting the \(x\)-axis perpendicularly and periodic orbits going through the origin. It does not explicitly classify them as primary or secondary, but it explicitly notes that these resonant families are exactly the kinds of orbit families that would commonly be regarded as secondary or family orbits associated with a resonance. The method is semi-numerical in the strict sense that it combines empirical formulas, harmonic-oscillator intuition, and calibration against direct numerical integration. The full equations are integrated with a **Bulirsh–Stöer** routine in double precision, with energy conserved to about the \(12\)th significant figure. The aim is not exact closed-form integration except in the symmetric \(1\!:\!1\) case, but accurate formulas for starting positions and periods. For periodic orbits starting perpendicularly from the \(x\)-axis, the \(1\!:\!1\) resonance is analytically tractable because axial symmetry yields a circular periodic orbit with exact period \[ T_c=\frac{2\pi}{\left[\omega^2+\frac{\mu}{(r_c^2+\alpha^2)^{3/2}}\right]^{1/2}}. \] The \(2\!:\!1\) case is less accurate in position than in period: the semi-numerical starting points differ from direct numerical integration by up to about \(9\%\), whereas the periods differ by less than \(1\%\). In the \(2\!:\!3\) and \(4\!:\!3\) resonances, the agreement is described as **very good**. For origin-crossing families, exact straight-line periodic solutions exist in the symmetric \(1\!:\!1\) case; good semi-numerical estimates are obtained for \(2\!:\!3\) and \(4\!:\!3\); and the authors **do not find** periodic orbits through the origin in the \(2\!:\!1\) case, even when energy is increased well beyond the tabulated range. The significance of this study is methodological as well as dynamical. It presents secondary-like resonant families through their geometric character—perpendicular intersection with the \(x\)-axis or passage through the origin—and through the resonance ratio constraining the initial data by the energy integral. In this formulation, secondary periodic orbits are not introduced as an abstract class; they are the concrete resonance-supported backbone orbits of the model. ## 3. Secondary bifurcation, invariant tori, and routes to chaos In forced internally resonant mechanical systems, secondary periodic orbits appear through a precise codimension-one mechanism. The full system considered in [2412.02100] is \[ \boldsymbol{M}\ddot{\boldsymbol{x}}+\boldsymbol{C}\dot{\boldsymbol{x}}+\boldsymbol{K}\boldsymbol{x}+\boldsymbol{f}(\boldsymbol{x},\dot{\boldsymbol{x}})=\epsilon \boldsymbol{f}^{\mathrm{ext}}(\Omega t), \qquad 0<\epsilon\ll1, \] with \(1\!:\!2\) internal resonance \[ \lambda_2 \approx 2\lambda_1,\qquad \bar{\lambda}_2 \approx 2\bar{\lambda}_1, \] and external resonance \[ \mathrm{i}\Omega \approx \lambda_1. \] Reduction on a spectral submanifold \(\mathcal{W}(\mathcal{E})\) yields a \(4\)D reduced-order model. The paper emphasizes a structural equivalence: a **fixed point** of the reduced vector field corresponds to a **periodic orbit** of the full system, whereas a **limit cycle** of the reduced model corresponds to a **2D invariant torus** of the full system. Thus, the secondary periodic orbit is periodic only in the reduced system; in the full system it is quasi-periodic. Within this framework, the forced response curve of periodic solutions is typically M-shaped and undergoes **secondary Hopf bifurcations**, also called **Neimark–Sacker bifurcations**. When the forcing amplitude \(\epsilon\) exceeds a critical value, the periodic response loses stability and a family of quasi-periodic responses appears. The paper continues these torus branches by continuing limit cycles of the reduced model and detects local bifurcations—**period-doubling** and **saddle-node** bifurcations of quasi-periodic orbits—as well as global events including **homoclinic bifurcations**, **isolas**, **simple bifurcations**, and **cusp bifurcations**. In the benchmark two-coupled-oscillator example, the quasi-periodic branch terminates in a homoclinic bifurcation at \(\epsilon=0.03\). In the shallow curved beam, the authors report **20 PD and 20 SN points on the FRC-QO at \(\epsilon=0.01\)**. In the high-dimensional shallow shell example, the finite-element model has **more than \(1300\) DOFs** and is reduced to a \(4\)D ROM. A different but related picture appears in plane Couette flow [1205.0347]. There, the primary coherent structure is the **Nagata–Busse–Clever state** created in a saddle-node at \[ Re = 163.8. \] The upper branch, initially stable in the symmetry-reduced subspace, loses stability in a **supercritical Hopf bifurcation** at \[ Re_H = 166.05, \] with limit-cycle amplitude scaling like \[ \text{amplitude} \sim \sqrt{Re - Re_H}. \] That first Hopf-born cycle is the first secondary periodic orbit of the flow. It then undergoes period doubling at \[ Re = 178.9, \] followed by an incomplete doubling cascade, a stable **period-3 window** from \[ Re = 183.8 \quad \text{to} \quad Re = 185.0, \] and a **boundary crisis** at \[ Re_c = 188.51. \] The paper identifies seven periodic orbits of symbolic period at most five, \[ \overline{1},\; \overline{01},\; \overline{001},\; \overline{011},\; \overline{0111},\; \overline{01101},\; \overline{01111}, \] and states that **all of these periodic orbits have exactly one unstable eigenvalue**. Their symbolic organization is “almost the same” as the universal Metropolis–Stein–Stein sequence for unimodal maps. Taken together, these studies show that secondary periodic orbits are often the dynamically relevant descendants of an initial periodic response. They can be born from torus bifurcations, period doublings, or saddle-node events, and they can organize transitions from regular motion to chaotic or transiently chaotic dynamics. ## 4. Continuation, symmetry, and family structure in celestial mechanics In spatial three-body dynamics, secondary periodic orbits often appear as continued families from singular or simplified limit problems. The spatial lunar problem studied in [1811.05897] starts from the Hamiltonian of the three-dimensional circular restricted three-body problem in a rotating frame and then applies the standard lunar scaling \[ \hat q=\mu^{-1/3}q,\qquad \hat p=\mu^{-1/3}p. \] At \(\mu=0\), corresponding to Hill’s lunar problem, the \(\hat q_3\)-axis is invariant, and one obtains a family of **consecutive collision orbits** on that axis. These orbits start at collision with the small primary, move out along the positive \(\hat q_3\)-axis, return to collision, and repeat. Their turning point \(d(h)\) satisfies \[ V(d)=h \quad\Longleftrightarrow\quad z^3-2hz-2=0, \] and the period is \[ T^*=2\int_0^{d(h)}\frac{dz}{\sqrt{2h+\frac2z-z^2}}, \] with the key estimate \[ T^*\le \pi. \] The main theorem asserts that, for \(\mu\) sufficiently small and for all but two energies where a non-degeneracy condition fails, there exists a **unique periodic orbit** \[ \hat\phi(t,h,\mu) \] continuing the collision orbit \(\hat\phi^*(t,h)\), with \(T(\mu)\to T^*\) as \(\mu\to0\). The family is denoted \(\mathcal F(h,\mu)\) and is symmetric with respect to the \(\hat q_2\hat q_3\)-plane. This is a canonical secondary family: it is not a primary Keplerian orbit, but a persistent periodic orbit obtained by continuation from a singular limit family. The same paper also gives numerical information on the family’s stability and geometry. In the rescaled lunar problem, the polar orbit is **elliptic** for \[ h<-1.03, \] and **complex hyperbolic** for \[ h>0.11. \] It undergoes a period-doubling bifurcation, two simple degeneracies, and a period-halving bifurcation. The reported intervals are \[ [-1.025245,-1.025225],\ [0.0909615,0.0909616] \] for period doubling and halving, \[ [-0.85556,-0.85555],\ [0.043843,0.043844] \] for simple degeneracies, and \[ [0.109989,0.109990] \] for eigenvalue collision. The oblate-secondary variant of the spatial Hill problem in [2001.00120] provides a different existence result. After passing to Poincaré-Delaunay variables \[ Q_1=l+g+h,\quad Q_2=-\sqrt{2(L-G)}\sin(g+h),\quad Q_3=l+g, \] \[ P_1=L-G+H,\quad P_2=\sqrt{2(L-G)}\cos(g+h),\quad P_3=G-H, \] and removing the short-period terms by two successive Lie-transform averagings, the paper proves the existence of a class of **doubly-symmetric, near-circular periodic solutions** around the oblate primary. The system is reversible with respect to two anti-symplectic involutions \(R_1\) and \(R_2\), and any orbit intersecting the corresponding Lagrangian planes \(L_1\) and \(L_2\) in the prescribed way is periodic with period \(4T\). The paper explicitly states that these are spatial analogues of the doubly symmetric secondary periodic orbits known in the spherical Hill problem. It also notes that stability and global bifurcations remain open. These celestial-mechanics examples make clear that secondary periodic orbits are frequently family-level objects. Their central features are persistence under perturbation, symmetry constraints, and the capacity to bridge singular limit dynamics and the full problem. ## 5. Relativistic geodesics and hierarchical orbital clocks In strong-field gravitation, secondary periodic-orbit structure appears both geometrically and observationally. For timelike geodesics around a Schwarzschild-like black hole immersed in a King-type dark matter halo [2511.10043], the motion is studied in a static spherically symmetric spacetime with effective potential \[ V_{\rm eff}=f(r)\left(1+\frac{L^{2}}{r^{2}}\right), \] and the periodic orbits are organized by the zoom-whirl triplet \((z,w,v)\). The rational classification parameter is \[ q=\frac{\omega_{\phi}}{\omega_r}-1=w+\frac{v}{z}. \] The authors emphasize that every periodic orbit can be labeled uniquely by \((z,w,v)\). The specific classes displayed are \[ (1,1,0),\ (1,2,0),\ (2,1,1),\ (2,2,1),\ (3,1,2),\ (3,2,2),\ (4,1,3),\ (4,2,3). \] In the paper’s interpretation, the “secondary” periodic orbits are the higher members of this hierarchy beyond the simplest class, such as \[ (2,1,1),\ (2,2,1),\ (3,1,2),\ (3,2,2),\ (4,1,3),\ (4,2,3). \] Increasing the halo parameters \(r_s\) and \(\rho_s\) moves both \(R_{\mathrm{MBO}}\) and \(R_{\mathrm{ISCO}}\) outward, increases \(L_{\mathrm{MBO}}\) and \(L_{\mathrm{ISCO}}\), and decreases \(E_{\mathrm{ISCO}}\). For periodic orbit classes, the halo causes the required energies to decrease slightly and the required angular momenta to increase. For example, at \(r_s=0.2\), the energy of \((1,1,0)\) changes from \[ 0.965425 \quad \text{at}\ \rho_s=0.0 \] to \[ 0.962949 \quad \text{at}\ \rho_s=0.3, \] while \(L_{(1,1,0)}\) changes from \[ 3.68358 \quad \text{to} \quad 4.00114. \] A complementary use of secondary orbital periodicity appears in the QPE timing analysis of eRO-QPE2 [2604.09788]. There the observed eruption times are modeled through an \(O-C\) analysis with a base period and super-periodic modulations. The recurrence inferred from the XMM1–4 baseline is \[ P = (8064 \pm 7)\,\mathrm{s} \approx 2.24~\mathrm{h}. \] The best-fit deterministic timing model adds two sinusoidal super-periods. The shorter modulation is measured at \[ P_{\mathrm{mod,1}} = (47.0\pm0.6)\,P \approx 4.4~\mathrm{d}, \] with amplitude \[ A_{\mathrm{mod,1}} \approx 289^{+30}_{-31}\,\mathrm{s}, \] and the longer modulation settles near \[ P_{\mathrm{mod,2}} \approx 1016\,P \approx 95~\mathrm{d}, \] with amplitude \[ A_{\mathrm{mod,2}} \approx 3521^{+803}_{-938}\,\mathrm{s}. \] The short modulation is interpreted as apsidal precession of an eccentric orbit, with \[ M_{\rm BH}\approx 1.5\times 10^5\,M_\odot,\qquad a \sim 141\,R_g,\qquad e \approx 0.1. \] The long modulation is inconsistent with EMRI nodal precession; disk precession is allowed only in a limited region of parameter space; and the paper identifies a stable hierarchical triple solution with \[ a_{\rm bin} \sim 0.4~\mathrm{mpc},\qquad M_2 \sim (0.1-1)\,M_{\rm BH}. \] The authors also report \[ |\dot{P}| \lesssim 2\times 10^{-6}\,\mathrm{s\,s^{-1}} \] at \(3\sigma\), disfavoring high-eccentricity white dwarfs and some IMBH secondaries via GW decay, and they state that the correlated odd/even \(O-C\) disfavors both disk crossings per orbit being observed. In this relativistic setting, secondary periodic orbits need not be separate branches in phase space. They can instead be secondary clocks: additional periodicities encoding precession or outer-orbit structure on top of a primary recurrence. ## 6. Symmetry, topology, and non-primary orbit generation Around irregular-shaped minor bodies, the language of secondary periodic orbits is most naturally translated into branch generation during continuation. In [1610.09440], the dynamics is formulated in the rotating body-fixed frame with gravitational potential obtained from the polyhedron model and effective potential \[ V = U - \frac{1}{2}\omega^2(x^2+y^2), \] while the Jacobi integral is \[ H = U + \frac{1}{2}( \dot{x}^2 + \dot{y}^2 + \dot{z}^2 ) - \frac{1}{2}\omega^2(x^2+y^2). \] The paper proves a conserved topological-degree quantity that restricts the number of periodic orbits on a fixed energy surface. Its consequences are that the number of non-degenerate periodic orbits varies in **pairs** during continuation and is always **odd** on the fixed energy surface. The numerical procedure fixes a Poincaré section, randomly searches for periodic orbits, uses a tolerance of \(1.0\times10^{-6}\), selects roughly \(60\) candidates from about \(10{,}000\) computed orbits, computes the monodromy matrix \(M\), and classifies the Floquet multipliers. The bifurcation structure is encoded by multiplier collisions. The paper states that a **period-doubling bifurcation** occurs when two multipliers cross \(-1\), a **tangent bifurcation** when two cross \(+1\), a **Neimark–Sacker bifurcation** when two collide on the unit circle away from \(\pm1\), and a **real saddle bifurcation** when two collide on the real axis away from \(\pm1\). It also stresses an important caveat: collision alone does not necessarily imply a true bifurcation. Around comet **1P/Halley**, the authors identify a stable approximately **1:1 resonant** family with period ratios in \([1.01193,1.01243]\) undergoing a **pseudo Neimark–Sacker bifurcation**, where the collision persists without genuine branch splitting. Around asteroid **216 Kleopatra**, a single continued family passes through \[ \text{N3} \to \text{DR1} \to \text{N1} \to \text{DR2} \to \text{N5} \to \text{PD5} \to \text{N7}, \] thereby exhibiting **two real saddle bifurcations** and **one period-doubling bifurcation**, with orbit periods varying from **2.4257 h to 2.3855 h**. In symmetric convex billiards [2504.03574], the corresponding non-primary objects are **non-Birkhoff periodic orbits**. A periodic billiard sequence is Birkhoff if cyclic order is preserved under time shifts; a non-Birkhoff periodic orbit fails this ordering property. For rational rotation number \(m/n\) with \(\gcd(m,n)=1\), a Birkhoff orbit has minimal period \(n\), whereas a non-Birkhoff orbit with the same rotation number may have minimal period \(p\gg n\). The principal existence criterion states that if \[ \kappa L < 2\sin\!\left(\frac{m\pi}{n}\right)\cos^2\!\left(\frac{N\pi}{p}\right), \] then the billiard has a non-Birkhoff periodic orbit with minimal period \(p=sn\), winding number \(q=sm\), rotation number \(q/p=m/n\), and spatiotemporal symmetry group \(H\). The paper further states that the lifts of the non-Birkhoff orbit and the reference Birkhoff orbit cross exactly \(2N\) times per period \(p\). If \[ \kappa L < 2\sin\!\left(\frac{m\pi}{n}\right), \] then there are **infinitely many non-Birkhoff periodic orbits** of rotation number \(m/n\) with arbitrarily long minimal periods. In the \(\mathbb D_2\)-symmetric case, the condition \[ \kappa L < 2 \] yields infinitely many non-Birkhoff periodic orbits of rotation number \(1/2\), with arbitrarily long minimal periods, of types I, II, and V. The paper also proves that any open neighborhood of the circular billiard in the analytic topology contains a billiard with infinitely many non-Birkhoff periodic orbits of any rational rotation number \(0