Third-Order Sextupole Resonance Conditions

• Third-order sextupole resonance conditions are tune relations where sextupole-driven nonlinearities become resonant, impacting beam stability and phase-space structures.
• They are fundamental in storage rings and synchrotrons, organizing resonance lines and governing phenomena like transverse resonance island buckets and slow extraction efficiency.
• Coupled normal-mode analysis reveals expanded resonance conditions and guides lattice optimization through resonance driving terms, offering improved tune flexibility and beam stability.

Third-order sextupole resonance conditions are tune relations at which sextupole-driven nonlinearities become resonant in the transverse dynamics of charged-particle beams. In storage rings and synchrotrons, they organize the appearance of resonance lines in tune space, the formation of separatrices and fixed points in phase space, the excitation of transverse resonance island buckets, and the efficiency of slow extraction. In uncoupled optics they appear as the familiar four third-order lines; in coupled normal-mode parametrization they expand to eight conditions; and in circular-mode operation most of these lines are naturally suppressed because the weak mode is intrinsically flat, with $J_2 \ll J_1$ (Gilanliogullari et al., 9 Sep 2025).

1. Uncoupled third-order resonance lines

For normal sextupoles in an uncoupled lattice, the third-order resonance conditions are the tune combinations

$$3\nu_x = p,\qquad 2\nu_x + \nu_y = p,\qquad \nu_x + 2\nu_y = p,\qquad 3\nu_y = p,$$

which are the third-order resonance lines in the tune diagram (Li et al., 25 Mar 2026). The same structure appears in machine-specific settings: slow extraction is conventionally formulated near

$\nu_x = \frac{k}{3} + \delta,$

with $\delta$ a small detuning from exact resonance (Nikitin, 2024), while the CESR study designed a $6$-GeV lattice with the horizontal tune near the third-order line $3\nu_x = 50$ (Wang et al., 2023).

Condition	Context in the data	Representative role
$3\nu_x = p$	Uncoupled third-order line	Horizontal sextupole resonance
$2\nu_x + \nu_y = p$	Uncoupled third-order line	Coupled-plane third-order line
$\nu_x + 2\nu_y = p$	Uncoupled third-order line	Coupled-plane third-order line
$3\nu_y = p$	Uncoupled third-order line	Vertical sextupole resonance

These conditions define where sextupole fields can excite unstable or weakly stable motion. The storage-ring literature in the data also makes clear that proximity to a resonance line is not, by itself, a sufficient criterion for a specific phase-space structure. At CESR, transverse resonance island buckets (TRIBs) are not always present near a resonance line; necessary conditions for their formation include a nonzero resonant driving term and a tune-dependent bifurcation constraint (Wang et al., 2023). This addresses a frequent simplification in which the presence of a line in the tune diagram is treated as equivalent to island formation.

2. Coupled normal-mode parametrization

When coupling is significant, the uncoupled Courant–Snyder description is replaced by normal-mode analysis. In the formulation summarized from the circular-mode study, a particle’s phase-space vector is written in terms of mode actions $$3\nu_x = p,\qquad 2\nu_x + \nu_y = p,\qquad \nu_x + 2\nu_y = p,\qquad 3\nu_y = p,$$0, phases $$3\nu_x = p,\qquad 2\nu_x + \nu_y = p,\qquad \nu_x + 2\nu_y = p,\qquad 3\nu_y = p,$$1, and eigenvectors $$3\nu_x = p,\qquad 2\nu_x + \nu_y = p,\qquad \nu_x + 2\nu_y = p,\qquad 3\nu_y = p,$$2, following the Lebedev–Bogacz parametrization (Gilanliogullari et al., 9 Sep 2025). In this setting, a normal sextupole introduces the potential

$$3\nu_x = p,\qquad 2\nu_x + \nu_y = p,\qquad \nu_x + 2\nu_y = p,\qquad 3\nu_y = p,$$3

with $$3\nu_x = p,\qquad 2\nu_x + \nu_y = p,\qquad \nu_x + 2\nu_y = p,\qquad 3\nu_y = p,$$4 the sextupole strength (Gilanliogullari et al., 9 Sep 2025).

After substitution of the coupled-mode coordinates into $$3\nu_x = p,\qquad 2\nu_x + \nu_y = p,\qquad \nu_x + 2\nu_y = p,\qquad 3\nu_y = p,$$5, the sextupole potential becomes a sum of oscillatory terms involving the harmonics $$3\nu_x = p,\qquad 2\nu_x + \nu_y = p,\qquad \nu_x + 2\nu_y = p,\qquad 3\nu_y = p,$$6 and $$3\nu_x = p,\qquad 2\nu_x + \nu_y = p,\qquad \nu_x + 2\nu_y = p,\qquad 3\nu_y = p,$$7. The summary gives the resulting structure as

