Spin-Induced Quadrupole Moments

• Spin-induced quadrupole moments are multipole deformations from spin that quantify deviations from spherical symmetry in compact objects.
• The κ parameter distinguishes Kerr black holes (κ = 1) from non-Kerr objects, affecting waveform phasing and binary parameter estimation.
• Advanced detectors and hierarchical methods are sharpening constraints on these deformations, aiding the separation of neutron stars, boson stars, and black holes.

Searching arXiv for recent and foundational papers on spin-induced quadrupole moments to ground the article. Spin-induced quadrupole moments are higher-multipole deformations generated by spin. In gravitational-wave theory for compact binaries, the standard parametrization is

$Q_i=-\kappa_i\,m_i^3\,\chi_i^2,$

with $m_i$ the mass, $\chi_i\equiv |S_i|/m_i^2$ the dimensionless spin, and $\kappa_i$ a dimensionless quadrupole parameter; for a Kerr black hole, general relativity fixes $\kappa_i=1$ exactly, while $\kappa_i\neq 1$ indicates non-Kerr structure (Saleem et al., 2021). Closely related forms, such as $Q_i=-\kappa_i S_i^2/m_i$, are also used (Krishnendu et al., 2018). The same phrase also appears in operator-based settings, where quadrupole tensors quantify anisotropic spin distributions in composite hadrons, nuclei, magnets, and quantum spin liquids (Krutov et al., 2018, Miura et al., 2022, Remund et al., 2022, Takahashi et al., 24 Mar 2026). Across these contexts, spin-induced quadrupole moments encode how spin degrees of freedom depart from spherical symmetry.

1. Kerr normalization, deformation parameters, and object classification

For compact objects, the central benchmark is the Kerr relation implied by the no-hair theorem: a Kerr black hole has $\kappa=1$, independent of mass or spin (Saleem et al., 2021). In the Hansen-Geroch-Thorne multipole language, the Kerr multipoles satisfy

$\mathcal{M}_\ell+i\mathcal{S}_\ell=M(ia)^\ell,$

so that the mass quadrupole is $Q_{\rm Kerr}=-a^2M=-\chi^2 M^3$ (Kong et al., 2024). A convenient deviation parameter is $m_i$0, defined by

$m_i$1

with $m_i$2 recovering the Kerr result (Kong et al., 2024).

The expected values of $m_i$3 are object-dependent. Neutron stars have $m_i$4 ranging roughly from $m_i$5 to $m_i$6, depending on the equation of state and rotation rate (Saleem et al., 2021). Boson stars can have $m_i$7, while gravastars can sit near $m_i$8 with $m_i$9 deviations in thin-shell models; some gravastar constructions even permit negative $\chi_i\equiv |S_i|/m_i^2$0 (Krishnendu et al., 2018, Saleem et al., 2021). This object dependence is what makes the parameter a direct probe of compact-object nature.

Binary analyses often replace $\chi_i\equiv |S_i|/m_i^2$1 by symmetric and antisymmetric combinations,

$\chi_i\equiv |S_i|/m_i^2$2

or, equivalently, $\chi_i\equiv |S_i|/m_i^2$3 and $\chi_i\equiv |S_i|/m_i^2$4 defined relative to the Kerr value (Krishnendu et al., 2018, Krishnendu et al., 2019). In the equal-deformation limit $\chi_i\equiv |S_i|/m_i^2$5, only $\chi_i\equiv |S_i|/m_i^2$6 survives (Krishnendu et al., 2018). This simplification is widely used because measuring the individual coefficients is strongly degenerate with other intrinsic parameters (Krishnendu et al., 2019).

A common misconception is that non-Kerr structure necessarily corresponds to a large positive quadrupole deformation. The literature does not support that simplification. Some models predict very large positive $\chi_i\equiv |S_i|/m_i^2$7, some allow negative $\chi_i\equiv |S_i|/m_i^2$8, and fully back-reacted string-theory constructions can even produce a positive four-dimensional mass quadrupole $\chi_i\equiv |S_i|/m_i^2$9, opposite in sign to Kerr (Saleem et al., 2021, Bena et al., 6 Oct 2025).

