Formation and X-ray emission from hot bubbles in planetary nebulae - III. The impact of [Wolf-Rayet]-type winds

Technical Summary: Formation and X-ray Emission from Hot Bubbles in Planetary Nebulae - III. The Impact of [Wolf-Rayet]-type Winds

Problem Statement
Planetary nebulae (PNe) exhibit hot, X-ray-emitting bubbles formed by the interaction of fast post-Asymptotic Giant Branch (post-AGB) winds with slower, dense AGB ejecta. While standard models assume hydrogen-rich (H-rich) central stars, a significant fraction of X-ray-emitting PNe host hydrogen-deficient, [Wolf-Rayet]-type ([WR]) central stars. Previous studies suggested that the high metallicity (C-rich) of [WR] winds enhances radiative cooling, potentially delaying bubble formation and altering thermal structures. However, existing models often relied on 1D simulations with artificial wind enhancements or assumed thermal conduction as the primary mechanism for cooling the hot plasma to observed soft X-ray temperatures ($T_X \approx 1-3 \times 10^6$ K). The specific impact of self-consistent [WR] wind parameters, derived from recent empirical mass-loss prescriptions, on the hydrodynamical evolution and synthetic X-ray emission of hot bubbles remains to be fully explored in a multidimensional context.

Methodology
The authors employ radiation-hydrodynamical simulations using the PLUTO code to model the circumstellar medium of PNe. The methodology involves a multi-stage approach:

1. Stellar Evolution: Stellar evolution models for progenitors of $1, 2,$ and $3 M_\odot$ are generated using MESA. Two evolutionary tracks are computed for each mass: a standard H-rich post-AGB track (using the Miller Bertolami 2016 mass-loss prescription) and a [WR] track utilizing a recent mass-loss rate prescription for H-deficient stars derived by Toalá, Todt & Sander (2024).
2. Hydrodynamical Setup:
   • 1D Simulations: Performed on radial grids to study global bubble formation and evolution, incorporating specific cooling curves for H-rich and H-deficient (C-rich) abundances generated via Cloudy.
   • 2D Simulations: Performed on spherical grids to capture hydrodynamical instabilities (Vishniac thin-shell and Rayleigh-Taylor) at the wind-wind interaction interface. These simulations track the mixing of AGB material and post-AGB wind using an advected scalar parameter.
3. X-ray Synthesis: Synthetic X-ray spectra (0.3–2.0 keV) are calculated post-processing using the CHIANTI atomic database. The Differential Emission Measure (DEM) is computed separately for H-rich and H-poor gas phases based on the advected scalar, allowing for the calculation of emissivity-weighted plasma temperatures ($T_X$) and luminosities ($L_X$) that account for the mixed composition of the bubble.

Key Contributions

• Self-Consistent [WR] Evolution: Unlike previous works that artificially enhanced wind parameters for a single stellar mass, this study utilizes MESA tracks with self-consistent [WR] mass-loss rates across a range of progenitor masses ($1-3 M_\odot$).
• Multidimensional Mixing: The study explicitly avoids thermal conduction, instead demonstrating that hydrodynamical instabilities in 2D simulations naturally generate turbulent mixing layers that regulate bubble temperatures.
• Dual-Component X-ray Analysis: The authors introduce a rigorous method to separate and recombine the X-ray contributions of H-rich and H-poor gas phases, acknowledging that [WR] bubbles contain a mixture of chemically distinct materials.

Results

• Delayed Bubble Formation: The enhanced radiative cooling in C-rich [WR] winds significantly delays the formation of hot bubbles compared to H-rich counterparts. For a $1 M_\odot$ model, bubble formation is delayed by approximately 6,000 years post-AGB onset, whereas higher mass models show less dramatic delays. This delay is attributed to the requirement that wind velocity must exceed a critical threshold ($v_{crit}$) to overcome rapid radiative cooling.
• Enhanced X-ray Luminosity: Models incorporating [WR] winds exhibit significantly higher X-ray luminosities ($L_X$) than H-rich models. This is driven by the higher mechanical wind luminosity ($L_{wind}$) resulting from the elevated mass-loss rates of [WR] stars.
• Convergence of Plasma Temperature: Despite the higher wind power and different chemical compositions, the emissivity-weighted plasma temperature ($T_X$) of the X-ray-emitting gas converges to values of $1-3 \times 10^6$ K for both [WR] and H-rich models once mixing is established.
• Role of Mixing vs. Thermal Conduction: The 2D simulations show that turbulent mixing efficiently transports cooler nebular material into the hot bubble, regulating the temperature to observed levels without requiring thermal conduction. Comparisons with 1D models reveal that 1D simulations systematically overestimate $T_X$ (by factors of 1.2–1.6, up to ~5 in massive models) due to the absence of this mixing mechanism.
• Surface Brightness: The X-ray surface brightness peaks early and declines rapidly (by $\sim 4$ orders of magnitude) as the bubble expands, suggesting that detectable diffuse X-ray emission is restricted to a specific evolutionary window where the bubble has expanded sufficiently but before it cools and dilutes too much.

Significance and Claims
The paper claims to reinforce the hypothesis that hydrodynamical mixing is a key mechanism for generating the observed soft X-ray emission in PNe, even for systems hosting [WR] central stars. The authors assert that their results:

1. Explain the absence of X-ray emission in young, compact [WR] PNe as a consequence of delayed bubble formation due to strong radiative cooling.
2. Demonstrate that the observed soft X-ray temperatures ($T_X \approx 1-3 \times 10^6$ K) are a natural outcome of multidimensional mixing, independent of whether the central star follows an H-rich or [WR] evolutionary path.
3. Provide a robust alternative to thermal conduction models, suggesting that soft X-ray emission can arise even in the presence of magnetic fields that might suppress electron thermal conduction.
4. Validate the use of empirical [WR] mass-loss prescriptions in hydrodynamical modeling, showing they reproduce observed stellar wind properties ($L, T_{eff}, \dot{M}, v_\infty$) of X-ray-emitting CSPNe.

The authors conclude that while their wind approximation is practical, future work should ideally model the "born-again" scenario explicitly to fully capture the transition to H-deficient winds. However, the current approach provides a physically motivated framework for understanding the X-ray properties of [WR] PNe.