Pith. sign in

REVIEW 4 major objections 7 minor 49 references

Numerous bow shocks in the outer Helix Nebula

T0 review · 4 major / 7 minor · reviewed 2026-08-15 · deepseek-v4-flash

Pith's one-line read Deep H-alpha imaging of the Helix Nebula's outer halo reveals twenty-two compact bow shocks driven by clumps of AGB ejecta, whose rapid outward shrinkage places the disruption of those clumps at roughly 10,000 years.

desk verdict Genuinely new imaging shows a striking factor-of-100 decline in bow-shock curvature across the Helix halo; the ~10^4 yr disruption time is plausible but hinges on a single-velocity age mapping that needs more care. read the letter →

arxiv 2608.11443 v1 pith:VCPKYV4Q submitted 2026-08-11 astro-ph.SR astro-ph.GA

classification astro-ph.SRastro-ph.GA
keywords planetarynebulaebowshocksasymptoticgiantbranchstarsinterstellarmediumrecyclingstellarmasslossHelixNebulaH-alphaimagingshock-cloudinteraction
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper reports the detection of twenty-two compact bow shocks in the eastern outer halo of the Helix Nebula, seen in H-alpha light with a partially completed array of telephoto lenses. The shocks are small, tied to individual clumps of gas rather than to the central star, and they become smaller and fuzzier with distance from the star. The authors read the outward sequence as a time sequence: the clumps are fragments of an older asymptotic-giant-branch shell, the envelope shed near the end of the star's life, and the shocks both mark the fragments and erode them as they plow through the surrounding interstellar medium. If this reading is correct, the e-folding decline of the bows' size implies that such fragments lose coherence and are mixed into the ISM on a timescale of about 10,000 years, giving galaxy-recycling models a direct empirical number to work with.

What carries the argument

The central object is the compact bow shock itself: a small H-alpha arc with no luminous source at its focus, interpreted as the working surface of a dense neutral fragment moving supersonically through the interstellar medium. The quantitative argument is carried by a single geometric measurement. Each bow is fitted with a parabola, chosen as the fiducial among four analytic fit families because parabolas keep their form under projection, and the apex radius of curvature R_c is read off the fit as a proxy for the size of the coherent fragment producing the shock. The paper deliberately attaches no dynamical significance to the fit family; it uses R_c only as a geometric scale, with the spread among fit families treated as systematic uncertainty. The relation log R_c = 0.34 - 1.59 r, with r the distance from the white dwarf, turns into a timescale because r is mapped to age through t = r / v_exp: the e-folding length of 0.27 pc becomes an e-folding time of about 7,000 years at v_exp = 40 km/s.

What would settle it

Measure radial velocities of the bow apexes or the fragments inferred at their foci across the eastern halo. If the outward motions deviate from the assumed ~40 km/s expansion field, or if fragments at nearly the same projected radius show a wide spread in R_c, then the R_c-r gradient is not a single age sequence and the inferred ~10,000-year disruption timescale is not uniquely determined.

Watch

Extended reading notes

Core claim

The paper's central claim is that the outer eastern halo of the Helix Nebula contains at least twenty-two compact H-alpha bow shocks, each produced not by the star's large-scale wind but by an individual, largely neutral fragment of an older shell ejected in the late asymptotic-giant-branch phase. Between 0.4 and 1.4 pc from the central white dwarf, the apex radius of curvature of these bows falls by a factor of roughly 100, and the bow morphology shifts from thin, well-defined arcs to broad, clumpy patches. The paper interprets this as progressive stripping and fragmentation: the same bow shock that makes a fragment visible is also the agent that ablates it and mixes its material into the flow. Reading each bow's radial distance as an age through a constant expansion velocity of about 40 km/s, the observed log-linear size decline corresponds to an e-folding time of about 7,000 years, which the paper rounds to a disruption timescale on the order of 10,000 years for fragmented stellar ejecta entering the ISM.

Load-bearing premise

The 10,000-year number rests on assuming that each bow's distance from the central star is a reliable clock: all fragments were launched at roughly the same time and have travelled outward at one constant speed of about 40 km/s, so that radius measures age and the shrinking curve measures destruction.