2. Post-Newtonian entry, precession dynamics, and waveform imprint

In the post-Newtonian expansion for compact binaries, the leading correction to the conservative dynamics from the rotation of each component appears at $\kappa_i$0PN order and is entirely captured by the second mass moment (Saleem et al., 2021). In frequency-domain inspiral waveforms, the spin-induced quadrupole enters through the phasing. One representative expression is

$\kappa_i$1

plus symmetric combinations of the two bodies’ contributions, with $\kappa_i$2 and $\kappa_i$3 (Saleem et al., 2021). In the aligned-spin decomposition used for third-generation forecasts, the waveform phase contains a $\kappa_i$4PN spin-spin term carrying the quadrupole dependence, while higher spin-spin corrections enter at $\kappa_i$5PN and the leading cubic-in-spin octupole term at $\kappa_i$6PN (Krishnendu et al., 2018).

The same dependence can be expressed in terms of $\kappa_i$7 and $\kappa_i$8. At $\kappa_i$9PN,

$\kappa_i=1$0

so analyses that set $\kappa_i=1$1 measure only the symmetric combination (Krishnendu et al., 2018). In TaylorF2-like parametrizations for large quadrupolar deviations, the $\kappa_i=1$2PN coefficient $\kappa_i=1$3 explicitly contains $\kappa_i=1$4 and $\kappa_i=1$5, while additional quadrupole terms appear in the $\kappa_i=1$6PN and $\kappa_i=1$7PN coefficients $\kappa_i=1$8 and $\kappa_i=1$9 (Chia et al., 2022).

For generic precessing binaries, spin-induced quadrupole moments also modify the precession dynamics. The precession-frequency vectors $\kappa_i\neq 1$0 include spin-orbit, spin-spin, and quadrupole-curvature couplings up to $\kappa_i\neq 1$1PN, and each object’s own $\kappa_i\neq 1$2 appears multiplying its spin $\kappa_i\neq 1$3 inside the precession equations (Lyu et al., 2023). This changes the Euler angles used in twisting-up procedures and therefore modulates higher-mode amplitudes and phases beyond the non-precessing $\kappa_i\neq 1$4 carrier (Lyu et al., 2023). The same work distinguishes the aligned-spin phase correction (“AI” effect) from the precession-induced amplitude and phase modulation (“PI” effect), and shows that for highly precessing systems about $\kappa_i\neq 1$5 of random orientations yield PI-effect mismatches larger than the AI effect (Lyu et al., 2023).

The extended-body formulation gives the same physics in a different language. In the Mathisson-Papapetrou-Dixon system, a spin-induced quadrupole is introduced through a quadrupole tensor $\kappa_i\neq 1$6 built from the STF part of the spin bilinear, with a coefficient $\kappa_i\neq 1$7 equal to unity for a Kerr black hole and of order $\kappa_i\neq 1$8 for neutron-star equations of state (Han et al., 2016). This formulation is used for orbital dynamics in Kerr spacetime rather than directly for comparable-mass waveform inference, but it encodes the same idea: spin induces a quadrupolar deformation that couples to curvature.

3. Current-detector inference: single events, hierarchical methods, and population tests

With second-generation detectors, single-event measurements are generally weak. The $\kappa_i\neq 1$9 parameter is weakly constrained unless the event has moderately high spins, $Q_i=-\kappa_i S_i^2/m_i$0, and good inspiral SNR, $Q_i=-\kappa_i S_i^2/m_i$1; otherwise $Q_i=-\kappa_i S_i^2/m_i$2 is effectively unmeasured, which dilutes population inference (Saleem et al., 2021). In dedicated simulations for Advanced LIGO-Virgo design sensitivity, $Q_i=-\kappa_i S_i^2/m_i$3 credible intervals remain very broad for low-spin systems, tighten to $Q_i=-\kappa_i S_i^2/m_i$4 for moderate effective spin, and to $Q_i=-\kappa_i S_i^2/m_i$5 for high $Q_i=-\kappa_i S_i^2/m_i$6 (Krishnendu et al., 2019).

