Spin-Axis Tossing: Cross-Disciplinary Dynamics

• Spin-Axis Tossing is defined as the phenomenon of significant reorientations in a body's spin, driven by external torques, resonance crossings, or state-transition frameworks.
• It has diverse applications, evidenced by planetary obliquity excitation during scattering, spin randomization in black-hole formation, stabilized tilting in circumbinary systems, and even rotational modeling in juggling.
• The concept unifies broad disciplinary insights while raising challenges in measurement, theoretical modeling, and the interpretation of spin dynamics across astrophysics and applied mathematics.

Spin-axis tossing denotes, in the literatures represented here, either a physical reorientation of a spin axis or a formal representation of rotational-state transitions. The phrase appears in planetary dynamics, where temporary spin-orbit resonances during planet-planet scattering can excite large obliquities; in circumbinary dynamics, where rapid orbital precession suppresses large obliquity excursions; in compact-object astrophysics, where black-hole spins may be randomized at birth and where numerical relativity must define spin axes despite gauge freedom; in small-body dynamics, where asteroid and tumbling-body spin orientations are measured or rapidly perturbed; and in mathematical juggling, where a “spin layer” extends toss state graphs to rotational configurations (Li, 2021, 2208.00018, Tauris, 2022, Owen et al., 2017, Aljbaae et al., 2021, Bowell et al., 2013, Varpanen, 2014).

1. Terminological scope across research areas

The phrase is not confined to a single field. In exoplanet dynamics it refers to excitation of planetary obliquity during chaotic orbital evolution; in black-hole formation it refers to randomization of a newborn black hole’s spin direction relative to the pre-collapse orbit; in asteroid and spacecraft dynamics it describes rapid, torque-driven reorientation of a rotating body; and in juggling mathematics it names an extension from timing-only state descriptions to spatial and rotational configurations (Li, 2021, Larsen et al., 5 Aug 2025, Aljbaae et al., 2021, Varpanen, 2014).

Domain	Meaning of spin-axis tossing	Representative mechanism or structure
Planet scattering	Planetary obliquity excitation	Temporary spin-orbit resonance
Circumbinary planets	Suppressed large obliquity variation	Binary quadrupole detuning
Binary black holes	Randomized BH spin orientation at formation	Core-collapse tossing
Numerical relativity	Apparent spin-axis motion contaminated by gauge	Spin-direction measure choice
Asteroids and spacecraft	Rapid reorientation under external torque	Earth-induced torques on Apophis
Juggling mathematics	Rotational/spatial state transitions	Spin layer on state graphs

This distribution of usages suggests a shared conceptual core: coupling between intrinsic spin dynamics and an external torque, evolving reference plane, or state-transition structure. A plausible implication is that “spin-axis tossing” functions as a cross-disciplinary descriptor for large or consequential changes in orientation, rather than a single standardized mechanism.

2. Planet-planet scattering and planetary obliquity excitation

In planetary dynamics, spin-axis tossing is developed as a mechanism by which planet-planet scattering can tilt a planet’s spin axis during close encounters and rapid orbital rearrangement (Li, 2021). The central result is that scattering does not merely reshape semi-major axes, eccentricities, and inclinations. It can also produce large spin-axis misalignment through temporary spin-orbit resonances, even in non-collisional events.

The relevant spin-axis precession coefficient is

$$\alpha = \frac{k_2}{2 \mathbb{C} n} \left[ \frac{\omega R_p^3}{\omega_b a^3 (1-e^2)^3} \right]^{1/2},$$

and the secular spin-axis evolution is written through the Hamiltonian

$$H(\xi, \psi, t) = \frac{1}{\alpha} \xi^2 + \sqrt{1-\xi^2}[A(t)\sin\psi + B(t)\cos\psi],$$

with $\xi=\cos\epsilon$. Resonance crossing occurs when the spin precession rate $\alpha$ becomes comparable to the orbital nodal precession rate $\dot{\Omega}$. During scattering, sudden changes in semi-major axis alter $\alpha$, so the system can pass through resonance quickly and receive a large obliquity “kick.”