$$3\nu_x = p,\qquad 2\nu_x + \nu_y = p,\qquad \nu_x + 2\nu_y = p,\qquad 3\nu_y = p,$$8

together with

$$3\nu_x = p,\qquad 2\nu_x + \nu_y = p,\qquad \nu_x + 2\nu_y = p,\qquad 3\nu_y = p,$$9

where the coefficients depend on the actions and optical functions (Gilanliogullari et al., 9 Sep 2025).

The stationary-phase condition for these terms produces eight resonance conditions in the coupled normal-mode parametrization: $\nu_x = \frac{k}{3} + \delta,$0

$\nu_x = \frac{k}{3} + \delta,$1

The paper explicitly labels $\nu_x = \frac{k}{3} + \delta,$2 and $\nu_x = \frac{k}{3} + \delta,$3 as single-mode resonances, and $\nu_x = \frac{k}{3} + \delta,$4 and $\nu_x = \frac{k}{3} + \delta,$5 as integer resonances (Gilanliogullari et al., 9 Sep 2025). In comparison, in an uncoupled lattice only four third-order resonance lines exist. The coupled description therefore enlarges the set of possible resonance relations, rather than merely rotating the uncoupled picture.

3. Circular modes and suppression of sextupole resonance lines

Circular modes are described in the data as round coupled beams with non zero angular momentum, providing an alternative beam motion and dynamics (Gilanliogullari et al., 9 Sep 2025). Their defining operational feature is strong flatness: one eigenmode emittance is much smaller than the other, $\nu_x = \frac{k}{3} + \delta,$6, equivalently $\nu_x = \frac{k}{3} + \delta,$7 (Gilanliogullari et al., 9 Sep 2025). The same summary states that the beam is “round” in coupled phase space but is dynamically one-dimensional.

In that regime the dynamics are dominated by mode $\nu_x = \frac{k}{3} + \delta,$8. Because $\nu_x = \frac{k}{3} + \delta,$9, all resonance terms with amplitudes proportional to $\delta$0 or mixed products vanish, and only the amplitudes $\delta$1 and $\delta$2 in the coupled sextupole expansion remain nonzero (Gilanliogullari et al., 9 Sep 2025). The surviving resonance conditions are therefore

$\delta$3

This is the central third-order result of the circular-mode paper: most third-order resonance lines are naturally suppressed in the tune diagram for a circular-mode beam (Gilanliogullari et al., 9 Sep 2025). The comparison given in the data is explicit. In coupled normal-mode lattices there are eight lines, but for circular-mode operation only two, those associated with the dominant eigenmode, contribute. The same source further states that simulations show that intentionally placing the weak mode on resonance does not lead to beam instability or loss, unlike the dominant mode. The immediate operational consequence is greater tune flexibility and a greatly expanded usable tune space, especially when space charge or other collective effects limit tune choice (Gilanliogullari et al., 9 Sep 2025).

4. Hamiltonian structure, separatrices, and adiabatic crossing

Near a third-integer resonance driven by sextupoles, the nonlinear dynamics are commonly represented by the Kobayashi Hamiltonian

$\delta$4

where the fractional tune is written as $\delta$5 with $\delta$6 and $\delta$7 (Niedermayer, 2024). In normalized phase space, the same work gives the coordinate transformation

$\delta$8

The phase-space geometry is triangular, and the separatrix boundary is located at

$\delta$9

with inscribed triangle radius

$6$0

(Niedermayer, 2024).

A central result of that analysis is that detuning near the resonance is not only amplitude-dependent but also phase-dependent. The three-turn phase advance is written as

$6$1

and the average detuning decreases as the trajectory approaches the separatrix, vanishing as $6$2 (Niedermayer, 2024). This explicitly refines amplitude-only descriptions of the resonance neighborhood.

For slow extraction, the adiabatic crossing theory summarized from VEPP-4M starts from

$6$3

and imposes the adiabaticity condition

$6$4

with invariant action integral

$6$5

until the particle reaches the separatrix (Nikitin, 2024). Chromaticity, momentum spread, and synchrotron oscillations modify the effective detuning, and with RF on the tune is modulated so that sideband resonances appear at

$6$6

(Nikitin, 2024). The same work states that introducing a small controlled acceleration during the crossing can monochromatize the extracted beam when

$6$7

