---
title: Spiral Wave-Driven Ejecta in Neutron-Star Mergers
url: https://www.emergentmind.com/topics/spiral-wave-driven-ejecta
type: topic
---

# Spiral Wave-Driven Ejecta in Neutron-Star Mergers

Spiral wave-driven ejecta denotes outflow whose unbinding is directly tied to spiral structure, most explicitly the “spiral-wave wind” or “spiral-wave-driven outflow” identified in long-lived neutron-star merger remnants, where non-axisymmetric spiral density waves launched by the surviving remnant into the disk carry angular momentum outward and drive a secular mass loss [2310.09934]. In broader astrophysical usage, closely related phenomena include spiral-organized, spiral-patterned, or spiral-shell ejecta in which spiral geometry is important but the immediate driver is instead jet pressure, reconnection, Parker-spiral current-sheet reconnection, or binary orbital motion; the literature therefore distinguishes direct spiral-wave launching from systems in which spiral structure primarily organizes or textures the ejecta [1907.04872][2501.09663][2503.11380][2104.06207][1309.4980][2308.15677].

## 1. Conceptual scope and taxonomy

The most precise use of the term applies to neutron-star mergers that leave a long-lived massive neutron star (MNS). In that setting, the remnant remains strongly non-axisymmetric, excites spiral density waves in the surrounding disk, and the waves “carry angular momentum outwards and drive the mass outflow” [2310.09934]. This is physically distinct from prompt tidal or shock dynamical ejecta, from later viscous secular outflows, and from high-latitude neutrino-driven winds.

A broader but qualified usage appears in several other subfields. In these cases, the spiral pattern is observationally or dynamically important, but the papers do not claim that a classical spiral density wave is itself the sole launching agent.

| Setting | Immediate driver | Role of spiral structure |
|---|---|---|
| Long-lived neutron-star merger remnant | Spiral density waves from remnant torques | Direct launching of a spiral-wave wind [2310.09934] |
| Blue-kilonova interpretation | Spiral density waves plus shock reprocessing | Direct hydrodynamic mass-ejection channel [1907.04872] |
| Milky Way \(^{26}\)Al ejecta | Forward superbubble blow-out | Spiral arms organize source offsets [1309.4980] |
| Common-envelope NS jets | Jet-inflated high-pressure regions and RT instability | Orbit-imprinted spiral ejecta morphology [2501.09663] |
| Solar spiral jet | External reconnection and twist transfer | Spiral jet with torsional Alfvén-wave-like propagation [2503.11380] |
| Parker-spiral current sheet | Pressure-curvature loss of equilibrium and reconnection | Spiral magnetic topology organizes plasmoid site [2104.06207] |
| Binary circumstellar shells | Orbital motion and hydrodynamic compression | Spiral-shell pattern records anisotropic ejecta kinematics [2308.15677] |

This taxonomy suggests that “spiral wave-driven ejecta” is exact in merger-disk hydrodynamics, but only approximate or metaphorical in several other environments.

## 2. Spiral-wave wind in long-lived neutron-star merger remnants

In the merger context, the spiral-wave-driven outflow is a post-merger ejecta channel generated by non-axisymmetric spiral density waves launched by the surviving neutron-star remnant into the disk [2310.09934]. The ejecta history shows an early pulse of unbound material of about \({\sim}3\times 10^{-3}\ M_\odot\) ejected during and shortly after merger, followed by a quasi-steady wind with \(\dot{M}\simeq 0.05\ M_\odot\ {\rm s}^{-1}\), which is stated to be “dominated by the spiral-wave wind.” The disk itself forms on a timescale of about \({\sim}20\ {\rm ms}\), direct mass ejection from the central remnant terminates after roughly \({\sim}20{-}30\ {\rm ms}\), and the spiral pattern forms shortly after merger and persists for the entire duration of the \({\sim}100\ {\rm ms}\) simulations.