Editorial extensions

If this is right

  • If the reading is right, the outer Helix is not a smooth wind-ISM boundary but a resolved population of individual fragments being destroyed in real time, so the final step of AGB mass recycling is observable directly.
  • The ~10,000-year coherence time is much shorter than the ~100,000-year lifetimes of classical star-centered bow shocks, so models of mass return to the ISM should treat fragmented ejecta as short-lived once exposed to the diffuse medium.
  • The clumps' inferred dynamical ages of 20,000-30,000 years at about 1 pc predate the planetary nebula itself (~12,000 years), meaning the bows trace the late-AGB circumstellar envelope rather than the nebula, and similar bows should appear around other fast-moving planetary nebulae.
  • Galaxy-formation simulations that currently rely on subgrid diffusion and mixing prescriptions for unresolved transport would have an empirical benchmark: exposed AGB fragments lose coherent identity in about 10,000 years.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: if the ~10,000-year coherence time is typical, then each AGB star returns its envelope to the ISM as a burst of short-lived dense parcels; combining the timescale with a mass-loss history would give a per-star fragment mass-return rate that chemical enrichment models could test.
  • Editorial inference: the west-side bows are a built-in control, since the same shock-velocity scaling predicts they should be 1-2 orders of magnitude fainter in H-alpha; observing them in a slower-shock tracer would test whether the disruption timescale depends on encounter speed or is set by the fragments themselves.
  • Editorial inference: the R_c-r slope should vary with the local ISM density if erosion is environmental; mapping the bows against the density structure of the eastern halo would show whether the e-folding length is a universal property of the fragments or a response to their surroundings.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 7 minor

Summary. The paper reports the detection of 22 compact bow shocks in the eastern outskirts of the Helix Nebula in deep Hα imaging with the partially completed MOTHRA telescope. The bows are argued to be shocks driven by dense, largely neutral fragments of AGB-shell ejecta as they move through the ISM. The authors fit parabolas to the bows, measure the radius of curvature R_c at the apex, and find a strong decrease of R_c with projected distance r from the central white dwarf (log R_c = 0.34 − 1.59 r, e-folding length 0.27 pc), accompanied by a morphological transition from thin, well-defined bows to fuzzy, patchy features. Interpreting r_i/v_exp as an age sequence with v_exp ≈ 40 km/s, they infer an e-folding disruption timescale for the fragments of τ_Rc ≈ 7×10^3 yr, i.e., of order 10^4 yr, and propose this as a direct constraint on the recycling timescale of fragmented stellar ejecta into the ISM.

Significance. The empirical detection of a forest of compact bow shocks in the outer halo of a PN, with no detected emission at the foci, is a striking and novel observational result. If the inferred ~10^4 yr disruption timescale holds, it would provide a valuable benchmark for models of AGB mass-loss recycling and turbulent mixing in galaxy simulations. The paper is commendably transparent about its assumptions, makes the reduced data publicly available, and tests the east-west brightness asymmetry with radiative shock models. However, the central quantitative claim—the disruption timescale—is currently supported only under a specific kinematic interpretation of the R_c–r relation, and alternative explanations based on velocity or density gradients are not ruled out. The significance is therefore conditional on the additional modeling and measurements requested below.