The process is described as analogous to secular spin-orbit resonances in migration scenarios, but often non-adiabatic because the orbital elements change over short timescales. Tilting is more likely when the planet’s semi-major axis is smaller, its mass is lower, and the system begins with higher $v_{orb}/v_{esc}$. The simulations report that about $80\%$ of non-collisional scattering cases undergo spin-orbit resonance crossing; planets with $a<1\,$AU do so at about $98\%$, compared with about $$H(\xi, \psi, t) = \frac{1}{\alpha} \xi^2 + \sqrt{1-\xi^2}[A(t)\sin\psi + B(t)\cos\psi],$$0 for $$H(\xi, \psi, t) = \frac{1}{\alpha} \xi^2 + \sqrt{1-\xi^2}[A(t)\sin\psi + B(t)\cos\psi],$$1AU. Large tilts $$H(\xi, \psi, t) = \frac{1}{\alpha} \xi^2 + \sqrt{1-\xi^2}[A(t)\sin\psi + B(t)\cos\psi],$$2 occur almost exclusively in systems where resonance crossing happens. The tilted planets retain nearly unchanged spin rates, distinguishing this channel from collisional tilting. Not all systems are equally susceptible: the massive, wide-orbit planet in 2M0122 is cited as unlikely to have acquired its large obliquity in this way.

An important misconception addressed by this framework is the earlier assumption that scattering leaves spin axes largely unaffected in the absence of direct collisions. The scattering calculations instead show that temporary frequency commensurabilities can efficiently transfer angular momentum and excite large obliquity even when inclination remains low.

3. Detuning and suppression in circumbinary planets

A complementary result appears for circumbinary planets: rather than amplifying spin-axis tossing, the stellar binary generally suppresses large obliquity variation by driving rapid orbital nodal precession (2208.00018). The paper identifies the binary quadrupole potential as the key ingredient. Because the close stellar pair acts as a strong quadrupole mass, the planetary orbital node precesses so rapidly that secular spin-orbit resonance is usually detuned.

The nodal precession frequency is given as

$$H(\xi, \psi, t) = \frac{1}{\alpha} \xi^2 + \sqrt{1-\xi^2}[A(t)\sin\psi + B(t)\cos\psi],$$3

while the spin-precession parameter for an oblate planet is

$$H(\xi, \psi, t) = \frac{1}{\alpha} \xi^2 + \sqrt{1-\xi^2}[A(t)\sin\psi + B(t)\cos\psi],$$4

with spin-axis precession rate $$H(\xi, \psi, t) = \frac{1}{\alpha} \xi^2 + \sqrt{1-\xi^2}[A(t)\sin\psi + B(t)\cos\psi],$$5. For typical circumbinary configurations, $$H(\xi, \psi, t) = \frac{1}{\alpha} \xi^2 + \sqrt{1-\xi^2}[A(t)\sin\psi + B(t)\cos\psi],$$6, so the resonance condition is not met. In the canonical secular theory, the coefficient $$H(\xi, \psi, t) = \frac{1}{\alpha} \xi^2 + \sqrt{1-\xi^2}[A(t)\sin\psi + B(t)\cos\psi],$$7, proportional to the binary’s quadrupole, dominates the dynamics, and the spin axis remains close to a nearly constant value.

The reported outcome is low spin-axis variation for observed circumbinary architectures, especially for planets in low-inclination orbits. Secular and numerical simulations find spin-axis variations $$H(\xi, \psi, t) = \frac{1}{\alpha} \xi^2 + \sqrt{1-\xi^2}[A(t)\sin\psi + B(t)\cos\psi],$$8, often much less. Large obliquity variation requires exceptional conditions, such as nearly orthogonal mutual inclinations or extremely tight binaries with periods below $$H(\xi, \psi, t) = \frac{1}{\alpha} \xi^2 + \sqrt{1-\xi^2}[A(t)\sin\psi + B(t)\cos\psi],$$9 day. The paper also notes that adding a massive close-in moon can destabilize circumbinary obliquity by increasing spin-axis precession sufficiently to re-enter the secular resonance regime.

This result corrects a common intuition imported from Solar-system discussions: a massive moon is not the only route to obliquity stabilization. In circumbinary systems, the central binary’s quadrupole can itself provide a stabilizing detuning mechanism.

4. Stellar-wind tilting of Sun-like stars

Spin-axis tossing also appears in the context of stellar obliquity evolution, although in this case the effect is modest rather than chaotic. Stellar winds can tilt the spin axes of Sun-like stars if the mean angular momentum vector carried by the wind is misaligned from the stellar spin axis (Spalding, 2019). The motivating datum is the Sun’s $\xi=\cos\epsilon$0 misalignment from the mean angular momentum plane of the Solar system.