The remnant origin of the outflow is diagnosed through persistent low-order azimuthal structure. The mode amplitudes are measured through
$$
C_m = \int_{z=0} e^{-{\rm i}\, m\, \varphi} \rho\, W\, \sqrt{\gamma}\, {\rm d}^2 x ,
$$
with \(m=1,2,3\) dominant, while the \(m=2\) deformation of the remnant, identified as the bar mode, generates the spiral-wave in the disk [2310.09934]. The associated disk-plane density contrast,
$$
\frac{\rho - \langle \rho \rangle}{\langle \rho \rangle},
$$
reveals a central \(m=2\) deformation and connected spiral arms extending through the disk. In this picture, the long-lived remnant is not a passive gravitating center but a continuously forcing non-axisymmetric source.

A related numerical-relativity study described the same mechanism as a new hydrodynamics-driven mass-ejection channel for the blue kilonova: spiral density waves launched by the nonaxisymmetric post-merger remnant propagate into the surrounding disk, transport angular momentum and energy outward, and eject matter from the outer disk layers [1907.04872]. That work characterized the spiral-wave wind by a mass \({\sim}10^{-2}~M_{\odot}\), velocity \({\sim}0.2c\), and electron fraction mostly distributed above \(0.25\). Its angular-momentum transport was diagnosed with
$$
J_r = T_{ra} (\partial_\phi)^a ,
$$
which showed strong outward angular momentum transport associated with the spiral pattern.

## 3. Geometry, composition, and nucleosynthesis

The spatial structure of the merger ejecta is anisotropic. The spiral-wave-driven outflow is predominantly located close to the disk orbital plane and has a broad distribution of electron fractions, while at higher latitudes a high electron-fraction wind is driven by neutrino radiation [2310.09934]. The angular distribution is summarized by the statement that “Most of the ejecta is equatorial and has a broad range of electron fractions. The polar ejecta is less massive and characterized by significantly larger electron fraction \(Y_e \gtrsim 0.4\).” In the gray M1 simulations, the ejecta have a broader distribution, with \(Y_e \lesssim 0.5\), whereas earlier M0 simulations produced \(Y_e<0.4\) at all latitudes.

Nucleosynthetically, the 2023 long-lived-MNS calculation concluded that the spiral-wave wind produces a full r-process yield and that the combined ejecta are in good agreement with Solar abundance measurements, with the main discrepancy around \(A\sim 140\) attributed to the adopted nuclear-mass model rather than the ejecta model itself [2310.09934]. By contrast, the 2019 blue-kilonova interpretation associated the spiral-wave wind with \(Y_e>0.25\), partial hydrodynamic shock reprocessing in the expanding arms, and predominantly production up to the second peak, \(A<130\), while the combination of dynamical ejecta and spiral-wave wind reproduced solar-system \(r\)-process abundances [1907.04872]. This suggests sensitivity to remnant properties, microphysics, and transport treatment.

A third merger study focused on very light elements and strontium in dynamical and spiral-wave wind ejecta [2009.08988]. There, spiral-wave wind ejection rates were reported as \(R_{\rm wind} \approx (1.57 \pm 0.42)\times 10^{-1}\,M_\odot\,{\rm s}^{-1}\) for BLh\_equal and \(R_{\rm wind} \approx (1.58 \pm 0.04)\times 10^{-1}\,M_\odot\,{\rm s}^{-1}\) for DD2\_equal, while the integrated wind strontium mass fraction was about \(1.5\%-4.5\%\). Hydrogen was negligible in the spiral-wave wind, helium was of order \(10^{-5}\) by mass, and strontium was synthesized efficiently enough that a long-lived remnant could produce an amount consistent with, or even larger than, that inferred for GW170817. The paper therefore argued that the remnant in GW170817 likely did not remain in the spiral-wave-wind-producing phase for much longer than about \(100\ {\rm ms}\).

## 4. Diagnostics and simulation frameworks

