-
Chemical tracers of a highly eccentric AGB-main sequence star binary
Authors:
T. Danilovich,
J. Malfait,
M. Van de Sande,
M. Montargès,
P. Kervella,
F. De Ceuster,
A. Coenegrachts,
T. J. Millar,
A. M. S. Richards,
L. Decin,
C. A. Gottlieb,
C. Pinte,
E. De Beck,
D. J. Price,
K. T. Wong,
J. Bolte,
K. M. Menten,
A. Baudry,
A. de Koter,
S. Etoka,
D. Gobrecht,
M. Gray,
F. Herpin,
M. Jeste,
E. Lagadec
, et al. (10 additional authors not shown)
Abstract:
Binary interactions have been proposed to explain a variety of circumstellar structures seen around evolved stars, including asymptotic giant branch (AGB) stars and planetary nebulae. Studies resolving the circumstellar envelopes of AGB stars have revealed spirals, discs and bipolar outflows, with shaping attributed to interactions with a companion. For the first time, we have used a combined chem…
▽ More
Binary interactions have been proposed to explain a variety of circumstellar structures seen around evolved stars, including asymptotic giant branch (AGB) stars and planetary nebulae. Studies resolving the circumstellar envelopes of AGB stars have revealed spirals, discs and bipolar outflows, with shaping attributed to interactions with a companion. For the first time, we have used a combined chemical and dynamical analysis to reveal a highly eccentric and long-period orbit for W Aquilae, a binary system containing an AGB star and a main sequence companion. Our results are based on anisotropic SiN emission, the first detections of NS and SiC towards an S-type star, and density structures observed in the CO emission. These features are all interpreted as having formed during periastron interactions. Our astrochemistry-based method can yield stringent constraints on the orbital parameters of long-period binaries containing AGB stars, and will be applicable to other systems.
△ Less
Submitted 23 July, 2024;
originally announced July 2024.
-
ATOMIUM: Molecular inventory of 17 oxygen-rich evolved stars observed with ALMA
Authors:
S. H. J. Wallstrom,
T. Danilovich,
H. S. P. Muller,
C. A. Gottlieb,
S. Maes,
M. Van de Sande,
L. Decin,
A. M. S. Richards,
A. Baudry,
J. Bolte,
T. Ceulemans,
F. De Ceuster,
A. de Koter,
I. El Mellah,
M. Esseldeurs,
S. Etoka,
D. Gobrecht,
E. Gottlieb,
M. Gray,
F. Herpin,
M. Jeste,
D. Kee,
P. Kervella,
T. Khouri,
E. Lagadec
, et al. (13 additional authors not shown)
Abstract:
The dusty winds of cool evolved stars are a major contributor of the newly synthesised material enriching the Galaxy and future generations of stars. However, the details of the physics and chemistry behind dust formation and wind launching have yet to be pinpointed. Recent spatially resolved observations show the importance of gaining a more comprehensive view of the circumstellar chemistry, but…
▽ More
The dusty winds of cool evolved stars are a major contributor of the newly synthesised material enriching the Galaxy and future generations of stars. However, the details of the physics and chemistry behind dust formation and wind launching have yet to be pinpointed. Recent spatially resolved observations show the importance of gaining a more comprehensive view of the circumstellar chemistry, but a comparative study of the intricate interplay between chemistry and physics is still difficult because observational details such as frequencies and angular resolutions are rarely comparable. Aiming to overcome these deficiencies, ATOMIUM is an ALMA Large Programme to study the physics and chemistry of the circumstellar envelopes of a diverse set of oxygen-rich evolved stars under homogeneous observing conditions at three angular resolutions between ~0.02"-1.4". Here we summarize the molecular inventory of these sources, and the correlations between stellar parameters and molecular content. Seventeen oxygen-rich or S-type asymptotic giant branch (AGB) and red supergiant (RSG) stars have been observed in several tunings with ALMA Band 6, targeting a range of molecules to probe the circumstellar envelope and especially the chemistry of dust formation close to the star. We systematically assigned the molecular carriers of the spectral lines and measured their spectroscopic parameters and the angular extent of the emission of each line from integrated intensity maps. Across the ATOMIUM sample, we detect 291 transitions of 24 different molecules and their isotopologues. This includes several first detections in oxygen-rich AGB/RSG stars: PO v=1, SO2 v1=1 and v2=2, and several high energy H2O transitions. We also find several first detections in S-type AGB stars: vibrationally excited HCN v2=2,3 and SiS v=4,5,6, as well as first detections of the molecules SiC, AlCl, and AlF in W Aql...
△ Less
Submitted 6 December, 2023;
originally announced December 2023.
-
3D simulations of AGB stellar winds -- II. Ray-tracer implementation and impact of radiation on the outflow morphology
Authors:
Mats Esseldeurs,
Lionel Siess,
Frederik De Ceuster,
Ward Homan,
Jolien Malfait,
Silke Maes,
Thomas Konings,
Thomas Ceulemans,
Leen Decin
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 tra…
▽ More
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.
△ Less
Submitted 19 April, 2023;
originally announced April 2023.
-
The VLT/SPHERE view of the ATOMIUM cool evolved star sample. I. Overview: Sample characterization through polarization analysis
Authors:
M. Montargès,
E. Cannon,
A. de Koter,
T. Khouri,
E. Lagadec,
P. Kervella,
L. Decin,
I. McDonald,
W. Homan,
L. B. F. M. Waters,
R. Sahai,
C. A. Gottlieb,
J. Malfait,
S. Maes,
B. Pimpanuwat,
M. Jeste,
T. Danilovich,
F. De Ceuster,
M. Van de Sande,
D. Gobrecht,
S. H. J. Wallström,
K. T. Wong,
I. El Mellah,
J. Bolte,
F. Herpin
, et al. (10 additional authors not shown)
Abstract:
Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars).
Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution dat…
▽ More
Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars).
Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution data. To help interpret the polarized signal, we produced synthetic maps of light scattering by dust, through 3D radiative transfer simulations with the RADMC3D code.
Results. The degree of linear polarization (DoLP) observed by ZIMPOL spreads across several optical filters. We infer that it primarily probes dust located just outside of the point spread function, and in or near the plane of the sky, with a total optical depth close to unity in the line of sight, representing only a fraction of the total circumstellar dust. The maximum DoLP ranges from 0.03-0.38 depending on the source, fractions that can be reproduced by our 3D pilot models for grains composed of common dust species. The spatial structure of the DoLP shows a diverse set of shapes. Only for three sources do we note a correlation between the ALMA CO and SiO lines, which trace the gas density, and the DoLP, which traces the dust.
Conclusion. The clumpiness of the DoLP and the lack of a consistent correlation between the gas and the dust location show that, in the inner circumstellar environment (CSE), dust formation occurs at very specific sites. This has potential consequences for the derived mass-loss rates and dust-to-gas ratio in the inner region of the CSE. Except for $π^1$~Gru and perhaps GY Aql, we do not detect interactions between the circumstellar wind and the hypothesized companions that shape the wind at larger scales. This suggests that the orbits of any other companions are tilted out of the plane of the sky.
△ Less
Submitted 5 January, 2023;
originally announced January 2023.
-
3D hydrodynamical survey of the impact of a companion on the morphology and dynamics of AGB outflow
Authors:
Jolien Malfait,
Silke Maes,
Ward Homan,
Jan Bolte,
Lionel Siess,
Frederik De Ceuster,
Leen Decin
Abstract:
With the use of high-resolution ALMA observations, complex structures that resemble those observed in post-AGB stars and planetary nebulae are detected in the circumstellar envelopes of low-mass evolved stars. These deviations from spherical symmetry are believed to be caused primarily by the interaction with a companion star or planet. With the use of three-dimensional hydrodynamic simulations, w…
▽ More
With the use of high-resolution ALMA observations, complex structures that resemble those observed in post-AGB stars and planetary nebulae are detected in the circumstellar envelopes of low-mass evolved stars. These deviations from spherical symmetry are believed to be caused primarily by the interaction with a companion star or planet. With the use of three-dimensional hydrodynamic simulations, we study the impact of a binary companion on the wind morphology and dynamics of an AGB outflow. We classifiy the wind structures and morphology that form in these simulations with the use of a classification parameter, constructed with characteristic parameters of the binary configuration. Finally we conclude that the companion alters the wind expansion velocity through the slingshot mechanism, if it is massive enough.
△ Less
Submitted 28 July, 2022;
originally announced July 2022.
-
Route towards complete 3D hydro-chemical simulations of companion-perturbed AGB outflows
Authors:
Silke Maes,
Lionel Siess,
Ward Homan,
Jolien Malfait,
Frederik De Ceuster,
Thomas Ceulemans,
Dion Donne,
Mats Esseldeurs,
Leen Decin
Abstract:
Low- and intermediate mass stars experience a significant mass loss during the last phases of their evolution, which obscures them in a vast, dusty envelope. Although it has long been thought this envelope is generally spherically symmetric in shape, recent high-resolution observations find that most of these stars exhibit complex and asymmetrical morphologies, most likely resulting from binary in…
▽ More
Low- and intermediate mass stars experience a significant mass loss during the last phases of their evolution, which obscures them in a vast, dusty envelope. Although it has long been thought this envelope is generally spherically symmetric in shape, recent high-resolution observations find that most of these stars exhibit complex and asymmetrical morphologies, most likely resulting from binary interaction. In order to improve our understanding about these systems, theoretical studies are needed in the form of numerical simulations. Currently, a handful of simulations exist, albeit they mainly focus on the hydrodynamics of the outflow. Hence, we here present the pathway to more detailed and accurate modelling of companion-perturbed outflows with Phantom, by discussing the missing but crucial physical and chemical processes. With these state-of-the-art simulations we aim to make a direct comparison with observations to unveil the true identity on the embedded systems.
△ Less
Submitted 24 June, 2022;
originally announced June 2022.
-
ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars: Motivation, sample, calibration, and initial results
Authors:
C. A. Gottlieb,
L. Decin,
A. M. S. Richards,
F. De Ceuster,
W. Homan,
S. H. J. Wallstrom,
T. Danilovich,
T. J. Millar,
M. Montarges,
K. T. Wong,
I. McDonald,
A. Baudry,
J. Bolte,
E. Cannon,
E. De Beck,
A. de Koter,
I. El Mellah,
S. Etoka,
D. Gobrecht,
M. Gray,
F. Herpin,
M. Jeste,
P. Kervella,
T. Khouri,
E. Lagadec
, et al. (11 additional authors not shown)
Abstract:
This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical e…
▽ More
This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. 17 oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (0.025-0.050 arcsec), medium (0.13-0.24 arcsec), and low (about 1 arcsec) angular resolution. The sensitivity per 1.3 km/s channel was 1.5-5 mJy/beam. 13 molecules were designated as primary molecules in the survey: CO, SiO, AlO, AlOH, TiO, TiO2, HCN, SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described; and we highlight one scientific result - the wind kinematics of the ATOMIUM sources. The ATOMIUM sources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion, and the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In 15 sources, some molecular transitions other than 12CO v=0 J=2-1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. Binary interaction with a (sub)stellar companion might (partly) explain the non-monotonic behaviour of the projected velocity field.
△ Less
Submitted 13 December, 2021; v1 submitted 8 December, 2021;
originally announced December 2021.
-
ATOMIUM: Halide molecules around the S-type AGB star W Aquilae
Authors:
T. Danilovich,
M. Van de Sande,
J. M. C. Plane,
T. J. Millar,
P. Royer,
M. A. Amor,
K. Hammami,
L. Decock,
C. A. Gottlieb,
L. Decin,
A. M. S. Richards,
E. De Beck,
A. Baudry,
J. Bolte,
E. Cannon,
F. De Ceuster,
A. de Koter,
S. Etoka,
D. Gobrecht,
M. Gray,
F. Herpin,
W. Homan,
M. Jeste,
P. Kervella,
T. Khouri
, et al. (14 additional authors not shown)
Abstract:
S-type asymptotic giant branch (AGB) stars are thought to be intermediates in the evolution of oxygen- to carbon-rich AGB stars. The chemical compositions of their circumstellar envelopes are also intermediate, but have not been studied in as much detail as their carbon- and oxygen-rich counterparts. We aim to determine the abundances of AlCl and AlF from rotational lines, which have been observed…
▽ More
S-type asymptotic giant branch (AGB) stars are thought to be intermediates in the evolution of oxygen- to carbon-rich AGB stars. The chemical compositions of their circumstellar envelopes are also intermediate, but have not been studied in as much detail as their carbon- and oxygen-rich counterparts. We aim to determine the abundances of AlCl and AlF from rotational lines, which have been observed for the first time towards an S-type AGB star, W Aql. In combination with models based on PACS observations, we aim to update our chemical kinetics network based on these results. We analyse ALMA observations towards W Aql of AlCl in the ground and first two vibrationally excited states and AlF in the ground vibrational state. Using radiative transfer models, we determine the abundances and spatial abundance distributions of Al$^{35}$Cl, Al$^{37}$Cl, and AlF. We also model HCl and HF emission and compare these models to PACS spectra to constrain the abundances of these species. AlCl is found in clumps very close to the star, with emission confined within 0.1$^{\prime\prime}$ of the star. AlF emission is more extended, with faint emission extending 0.2$^{\prime\prime}$ to 0.6$^{\prime\prime}$ from the continuum peak. We find peak abundances, relative to H$_2$, of $1.7\times 10^{-7}$ for Al$^{35}$Cl, $7\times 10^{-8}$ for Al$^{37}$Cl and $1\times 10^{-7}$ for AlF. From the PACS spectra, we find abundances of $9.7\times 10^{-8}$ and $\leq 10^{-8}$, relative to H$_2$, for HCl and HF, respectively. The AlF abundance exceeds the solar F abundance, indicating that fluorine synthesised in the AGB star has already been dredged up to the surface of the star and ejected into the circumstellar envelope. From our analysis of chemical reactions in the wind, we conclude that AlF may participate in the dust formation process, but we cannot fully explain the rapid depletion of AlCl seen in the wind.
△ Less
Submitted 10 September, 2021;
originally announced September 2021.
-
SPH modelling of wind-companion interactions in eccentric AGB binary systems
Authors:
J. Malfait,
W. Homan,
S. Maes,
J. Bolte,
L. Siess,
F. De Ceuster,
L. Decin
Abstract:
The late evolutionary stages of low- and intermediate-mass stars are characterised by mass loss through a dust-driven stellar wind. Recent observations reveal complex structures within these winds, that are believed to be formed primarily via interaction with a companion. How these complexities arise, and which structures are formed in which type of systems, is still poorly understood. Particularl…
▽ More
The late evolutionary stages of low- and intermediate-mass stars are characterised by mass loss through a dust-driven stellar wind. Recent observations reveal complex structures within these winds, that are believed to be formed primarily via interaction with a companion. How these complexities arise, and which structures are formed in which type of systems, is still poorly understood. Particularly, there is a lack of studies investigating the structure formation in eccentric systems. We aim to improve our understanding of the wind morphology of eccentric AGB binary systems by investigating the mechanism responsible for the different small-scale structures and global morphologies that arise in a polytropic wind with different velocities. Using the smoothed particle hydrodynamics (SPH) code Phantom, we generate nine different high-resolution, 3D simulations of an AGB star with a solar-mass companion with various wind velocity and eccentricity combinations. The models assume a polytropic gas, with no additional cooling. We conclude that for models with a high wind velocity, the short interaction with the companion results in a regular spiral morphology, that is flattened. In the case of a lower wind velocity, the stronger interaction results in the formation of a high-energy region and bow-shock structure that can shape the wind into an irregular morphology if instabilities arise. High-eccentricity models show a complex, phase-dependent interaction leading to wind structures that are irregular in three dimensions. However, the significant interaction with the companion compresses matter into an equatorial density enhancement, irrespective of eccentricity.
△ Less
Submitted 2 July, 2021;
originally announced July 2021.
-
SPH modelling of companion-perturbed AGB outflows including a new morphology classification scheme
Authors:
S. Maes,
W. Homan,
J. Malfait,
L. Siess,
J. Bolte,
F. De Ceuster,
L. Decin
Abstract:
Asymptotic giant branch (AGB) stars are known to lose a significant amount of mass by a stellar wind, which controls the remainder of their stellar lifetime. High angular-resolution observations show that the winds of these cool stars typically exhibit mid- to small-scale density perturbations such as spirals and arcs, believed to be caused by the gravitational interaction with a (sub-)stellar com…
▽ More
Asymptotic giant branch (AGB) stars are known to lose a significant amount of mass by a stellar wind, which controls the remainder of their stellar lifetime. High angular-resolution observations show that the winds of these cool stars typically exhibit mid- to small-scale density perturbations such as spirals and arcs, believed to be caused by the gravitational interaction with a (sub-)stellar companion. We aim to explore the effects of the wind-companion interaction on the 3D density and velocity distribution of the wind, as a function of three key parameters: wind velocity, binary separation and companion mass. For the first time, we compare the impact on the outflow of a planetary companion to that of a stellar companion. We intend to devise a morphology classification scheme based on a singular parameter. With our grid of models we cover the prominent morphology changes in a companion-perturbed AGB outflow: slow winds with a close, massive binary companion show a more complex morphology. Additionally, we prove that massive planets are able to significantly impact the density structure of an AGB wind. We find that the interaction with a companion affects the terminal velocity of the wind, which can be explained by the gravitational slingshot mechanism. We distinguish between two types of wind focussing to the orbital plane resulting from distinct mechanisms: global flattening of the outflow as a result of the AGB star's orbital motion and the formation of an EDE as a consequence of the companion's gravitational pull. We investigate different morphology classification schemes and uncover that the ratio of the gravitational potential energy density of the companion to the kinetic energy density of the AGB outflow yields a robust classification parameter for the models presented in this paper.
△ Less
Submitted 1 July, 2021;
originally announced July 2021.
-
Atomium: The astounding complexity of the near circumstellar environment of the M-type AGB star R Hydrae. I. Morpho-kinematical interpretation of CO and SiO emission
Authors:
Ward Homan,
Bannawit Pimpanuwat,
Fabrice Herpin,
Taissa Danilovich,
Iain McDonald,
Sofia H. J. Wallström,
Anita M. S. Richards,
Alain Baudry,
Raghvendra Sahai,
Tom J. Millar,
Alex de Koter,
C. A. Gottlieb,
Pierre Kervella,
Miguel Montargès,
Marie Van de Sande,
Leen Decin,
Albert Zijlstra,
Sandra Etoka,
Manali Jeste,
Holger S. P. Müller,
Silke Maes,
Jolien Malfait,
Karl Menten,
John Plane,
Kelvin Lee
, et al. (14 additional authors not shown)
Abstract:
Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we…
▽ More
Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we observed the M-type asymptotic giant branch (AGB) star R Hydrae with ALMA. The morphology of the inner wind of R Hya, which has a known companion at ~3500 au, was determined from maps of CO and SiO obtained at high angular resolution. A map of the CO emission reveals a multi-layered structure consisting of a large elliptical feature at an angular scale of ~10'' that is oriented along the north-south axis. The wind morphology within the elliptical feature is dominated by two hollow bubbles. The bubbles are on opposite sides of the AGB star and lie along an axis with a position angle of ~115 deg. Both bubbles are offset from the central star, and their appearance in the SiO channel maps indicates that they might be shock waves travelling through the AGB wind. An estimate of the dynamical age of the bubbles yields an age of the order of 100 yr, which is in agreement with the previously proposed elapsed time since the star last underwent a thermal pulse. When the CO and SiO emission is examined on subarcsecond angular scales, there is evidence for an inclined, differentially rotating equatorial density enhancement, strongly suggesting the presence of a second nearby companion. The position angle of the major axis of this disc is ~70 deg in the plane of the sky. We tentatively estimate that a lower limit on the mass of the nearby companion is ~0.65 Msol on the basis of the highest measured speeds in the disc and the location of its inner rim at ~6 au from the AGB star.
△ Less
Submitted 15 April, 2021;
originally announced April 2021.
-
Magritte, a modern software library for 3D radiative transfer: II. Adaptive ray-tracing, mesh construction and reduction
Authors:
Frederik De Ceuster,
Jan Bolte,
Ward Homan,
Silke Maes,
Jolien Malfait,
Leen Decin,
Jeremy Yates,
Peter Boyle,
James Hetherington
Abstract:
Radiative transfer is a notoriously difficult and computationally demanding problem. Yet, it is an indispensable ingredient in nearly all astrophysical and cosmological simulations. Choosing an appropriate discretization scheme is a crucial part of the simulation, since it not only determines the direct memory cost of the model but also largely determines the computational cost and the achievable…
▽ More
Radiative transfer is a notoriously difficult and computationally demanding problem. Yet, it is an indispensable ingredient in nearly all astrophysical and cosmological simulations. Choosing an appropriate discretization scheme is a crucial part of the simulation, since it not only determines the direct memory cost of the model but also largely determines the computational cost and the achievable accuracy. In this paper, we show how an appropriate choice of directional discretization scheme as well as spatial model mesh can help alleviate the computational cost, while largely retaining the accuracy. First, we discuss the adaptive ray-tracing scheme implemented in our 3D radiative transfer library Magritte, that adapts the rays to the spatial mesh and uses a hierarchical directional discretization based on HEALPix. Second, we demonstrate how the free and open-source software library Gmsh can be used to generate high quality meshes that can be easily tailored for Magritte. In particular, we show how the local element size distribution of the mesh can be used to optimise the sampling of both analytically and numerically defined models. Furthermore, we show that when using the output of hydrodynamics simulations as input for a radiative transfer simulation, the number of elements in the input model can often be reduced by an order of magnitude, without significant loss of accuracy in the radiation field. We demonstrate this for two models based on a hierarchical octree mesh resulting from adaptive mesh refinement (AMR), as well as two models based on smoothed-particle hydrodynamics (SPH) data.
△ Less
Submitted 30 November, 2020;
originally announced November 2020.