A Bayesian single-event framework measures the symmetric deformation parameter under the assumption $Q_i=-\kappa_i S_i^2/m_i$7, inserts the inspiral corrections into IMRPhenomPv2, and truncates the waveform to avoid unmodeled post-inspiral effects (Krishnendu et al., 2019). Applied to GW151226 and GW170608, the posteriors remain consistent with the binary-black-hole hypothesis. The reported $Q_i=-\kappa_i S_i^2/m_i$8 intervals and Bayes factors were:

Event	$Q_i=-\kappa_i S_i^2/m_i$9 interval on $\kappa=1$0	$\kappa=1$1
GW151226	$\kappa=1$2	$\kappa=1$3
GW170608	$\kappa=1$4	$\kappa=1$5

These results place $\kappa=1$6 well inside the inferred intervals and do not favor the non-BH model (Krishnendu et al., 2019).

Population analyses were introduced to overcome the weakness of single-event constraints. One approach posits that the true $\kappa=1$7 values follow a Gaussian hyperdistribution $\kappa=1$8, and combines event-level posteriors through

$\kappa=1$9

A narrow $\mathcal{M}_\ell+i\mathcal{S}_\ell=M(ia)^\ell,$0, $\mathcal{M}_\ell+i\mathcal{S}_\ell=M(ia)^\ell,$1 is consistent with an all-Kerr population (Saleem et al., 2021). The second approach is a hierarchical mixture likelihood in which each event is either a BBH with $\mathcal{M}_\ell+i\mathcal{S}_\ell=M(ia)^\ell,$2 exactly or a non-BBH with $\mathcal{M}_\ell+i\mathcal{S}_\ell=M(ia)^\ell,$3, with mixture fraction $\mathcal{M}_\ell+i\mathcal{S}_\ell=M(ia)^\ell,$4. The single-event likelihood is

$\mathcal{M}_\ell+i\mathcal{S}_\ell=M(ia)^\ell,$5

and the $\mathcal{M}_\ell+i\mathcal{S}_\ell=M(ia)^\ell,$6-event likelihood is the product over events (Saleem et al., 2021).

Both methods are equally effective to hint at inhomogeneous populations, but the mixture-likelihood approach was found to be more natural for mixture populations comprising compact objects of diverse classes (Saleem et al., 2021). On simulated high-quality injections satisfying $\mathcal{M}_\ell+i\mathcal{S}_\ell=M(ia)^\ell,$7 and inspiral SNR $\mathcal{M}_\ell+i\mathcal{S}_\ell=M(ia)^\ell,$8, the mixture method recovered $\mathcal{M}_\ell+i\mathcal{S}_\ell=M(ia)^\ell,$9 for a pure non-BH set and $Q_{\rm Kerr}=-a^2M=-\chi^2 M^3$0 for a $Q_{\rm Kerr}=-a^2M=-\chi^2 M^3$1 mixture (Saleem et al., 2021). Applied to the LIGO-Virgo detections published in GWTC-2, the posterior on $Q_{\rm Kerr}=-a^2M=-\chi^2 M^3$2 peaks at zero, $Q_{\rm Kerr}=-a^2M=-\chi^2 M^3$3 is strongly disfavored, and the hierarchical-Gaussian analysis likewise finds $Q_{\rm Kerr}=-a^2M=-\chi^2 M^3$4 and $Q_{\rm Kerr}=-a^2M=-\chi^2 M^3$5, all consistent with a pure BH population within present statistical precision (Saleem et al., 2021).

4. Systematics, degeneracies, and reduced waveform constructions