major comments (4)
  1. [Main text, 'Interpreting the radial locations...'] The mapping r_i → t_i = r_i/v_exp requires a common ejection epoch and a single constant expansion velocity for all 22 fragments. This assumption is in tension with the paper's own discussion of bipolar outflows: the text states that late-AGB 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/northeast of the Helix to the weak bows in the west/southwest.' If the outer fragments were ejected with somewhat higher velocities along that axis, their larger r would reflect velocity rather than age, and the e-folding of R_c would no longer measure a disruption time. Please (i) provide a quantitative argument or model that justifies treating the radial sequence as a single-age sequence, or (ii) test the assumption directly, e.g., by measuring proper motions of the bows (the Gaia epoch and MOTHRA epoch could provide a baseline for the brightest arcs) or by restricting the fit to bows that share a common position angle and comparing the inferred e-folding length.
  2. [Methods, 'Shock velocity' and main text 'The expansion velocity...'] The expansion velocity v_exp = 35–45 km/s is derived from MAPPINGS V shock models fit to line ratios measured only for the brightest region of arc 14, and is then applied to all 22 bows. The inferred shock velocity 80–90 km/s depends on the assumed magnetic-field parameter α (α=1 vs α=2; see Extended Data Fig. 2), so the systematic uncertainty in v_exp is not fully captured by the quoted 35–45 km/s range. Moreover, if different fragments have different velocities—as suggested by the bipolar geometry—the timescale τ_Rc = 0.27 pc / v_exp would be biased. Please propagate the full uncertainty (including the α dependence and any bow-to-bow scatter) into the quoted τ_Rc ~ 7×10^3 yr, or state explicitly that the conclusion is insensitive to a plausible spread in v_exp.
  3. [Main text, 'The relation between R_c and the distance from the white dwarf' (Fig. 4)] The log-linear fit log R_c = 0.34 − 1.59 r is quoted without uncertainties on the slope, intercept, or scatter. The Methods state that rms ranges from the four fit families are used as uncertainties, but these are not reported in the text. Without the slope uncertainty, the reader cannot assess the significance of the trend or the uncertainty on the e-folding length (0.27 pc) and on τ_Rc. Please report the fit parameters with uncertainties, the scatter of individual points, and propagate the scatter and the four-family rms into the derived timescale. Also specify the base of the logarithm used.
  4. [Main text, 'Because R_c is a purely geometric quantity...'] The paper uses R_c as a direct proxy for the coherent fragment size, but the curvature of a bow shock also depends on the relative velocity and the ambient density. For a wind-driven bow, R_0 ∝ (Mdot v_w / ρ_ISM v_rel^2)^{1/2}; for a fragment-driven bow, an analogous dependence can enter through the mass-loading and the local ram pressure. A radial gradient in the ambient density or in the fragment velocities across the outer halo could therefore produce a decreasing R_c with r even for identical fragments. The text acknowledges that R_c 'does not by itself specify the detailed momentum balance within the flow,' but then uses the R_c–r trend to infer physical disruption. Please provide a scaling argument or simple model that connects R_c to the fragment size for the relevant flow geometry, and demonstrate that plausible density/velocity gradients (e.g., from the known ISM interaction) cannot reproduce the observed factor ~40 decrease over ~1 pc.
minor comments (7)
  1. [Page 1, main text] The word 'undected' in 'The arcs are undected in [O III]' should be 'undetected'.
  2. [Abstract and Fig. 4] The abstract states that R_c decreases by a factor of ~100 over r=0.4–1.4 pc, but the quoted log-linear fit implies a factor 10^1.59 ≈ 39 over the same range; please reconcile the numbers.
  3. [Main text, Fig. 4 caption] The base of the logarithm in the fit log R_c = 0.34 − 1.59 r should be stated (presumably base 10, consistent with the e-folding length of 0.27 pc).
  4. [References] References 14 and 38 are the same paper (Hora et al. 2006, ApJ 652, 426); please remove the duplicate.
  5. [Methods, 'Shock velocity'] The line ratios for arc 14 are quoted as [N II]/Hα = 0.06±0.01 and [O III]/Hα <0.015 (2σ); please specify the extraction aperture and the noise estimate used for the upper limit.
  6. [Extended Data Figure 3] The right panel should label the axes (Hα luminosity and shock velocity) directly on the figure; as rendered, the axis labels are not visible.
  7. [Main text, 'Data availability'] The phrase 'Figure Extended Data Figure 1' should be 'Extended Data Figure 1'.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the R_c–r trend is directly measured and the disruption timescale is obtained by dividing its e-folding length by an independently derived expansion velocity.

full rationale

