The Helix Nebula (NGC 7293) is one of the closest and brightest planetary nebulae, and therefore a benchmark for resolving how the late-stage ejecta of stars couple to their surroundings. Using the Gaia EDR3/DR3 astrometric solution for the central star (WD 2226-210), we adopt a distance of d = 198.6 (+1.6/−1.8) pc (ref. ¹⁰). Imaging of the Helix has shown it to be highly complex. Its bright main nebula comprises an inner disc and a surrounding outer torus, embedded within a larger structure whose upstream side is truncated, consistent with interaction between the expanding asymptotic giant branch (AGB) ejecta and the ambient interstellar medium (ISM)^11,12,13. The ionized nebula is threaded by thousands of dense cometary knots and associated molecular material, indicating that much of the ejected material remains in a clumpy, only partially processed phase^11,14. Deep imaging and spectroscopy have also revealed a bow-shock feature in the faint outer halo, in the direction of the nebula’s motion through the local ISM^12,15. Together, these properties make the Helix uniquely suited to place direct constraints on how fragmented stellar ejecta are dispersed and mixed into the ISM: an essential step in the recycling of mass, dust and newly synthesized elements in galaxies^4,9.

Thank you for reading this post, don't forget to subscribe!

The Helix Nebula was observed in the light of Hα, [N ii] and [O iii] with the partially built Modular Optical Telephoto Hyperspectral Robotic Array (MOTHRA). MOTHRA is an array of high-end telephoto lenses equipped with tiltable ultra-narrow interference filters, located at the El Sauce Observatory in Chile. Its design evolved from the Dragonfly Spectral Line Mapper at New Mexico Skies Observatory^16,17. When completed, MOTHRA will have 1,140 lenses distributed over 30 mounts, and be optically equivalent to a 4.8 m f/0.08 refractor. The data described here are equivalent to roughly 20 minutes of on-source exposure time with the completed array.

The MOTHRA Hα image of the Helix is shown in Fig. 1. It shows many features that have not been seen in ionized gas before, such as extensions of the plumes in the northwest and southeast^11,13 and turbulent and complex low surface brightness Hα emission in the southwest¹³. The most striking feature in the Hα image is a forest of arcs and partial arcs on the eastern side of the nebula. We identify at least 22 arcs on the eastern side, labelled 1–22 in Fig. 3 in order of increasing distance from the central white dwarf. Although most of the features appear to be new discoveries, several can be seen in previous GALEX and Hα images. Besides the large and complex feature 14 this includes arcs 3, 10 and 13, among others^11,13,15. The arcs are undetected in [O iii] and faint in [N ii]; from the brightest region of arc 14 we measure [N ii]/Hα = 0.06 ± 0.01 and [O iii]/Hα < 0.015 (2σ). We also find faint arc-like features on the western side, at a similar distance from the white dwarf as the much brighter ones in the east.

Following earlier studies we interpret the eastern features as bow shocks in which expanding stellar ejecta encounter the ISM at supersonic velocities^13,15. The velocity of the Helix with respect to the ISM is roughly 36 km s⁻¹ towards 95° east of north in the plane of the sky and roughly 27 km s⁻¹ along the line of sight^15,18, for a combined ISM velocity of v_ISM ≈ 45 km s⁻¹ with respect to the systemic velocity of the nebula. On the eastern side the shock velocity is the sum of v_ISM and the expansion velocity of the ejecta, v_{shock, e} ≈ |v_ISM + v_exp|. On the western side the ejecta encounter a turbulent postshock wake that has passed through, and mixed with, the planetary nebula, v_{shock, w} ≈ |v_wake − v_exp|. The wake velocity is expected to be small with respect to the systemic velocity of the nebula. In the case of the well-studied AGB star Mira, the processed gas being shed from the bow shock into the immediate downstream tail lags the star by only roughly 10% (ref. ¹⁹); applying the same scaling to the Helix gives v_wake = 0−10 km s⁻¹ for the flow on the western side.

Using the MAPPINGS V code²⁰ we derive shock velocities on the eastern side of v_{shock, e} = 80−90 km s⁻¹ from the [O ii]/Hα and [N ii]/Hα line ratios (Methods). Subtracting the 45 km s⁻¹ ISM velocity gives an expansion velocity of the ejecta of v_exp = 35−45 km s⁻¹, consistent with previously measured Hα kinematics in the region of the brightest bow¹⁵. The implied velocity field around the Helix is shown in Fig. 2. Near the east–west axis the shock velocities are roughly 80 km s⁻¹ in the east and roughly 35 km s⁻¹ in the west. For these velocities, shock models predict Hα luminosities that are 1 to 2 orders of magnitude fainter in the west than in the east, consistent with the appearance of the bows on the two sides of the Helix (Methods).

The expansion velocity of v_exp = 35–45 km s⁻¹ indicates a dynamical age of the clumps of 20,000–30,000 years at r ≈ 1 pc, predating the formation of the planetary nebula roughly 12,000 years ago²¹. The material therefore most likely belongs to an older circumstellar envelope that was ejected during the late-AGB phase, and has now fragmented into many individual clumps. This interpretation is strengthened by the fact that the bows lie at approximately the same radius as the roughly 40′ outer WISE 12-μm halo, which has been associated with dust from an AGB wind¹³. Fast winds and outflows, with velocities that can exceed the canonical roughly 5−20 km s⁻¹ range of the main AGB phase²², are commonly observed during the late-AGB and early post-AGB phases^23,24. These flows are often bipolar rather than isotropic²⁵, which may explain why there appears to be a preferred axis connecting the strong bows in the east to northeast of the Helix to the weak bows in the west to southwest.