The main limitations of present-day inference are not only statistical. Strong degeneracies exist between $Q_{\rm Kerr}=-a^2M=-\chi^2 M^3$6, effective spin, and mass ratio, and overlap calculations show that for a given injected $Q_{\rm Kerr}=-a^2M=-\chi^2 M^3$7 there can exist black-hole waveform parameters $Q_{\rm Kerr}=-a^2M=-\chi^2 M^3$8 yielding overlaps above $Q_{\rm Kerr}=-a^2M=-\chi^2 M^3$9 (Krishnendu et al., 2019). This explains skewed posteriors and the slow growth of evidence for non-Kerr structure in second-generation data.

Prior choices are also important. In one population study the prior range on $m_i$00 in single-event runs was uniform in $m_i$01; unless the true non-BBH $m_i$02 lie in similarly wide ranges, Occam’s penalty disfavors non-BH interpretations (Saleem et al., 2021). The same analysis imposed $m_i$03, i.e. $m_i$04, so events with mismatched $m_i$05 would be mis-modeled and could pull $m_i$06 downward (Saleem et al., 2021). Selection biases as a function of $m_i$07 were not included there and were identified as a necessary ingredient for a fully astrophysical analysis (Saleem et al., 2021).

Waveform-model dependence introduces an additional systematic. Re-analyses with IMRPhenomXPHM plus the new precession module showed that IMRPhenomPv2 can artificially tighten $m_i$08 bounds when its internal phase differences relative to the true signal are larger than the true-signal phase shifts due to $m_i$09 (Lyu et al., 2023). In those injected-signal studies, mismatch-versus-$m_i$10 curves confirmed that IMRPhenomPv2 overestimates constraining power (Lyu et al., 2023).

Large deviations from Kerr require dedicated search models. A six-dimensional post-Newtonian waveform with $m_i$11 dependence was reduced to five- and four-dimensional models by reparameterizing the dominant quadrupole terms into effective parameters and truncating subdominant $m_i$12-dependent contributions (Chia et al., 2022). The same work showed that large positive $m_i$13 lowers the minimum-binding-energy frequency substantially, providing a natural inspiral cutoff; the analytic approximation to that cutoff reproduces the true PN ISCO to about $m_i$14 even for $m_i$15 up to $m_i$16 (Chia et al., 2022). In injection studies drawn from $m_i$17, $m_i$18, and $m_i$19, both reduced models achieved mean effectualness $m_i$20, more than $m_i$21 of injections had $m_i$22, and even the worst cases remained above $m_i$23; by contrast, a standard BBH template bank had mean $m_i$24 and failed for $m_i$25 (Chia et al., 2022).

5. Measurement prospects with third-generation and space-borne detectors

Third-generation ground-based detectors substantially improve the measurement of spin-induced quadrupoles. For a single optimally oriented binary at $m_i$26 Mpc, mass ratio $m_i$27, $m_i$28, and $m_i$29, forecasts for Cosmic Explorer and Einstein Telescope ET-D give

$m_i$30

across total masses $m_i$31, with the tightest bound reaching $m_i$32 (Krishnendu et al., 2018). The error improves for larger total mass, higher mass ratio, and spins aligned with $m_i$33; relative to Advanced LIGO, the error shrinks by about $m_i$34 for CE (Krishnendu et al., 2018). For simulated populations, roughly $m_i$35 of detected binaries satisfy $m_i$36, and about $m_i$37 satisfy $m_i$38 (Krishnendu et al., 2018).

Precession adds further information for mass-gap systems. In representative A#-network studies at $m_i$39 Mpc, a GW190814-like configuration with $m_i$40, $m_i$41, and injected $m_i$42 yielded

$m_i$43

for A#, while Cosmic Explorer gave

$m_i$44

The same analysis found that, in a high-SNR O4 injection with moderate precession, including the PI effect shrinks $m_i$45 upper limits by about $m_i$46 relative to AI-only recovery (Lyu et al., 2023).

