Quadrupole Flow: Dynamic Multipolar Phenomena

• Quadrupole flow is a context-dependent descriptor of m=2 angular dynamics found in systems ranging from celestial mechanics and nuclear collisions to atomic interactions and photonic topologies.
• Research shows that quadrupole-driven effects organize potential deformations, angular correlations, and chaotic behaviors, with quantitative insights from Melnikov analysis and flow coefficient measurements.
• Across classical, nuclear, and quantum systems, quadrupole flow underpins phenomena such as symmetry breaking, guided excitation, and transport, offering actionable metrics in both theory and experiments.

In current research usage, quadrupole flow is not a single phenomenon but a family of domain-specific concepts tied to quadrupole moments, quadrupole deformations, or $m=2$ angular structures. In celestial mechanics it denotes Hamiltonian flow in a Kepler problem with a quadrupole perturbation; in high-energy nuclear collisions it denotes the azimuthal $\cos(2\Delta\phi)$ component conventionally associated with $v_2$; in galactic dynamics it denotes gas motion driven by a bar-like quadrupole potential; in atomic, molecular, and Rydberg systems it denotes transport and binding controlled by $1/R^5$ quadrupole-quadrupole interactions; and in topological photonics it denotes light confinement governed by a quantized quadrupole moment $q_{xy}$ (Depetri et al., 2013, Trainor, 2015, Sormani et al., 2015, Han, 2021, Guo et al., 8 May 2025). The shared element is a quadrupolar angular structure, but the relevant dynamical variables, observables, and mechanisms differ sharply across fields.

1. Terminological scope and mathematical motifs

Several literatures use the term through explicit quadrupolar objects. In nuclear-structure and heavy-ion work, the intrinsic nuclear shape is parameterized by the quadrupole deformation $\beta_2$ in a deformed Woods–Saxon radius $R(\theta,\phi)$ (Pandey et al., 2023, Jia et al., 2021). In correlation analyses of nuclear collisions, the quadrupole is the $\cos(2\Delta\phi)$ harmonic of the two-particle angular correlation function, often written through $A_Q=\bar\rho_0 v_2^2$ (Trainor, 2015). In barred-galaxy gas dynamics, the quadrupole is the $m=2$ term in the planar potential,

$\cos(2\Delta\phi)$0

which drives non-circular streaming, shocks, and spiral structure (Sormani et al., 2015). In Rydberg and cold-molecule problems, the interaction scales as $\cos(2\Delta\phi)$1, the characteristic quadrupole-quadrupole law (Han, 2021, Lepers et al., 2010). In photonics, a quadrupole topological insulator is defined by vanishing dipole polarization $\cos(2\Delta\phi)$2 and a quantized bulk quadrupole moment $\cos(2\Delta\phi)$3 in suitable units (Guo et al., 8 May 2025).

These usages are related by symmetry rather than by a common constitutive theory. In some settings the quadrupole organizes a potential; in others it organizes a correlation harmonic, a deformation field, or a topological invariant. A plausible implication is that the phrase is best read locally, with its meaning fixed by the surrounding formalism rather than by a universal definition.

2. Hamiltonian and gravitational quadrupole flow

In the Kepler problem with a quadrupole correction, the Hamiltonian can be written in polar coordinates in the meridional plane as

$\cos(2\Delta\phi)$4

with

$\cos(2\Delta\phi)$5

The sign of $\cos(2\Delta\phi)$6 distinguishes prolate ($\cos(2\Delta\phi)$7) from oblate ($\cos(2\Delta\phi)$8) deformations, while $\cos(2\Delta\phi)$9 supplies the angular dependence that destroys spherical symmetry (Depetri et al., 2013).

The unperturbed zero-energy manifold $v_2$0 corresponds to parabolic Kepler orbits. In that setting the effective total angular momentum is

$v_2$1

and the parabolic orbit can be parametrized by

$v_2$2

That explicit form makes the zero-energy manifold suitable for Melnikov analysis (Depetri et al., 2013).