The angular momentum loss is written as

$\xi=\cos\epsilon$1

and the stellar tilt evolution includes

$\xi=\cos\epsilon$2

The associated spin-down time is

$\xi=\cos\epsilon$3

Time variability in the wind offset is modeled as an Ornstein-Uhlenbeck process with autocorrelation time $\xi=\cos\epsilon$4, and the buildup of obliquity depends on the ratio of the early spin-down time to that coherence time.

The model shows that Solar-like tilts of $\xi=\cos\epsilon$5 naturally arise during the first $\xi=\cos\epsilon$6-$\xi=\cos\epsilon$7 Myr after planet formation if the wind deviates by about $\xi=\cos\epsilon$8 from the spin axis and remains coherent over comparable timescales. Most tilting occurs early, when the star is spinning down rapidly. The paper further hypothesizes that the gaseous environments of young open clusters may provide forcing over sufficient timescales to tilt the astrospheres of young stars. At the same time, it states that the more extreme, retrograde stellar obliquities of extrasolar planetary systems likely arise through separate mechanisms.

This suggests a useful distinction within the broader concept: some forms of spin-axis tossing are abrupt and resonance-driven, whereas others are cumulative, stochastic, and torque-mediated over Myr timescales.

5. Black-hole formation, population inference, and spin-axis definition

In compact-object astrophysics, spin-axis tossing refers most directly to randomization of the spin orientation of the second-born black hole during or after core collapse (Tauris, 2022). The observable most strongly affected is the effective inspiral spin,

$\xi=\cos\epsilon$9

or equivalently in the later study,

$\alpha$0

Without tossing, isolated binary evolution predicts closely aligned spins and has difficulty producing the observed negative-$\alpha$1 tail. With tossing, the second-born black hole is drawn from an isotropic distribution,

$\alpha$2

so anti-aligned and weakly aligned configurations become natural.

The 2022 study argues that isolated binary evolution can explain the observed data if black holes have their spin axis tossed during formation, whereas black-hole formation without tossing has difficulty reproducing the data even if alignment prior to the second core collapse is disregarded (Tauris, 2022). The 2025 study sharpens this conclusion with Monte Carlo simulations over a 9-dimensional phase space, empirical kernel density estimates, functional data analysis, and Kolmogorov-Smirnov, Cramer-von Mises, and Anderson-Darling tests (Larsen et al., 5 Aug 2025). For isolated binaries with isotropic tossing, the optimal models achieve p-values up to $\alpha$3 for KS, $\alpha$4 for CvM, and $\alpha$5 for AD. Without tossing, all p-values fall below $\alpha$6 for isolated binaries. A statistically acceptable no-tossing solution emerges only if about $\alpha$7 of detected binary black-hole mergers arise from dynamical channels that randomize spin directions, and for an isolated binary origin the simulations prefer mass reversal in about $\alpha$8 of progenitors.

A separate but related problem concerns how a black-hole spin axis is defined in numerical relativity (Owen et al., 2017). There the difficulty is not birth randomization but gauge freedom: the spin axis is “seriously undermined by the coordinate freedom of general relativity.” Traditional measures based on the Euclidean line between horizon poles or coordinate moments of horizon scalars can show large, unphysical nutations. The recommended alternative is a boost-fixed, centroid-aware quasilocal angular momentum measure,

$\alpha$9

which gives much better agreement with post-Newtonian expectations and removes much of the spurious in-plane nutation. This is relevant to “spin-axis tossing” because apparent rapid changes in the spin direction can be partly diagnostic artifact if the direction measure is contaminated by coordinate, gauge, or tidal effects.

Two controversies are therefore separated by this literature. The first is astrophysical: whether the LVK $\dot{\Omega}$0 distribution requires spin-axis tossing in isolated binaries or instead a large dynamical formation fraction. The second is definitional: whether an inferred spin-axis motion is physical or partly a consequence of how spin direction is measured in numerical relativity.

6. Small bodies: asteroid populations and the Apophis flyby