Space-borne detectors target a different regime. For a canonical massive black hole binary with $m_i$47, $m_i$48, true $m_i$49, and $m_i$50 Gpc, TianQin forecasts a total one-year SNR of about $m_i$51, with the final day alone carrying $m_i$52 of that SNR (Kong et al., 2024). The corresponding Fisher errors are $m_i$53 for a full year or one month of data and $m_i$54 for one day, indicating constraints at the $m_i$55 level (Kong et al., 2024). The optimal mass range is around $m_i$56, larger mass ratios improve the constraints, and including higher modes tightens $m_i$57 by roughly a factor of $m_i$58 in the Bayesian study (Kong et al., 2024).

These forecasted accuracies imply progressively different tests. Third-generation ground networks can directly probe the inspiral quadrupole coefficient $m_i$59 at the level needed to separate Kerr binaries from neutron-star or boson-star mimickers in favorable events (Krishnendu et al., 2018). Space-borne observations of massive black hole inspirals can test the Kerr relation $m_i$60 to about $m_i$61 and use Bayes factors to distinguish black-hole and boson-star-like injections decisively when higher modes are included (Kong et al., 2024).

6. Extended-body dynamics, sign structure, and broader spin-system realizations

In extreme-mass-ratio motion, spin-induced quadrupoles can be incorporated directly in the Mathisson-Papapetrou-Dixon equations. Numerical studies of circular and generic Kerr orbits found that the orbital-frequency shift scales as $m_i$62, peaks at an intermediate radius, and is only about $m_i$63 for typical EMRI parameters $m_i$64, $m_i$65, and $m_i$66 (Han et al., 2016). The same analysis concluded that one may safely neglect spin-induced quadrupoles in EMRI waveform templates for space-based detectors, although periodic variations in the kinematical spin $m_i$67 at the level of $m_i$68 could in principle be relevant for pulsar timing around intermediate-mass black holes (Han et al., 2016).

The sign of the quadrupole is not universal. In the Kerr family one has $m_i$69, or equivalently $m_i$70 (Bena et al., 6 Oct 2025). By contrast, singularity-free Running-Kerr-Taub-Bolt solutions descending from eleven-dimensional supergravity can have strictly positive $m_i$71 for a range of charges. At $m_i$72, the explicit formula given for $m_i$73 is positive, and the sign changes as the bolt’s running speed increases (Bena et al., 6 Oct 2025). The paper interprets positive $m_i$74 as elongation along the spin axis rather than Kerr-like pancaking (Bena et al., 6 Oct 2025). A plausible implication is that future measurements sensitive to the sign structure of the spin-induced quadrupole would test not only the magnitude of non-Kerr deviations but also the class of microphysics supporting the compact object.

Outside relativistic compact-object astrophysics, spin-induced quadrupole moments appear as operator-valued observables. In spin-1 magnets, the on-site traceless symmetric tensor

$m_i$75

quantifies quadrupolar or spin-nematic order, and the classical limit of spin-1 is not a single $m_i$76 vector but a point on $m_i$77 carrying both dipole and quadrupole components (Remund et al., 2022). In magnetic materials with transition metals, the quadrupole tensor of the spin density contributes to magnetocrystalline anisotropy through the spin-flip term

$m_i$78

with $m_i$79 the intra-atomic magnetic dipole or spin-quadrupole moment (Miura et al., 2022). In relativistic two-body $m_i$80-wave spin-1 composites, a nonzero static quadrupole moment arises even when $m_i$81, because Wigner (Melosh) rotations break the non-relativistic cancellation (Krutov et al., 2018). In the Kitaev spin liquid with spin-$m_i$82 impurities, the impurity quadrupole moment exhibits discontinuous jumps at flux-sector transitions and thus acts as a local probe of the underlying $m_i$83 gauge sector (Takahashi et al., 24 Mar 2026).

Across these realizations, the shared structure is precise but not identical. In compact-binary relativity, spin-induced quadrupole moments are usually effective multipole coefficients in orbital dynamics and waveform phasing. In composite and many-body spin systems, they are often expectation values of traceless rank-2 operators that diagnose anisotropic spin structure. The common theme is that spin generates quadrupolar information beyond the dipole level; the operational meaning depends on the dynamical framework in which the quadrupole is defined.