---
title: 'Quadrupole Flow: Dynamic Multipolar Phenomena'
url: https://www.emergentmind.com/topics/quadrupole-flow
type: topic
---

# Quadrupole Flow: Dynamic Multipolar Phenomena

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}\) [1306.6021], [1512.01857], [1507.03078], [2106.01479], [2505.05320]. 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)\) [2312.05853], [2105.05713]. 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\) [1512.01857]. In barred-galaxy gas dynamics, the quadrupole is the \(m=2\) term in the planar potential,
\[
\Phi(R,\phi)=\Phi_0(R)+\Phi_2(R)\cos(2\phi),
\]
which drives non-circular streaming, shocks, and spiral structure [1507.03078]. In Rydberg and cold-molecule problems, the interaction scales as \(1/R^5\), the characteristic quadrupole-quadrupole law [2106.01479], [1006.5314]. In photonics, a quadrupole topological insulator is defined by vanishing dipole polarization \(\mathbf P=(0,0)\) and a quantized bulk quadrupole moment \(q_{xy}=1/2\) in suitable units [2505.05320].

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
\[
H=H_0+qW_1(R)+qW_2(R,\theta),
\]
with
\[
W_1(R)=-\frac{1}{R^3}, \qquad W_2(R,\theta)=\frac{3\cos^2\theta}{2R^3}.
\]
The sign of \(q\) distinguishes prolate (\(q>0\)) from oblate (\(q<0\)) deformations, while \(W_2\) supplies the angular dependence that destroys spherical symmetry [1306.6021].

The unperturbed zero-energy manifold \(H_0=0\) corresponds to parabolic Kepler orbits. In that setting the effective total angular momentum is
\[
G^2=R^4\dot\theta^2+\frac{L_z^2}{\cos^2\theta},
\]
and the parabolic orbit can be parametrized by
\[
R(\theta)=R_{\min}\,\sec^2\!\left[\frac{1}{4A}\ln\left|\frac{A+\sin\theta}{A-\sin\theta}\right|\right],
\qquad
A=\sqrt{1-\frac{L_z^2}{G^2}}, \quad 0<A\le 1.
\]
That explicit form makes the zero-energy manifold suitable for Melnikov analysis [1306.6021].

The key result is that the first Melnikov function vanishes,
\[
M_1(\theta_0)=0,
\]
whereas the second Melnikov function has simple zeros whenever a numerically evaluated coefficient \(I_2\neq 0\). The analysis shows \(I_2\neq 0\) for \(0<A<1\), with \(I_2=0\) only in the limiting cases \(A=0\) and \(A=1\). 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 [1306.6021].

This result is notable because it removes a sign-based distinction sometimes inferred from numerics. In this formulation, the sign of \(q\) rescales the perturbation but does not eliminate the transverse intersections responsible for chaos. The paper therefore identifies the angular structure of \(W_2\), not the prolate-versus-oblate sign alone, as the decisive ingredient in the onset of quadrupole-driven nonintegrability [1306.6021].

## 3. Azimuthal quadrupole flow in collider phenomenology