The key result is that the first Melnikov function vanishes,

$v_2$3

whereas the second Melnikov function has simple zeros whenever a numerically evaluated coefficient $v_2$4. The analysis shows $v_2$5 for $v_2$6, with $v_2$7 only in the limiting cases $v_2$8 and $v_2$9. As a consequence, the stable and unstable manifolds of the parabolic invariant set intersect transversely, implying Smale horseshoes, symbolic dynamics, and chaotic flow on the zero-energy manifold for both prolate and oblate perturbations (Depetri et al., 2013).

This result is notable because it removes a sign-based distinction sometimes inferred from numerics. In this formulation, the sign of $1/R^5$0 rescales the perturbation but does not eliminate the transverse intersections responsible for chaos. The paper therefore identifies the angular structure of $1/R^5$1, not the prolate-versus-oblate sign alone, as the decisive ingredient in the onset of quadrupole-driven nonintegrability (Depetri et al., 2013).

3. Azimuthal quadrupole flow in collider phenomenology

In high-energy nuclear collisions, quadrupole flow most commonly denotes the azimuthal second harmonic conventionally written as $1/R^5$2. A full two-dimensional angular-correlation analysis separates a soft component, a dijet component consisting of a same-side 2D peak and an away-side 1D peak, and a nonjet (NJ) quadrupole, a pure azimuthal $1/R^5$3 term that is not localized in $1/R^5$4. Its amplitude is written

$1/R^5$5

with $1/R^5$6 (Trainor, 2015).

In this framework, soft, dijet, and NJ quadrupole components are distinct. For $1/R^5$7 GeV $1/R^5$8-$1/R^5$9 collisions, the hard yield obeys

$q_{xy}$0

or equivalently $q_{xy}$1. At the level of correlated pair yields, soft pairs scale as $q_{xy}$2, dijet pairs as $q_{xy}$3, and NJ quadrupole pairs as $q_{xy}$4. The NJ quadrupole therefore appears as a third component rather than as an extension of either the soft or hard contribution (Trainor, 2015).

The same body of work argues that conventional $q_{xy}$5 methods mix the NJ quadrupole with jet-related Fourier content from the same-side peak. In the 2D-correlation language of Trainor and related analyses, the extracted nonjet quadrupole is denoted $q_{xy}$6, and the total quadrupole inferred by conventional methods is approximated by

$q_{xy}$7

with $q_{xy}$8 the jet-related contribution from the same-side peak (Trainor et al., 2013, Trainor, 2016). For $q_{xy}$9 and $\beta_2$0 GeV Au–Au collisions, the nonjet quadrupole amplitude is reported to factorize as

$\beta_2$1

with $\beta_2$2 a logarithmic energy factor above about $\beta_2$3 GeV (Kettler et al., 2015).

A second line of analysis reconstructs quadrupole spectra from identified-hadron $\beta_2$4 data. There the inferred quadrupole source has a common monopole boost $\beta_2$5, and in the boost frame the quadrupole spectra of several hadron species collapse onto a universal Lévy distribution. At RHIC the reported parameters are $\beta_2$6 MeV and $\beta_2$7; at $\beta_2$8 TeV Pb–Pb they are $\beta_2$9 MeV and $R(\theta,\phi)$0 (Trainor, 2016, Trainor, 2016, Trainor, 8 Dec 2025). These studies further state that only a small fraction of hadrons, less than about $R(\theta,\phi)$1, carry the NJ quadrupole signal in A–A collisions (Trainor, 2015).

The interpretation of this collider quadrupole is contested. Trainor’s analyses, together with related correlation-based studies, argue that the NJ quadrupole is inconsistent with a hydrodynamic bulk-flow interpretation and instead indicates a distinct nonflow QCD mechanism (Trainor, 2015, Trainor, 2016, Trainor, 2016, Trainor et al., 2013). By contrast, deformation-centered hydrodynamic studies treat $R(\theta,\phi)$2 precisely as a collective response to the initial geometry; that usage is discussed separately below.