Small-body studies approach spin-axis tossing through both population statistics and explicitly time-dependent torque modeling. Using the Lowell Observatory photometric database, spin-axis longitudes were estimated for more than $\dot{\Omega}$1 asteroids, a large increase over earlier samples (Bowell et al., 2013). The magnitude method fits brightness as a function of ecliptic longitude, with the peak phase identifying the spin-axis longitude up to a $\dot{\Omega}$2 ambiguity. The main-belt distribution is reported to be clearly non-random, with an excess of longitudes from $\dot{\Omega}$3 to $\dot{\Omega}$4 and a paucity between $\dot{\Omega}$5 and $\dot{\Omega}$6. Jupiter Trojans show a similar pattern, whereas near-Earth asteroids and Mars crossers do not show statistically significant anisotropy. The paper does not provide a definitive explanation and explicitly leaves open whether the pattern reflects observational bias, selection effects, a real physical process, or another mechanism.

For the tumbling asteroid Apophis, the issue is not long-term anisotropy but rapid reorientation during the 2029 close approach to Earth (Aljbaae et al., 2021). Apophis is described as a $\dot{\Omega}$7 m tumbling asteroid passing within $\dot{\Omega}$8 km of Earth’s center on 13 April 2029. Its main rotation period is $\dot{\Omega}$9 h, its secondary wobble period is $\alpha$0 h, and the wobble amplitude is $\alpha$1-$\alpha$2. Earth’s tidal and gravitational torques induce significant and rapid changes in obliquity and precession angle. The spacecraft-dynamics model uses a polyhedral body composed of $\alpha$3 point masses and studies extreme initial spin-axis cases, with the minimum-variation case at $\alpha$4 and the maximum-variation case at $\alpha$5.

The orbiter consequences are severe. Three diagnostic tools are used: MEGNO, PMap, and Time-Series prediction. For circular retrograde equatorial orbits with initial $\alpha$6 and $\alpha$7 km, peak orbital variations in $\alpha$8, $\alpha$9, and $v_{orb}/v_{esc}$0 are tabulated for the two extreme spin-axis scenarios. Nearly all spacecraft orbits become unstable, and about $v_{orb}/v_{esc}$1 do so within $v_{orb}/v_{esc}$2 days. The central conclusion is that no spacecraft with natural orbits could survive the high perturbations caused by the close encounter with Earth.

Taken together, these papers show two distinct small-body meanings of spin-axis tossing: unresolved non-isotropic orientation statistics in asteroid populations, and explicit, fast, externally forced axis motion in a tumbling near-Earth asteroid.

7. Graph-theoretic formalization in juggling

In mathematical juggling, spin-axis tossing is not a torque-driven astrophysical process but an extension of the state approach from toss juggling to spin juggling (Varpanen, 2014). The classical toss-juggling state model represents a pattern by a bijection $v_{orb}/v_{esc}$3 with $v_{orb}/v_{esc}$4, where the throw at beat $v_{orb}/v_{esc}$5 is $v_{orb}/v_{esc}$6. States are $v_{orb}/v_{esc}$7-tuples or $v_{orb}/v_{esc}$8-subsets of $v_{orb}/v_{esc}$9, and the state graph encodes legal throw transitions.

Spin juggling generalizes this particle-centric description by adding spatial and rotational degrees of freedom. The paper introduces a “spin layer” that can be added to toss state graphs, tracking how hand and tether configurations, crossings, and orientations evolve. In the poi state graph, each arrow corresponds to a $80\%$0 rotation per beat; blue arrows denote non-crossing motions and red arrows denote motions that cross the body axis. The encoding yields named patterns such as the ground pattern $80\%$1, the $80\%$2-weave $80\%$3, and the $80\%$4-weave $80\%$5. The framework can also be overlaid with toss states, as in “522 siteswap with 3-weave.”

Unlike the toss-juggling case, where explicit steady-state probabilities for random walks on state graphs are given, the spin-juggling formulation remains primarily graph-theoretic. The paper explicitly identifies rigorous analytical formulas for spin-state transition probabilities or combinatorics as an open area of research. This use of the term therefore emphasizes state-space representation of rotational structure rather than physical reorientation under torque. A plausible implication is that “spin-axis tossing” in this setting names a mathematical language for describing and generating coupled temporal and spatial patterns.

Across these literatures, spin-axis tossing spans excitation, suppression, randomization, measurement, and state representation. The unifying theme is sensitivity of orientation dynamics to external forcing, resonance structure, stochastic accretion, coordinate choice, or graph-constrained transitions. The main points of disagreement are equally clear: whether observed black-hole spin populations require tossing or a dominant dynamical channel; whether asteroid longitude anisotropy is physical or observational; and whether an apparent spin-axis motion in numerical relativity is intrinsic or gauge-contaminated.