In high-energy nuclear collisions, quadrupole flow most commonly denotes the azimuthal second harmonic conventionally written as \(v_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 \(\cos(2\Delta\phi)\) term that is not localized in \(\Delta\eta\). Its amplitude is written
\[
A_Q=\bar\rho_0\,v_2^2,
\]
with \(\bar\rho_0=n_{ch}/\Delta\eta\) [1512.01857].

In this framework, soft, dijet, and NJ quadrupole components are distinct. For \(200\) GeV \(p\)-\(p\) collisions, the hard yield obeys
\[
n_h \approx 0.01\,n_s^2 \qquad (\Delta\eta=1),
\]
or equivalently \(\bar\rho_h\propto \bar\rho_s^2\). At the level of correlated pair yields, soft pairs scale as \(\bar\rho_s\), dijet pairs as \(\bar\rho_s^2\), and NJ quadrupole pairs as \(\bar\rho_s^3\). The NJ quadrupole therefore appears as a third component rather than as an extension of either the soft or hard contribution [1512.01857].

The same body of work argues that conventional \(v_2\) 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 \(v_2\{2D\}\), and the total quadrupole inferred by conventional methods is approximated by
\[
A_Q\{2\}\approx A_Q\{2D\}+A_Q\{SS\},
\]
with \(A_Q\{SS\}\) the jet-related contribution from the same-side peak [1302.0300], [1610.06256]. For \(62\) and \(200\) GeV Au–Au collisions, the nonjet quadrupole amplitude is reported to factorize as
\[
A_Q\{2D\}(b,\sqrt{s_{NN}})=C\,R(\sqrt{s_{NN}})\,N_{bin}(b)\,\epsilon_{2,\mathrm{opt}}^2(b),
\]
with \(R(\sqrt{s_{NN}})\) a logarithmic energy factor above about \(13.5\) GeV [1504.02741].

A second line of analysis reconstructs **quadrupole spectra** from identified-hadron \(v_2(p_t)\) data. There the inferred quadrupole source has a common monopole boost \(\Delta y_{t0}\approx 0.6\), 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 \(T_2\approx 92\) MeV and \(n_2=14\); at \(2.76\) TeV Pb–Pb they are \(T_2=94\) MeV and \(n_2=12\) [1609.07693], [1610.06256], [2512.08003]. These studies further state that only a small fraction of hadrons, less than about \(5\%\), carry the NJ quadrupole signal in A–A collisions [1512.01857].

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 [1512.01857], [1609.07693], [1610.06256], [1302.0300]. By contrast, deformation-centered hydrodynamic studies treat \(v_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 \(\sqrt{s_{NN}}=200\) GeV, AMPT-based calculations find a dipolar distribution of \(\mathbf E\cdot\mathbf B\) in non-central collisions. Coupling that dipole to the magnetic field \(\mathbf B\) 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 \(v_2\) without requiring a chiral magnetic wave [1901.04156].

## 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 \(\beta_2\). In isotopic Xe–Xe collisions at \(\sqrt{s_{NN}}=5.44\) TeV, deformation is implemented through a modified Woods–Saxon density,
\[
\rho(r,\theta,\phi)=\frac{\rho_0}{1+\exp[(r-R'(\theta,\phi))/a]},
\qquad
R'(\theta,\phi)=R\big[1+\beta_2 Y_2^0(\theta,\phi)+\beta_3Y_3^0+\beta_4Y_4^0+\dots\big],
\]
with only \(\beta_2\) retained in that study [2312.05853]. In HYDJET++, elliptic flow is parameterized through spatial and flow-velocity anisotropies, and the resulting \(v_2\), \(v_3\), and \(\langle p_T\rangle\) show systematic sensitivity to \(\beta_2\) [2312.05853].

For the Xe isotopic chain, the reported behavior differs between even-\(A\) and odd-\(A\) nuclei. In even-\(A\) isotopes, \(v_n\) decreases as \(A\) increases and \(\beta_2\) decreases, corresponding to a direct correlation \(v_n\propto\beta_2\) in the studied range. In odd-\(A\) isotopes, \(v_n\) increases as \(A\) increases and \(\beta_2\) decreases, corresponding to an inverse correlation in that range, and flow in odd-\(A\) nuclei is suppressed relative to even-\(A\) nuclei [2312.05853]. The same study reports an approximately linear dependence
\[
\frac{(v_2^2)_{\rm Xe}}{(v_2^2)_{\rm Pb}} \approx a_2+b_2\,\beta_2^2,
\]
a positive correlation of \(\langle p_T\rangle\) with deformation, and a negative correlation of \(\langle p_T\rangle\) with system size. In body–body and tip–tip orientations,
\[
\langle p_T\rangle_{\text{tip-tip}}>\langle p_T\rangle_{\text{min-bias}}>\langle p_T\rangle_{\text{body-body}},
\]
which is interpreted as a difference in effective fireball compactness and radial acceleration [2312.05853].

Flow-distribution analyses make the same point at the level of event-by-event \(v_2\) statistics. For spherical nuclei in central collisions, the radial distribution \(P_r(v_n)\) is well approximated by a Bessel–Gaussian. For deformed nuclei, the distribution requires a shifted Gram–Charlier expansion with a non-zero \(\bar v_n\) and higher radial cumulants \(R_n\{4\}\) and \(R_n\{6\}\). In simulations of deformed U+U with \(\beta_2=0.265\), the Bessel–Gaussian alone fails, while the shifted radial distribution reproduces the broader and more skewed \(v_2\) distribution [2301.07770].

Correlators between anisotropic and radial flow provide another deformation-sensitive observable. Using AMPT and Glauber-based geometry, the Pearson coefficient \(\rho_2\) between \(v_2^2\) and event-wise mean transverse momentum \([p_T]\) is found to be particularly sensitive to quadrupole deformation. Prolate deformation \(\beta>0\) reduces \(\rho_2\) in ultra-central collisions, whereas oblate deformation \(\beta<0\) enhances it. Because \(\beta>0\) and \(\beta<0\) are the two extremes of triaxiality in the chosen convention, the sign and magnitude of the \(v_2^2\)–\([p_T]\) correlation can be used to probe triaxiality [2105.05713].

Recent viscous-hydrodynamic work extends the deformation program to the isotropic sector of the flow. There the radial-flow fluctuation amplitude is defined by
\[
v_0 \equiv \frac{\sigma_{[p_T]}}{\langle [p_T]\rangle},
\]
and the differential observable \(v_0(p_T)\) is constructed from correlations between \(\delta n(p_T)\) and \(\delta [p_T]\). In Trento-3D plus CLVisc calculations, both \(v_0\) and \(v_0(p_T)\) increase with \(\beta_2\) in central collisions, while the Pearson coefficient \(\rho(n(p_T),[p_T])\) exhibits a universal step-like behavior across collision systems and centralities. The same analysis reports that large \(\beta_2\) suppresses longitudinal decorrelation of radial flow in central collisions, whereas \(\beta_4\) enhances it [2602.04148].

Taken together, these deformation-centered studies treat quadrupole flow as a mapping
\[
\beta_2 \;\to\; \text{initial geometry} \;\to\; v_2,\; v_3,\; \langle p_T\rangle,\; P(v_2),\; \rho_2,\; v_0(p_T),
\]
and they regard that mapping as sufficiently systematic to constrain nuclear structure parameters directly from heavy-ion data [2312.05853], [2301.07770], [2105.05713], [2602.04148].

## 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
\[
\Phi(R,\phi)=\Phi_0(R)+\Phi_2(R)\cos(2\phi),
\]
with the quadrupole generated by a three-dimensional density
\[
\rho_2(r,\phi,\theta)=\frac{K A}{L^2}\exp\!\left(-\frac{2r}{L}\right)\sin^2\theta\cos(2\phi),
\]
where \(A\) is the bar strength and \(L\) its exponential scale length [1507.03078]. In 2D isothermal, non-self-gravitating, inviscid simulations, this quadrupole organizes \(x_1\) and \(x_2\) orbit families, offset shocks, a central \(x_2\) disc, and bar-driven spiral arms. Comparison with longitude–velocity data leads to the constraints \(L\gtrsim1.5\) kpc, \(A\gtrsim0.4\), and a preferred pattern speed \(\Omega_p\simeq40\,\mathrm{km\,s^{-1}\,kpc^{-1}}\), although no single parameter set reproduces all observed features simultaneously [1507.03078].

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,
\[
X(t)=\frac{A}{2}\big[\sin(\Omega t)+\sin(\alpha\Omega t)\big],
\]
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 \(\alpha=2\), 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 \(a+b\) when the frequencies are \(a\) and \(b\) with \(\gcd(a,b)=1\) [2510.01344].

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 \(\mathbf B_2(\theta)\), built from \(P_2(\cos\theta)\). 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 [1803.11037]. 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
\[
Q=e r^2,
\]
and the two-body interaction has the form
\[
V_{12}=\frac{Q_1Q_2}{4\pi\epsilon_0 R^5}\times \text{(angular/tensor structure)},
\]
so the interaction scales as \(1/R^5\) and is strongly angle dependent [2106.01479]. Because \(r\sim n^2\) for Rydberg states, \(Q\sim e n^4\), and at fixed separation the maximum repulsive frequency shift obeys
\[
\Delta\nu_{\max}(R=0.5\,\mu\mathrm m)\propto n^8
\]
[2106.01479].

For three atoms at the vertices of an equilateral triangle, the Hamiltonian is written with additive pairwise interactions,
\[
V_{qq}=V_{12}+V_{23}+V_{31}.
\]
The resulting spectra show densely packed levels, avoided crossings, attractive and repulsive branches, potential wells, and repulsive peaks of order \(\sim100\) THz or higher at small \(R\). The paper interprets these features as the microscopic basis for quadrupole blockade, bound few-body complexes, and geometry-dependent excitation transport in two dimensions [2106.01479].

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
\[
\hat V^{qq}_{AB}(R)=\frac{24}{R^5}\sum_{M=-2}^{2}\frac{\hat Q_2^M(\hat r_A)\,\hat Q_2^{-M}(\hat r_B)}{(2+M)!(2-M)!},
\]
and for Cs\(_2\)+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 [1006.5314].

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,
\[
\mathbf P=(P_x,P_y)=(0,0),
\]
and a quantized quadrupole moment
\[
q_{xy}=2p_x p_y=\frac12
\]
in the relevant Wannier sector [2505.05320]. 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 \(Q\gtrsim10^4\), experimental lasing near \(1567\) nm in the telecom C-band, a threshold pump power \(P_{\rm th}\approx0.5\,\mu\mathrm W\), an experimental quality factor \(Q_{\rm exp}\approx2300\), and wavelength tunability from about \(1519\) nm to \(1543\) nm [2505.05320].

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 \(1/R^5\) interactions or by a quantized quadrupole invariant. The recurring feature is not a single mechanism but the dynamical role of quadrupole symmetry.

Source: https://www.emergentmind.com/topics/quadrupole-flow