A further variant appears in charge-dependent flow studies. In Au+Au collisions at $R(\theta,\phi)$3 GeV, AMPT-based calculations find a dipolar distribution of $R(\theta,\phi)$4 in non-central collisions. Coupling that dipole to the magnetic field $R(\theta,\phi)$5 induces an electric quadrupole moment, which can then lead to a difference in elliptic flows between positive and negative particles, providing an alternative interpretation of the observed charge-dependent pion $R(\theta,\phi)$6 without requiring a chiral magnetic wave (Zhao et al., 2019).

4. Nuclear deformation, flow distributions, and radial-flow probes

A distinct and widely used meaning of quadrupole flow treats flow observables as precision probes of the intrinsic nuclear quadrupole deformation $R(\theta,\phi)$7. In isotopic Xe–Xe collisions at $R(\theta,\phi)$8 TeV, deformation is implemented through a modified Woods–Saxon density,

$R(\theta,\phi)$9

with only $\cos(2\Delta\phi)$0 retained in that study (Pandey et al., 2023). In HYDJET++, elliptic flow is parameterized through spatial and flow-velocity anisotropies, and the resulting $\cos(2\Delta\phi)$1, $\cos(2\Delta\phi)$2, and $\cos(2\Delta\phi)$3 show systematic sensitivity to $\cos(2\Delta\phi)$4 (Pandey et al., 2023).

For the Xe isotopic chain, the reported behavior differs between even-$\cos(2\Delta\phi)$5 and odd-$\cos(2\Delta\phi)$6 nuclei. In even-$\cos(2\Delta\phi)$7 isotopes, $\cos(2\Delta\phi)$8 decreases as $\cos(2\Delta\phi)$9 increases and $A_Q=\bar\rho_0 v_2^2$0 decreases, corresponding to a direct correlation $A_Q=\bar\rho_0 v_2^2$1 in the studied range. In odd-$A_Q=\bar\rho_0 v_2^2$2 isotopes, $A_Q=\bar\rho_0 v_2^2$3 increases as $A_Q=\bar\rho_0 v_2^2$4 increases and $A_Q=\bar\rho_0 v_2^2$5 decreases, corresponding to an inverse correlation in that range, and flow in odd-$A_Q=\bar\rho_0 v_2^2$6 nuclei is suppressed relative to even-$A_Q=\bar\rho_0 v_2^2$7 nuclei (Pandey et al., 2023). The same study reports an approximately linear dependence

$A_Q=\bar\rho_0 v_2^2$8

a positive correlation of $A_Q=\bar\rho_0 v_2^2$9 with deformation, and a negative correlation of $m=2$0 with system size. In body–body and tip–tip orientations,

$m=2$1

which is interpreted as a difference in effective fireball compactness and radial acceleration (Pandey et al., 2023).

Flow-distribution analyses make the same point at the level of event-by-event $m=2$2 statistics. For spherical nuclei in central collisions, the radial distribution $m=2$3 is well approximated by a Bessel–Gaussian. For deformed nuclei, the distribution requires a shifted Gram–Charlier expansion with a non-zero $m=2$4 and higher radial cumulants $m=2$5 and $m=2$6. In simulations of deformed U+U with $m=2$7, the Bessel–Gaussian alone fails, while the shifted radial distribution reproduces the broader and more skewed $m=2$8 distribution (Mehrabpour et al., 2023).

Correlators between anisotropic and radial flow provide another deformation-sensitive observable. Using AMPT and Glauber-based geometry, the Pearson coefficient $m=2$9 between $\cos(2\Delta\phi)$00 and event-wise mean transverse momentum $\cos(2\Delta\phi)$01 is found to be particularly sensitive to quadrupole deformation. Prolate deformation $\cos(2\Delta\phi)$02 reduces $\cos(2\Delta\phi)$03 in ultra-central collisions, whereas oblate deformation $\cos(2\Delta\phi)$04 enhances it. Because $\cos(2\Delta\phi)$05 and $\cos(2\Delta\phi)$06 are the two extremes of triaxiality in the chosen convention, the sign and magnitude of the $\cos(2\Delta\phi)$07–$\cos(2\Delta\phi)$08 correlation can be used to probe triaxiality (Jia et al., 2021).