The experimental comparison at the VEPP-4M storage ring used the third-order resonance $6$8 and low chromaticity $6$9. The measured time profiles of losses were asymmetric, with a slow rise and a sharper leading edge, matching predictions of the adiabatic theory for low chromaticity (Nikitin, 2024). The same source notes that non-total beam loss is attributed to effects not included in the adiabatic theory, such as radiative damping and/or higher-order nonlinearities.

5. Resonance driving terms and beam-based diagnosis

The resonance-driving-term formalism provides a lattice-based description of third-order sextupole resonances. In the storage-ring optimization study, the third-order RDT at observation point $3\nu_x = 50$0 is written as

$3\nu_x = 50$1

where the indices satisfy $3\nu_x = 50$2 and the coefficient for a normal sextupole is

$3\nu_x = 50$3

(Li et al., 25 Mar 2026). The denominator becomes small when

$3\nu_x = 50$4

which reproduces the third-order tune combinations $3\nu_x = 50$5, $3\nu_x = 50$6, $3\nu_x = 50$7, and $3\nu_x = 50$8 (Li et al., 25 Mar 2026).

The turn-by-turn beam-position-monitor measurement framework gives a complementary experimental realization of the same objects. The first-order sextupole RDTs are written in the same general form, and single-BPM data provide combined RDTs rather than individual ones (Franchi et al., 2014). The data explicitly list the sextupolar CRDT combinations

$3\nu_x = 50$9

$3\nu_x = p$0

associated with the spectral lines $3\nu_x = p$1, $3\nu_x = p$2, $3\nu_x = p$3, and $3\nu_x = p$4, respectively (Franchi et al., 2014). These observables are extracted from FFT amplitudes and phases after normalization to Courant–Snyder coordinates.

This measurement formalism is not only diagnostic but corrective. The ESRF study states that effective sextupole magnetic errors and tilts were evaluated and corrected when possible; octupolar RDTs were also measured; and most of the deviations from the model observed in the sextupolar RDTs turned out to be generated by focusing errors rather than by sextupole errors (Franchi et al., 2014). That point is important because it distinguishes resonance diagnosis from magnet attribution: a measured sextupolar spectral signature does not uniquely imply a sextupole fault.

6. TRIBs, lattice design, and dynamic-aperture optimization

Near the resonance $3\nu_x = p$5, the CESR study describes the horizontal dynamics with the normalized Hamiltonian

$3\nu_x = p$6

where $3\nu_x = p$7, $3\nu_x = p$8 are amplitude-dependent tune shift coefficients, and $3\nu_x = p$9 is the complex resonance-strength coefficient (Wang et al., 2023). Stable resonance islands form only if $2\nu_x + \nu_y = p$0. The bifurcation condition for the existence of stable fixed points is

$2\nu_x + \nu_y = p$1

and only solutions with $2\nu_x + \nu_y = p$2 are physical (Wang et al., 2023).

The same work distinguishes two classes of TRIBs. In the first type, three stable fixed points are intercalated with three unstable fixed points in phase space. In the second type, stable and unstable fixed points coincide in phase, and the operational window is narrower (Wang et al., 2023). At CESR, particle tracking, PTC map-based analysis, and visible-light beam-size-monitor observations all showed islands near the designed third-order resonance, and a sextupole knob changed the phase $2\nu_x + \nu_y = p$3 so that the islands rotated in phase space, in good agreement with the theoretical calculation (Wang et al., 2023).

For broader lattice optimization, the longitudinal variation of third-order RDTs is quantified by

$2\nu_x + \nu_y = p$4

and the same study proves that

$2\nu_x + \nu_y = p$5

is a convex quadratic function of the sextupole strengths (Li et al., 25 Mar 2026). The iso-surfaces of $2\nu_x + \nu_y = p$6 are a series of concentric and coaxial ellipsoidal surfaces, with the central position possessing minimum $2\nu_x + \nu_y = p$7. Global scanning in FODO, HLS-III, and SSRF shows a strong consistency between the distributions of $2\nu_x + \nu_y = p$8 and dynamic aperture, so that dynamic-aperture optimization can be regarded as a roughly approximate convex optimization problem (Li et al., 25 Mar 2026).

Taken together, these results place third-order sextupole resonance conditions at the intersection of tune selection, nonlinear diagnostics, beam extraction, and lattice optimization. In uncoupled optics they delimit the standard third-order lines; in coupled optics they acquire a larger normal-mode structure; in circular modes most of that structure becomes dynamically irrelevant; and in RDT-based analysis they become measurable, optimizable, and, in specific machines, directly visible through resonance islands and spill profiles (Gilanliogullari et al., 9 Sep 2025).