The spiral-wave wind has been studied in fully three-dimensional relativistic merger calculations. One implementation used 3D general-relativistic neutrino-radiation hydrodynamics with the gray moment-based M1 code THC\_M1, the Z4c spacetime formulation, and 7 AMR levels, with finest spacing \(h = 0.125\ G M_\odot/c^2 \simeq 185\ {\rm m}\); magnetic fields were neglected [2310.09934]. Ejecta were identified on a sphere at coordinate radius \(R = 300\ G M_\odot/c^2 \simeq 443\ {\rm km}\), including only material with outgoing radial velocity and unbound according to the Bernoulli criterion
$$
-h u_t \geq h_{\rm min}.
$$
The remnant and disk were separated by a density threshold \(\rho = 10^{13}\ {\rm g}\ {\rm cm}^{-3}\).

The earlier spiral-wave-wind kilonova calculations used full \(3+1\) numerical relativity with WhiskyTHC, microphysical equations of state, approximate neutrino transport, and GRLES turbulent viscosity [1907.04872]. Equal-mass binaries with \(M=(1.364+1.364)\,M_\odot\) were evolved with HS(DD2) and LS220, through merger and remnant phases for at least \(30\ {\rm ms}\), up to \(100\ {\rm ms}\) depending on the model. Seven refinement levels were used, with finest spacings \(246\), \(185\), and \(123\ {\rm m}\), and ejecta were measured on coordinate spheres at \(r=294\ {\rm km}\). Dynamical ejecta were computed with the geodesic criterion, whereas wind ejecta were computed with the Bernoulli criterion.

The light-element and strontium study used the same broad merger framework—WhiskyTHC, 3+1 Z4c evolution, finite-temperature composition-dependent EOSs, approximate neutrino transport with absorption in optically thin conditions, and GRLES viscosity—and extended the model set to BLh\_equal, BLh\_unequal, and DD2\_equal binaries targeted to GW170817 [2009.08988]. For the wind component, the asymptotic speed was defined as
$$
v_{\infty} = c \left( 1-\gamma^{-2} \right)^{1/2},
$$
with \(\gamma = u_t\) for geodetic extraction and \(\gamma = u_t h\) for Bernoulli extraction. The nucleosynthesis trajectories were evolved with SkyNet.

## 5. Related spiral-organized ejecta in other astrophysical settings

Outside merger disks, spiral structure often remains central to ejecta morphology or kinematics, but the causal chain differs. In the Milky Way \(^{26}\)Al study, radial velocities vary systematically with Galactic longitude, reaching approximately \(+\SI{200}{\km\per\s}\) in the first quadrant and \(-\SI{200}{\km\per\s}\) in the fourth, with the average bulk velocity being \({\sim}200\ {\rm km\,s^{-1}}\) larger than the Galactic-rotation prediction; the preferred explanation is spiral-arm source geometry plus a forward superbubble blow-out with \(v_\mathrm{bo} = \SI{225 \pm 50}{\km\per\s}\), not direct acceleration by a spiral density wave [1309.4980]. In that interpretation, spiral structure organizes where massive stars form relative to dense gas, and the offset toward the arm leading edge shapes the direction of ejecta blow-out.

In jetting common-envelope evolution, three-dimensional hydrodynamical simulations of a neutron star spiraling inside a rotating red-supergiant envelope found that envelope rotation leads to more prominent spiral structures of the ejecta than in the non-rotating case, while the ejecta themselves are accelerated by jet-inflated high-pressure volumes and fragmented by Rayleigh–Taylor instability [2501.09663]. The most obvious difference among models with \(\alpha=0\), \(0.25\), and \(0.5\) is that the faster the rotation is, the more pronounced the spiral structure in planes parallel to the orbital plane is. The paper explicitly states that the jet-inflated high-pressure volumes around the neutron star accelerate the envelope, a process prone to Rayleigh–Taylor instability.