Recent viscous-hydrodynamic work extends the deformation program to the isotropic sector of the flow. There the radial-flow fluctuation amplitude is defined by

$\cos(2\Delta\phi)$09

and the differential observable $\cos(2\Delta\phi)$10 is constructed from correlations between $\cos(2\Delta\phi)$11 and $\cos(2\Delta\phi)$12. In Trento-3D plus CLVisc calculations, both $\cos(2\Delta\phi)$13 and $\cos(2\Delta\phi)$14 increase with $\cos(2\Delta\phi)$15 in central collisions, while the Pearson coefficient $\cos(2\Delta\phi)$16 exhibits a universal step-like behavior across collision systems and centralities. The same analysis reports that large $\cos(2\Delta\phi)$17 suppresses longitudinal decorrelation of radial flow in central collisions, whereas $\cos(2\Delta\phi)$18 enhances it (Zhu et al., 4 Feb 2026).

Taken together, these deformation-centered studies treat quadrupole flow as a mapping

$\cos(2\Delta\phi)$19

and they regard that mapping as sufficiently systematic to constrain nuclear structure parameters directly from heavy-ion data (Pandey et al., 2023, Mehrabpour et al., 2023, Jia et al., 2021, Zhu et al., 4 Feb 2026).

5. Quadrupole-driven flows in classical fluids and astrophysical gas

In galactic gas dynamics, quadrupole flow refers to gas motion driven by the non-axisymmetric part of a barred gravitational potential. In the Milky Way modeling of Sormani, Binney, and Magorrian, the planar potential is

$\cos(2\Delta\phi)$20

with the quadrupole generated by a three-dimensional density

$\cos(2\Delta\phi)$21

where $\cos(2\Delta\phi)$22 is the bar strength and $\cos(2\Delta\phi)$23 its exponential scale length (Sormani et al., 2015). In 2D isothermal, non-self-gravitating, inviscid simulations, this quadrupole organizes $\cos(2\Delta\phi)$24 and $\cos(2\Delta\phi)$25 orbit families, offset shocks, a central $\cos(2\Delta\phi)$26 disc, and bar-driven spiral arms. Comparison with longitude–velocity data leads to the constraints $\cos(2\Delta\phi)$27 kpc, $\cos(2\Delta\phi)$28, and a preferred pattern speed $\cos(2\Delta\phi)$29, although no single parameter set reproduces all observed features simultaneously (Sormani et al., 2015).

In classical viscous-fluid mechanics, quadrupole flow appears as the standard steady-streaming pattern around a singly oscillating cylinder. For single-frequency forcing, steady streaming is second order in amplitude and produces symmetric quadrupole-like flows with no net pumping. For dual-frequency oscillation,

$\cos(2\Delta\phi)$30

the symmetry is broken: asymmetrical streaming and a non-zero mean flux appear, with the direction set by the polarity of the oscillation. For the case $\cos(2\Delta\phi)$31, the analysis and simulations show that pumping occurs at third order in amplitude; for general rational frequency ratios, pumping requires one integer frequency to be even and the other odd, and the minimum order is $\cos(2\Delta\phi)$32 when the frequencies are $\cos(2\Delta\phi)$33 and $\cos(2\Delta\phi)$34 with $\cos(2\Delta\phi)$35 (Lee et al., 1 Oct 2025).

A related multipolar fluid problem appears in Felderhof’s analysis of a sphere with axisymmetric surface distortions. There the quadrupole is encoded in the vector spherical harmonic $\cos(2\Delta\phi)$36, built from $\cos(2\Delta\phi)$37. In both the dipole–quadrupole swimmer and the quadrupole–octupole swimmer, the quadrupole component shapes the first-order irrotational flow, contributes to the steady vortex-ring structure of the mean second-order flow, and participates in second-harmonic vortex shedding (Felderhof et al., 2018). The quadrupole therefore acts as an organizing mode for both streaming topology and propulsion when coupled to another multipole.