The paper's central chain is: measure radii of curvature R_c of 22 bow shocks, fit log R_c = 0.34 − 1.59 r (e-folding length 0.27 pc), and convert to a timescale τ_Rc ≈ 7×10^3 yr by dividing by v_exp ≈ 40 km/s. The expansion velocity is not fitted to the R_c–r data; it is derived from MAPPINGS V shock models applied to the [O III]/Hα and [N II]/Hα ratios of the brightest bow, combined with the Gaia-based ISM velocity. The Hα luminosity asymmetry between east and west is checked against the same shock models as a consistency test, not used to force the R_c trend. The assumption that radial distance maps to age t_i ≈ r_i/v_exp is an interpretive kinematic assumption, which could be incorrect if ejection times or velocities vary, but it is not a mathematical identity and does not reduce the timescale to the input of the fit. Self-citations concern data-reduction tools (Dragonfly Spectral Line Mapper, maskfill, dfreproject) and are not load-bearing for the physical conclusion. The lack of propagated uncertainties in the log-linear fit is a statistical reporting concern, not circularity. The paper explicitly cautions that the timescale is a survival time of the coherent fragment/bow system rather than a direct measurement of mass-loss rate, further indicating that the interpretation is not disguised as an input.

Assumptions & free parameters 5 free parameters · 5 assumptions · 1 invented entities

The central claim rests on the interpretation of the arcs as bow shocks from unresolved neutral fragments, the geometric R_c trend, and a constant-velocity clock. The MAPPINGS V model inputs (n_H, alpha) bracket the shock velocity but are not fitted to the R_c-r data. No new physical entities are required beyond the inferred dense neutral fragments, which have a falsifiable handle in predicted CO and H2 emission.

free parameters (5)
  • Expansion velocity v_exp = 35-45 km/s (fiducial 40 km/s)
    Derived from MAPPINGS V shock velocity (80-90 km/s) minus the Helix ISM velocity (45 km/s); used to convert radial distances to ages and the R_c e-folding length to a disruption timescale. Not fitted to the R_c-r data but the central value sets the timescale.
  • Log-linear R_c-r fit (intercept and slope) = intercept 0.34, slope -1.59 per pc
    Descriptive fit to the 22 measured curvatures; the slope is the observational basis for the e-folding length and timescale. No error bars quoted.
  • Post-shock wake velocity v_wake = 0-10 km/s (fiducial 5 km/s)
    Adopted from the Mira bow-shock scaling; used only for the western-side velocity estimate and the brightness asymmetry argument, not for the timescale.
  • Pre-shock hydrogen density n_H = 5 cm^-3
    Input to MAPPINGS V shock grid; the paper states the line ratios are insensitive to this choice.
  • Magnetic field parameter alpha = B / sqrt(n_H) = 1 and 2 (two regimes)
    Chosen to bracket magnetic compression effects in the MAPPINGS V models; affects the inferred shock velocity (80-90 km/s) and therefore v_exp and the timescale.
assumptions (5)
  • domain assumption The Helix Nebula is at distance d = 198.6 (+1.6/-1.8) pc from Gaia EDR3/DR3 astrometry.
    All physical scales (r, R_c in pc) and the linear sizes of the bows scale with this distance. Adopted from Bailer-Jones et al. (2021), reference 10.
  • domain assumption The eastern arcs are bow shocks driven by dense, largely neutral fragments moving supersonically through the ISM.
    Central interpretation established in the main text. The absence of H-alpha, [O III], and [N II] emission at the parabolic foci supports but does not prove that the fragments are neutral and unresolved.
  • domain assumption Projected radius r_i maps to ejection time as t_i = r_i / v_exp with one constant expansion velocity for all fragments.
    Used in the main text to turn the R_c-r slope into a timescale. Ignores spread in 3D orientation, ejection epoch, and velocity.
  • domain assumption The apex radius of curvature R_c is a monotonic proxy for the coherent size of the bow-forming fragment.
    The paper explicitly uses R_c as a size indicator without adopting a specific bow-shock solution; the monotonic mapping is assumed.
  • domain assumption MAPPINGS V shock models with solar abundances, n_H = 5 cm^-3, and alpha = 1 or 2 reproduce the observed line ratios.
    Used to infer v_shock = 80-90 km/s from [N II]/H-alpha and [O III]/H-alpha. The choice of model and the magnetic-field parameterization affects the inferred velocity and hence v_exp and the timescale.