We fit the bow morphologies with a range of functional forms: parabolas, hyperbolas, ellipses and Wilkinoids²⁶. We adopt the parabolic fits as our fiducial model, as they provide a reasonable description of the data and retain the same functional form under projection²⁷. As many of the structures do not show a complete bow, and some are broader and less sharply bounded than ideal thin-shell bow shocks, we do not attach direct dynamical meaning to the adopted fit family itself. Instead, we use the fits to extract geometric quantities (in particular the characteristic curvature scale), and use the variation between fit families as an estimate of the systematic uncertainty. The fitting procedure is detailed in the Methods, with the results shown in Fig. 3. Most of the 22 bows are reasonably well fit by a parabola, although the wings are often closer to hyperbolic²⁷. The brightest bow (14) is more sharply peaked than the model curve; inspection of the region near the apex shows that it is broken up in a complex network of smaller shocks.

The foci of the parabolic fits are indicated with red dots in Fig. 3. These are the approximate locations of the shell fragments that produce the shocks, although the exact location depends on the three-dimensional orientation and shape of the bows²⁷. There is generally nothing visible in Hα, [O iii] or [N ii] at or near these locations. The lack of detected emission counterparts at most of the inferred obstacle locations is consistent with the fragments being largely neutral. Our observations thus represent a new observational window on the mixing of AGB and planetary nebula ejecta into the ISM, in which the fragments are identified not by their intrinsic emission but by the shocks that they drive.

The nature of the bows changes systematically with distance from the white dwarf: close to the star they are large, thin and well-defined, whereas in the outskirts they are smaller and fuzzier. To quantify this trend without assuming a particular steady-state bow-shock solution, we characterize each structure by the radius of curvature at its apex R_c, as measured from the best-fitting parabola. The relation between R_c and the distance from the white dwarf, r, is shown in Fig. 4. The characteristic size of the bows decreases by two orders of magnitude over the radial range 0.4 pc ≲ r ≲ 1.4 pc. A log-linear fit gives log R_c = 0.34−1.59r, corresponding to an e-folding length of 0.27 pc.

Because R_c is a purely geometric quantity, it does not by itself specify the detailed momentum balance within the flow. It does, however, show that the spatial scale of the coherent bow-forming obstacle decreases strongly with outward distance. This geometric trend is accompanied by a systematic morphological transition: the inner bows are thin and sharply bounded, whereas the outer structures are broader, more irregular, and increasingly clumpy. Taken together, these changes suggest progressive stripping and fragmentation of the dense AGB-shell remnants as they interact with the ambient medium. As material is ablated from the fragments and mixed into the surrounding flow, the surviving dense heads become smaller and more porous, and a larger fraction of the Hα emission probably arises in fragment-associated, mass-loaded mixed gas rather than in a geometrically thin, well-defined forward shock^8,28,29.

In this interpretation, the bows are both signposts of the fragments and agents of their destruction. The shocks are powered by the relative kinetic energy of the fragments and the ambient flow, and their presence implies ongoing momentum transfer, ablation and mixing^8,28. Interpreting the radial locations of the bows, r_i, as a time sequence for outward evolution, t_i ≈ r_i/v_exp, the slope of the R_c–r relation implies a characteristic timescale for the loss of coherent bow structure. For the expansion velocity derived from the shock velocities in the east, v_exp ≈ 40 km s⁻¹, the curvature scale declines with an e-folding time τ_Rc ≈ 7 × 10³ years. We therefore infer that the bow-forming AGB-shell fragments are disrupted on a timescale of order 10⁴ years. This estimate should be interpreted as the survival time of the coherent dense fragment–bow system, rather than as a direct measurement of a specific momentum or mass-loss rate.

Classical AGB–ISM bow shocks such as Mira trace a single, wind-driven stand-off interaction centred on the mass-losing star⁶, and far-infrared surveys show that such global wind–ISM interaction structures are common around evolved stars⁷. However, these long-lived⁶ (roughly 10⁵ years) large-scale bows primarily map where the wind meets the ISM; they do not directly constrain how fragmented ejecta are ultimately assimilated. In the Helix outer halo we instead resolve numerous compact bow shocks with no luminous source at their foci, indicating that the AGB–ISM interaction has fragmented into dense, line-dark obstacles that are progressively ablated and entrained. The roughly 10⁴ years characteristic disruption time for the fragments can therefore be interpreted as the relevant timescale for recycling of late-AGB ejecta into the ISM.

More broadly, stellar mass loss is a major channel by which galaxies recycle gas, metals and dust back into the ISM³⁰, yet the efficiency and duration of the final, fragment-driven AGB–ISM assimilation step remain poorly constrained observationally^6,31. Galaxy formation simulations therefore rely on subgrid turbulent mixing and diffusion prescriptions to represent unresolved transport^32,33,34. Our empirically inferred roughly 10⁴ years disruption time implies that once AGB ejecta are fragmented and exposed to the diffuse medium they lose their coherent identity rapidly, providing a benchmark for models of recycling and feedback.

Our conclusions can be tested and extended in several ways. The Helix is a fairly typical planetary nebula, and we should see similar fragment-driven bow shocks in the outskirts of other planetary nebulae. Because of the steep relation between Hα luminosity and shock velocity such features will be most readily detected when the bulk motion of the nebula with respect to the ISM exceeds roughly 40 km s⁻¹. When more examples are found, it will be interesting to see if the mixing timescale depends on the shock velocities. Such observations will be within easy reach of the completed MOTHRA. Furthermore, by analogy with the molecular cometary knots in the Helix^14,35,36,37, the Hα-dark clumps whose presence is inferred from the shocks could be detectable in CO rotational lines and, particularly, in the H₂ 1–0 S(1) 2.12-μm line.