Across these fluid examples, the quadrupole does not merely label an angular pattern. It acts as a symmetry-breaking driver of shocks, steady-streaming cells, vortex shedding, and net transport, depending on how it is coupled to background rotation, additional frequencies, or higher multipoles.

6. Quadrupole-guided transport and confinement in atomic, molecular, and photonic systems

In Rydberg-atom and cold-molecule physics, quadrupole flow is tied to long-range electric quadrupole-quadrupole forces. For two Rydberg atoms, the quadrupole moment is modeled as

$\cos(2\Delta\phi)$38

and the two-body interaction has the form

$\cos(2\Delta\phi)$39

so the interaction scales as $\cos(2\Delta\phi)$40 and is strongly angle dependent (Han, 2021). Because $\cos(2\Delta\phi)$41 for Rydberg states, $\cos(2\Delta\phi)$42, and at fixed separation the maximum repulsive frequency shift obeys

$\cos(2\Delta\phi)$43

(Han, 2021).

For three atoms at the vertices of an equilateral triangle, the Hamiltonian is written with additive pairwise interactions,

$\cos(2\Delta\phi)$44

The resulting spectra show densely packed levels, avoided crossings, attractive and repulsive branches, potential wells, and repulsive peaks of order $\cos(2\Delta\phi)$45 THz or higher at small $\cos(2\Delta\phi)$46. The paper interprets these features as the microscopic basis for quadrupole blockade, bound few-body complexes, and geometry-dependent excitation transport in two dimensions (Han, 2021).

An analogous long-range mechanism appears in cold atom–molecule photoassociation. For a ground-state diatomic molecule and an excited atom, the quadrupole-quadrupole interaction is

$\cos(2\Delta\phi)$47

and for Cs$\cos(2\Delta\phi)$48+Cs(6P) it is attractive enough to bind trimers. The paper therefore identifies photoassociation of a cold atom–molecule pair into a long-range trimer as a promising application, and notes that the same multipolar formalism can be generalized to tetramer formation in molecule–molecule systems (Lepers et al., 2010).

In topological photonics, the term enters through higher-order topology rather than through electrostatics. A photonic quadrupole topological insulator is characterized by vanishing dipole polarization,

$\cos(2\Delta\phi)$49

and a quantized quadrupole moment

$\cos(2\Delta\phi)$50

in the relevant Wannier sector (Guo et al., 8 May 2025). In defect-evolved photonic crystal slabs, clockwise and counterclockwise defect evolution generate two topologically distinct quadrupole phases. Their interface hosts a corner state that functions as a nanocavity. The reported device exhibits simulated $\cos(2\Delta\phi)$51, experimental lasing near $\cos(2\Delta\phi)$52 nm in the telecom C-band, a threshold pump power $\cos(2\Delta\phi)$53, an experimental quality factor $\cos(2\Delta\phi)$54, and wavelength tunability from about $\cos(2\Delta\phi)$55 nm to $\cos(2\Delta\phi)$56 nm (Guo et al., 8 May 2025).

These atomic, molecular, and photonic usages replace the language of hydrodynamic flow with the language of guided excitation, capture, or confinement. In each case, the quadrupole structure creates preferred channels, wells, or corner-localized sinks that organize the motion of atoms, excitations, or light.

Quadrupole flow is therefore best understood as a context-dependent descriptor of quadrupolar dynamics. In some settings it means chaotic Hamiltonian flow near a quadrupole perturbation; in others, an azimuthal second harmonic in collision data, a deformation-sensitive hydrodynamic response, a bar-driven gas pattern, a quadrupole-streaming cell structure, or a transport regime controlled by $\cos(2\Delta\phi)$57 interactions or by a quantized quadrupole invariant. The recurring feature is not a single mechanism but the dynamical role of quadrupole symmetry.