invented entities (1)
  • Dense neutral AGB shell fragments at the bow foci independent evidence
    purpose: To explain the 22 bow shocks and their shrinking curvature with radius; the fragments are the obstacles whose stripping produces the morphological sequence.
    The fragments are not directly detected (no H-alpha, [O III], or [N II] at the foci), but the paper predicts they should be detectable in CO rotational lines and H2 1-0 S(1) 2.12 micron emission, which is a falsifiable handle outside this data set.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Numerous bow shocks in the outer Helix Nebula." pith.science (2026). https://pith.science/paper/VCPKYV4Q

@misc{pith2026260811443,
  author       = {Pith},
  title        = {Pith review of: Numerous bow shocks in the outer Helix Nebula},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/VCPKYV4Q}},
  note         = {Machine review of arXiv:2608.11443}
}
abstract

Near the end of their lives, low- and intermediate-mass stars expel metal-enriched material in winds and outflows, ultimately producing planetary nebulae (PNe). The ejected material is expected to fragment and mix into the interstellar medium (ISM), but this final assimilation step has been difficult to observe directly. Here we report evidence for this process in the form of twenty-two bow shocks in the eastern outskirts of the Helix Nebula, detected in H$\alpha$ emission with the partially completed MOTHRA telescope. Unlike the large-scale wind-ISM bow shocks commonly observed around evolved stars, the shocks are compact and associated with individual clumps of gas. Going outward from the central star, the radius of curvature $R_c$ decreases by a factor of ~100 over the radial range r=0.4-1.4 pc. This is accompanied by a morphological transition from thin, well-defined bows to fuzzy, patchy structures. We interpret these changes as progressive stripping and fragmentation of AGB-shell remnants as they interact with the ISM. The slope of the observed $R_c - r$ relation implies a loss of fragment coherence on a timescale of ~10,000 yr, providing a rare direct constraint on the time-scale for disruption and entrainment of fragmented stellar ejecta into the ISM.

Figures

Figures reproduced from arXiv: 2608.11443 by the authors.