In solar-coronal spiral jets, a data-constrained 3D MHD model found that a tiny spiral jet formed when an unstable mini-flux rope erupted in a fan–spine magnetic topology, the pre-existing magnetic null collapsed into a curved 3D current sheet, and external reconnection transferred both twist and cool material onto the outer spine [2503.11380]. The twist then propagated upward at about \(230\ {\rm km\,s^{-1}}\), close to a local estimated Alfvén speed of about \(300\ {\rm km\,s^{-1}}\), which the authors interpreted as consistent with a nonlinear torsional Alfvén wave. Here again, the spiral character is real, but the eruption is primarily reconnection-driven.

A laboratory Parker-spiral current-sheet analogue provided a different near-analogue: quasi-periodic plasmoids form at the tip of a streamer structure when a high-\(\beta\) region undergoes pressure-curvature-driven dynamic loss of equilibrium, and a simple heuristic model extrapolated to solar streamers yielded plasmoid periods of about \(90\) minutes at \(3R_\odot\) [2104.06207]. The “spiral” in that case is the Parker-spiral background magnetic topology rather than a spiral wave that directly launches ejecta.

Binary circumstellar outflows offer a kinematic analogue. Reanalysis of 3D hydrodynamic models for spiral-shell patterns induced by the orbital motion of a mass-losing post-main-sequence star found that the variation of the transverse wind velocity is as large as half the average wind velocity over the entire three dimensional domain in the simulated models [2308.15677]. The paper argued that the often-adopted isotropic wind assumption and the binary hypothesis as the underlying origin for circumstellar multilayered shells are mutually incompatible, so the spiral-shell pattern should be read as a direction-dependent record of ejecta kinematics.

## 6. Dependencies, caveats, and conceptual limits

Within neutron-star mergers, the prominence of spiral-wave-driven ejecta depends on remnant survival. The mechanism is a hallmark of long-lived MNS remnants and is strongly suppressed if the remnant collapses promptly to a black hole, because the non-axisymmetric stellar core that launches the waves disappears [2310.09934]. The calculations also remain limited to the early secular phase, typically the first \({\sim}100\ {\rm ms}\), so quoted wind rates characterize that interval rather than the full seconds-long evolution. Additional caveats include omitted magnetic fields in the THC\_M1 calculation, approximate neutrino transport in both the M1 and earlier leakage+M0 schemes, model dependence on EOS, compactness, total mass, and mass ratio, and the explicit caution that the Cartesian grid may artificially enhance the \(m=4\) mode [2310.09934][1907.04872][2009.08988].

Across the broader literature, the phrase itself requires care. The \(^{26}\)Al paper does not show spiral density waves directly accelerating ejecta by \(200\ {\rm km\,s^{-1}}\); it supports a mediated claim in which spiral/bar geometry and source offsets create preferential superbubble breakout [1309.4980]. The common-envelope jet simulations do not identify a coherent spiral wave mode; they instead show jet-driven, pressure-driven ejecta with a spiral morphology imposed by inspiral and enhanced by envelope rotation [2501.09663]. The Parker-spiral current-sheet work attributes ejecta to pressure buildup, bad magnetic curvature, and reconnection/plasmoid formation rather than to a propagating spiral wave [2104.06207]. The circumstellar-shell analysis shows that even when a spiral pattern is the dominant observable, pattern propagation in the sky plane should not automatically be interpreted as a unique bulk ejecta expansion speed [2308.15677]. In the solar-jet case, the torsional Alfvén-wave-like disturbance is embedded within a reconnection-driven eruption rather than replacing it as the primary trigger [2503.11380].

Taken together, these studies establish two linked meanings of the subject. In the strict sense, spiral wave-driven ejecta refers to the equatorially concentrated, broad-\(Y_e\), post-merger wind unbound when persistent non-axisymmetric spiral density waves torque a merger disk [2310.09934][1907.04872]. In the broader comparative sense, spiral structure can also organize where ejecta are launched, how they are collimated, and how their kinematics are recorded, even when the immediate driver is jet pressure, reconnection, or orbital-motion-induced hydrodynamic shaping.

Source: https://www.emergentmind.com/topics/spiral-wave-driven-ejecta