3D simulations of AGB stellar winds -- II. Ray-tracer implementation and impact of radiation on the outflow morphology

Mats Esseldeurs,Lionel Siess,Frederik De Ceuster,Ward Homan,Jolien Malfait,Silke Maes,Thomas Konings,Thomas Ceulemans,Leen Decin

DOI: https://doi.org/10.1051/0004-6361/202346282

2023-04-20

Abstract:Stars with an initial mass below ~ 8 Msun evolve through the asymptotic giant branch (AGB) phase, during which they develop strong stellar winds. Recent observations have revealed significant morphological complexities in their outflows, most likely caused by a companion. We study the impact of the radiation force on such companion-perturbed AGB outflows. We present the implementation of a ray tracer for radiative transfer in smoothed particle hydrodynamics (SPH) and compared four different descriptions of radiative transfer: the free-wind, the geometrical, the Lucy, and the attenuation approximation. For both low and high mass-loss rates, the velocity profile of the outflow is modified when going from the free-wind to the geometrical approximation, also resulting in a different morphology. In the case of a low mass-loss rate, the effect of the Lucy and attenuation approximation is negligible due to the low densities but morphological differences appear in the high mass-loss rate regime. By comparing the radiative equilibrium temperature and radiation force to full 3D radiative transfer, we show that the Lucy approximation works best. Although, close to the companion, artificial heating occurs and it fails to simulate the shadow cast by the companion. The attenuation approximation produces a lower equilibrium temperature and weaker radiation force, but it produces the shadow cast by the companion. From the predictions of the 3D radiative transfer, we also conclude that a radially directed radiation force is a reasonable assumption. The radiation force thus plays a critical role in dust-driven AGB winds, impacting the velocity profile and morphological structures. For low mass-loss rates, the geometrical approximation suffices, while high mass-loss rates require a more rigorous method, where the Lucy approximation provides the most accurate results although not accounting for all effects.

Solar and Stellar Astrophysics,Astrophysics of Galaxies

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: How to describe the morphology of asymptotic giant branch (AGB) stellar winds more accurately, especially when considering the influence of radiation forces on dust - driven winds and the influence of companion stars on the morphology of AGB stellar winds. Specifically, the research aims to simulate radiation transfer by implementing a ray tracer in the smoothed - particle hydrodynamics (SPH) code Phantom to create a three - dimensional optical depth map around AGB stars and analyze the effects of different radiation transfer description methods.

### Research Background and Objectives

**Background**:
- **AGB Stars**: Stars with an initial mass lower than about 8 solar masses will experience the AGB phase, during which they will produce strong stellar winds due to the radiation pressure on newly - formed dust grains.
- **Observational Findings**: Recent observations have revealed significant morphological complexity in AGB outflows, which are most likely caused by the interaction of companion stars.

**Objectives**:
- **Describe the Morphology of AGB Winds More Accurately**: The research aims to provide a more accurate description by considering the role of radiation forces in dust - driven winds and the influence of companion stars on the morphology of AGB winds.
- **Method**: Implement a ray tracer in the SPH code Phantom for the calculation of radiation transfer and compare four radiation transfer description methods with different levels of complexity: free - wind approximation, geometric approximation, Lucy approximation, and attenuation approximation.

### Methods

**Implementation of the Ray Tracer**:
- **Calculation of Optical Depth**: Calculate the optical depth of each SPH particle to the AGB star through the ray tracer to create a three - dimensional optical depth map.
- **Approximation Methods for Radiation Transfer**:
- **Free - Wind Approximation**: Do not consider radiation transfer and assume that Γ is equal to 1.
- **Geometric Approximation**: Consider the geometric dilution effect.
- **Lucy Approximation**: Consider that the radiation comes not only from the AGB star but also from the re - emission of the surrounding medium.
- **Attenuation Approximation**: Consider the attenuation of stellar radiation.

**Numerical Implementation**:
- **SPH Code**: Use the Phantom code for 3D hydrodynamics simulations.
- **Interpolation of Optical Depth**: Interpolate through a predefined number of rays to reduce the computational cost.

### Results

**Effects of Different Radiation Transfer Treatments**:
- **Low and High Mass - Loss Rates**: The velocity profiles and wind morphologies at low mass - loss rates (10^(-8) solar masses per year) and high mass - loss rates (3×10^(-6) solar masses per year) were studied.
- **Single AGB Star**: Under the free - wind approximation, the material accelerates immediately and propagates at an almost constant velocity; under the geometric approximation, the material starts to accelerate when it is close to the dust condensation radius.
- **AGB Binary Systems**: The presence of a companion star leads to different wind morphologies, especially at high mass - loss rates, and the effects of the Lucy approximation and the attenuation approximation are more obvious.

### Conclusions

- **The Crucial Role of Radiation Forces**: Radiation forces play a crucial role in dust - driven AGB winds, affecting the velocity profile and morphological structure.
- **Applicability**: For low mass - loss rates, the geometric approximation is sufficient; for high mass - loss rates, a more rigorous method is required. Among the methods studied, the Lucy approximation provides the most accurate results, although it cannot fully simulate the shadow effect of the companion star.

Through these studies, the paper provides important insights into understanding the formation mechanism and morphological complexity of AGB stellar winds.