Figure 1
Figure 1. MOTHRA Hα imaging of the Helix nebula. The MOTHRA continuum-subtracted Hα image is shown with an inverted grey scale, emphasizing faint outer features. In the bright central regions a combined Hubble Space Telescope (HST) and Kitt Peak 4m image is superposed on the MOTHRA data.11 The HST image was generated from data in the Advanced Camera for Surveys F502N ([O III]) and F658N (Hα) filters. The arrow indicates the G… view at source ↗
Figure 2
Figure 2. Velocities experienced by fragments. Schematic veloc￾ity field around the Helix, for a bulk velocity with respect to the ISM of vISM = 45 km s−1 eastward, a post-shock flow velocity of vwake = 5 km s−1 , and a radial expansion velocity of the ejecta of vexp = 40 km s−1 . The highest velocities are found on the east side, where we see the strong bow shocks. The scale of the image matches that of [PITH_FULL_IMAGE:fig… view at source ↗
Figure 3
Figure 3. [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Relation between the size and morphology of bow shocks and their distance from the central star. Main panel: Radius of curvature Rc, determined from fitting parabolas to the bows, versus distance from the central white dwarf r. Grey points indicate the location of the …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

49 extracted references · 44 canonical work pages

  1. [1]

    Numerous bow shocks in the outer Helix Nebula

    It shows many features that have not been seen in ionized gas before, such as extensions of the plumes in the north- west and southeast 11, 13 and turbulent and complex low sur- face brightness Hαemission in the southwest. 13 The most striking feature in the Hαimage is a forest of arcs and par- tial arcs on the eastern side of the nebula. We identify at l...

  2. [2]

    Kwok, S., Purton, C. R. & FitzGerald, P. M. On the origin of planetary nebulae.The Astrophysical Journal Letters219, L125–L127 (1978)

  3. [3]

    & Renzini, A

    Iben, I. & Renzini, A. Asymptotic giant branch evo- lution and beyond.Annual Review of Astronomy and Astrophysics21, 271–342 (1983)

  4. [4]

    J., Sarazin, C

    Borkowski, K. J., Sarazin, C. L. & Soker, N. Interaction of planetary nebulae with the interstellar medium.The Astrophysical Journal360, 173–183 (1990)

  5. [5]

    J., Zijlstra, A

    Wareing, C. J., Zijlstra, A. A. & O’Brien, T. J. V ortices in the wakes of asymptotic giant branch stars.The As- trophysical Journal Letters660, L129–L132 (2007)

  6. [6]

    J.et al.Detached shells as tracers of asymp- totic giant branch–interstellar medium bow shocks

    Wareing, C. J.et al.Detached shells as tracers of asymp- totic giant branch–interstellar medium bow shocks. Monthly Notices of the Royal Astronomical Society: Let- ters372, L63–L67 (2006)

  7. [7]

    C.et al.A turbulent wake as a tracer of 30,000 years of mira’s mass loss history.Nature448, 780–783 (2007)

    Martin, D. C.et al.A turbulent wake as a tracer of 30,000 years of mira’s mass loss history.Nature448, 780–783 (2007)

  8. [8]

    Cox, N. L. J.et al.A far-infrared survey of bow shocks and detached shells around agb stars and red super- giants.Astronomy & Astrophysics537, A35 (2012). 1110.5486

Show all 49 references
  1. [9]

    I., McKee, C

    Klein, R. I., McKee, C. F. & Colella, P. On the hy- drodynamic interaction of shock waves with interstellar clouds. 1: Nonradiative shocks in small clouds.Astro- physical Journal420, 213–236 (1994)

  2. [10]

    Karakas, A. I. & Lattanzio, J. C. The dawes re- view 2: Nucleosynthesis and stellar yields of low- and intermediate-mass single stars.Publications of the As- tronomical Society of Australia31, e030 (2014)

  3. [11]

    Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M. & Andrae, R. Estimating distances from parallaxes. v. geometric and photogeometric distances to 1.47 billion stars in gaia early data release 3.The Astro- nomical Journal161, 147 (2021)

  4. [12]

    R., McCullough, P

    O’Dell, C. R., McCullough, P. R. & Meixner, M. Un- raveling the Helix Nebula: Its Structure and Knots.AJ 128, 2339–2356 (2004).astro-ph/0407556

  5. [13]

    Meaburn, J.et al.The creation of the helix planetary nebula (ngc 7293) by multiple events.Monthly Notices of the Royal Astronomical Society360, 963–976 (2005)

  6. [14]

    & Kwok, S

    Zhang, Y ., Hsia, C.-H. & Kwok, S. Discovery of a Halo around the Helix Nebula NGC 7293 in the WISE All- sky Survey.ApJ755, 53 (2012).1207.4606

  7. [15]

    L., Latter, W

    Hora, J. L., Latter, W. B., Smith, H. A. & Marengo, M. Infrared observations of the helix planetary nebula.The Astrophysical Journal652, 426–441 (2006)

  8. [16]

    & Akras, S

    Meaburn, J., Boumis, P. & Akras, S. The bow-shock and high-speed jet in the faint, 40 arcmin diameter, outer halo of the evolved Helix planetary nebula (NGC 7293). MNRAS435, 3462–3468 (2013).1308.5460

  9. [17]

    M.et al.Realizing the potential of the Dragonfly Spectral Line Mapper: calibration methods and on-sky performance

    Lokhorst, D. M.et al.Realizing the potential of the Dragonfly Spectral Line Mapper: calibration methods and on-sky performance. In Marshall, H. K., Spy- romilio, J. & Usuda, T. (eds.)Ground-based and Air- borne Telescopes X, vol. 13094 ofSociety of Photo- Optical Instrumentati...

  10. [18]

    Chen, S.et al.First Light with the 120-lens Dragonfly Spectral Line Mapper.PASP137, 084103 (2025)

  11. [19]

    Gaia Collaboration, Brown, A. G. A.et al.Gaia early data release 3: Summary of the contents and survey properties.Astronomy & Astrophysics649, A1 (2021)

  12. [20]

    D., Libert, Y ., G´erard, E., Le Bertre, T

    Matthews, L. D., Libert, Y ., G´erard, E., Le Bertre, T. & Reid, M. J. The discovery of an h i counterpart to mira’s far-ultraviolet-emitting tail.The Astrophysical Journal 684, 603–615 (2008)

  13. [21]

    Sutherland, R. S. & Dopita, M. A. Mappings v: A new shock and photoionization code.Astrophys. J. Suppl. 237, 12 (2018). 5

  14. [22]

    Meaburn, J., L ´opez, J. A. & Richer, M. G. Optical line profiles of the helix planetary nebula (ngc 7293) to large radii.Monthly Notices of the Royal Astronomical Soci- ety384, 497–503 (2008)

  15. [23]

    & Olofsson, H

    H ¨ofner, S. & Olofsson, H. Mass loss of stars on the asymptotic giant branch. mechanisms, models and mea- surements.The Astronomy and Astrophysics Review26, 1 (2018)

  16. [24]

    & Frank, A

    Balick, B. & Frank, A. Shapes and Shaping of Planetary Nebulae.ARA&A40, 439–486 (2002)

  17. [25]

    Hst imaging of proto-planetary nebulae and very young planetary nebulae—towards a new under- standing of their formation

    Sahai, R. Hst imaging of proto-planetary nebulae and very young planetary nebulae—towards a new under- standing of their formation. InPost-AGB Objects as a Phase of Stellar Evolution, vol. 265 ofAstrophysics and Space Science Library, 53–63 (Springer, 2001)

  18. [26]

    & S´anchez Contreras, C

    Bujarrabal, V ., Castro-Carrizo, A., Alcolea, J. & S´anchez Contreras, C. Mass, linear momentum and ki- netic energy of bipolar flows in protoplanetary nebulae. Astronomy & Astrophysics377, 868–897 (2001)

  19. [27]

    Wilkin, F. P. Exact analytic solutions for stellar wind bow shocks.Astrophys. J. Lett.459, L31–L34 (1996)

  20. [28]

    Tarango-Yong, J. A. & Henney, W. J. True versus ap- parent shapes of bow shocks.MNRAS477, 2431–2454 (2018).1712.02300

  21. [29]

    W., Dyson, J

    Hartquist, T. W., Dyson, J. E., Pettini, M. & Smith, L. J. Mass-loaded astronomical flows – i. general principles and their application to rcw 58.Monthly Notices of the Royal Astronomical Society221, 715–726 (1986)

  22. [30]

    M., Hartquist, T

    Pittard, J. M., Hartquist, T. W. & Dyson, J. E. Self- similar evolution of wind-blown bubbles with mass loading by hydrodynamic ablation.Astronomy & As- trophysics373, 1043–1055 (2001)

  23. [31]

    Leitner, S. N. & Kravtsov, A. V . Fuel efficient galaxies: Sustaining star formation with stellar mass loss.The Astrophysical Journal734, 48 (2011).1011.1252

  24. [32]

    D., Le Bertre, T., G ´erard, E

    Matthews, L. D., Le Bertre, T., G ´erard, E. & John- son, M. C. An h i imaging survey of asymptotic giant branch stars.The Astronomical Journal145, 97 (2013). 1301.7429

  25. [33]

    W., Veeravalli, G

    Wadsley, J. W., Veeravalli, G. & Couchman, H. M. P. On the treatment of entropy mixing in numerical cosmol- ogy.Monthly Notices of the Royal Astronomical Society 387, 427–438 (2008)

  26. [34]

    & Stinson, G

    Shen, S., Wadsley, J. & Stinson, G. The enrichment of the intergalactic medium with adiabatic feedback – i. metal cooling and metal diffusion.Monthly Notices of the Royal Astronomical Society407, 1581–1596 (2010). 0910.5956

  27. [35]

    F., Dav ´e, R

    Rennehan, D., Babul, A., Hopkins, P. F., Dav ´e, R. & Moa, B. Dynamic localized turbulent diffusion and its impact on the galactic ecosystem.Monthly Notices of the Royal Astronomical Society483, 3810–3831 (2019)

  28. [36]

    J., Bachiller, R., Cox, P

    Huggins, P. J., Bachiller, R., Cox, P. & Forveille, T. Co in the cometary globules of the helix nebula.The Astro- physical Journal Letters401, L43–L46 (1992)

  29. [37]

    J.et al.High-resolution co and h 2 molecu- lar line imaging of a cometary globule in the helix neb- ula.The Astrophysical Journal Letters573, L55–L58 (2002).astro-ph/0205516

    Huggins, P. J.et al.High-resolution co and h 2 molecu- lar line imaging of a cometary globule in the helix neb- ula.The Astrophysical Journal Letters573, L55–L58 (2002).astro-ph/0205516

  30. [38]

    firework

    Matsuura, M.et al.A “firework” of h 2 knots in the plan- etary nebula ngc 7293 (the helix nebula).The Astrophys- ical Journal700, 1067–1077 (2009).0906.2870

  31. [39]

    L., Latter, W

    Hora, J. L., Latter, W. B., Smith, H. A. & Marengo, M. Infrared observations of the helix planetary neb- ula.The Astrophysical Journal652, 426–441 (2006). astro-ph/0607541

  32. [40]

    Evolution of planetary nebulae

    Balick, B. Evolution of planetary nebulae. i. structures, ionizations, and morphological sequences.The Astro- nomical Journal94, 671–678 (1987)

  33. [41]

    Corradi, R. L. M., Manso, R., Mampaso, A. & Schwarz, H. E. Unveiling low-ionization microstructures in plan- etary nebulae.Astronomy & Astrophysics313, 913–923 (1996)

  34. [42]

    R., Corradi, R

    Gonc ¸alves, D. R., Corradi, R. L. M. & Mampaso, A. Low-ionization structures in planetary nebulae: Con- fronting models with observations.The Astrophysical Journal547, 302–310 (2001)

  35. [43]

    Richer, M. G. Optical line profiles of the helix planetary nebula (ngc 7293) to large radii.Monthly Notices of the Royal Astronomical Society384, 497–503 (2008)

  36. [44]

    C., van Hoof, P

    Van de Steene, G. C., van Hoof, P. A. M., Exter, K. M., Barlow, M. J.et al.Herschel imaging of the dust in the helix nebula (ngc 7293).Astronomy & Astrophysics 574, A134 (2015)

  37. [45]

    & Pasha, I

    van Dokkum, P. & Pasha, I. A Robust and Simple Method for Filling in Masked Data in Astronomical Im- ages.PASP136, 034503 (2024).2312.03064

  38. [46]

    A.et al.Tidal Distortions in NGC1052-DF2 and NGC1052-DF4: Independent Evidence for a Lack of Dark Matter.ApJ935, 160 (2022).2109.09778

    Keim, M. A.et al.Tidal Distortions in NGC1052-DF2 and NGC1052-DF4: Independent Evidence for a Lack of Dark Matter.ApJ935, 160 (2022).2109.09778

  39. [47]

    & Raymond, J

    Jin, Y . & Raymond, J. Dialog Concerning the Two Shock Codes.ApJ989, 203 (2025).2507.03225

  40. [48]

    6 Figure Extended Data Figure 1:Color representation of the Helix Nebula.Combination of MOTHRA [O III], Hα, and [N II] imaging of the Helix

    Rhea, C.et al.dfreproject: A Python package for as- tronomical reprojection.The Journal of Open Source Software10, 8525 (2025).2505.03932. 6 Figure Extended Data Figure 1:Color representation of the Helix Nebula.Combination of MOTHRA [O III], Hα, and [N II] imaging of the Heli...

  41. [49]

    As discussed in the main text, we use the parabola as the fiducial profile, and take the rms ranges inR c andR 0 from the four fits as uncertainties

    The fits are stable for the 22 bows that we identify in the image. As discussed in the main text, we use the parabola as the fiducial profile, and take the rms ranges inR c andR 0 from the four fits as uncertainties. Data availability.The reduced, continuum- subtracted Hαimage...

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.