Pith. sign in

REVIEW 4 major objections 4 minor 69 references

Impact of a binary companion in AGB outflows on CO spectral lines

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

Pith's one-line read Low-J CO lines from an AGB wind carry a binary companion's imprint — double peaks face-on, a central bump edge-on — that can act as a binary diagnostic, though beams and noise can hide them.

desk verdict Competent, honest 3D SPH+NLTE study of binary AGB line profiles; the qualitative binary diagnostics are plausible, but the specific line shapes rest on the free-wind velocity field and should be treated as conditional. read the letter →

arxiv 2507.03583 v2 pith:SQFMQ4D2 submitted 2025-07-04 astro-ph.SR

classification astro-ph.SR
keywords AGBstarsbinarycompanionsCOspectrallinescircumstellaroutflowssyntheticlineprofilesNLTEradiativetransfersmoothedparticlehydrodynamicsphotodissociation
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

This paper claims that a binary companion orbiting an AGB star leaves characteristic, identifiable marks on the low-J CO spectral lines of the star's wind: instead of the smooth parabolic or flat-topped profiles expected from a spherical outflow, the lines show two peaks near the terminal velocity when viewed face-on and a pronounced bump near the central velocity when viewed edge-on. The claim is built on a grid of nine 3D hydrodynamic models, spanning outflow velocities of 5, 10, and 20 km/s and orbital separations of 9, 15, and 25 au, post-processed with full non-local thermodynamic equilibrium radiative transfer. If the claim is right, molecular line profiles become a binary diagnostic for the many AGB outflows that are observed only through spectra, not resolved images. The paper adds a caution: the same features can be concealed by the telescope beam and by noise, so binaries can be misclassified as single stars, and standard single-star modelling of their lines yields systematically wrong mass-loss rates, around a factor of two in the cases studied.

What carries the argument

The load-bearing mechanism is the geometry of velocity bins: in a spherical outflow the projected-velocity bins are conical and roughly equal in volume, producing flat-topped or parabolic lines, whereas the companion's gravity and the star's orbital motion distort the velocity field so that certain bins become arc-shaped, enlarged, or double-valued in projected velocity. Those distortions concentrate or deplete emission at specific velocities, creating the two peaks near the terminal velocity and the central bump. The machinery has three parts: the smoothed particle hydrodynamics code Phantom, which evolves a free wind (launched at the dust-condensation radius with fixed speed) in the binary potential with a polytropic equation of state and H I cooling; a new ray-tracing CO photodissociation calculation that generalises the spherical method of Groenewegen (2017) to 3D using HEALPix-distributed rays and nearest-ray interpolation, setting the size and shape of the emitting envelope; and the 3D NLTE line radiative transfer code Magritte, which solves for the level populations of the first 40 CO rotational transitions and produces the synthetic lines at chosen inclinations and position angles.

What would settle it

Observe a confirmed binary AGB system whose orbital parameters are known from astrometry or resolved imaging, with a slow wind (≈5 km/s) and a close orbit (≈9 au), and measure its low-J CO lines at high signal-to-noise with beam effects removed: if the lines do not show the predicted double-peaked (face-on) or central-bump (edge-on) structure across several transitions, the velocity-bin mechanism is falsified. A positive control would be a known face-on narrow-spiral system, which the models say should show the two-peaked profile rather than a smooth parabola.

Watch

Extended reading notes

Core claim

The central discovery is that the companion's gravitational interaction reshapes both the density and the velocity field of the outflow, and the reshaped velocity field translates directly into specific low-J CO line morphologies. All nine models develop a pronounced spiral in the orbital plane, with the meridional structure, the degree of global flattening, and the emission shapes varying with orbital separation and wind speed; a novel 3D photodissociation treatment, generalising the spherical scheme of Groenewegen (2017), sets a non-spherical emitting-envelope size that closely follows the morphology. The synthetic lines deviate strongly from parabolic or flat-topped profiles, with the two most pronounced features being two peaks near the terminal velocity (viewed face-on) and a bump near the central velocity (viewed edge-on), whose presence, strength, and position depend on inclination, position angle, and the CO transition. In specific cases the profiles appear quasi-parabolic, hiding the companion entirely. The authors' conclusion is that molecular line profiles can serve as a binary diagnostic, but the companion's influence can go unnoticed because the features are easily concealed by the beam profile and observational noise, making it easy to misclassify such systems as single stars and to derive systematically biased mass-loss rates.

Load-bearing premise

The central claim rests on the free-wind approximation, in which the primary's gravity is artificially balanced by radiation force and wind particles are launched at the dust-condensation radius with a fixed velocity; the paper itself states (Sect. 4.5) that the resulting inner-wind velocity profile is unrealistic, and the wind speed at the companion is exactly what controls the morphology that creates the predicted line features.

Editorial extensions

If this is right

  • Observed low-J CO lines that deviate from parabolic or flat-topped profiles — two peaks near the terminal velocity or a bump near the central velocity — can be read as indirect evidence of a hidden companion even when the outflow is not spatially resolved.
  • Mass-loss rates derived by fitting single-star spherical models to such lines can be off by around a factor of two, with the error depending on inclination, so surveys that ignore binarity carry a systematic bias.
  • Observing several CO transitions together helps lift the degeneracy between orbital separation and wind velocity, because the line shape and strength evolve differently with J for different binary parameters.
  • In wide-separation, fast-wind models the higher-J lines probe the spiral structure directly, and the line shape's dependence on position angle becomes a potential way to locate the companion.
  • The terminal velocity inferred from the lines depends on orbital parameters and inclination rather than simply the wind speed, so companion shaping can masquerade as a different wind speed in single-star fits.

Reading between the lines

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

  • Re-examining existing single-dish CO spectra of 'single' AGB stars for the specific double-peak and central-bump morphologies could reveal a population of hidden companions; the grid predicts the strongest signatures for slow winds (5 km/s) and close orbits (9-15 au).
  • The same velocity-bin argument implies that other velocity-field distortions — pulsations, convection cells, eccentric orbits — could produce line shapes resembling the binary signatures, so the diagnostic is most secure when combined with resolved imaging or multi-transition data.
  • The 3D photodissociation scheme could be extended to other molecules (e.g., 13CO, HCN) or to higher mass-loss rates, where non-spherical UV shielding should be even more pronounced and the line-shape effects correspondingly stronger.
  • If the diagnostic holds, comparing the inferred incidence of hidden companions in CO line surveys with the known binarity fraction of AGB progenitors would offer an observational test of whether companions are indeed the dominant shaping mechanism of aspherical mass loss.
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 / 4 minor

Summary. The paper presents a grid of nine 3D SPH models of an AGB star with a 1 M_sun companion at orbital separations of 9, 15, and 25 au, using free-wind injection velocities of 5, 10, and 20 km/s. The models are post-processed with the NLTE line radiative transfer code Magritte to compute synthetic low-J CO lines at multiple inclinations and position angles, including a 3D treatment of CO photodissociation. The authors report non-spherical morphologies with spirals and global flattening, and show that the synthetic line profiles can deviate strongly from the parabolic or flat-topped shapes expected from spherical winds, with the most pronounced features being two peaks near the terminal velocity in face-on views and a central bump near the central velocity in edge-on views. They also convolve the lines with Gaussian beams, add noise, and use the De Beck et al. (2010) fitting formula to estimate retrieved mass-loss rates. The paper concludes that molecular line profiles can serve as a binary diagnostic, but that the characteristic features can be hidden by beam and noise effects, potentially leading to misclassification and systematic errors in derived mass-loss rates.

Significance. If the qualitative conclusions hold, this is a useful step toward connecting binary-induced outflow morphologies with the most commonly observed molecular line observations of AGB stars. The paper's strengths include the use of open-source, reproducible codes (Phantom, Magritte), a validation of the radiative transfer against SKIRT in Appendix A (roughly 1% agreement for J=3-2 and 5-4, 5-10% for J=1-0), a validation of the 3D CO photodissociation interpolation in Appendix B, and an explicit treatment of beam convolution and noise. The authors are also commendably transparent about the limitations of their wind-driving prescription in Sect. 4.5. However, the central diagnostic signatures are tightly linked to a simplified free-wind velocity field, the resolution study is only qualitative, and the mass-loss-bias argument is not yet controlled. The manuscript therefore establishes a promising framework and a set of well-documented model predictions, but the specific claim that low-J CO line profiles can serve as a binary diagnostic still needs additional quantitative support.

major comments (4)
  1. [§2.1, §3.2.1, §4.5]
  2. [§2.2, Table 1]
  3. [§4.1, Table 3]
  4. [§3.2.5, Fig. 13]
minor comments (4)
  1. [§3.2.5]
  2. [§4.2]
  3. [§3.2.1]
  4. [§3.2.5]

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the synthetic lines are forward-modeled from binary hydrodynamics, with a single-star control and an independent SKIRT benchmark; no prediction reduces to a fitted input.

full rationale

The central claim, that a binary companion can imprint double peaks and central bumps on low-J CO lines, is obtained by running Phantom hydrodynamics with a binary companion, post-processing with Magritte, and comparing against a single-star Phantom control (Sect. 3.2.1, Fig. 5). The spectral features arise from velocity-bin geometries created by the companion's gravitational interaction; they are emergent outputs, not quantities fitted to the model inputs. The radiative transfer is benchmarked against the independent code SKIRT in Appendix A, with about 1% agreement for the J=3-2 and J=5-4 lines, so the line-formation step is externally checked. The mass-loss retrieval in Sect. 4.1 uses the external spherical fitting function of De Beck et al. (2010) only to compare relative trends between models and inclinations; it is not used to construct the predicted line shapes, so there is no fitted-input-called-prediction pattern. The paper does contain many self-citations, including Maes et al. (2021), Malfait et al. (2021, 2024a,b), Siess et al. (2022), and Esseldeurs et al. (2023), but these support simulation setup, cooling prescriptions, and prior morphology classifications rather than the load-bearing spectral-line claim itself, so they are not circular. The explicit limitation in Sect. 4.5, 'we neglect the radiative acceleration in our models... the velocity profile in our simulations is also unrealistic,' is weighed here: it is a genuine physical-modeling risk that could alter the predicted double peaks and central bumps, but it is an acknowledged approximation rather than a circular definition or a renamed input. No equation in the paper reduces to its own inputs by construction, so the appropriate circularity score is 0.

Assumptions & free parameters 7 free parameters · 8 assumptions · 0 invented entities

The central predictions rest on a chain of modeling choices: free wind driving, no dust, simplified cooling, assumed stellar parameters and CO abundance, and a fixed low mass-loss rate. None of these are fitted to observations, so the circularity burden is low, but they are unverified inputs that could shift the quantitative line shapes and flux levels. The grid does include the intended parameter-space variables (wind velocity and orbital separation), so those are not counted as free parameters.

free parameters (7)
  • Stellar effective temperature T* = 2750 K
    Chosen arbitrarily (Sect. 2.2) as representative of AGB stars; sets the temperature power law and the dust condensation radius.
  • Stellar radius R* = 1.267 au
    Chosen arbitrarily (Sect. 2.2); sets the inner boundary and the wind injection radius together with T*.
  • Mass-loss rate Mdot = 1e-7 Msun/yr
    Fixed for all models (Sect. 2.2); controls density and the CO photodissociation radius and is key to line fluxes.
  • CO fractional abundance = 1.5e-4 relative to H nuclei (3e-4 relative to H2)
    Assumed constant for oxygen-rich AGB stars (Sect. 2.3); directly controls CO column densities and line optical depths.
  • Secondary accretion radius = 0.2 au
    Enlarged beyond physical expectation to suppress numerical instabilities at the adopted resolution (Sect. 2.2); authors argue it does not affect line profiles.
  • Turbulent velocity = 1 km/s
    Constant microturbulence added in Magritte (Sect. 2.3); broadens the synthetic lines.
  • Interstellar radiation field strength chi = 1 (Draine field)
    Assumed in the CO photodissociation calculation (Sect. 2.4); sets the destruction rate and the size of the CO envelope.
assumptions (8)
  • domain assumption Free-wind approximation: radiation pressure is replaced by an artificial force balancing gravity, and the wind is injected at a fixed velocity from the dust condensation radius.
    Invoked in Sect. 2.2; determines the velocity field that shapes the line profiles.
  • domain assumption Polytropic equation of state P=(gamma-1)rho u with computed gamma, plus Hi cooling only.
    Sect. 2.1; simplified thermodynamics ignores many molecular and dust cooling channels that affect temperature and level populations.
  • domain assumption No dust is included in the hydrodynamic models or radiative transfer.
    Sect. 2.2 and 4.5; dust affects wind driving, IR pumping, and attenuation of CO lines.
  • domain assumption CO is destroyed only by external UV photodissociation; no molecule formation is included.
    Sect. 2.4, Eq. 17; sets the CO abundance distribution in the envelope.
  • domain assumption H and H2 are in equilibrium chemistry at the local gas temperature.
    Sect. 2.3, Eqs. 12-14; determines H2 abundance and shielding.
  • domain assumption Circular orbit with zero eccentricity for the binary.
    Sect. 2.2; eccentricity would add asymmetry and is listed as future work.
  • ad hoc to paper Low-resolution SPH runs are sufficient to capture inner-wind morphology.
    Sect. 2.2 asserts a test with higher resolution shows no significant differences, but no quantitative convergence evidence is presented.
  • domain assumption Spherically-symmetric power-law temperature and beta-law velocity are assumed in the region between the stellar surface and the injection radius.
    Sect. 2.3, Eqs. 15-16; this region contributes to the lines, especially higher J transitions.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Impact of a binary companion in AGB outflows on CO spectral lines." pith.science (2026). https://pith.science/paper/SQFMQ4D2

@misc{pith2026250703583,
  author       = {Pith},
  title        = {Pith review of: Impact of a binary companion in AGB outflows on CO spectral lines},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SQFMQ4D2}},
  note         = {Machine review of arXiv:2507.03583}
}
read the original abstract

In the late stage of their evolution, low- to intermediate-mass stars pass through the asymptotic giant branch (AGB) phase, characterised by strong mass loss through dust driven winds. High angular resolution observations reveal that these winds harbour strong deviations from spherical symmetry, such as spirals and arcs, believed to be caused by hidden (sub-)stellar companions. Much more often, one observes spectral lines, where the presence of a companion is less clear. We study the impact of a binary companion on low-J CO spectral lines of AGB star outflows. By varying the orbital separation and wind velocity, we aim to find line shapes characteristic of more complex binary-induced morphologies. We generated a grid of nine 3D models of a mass-losing AGB star using the smoothed particle hydrodynamics code Phantom, with three values for both the outflow velocity and orbital separation. Utilising the radiative transfer code Magritte, we created synthetic spectral lines for the low rotational transitions of CO at different inclinations and position angles. Our simulations show a variety of morphologies, always with a pronounced spiral structure arising in the orbital plane, but with varying shapes in the meridional plane, and different degrees of global flattening. We find that the CO line profiles can deviate strongly from the parabolic or flat-topped profiles expected from spherically symmetric outflows. A variety of line shapes emerge, with two peaks near the terminal velocity, and a central bump near the central velocity being the most pronounced. In specific cases, the spectral lines can appear parabolic, hiding the presence of a binary companion. We find the CO spectral lines can serve as a binary diagnostic. The influence of the companion on the line can however also go easily unnoticed, as the features can be concealed by the beam profile and the noise of the observations.

Figures

Figures reproduced from arXiv: 2507.03583 by the authors.

Figure 1
Figure 1. Density distributions in slices through the orbital plane (upper row) and the meridional plane (lower row) for the models with initial wind velocities of 10 km s−1 . From left to right, the orbital separation is 9, 15, and 25 au. models in slices through both the orbital (x − y) and meridional (x−z) plane for the inner regions in [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Number density of CO in slices through the meridional plane after accounting for photodissociation. The solid and dotted lines represent the location where x(r) is 0.5 and 0.01, respectively. abundance of CO with respect to H2, compared to the initial frac￾tional abundance, is 0.5 and 0.01, respectively. We always show a slice through the x−z plane, as in the orbital plane the CO map is circular in all models. In th… view at source ↗
Figure 3
Figure 3. Synthetic CO J = 2→1 spectral lines for all models. The columns correspond, from left to right, to initial outflow velocities of 5, 10, and 20 km s−1 . The solid, dashed, and dotted lines represent orbital separations of 9, 15, and 25 au. The spectral lines are creating viewing face-on at an inclination of 0° (upper row) and edge-on at 90° (lower row). sitions as they are only excited within this region, and the out… view at source ↗
Figures from the paper (9 more)
Figure 4
Figure 4. Figure 4: Velocities in the z direction (upper row) and in the x direction (lower row) in slices through the meridional plane for the models with initial wind velocities of 10 km s−1 . Contours are shown on the plots around velocities interesting for the line profiles, indicated…
Figure 5
Figure 5. Figure 5: Synthetic CO spectral lines of the first six transitions from a single-star Phantom model, with a mass-loss rate of 10−7 M⊙ yr−1 and an initial wind velocity of 10 km s−1 . When comparing the different spectral lines more closely, we see that the flux changes significa…
Figure 6
Figure 6. Figure 6: Synthetic CO J = 2→1 spectral lines for the models with initial wind velocities of 10 km s−1 , viewed under inclinations of 18°, 36°, 54°, and 72°. The solid, dashed, and dotted lines represent orbital separations of 9, 15, and 25 au. 3.2.3. Different CO transitions Th…
Figure 7
Figure 7. Figure 7: Synthetic CO spectral lines of the first six rotational transitions, for the models with initial outflow velocities of 5 km s−1 , viewed under an inclination of 0°. The solid, dashed, and dotted lines represent orbital separations of 9, 15, and 25 au. and meridional pl…
Figure 11
Figure 11. Figure 11: Temperature distributions in slices through the orbital plane (upper row) and the meridional plane (lower row) for the v05a09 (left) and v20a25 (right) models. −10 −5 0 5 10 Velocity [km s−1 ] 0 1 2 3 4 Flux [Jy] v05a09 −20 0 20 Velocity [km s−1 ] 0.0 0.1 0.2 0.3 0.4 …
Figure 12
Figure 12. Figure 12: Synthetic CO J = 6→5 spectral lines for the v05a09 (left) and v20a25 (right) models, viewed edge-on and at four different values of the PA (ϕ). The different colours correspond to the different values of the PA. 3.2.5. Single dish telescopes When observing AGB outflow…
Figure 10
Figure 10. Figure 10: Synthetic CO J = 1→0 (top) and J = 6→5 (bottom) spec￾tral lines for the v05 (left) and v10 (right) models viewed under an in￾clination of 90°. The solid, dashed, and dotted lines represent orbital separations of 9, 15, and 25 au. tion and the strength of the stronger …
Figure 13
Figure 13. Figure 13: Synthetic CO J = 2→1 spectral lines for the models with initial wind velocities of 10 km s−1 , where a beam with five different values of θb has been applied. The spectral lines are viewed face-on at an inclination of 0° and edge-on at 90° in the upper and lower row, …
Figure 14
Figure 14. Figure 14: Synthetic CO spectral lines for three selected models. From left to right these are v05a09, v10a25, and v20a25 with the J = 2→1, 4→3, 6→5 transitions, respectively. The spectral lines are created viewing face-on (upper row) and edge-on (lower row). The black dashed li…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

69 extracted references · 57 canonical work pages

  1. [1]

    2023, A&A, 669, A49

    Ahmad, A., Freytag, B., & Höfner, S. 2023, A&A, 669, A49

  2. [2]

    Andriantsaralaza, M., Ramstedt, S., Vlemmings, W. H. T., & De Beck, E. 2022, A&A, 667, A74

  3. [3]

    & Mohamed, S

    Aydi, E. & Mohamed, S. 2022, MNRAS, 513, 4405

  4. [4]

    2022, A&A, 666, A101

    Baes, M., Camps, Peter, & Matsumoto, Kosei. 2022, A&A, 666, A101

  5. [5]

    I., Dejonghe, H., et al

    Baes, M., Davies, J. I., Dejonghe, H., et al. 2003, MNRAS, 343, 1081

  6. [6]

    2022, A&A, 666, A40

    Bollen, D., Kamath, D., Van Winckel, H., et al. 2022, A&A, 666, A40

  7. [7]

    2024, Astronomy and Com- puting, 49, 100889

    Ceulemans, T., De Ceuster, F., Decin, L., & Yates, J. 2024, Astronomy and Com- puting, 49, 100889

  8. [8]

    G., Nordhaus, J., & Carroll-Nellenback, J

    Chen, Z., Frank, A., Blackman, E. G., Nordhaus, J., & Carroll-Nellenback, J. 2017, MNRAS, 468, 4465–4477

Show all 69 references
  1. [9]

    2020, AJ, 892, 110

    Chen, Z., Ivanova, N., & Carroll-Nellenback, J. 2020, AJ, 892, 110

  2. [10]

    2015, A&A, 581, A60 De Beck, E., Decin, L., de Koter, A., et al

    Danilovich, T., Teyssier, D., Justtanont, K., et al. 2015, A&A, 581, A60 De Beck, E., Decin, L., de Koter, A., et al. 2010, A&A, 523, A18 De Beck, E., Lombaert, R., Agúndez, M., et al. 2012, A&A, 539, A108 De Ceuster, F., Bolte, J., Homan, W., et al. 2020, MNRAS, 499, 5194 De ...

  3. [11]

    2021, Annual Review of Astronomy and Astrophysics, 59, 337

    Decin, L. 2021, Annual Review of Astronomy and Astrophysics, 59, 337

  4. [12]

    2019, Nature Astronomy, 3, 408–415

    Decin, L., Homan, W., Danilovich, T., et al. 2019, Nature Astronomy, 3, 408–415

  5. [13]

    2006, A&A, 456, 549

    Decin, L., Hony, S., de Koter, A., et al. 2006, A&A, 456, 549

  6. [14]

    Decin, L., Montargès, M., Richards, A. M. S., et al. 2020, Science, 369, 1497–1500 Díaz-Luis, J. J., Alcolea, J., Bujarrabal, V ., et al. 2019, A&A, 629, A94

  7. [15]

    Draine, B. T. 1978, ApJS, 36, 595 El Mellah, I., Bolte, J., Decin, L., Homan, W., & Keppens, R. 2020, A&A, 637, A91

  8. [16]

    2023, A&A, 674, A122

    Esseldeurs, M., Siess, L., De Ceuster, F., et al. 2023, A&A, 674, A122

  9. [17]

    & Höfner, S

    Freytag, B. & Höfner, S. 2023, A&A, 669, A155

  10. [18]

    Fulton, B. J. & Petigura, E. A. 2018, AJ, 156, 264

  11. [19]

    2006, Proceedings of the International Astronomical Union Symposium Górski, K

    Giesen, T., Brünken, S., Caris, M., et al. 2006, Proceedings of the International Astronomical Union Symposium Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759

  12. [20]

    Groenewegen, M. A. T. 2017, A&A, 606, A67

  13. [21]

    & Olofsson, H

    Habing, H. & Olofsson, H. 2004, Asymptotic Giant Branch Stars (Springer New

  14. [22]

    & Olofsson, H

    York, NY) Höfner, S. & Olofsson, H. 2018, A&A Rev., 26, 1

  15. [23]

    2016, A&A, 596, A91

    Homan, W., Boulangier, J., Decin, L., & de Koter, A. 2016, A&A, 596, A91

  16. [24]

    2020, A&A, 642, A93

    Homan, W., Cannon, E., Montargès, M., et al. 2020, A&A, 642, A93

  17. [25]

    2015, A&A, 579, A118

    Homan, W., Decin, L., de Koter, A., et al. 2015, A&A, 579, A118

  18. [26]

    2018, A&A, 616, A34

    Homan, W., Richards, Anita, Decin, Leen, de Koter, Alex, & Kervella, Pierre. 2018, A&A, 616, A34

  19. [27]

    & Morris, M

    Jura, M. & Morris, M. 1981, ApJ, 251, 181

  20. [28]

    & Olofsson, H

    Kerschbaum, F. & Olofsson, H. 1999, Astron. Astrophys. Suppl. Ser., 138, 299

  21. [29]

    Kervella, P., Homan, W., Richards, A. M. S., et al. 2016, A&A, 596, A92

  22. [30]

    Kim, H., Liu, S.-Y ., & Taam, R. E. 2019, The Astrophysical Journal Supplement Series, 243, 35

  23. [31]

    2013, Stellar Structure and Evolution (Springer Berlin, Heidelberg)

    Kippenhahn, R., Weigert, A., & Weiss, A. 2013, Stellar Structure and Evolution (Springer Berlin, Heidelberg)

  24. [32]

    I., & Lugaro, M

    Kobayashi, C., Karakas, A. I., & Lugaro, M. 2020, ApJ, 900, 179

  25. [33]

    Lamers, H. J. G. L. M. & Cassinelli, J. P. 1999, Introduction to Stellar Winds (Cambridge University Press)

  26. [34]

    2022, ApJ, 931, 142

    Lee, Y .-M., Kim, H., & Lee, H.-W. 2022, ApJ, 931, 142

  27. [35]

    2014, A&A, 568

    Li, X., Millar, T., Walsh, C., Heays, A., & Dishoeck, E. 2014, A&A, 568

  28. [36]

    2024, A&A, 687, A112

    Maercker, M., De Beck, E., Khouri, T., et al. 2024, A&A, 687, A112

  29. [37]

    Maercker, M., Mohamed, S., Vlemmings, W. H. T., et al. 2012, Nature, 490, 232

  30. [38]

    Maercker, M., Vlemmings, W. H. T., Brunner, M., et al. 2016, A&A, 586, A5

  31. [39]

    2024, AJ, 969, 79

    Maes, S., De Ceuster, F., Van de Sande, M., & Decin, L. 2024, AJ, 969, 79

  32. [40]

    2021, A&A, 653

    Maes, S., Homan, W., Malfait, J., et al. 2021, A&A, 653

  33. [41]

    2021, A&A, 652

    Malfait, J., Homan, W., Maes, S., et al. 2021, A&A, 652

  34. [42]

    A., Glassgold, A

    Mamon, G. A., Glassgold, A. E., & Huggins, P. J. 1988, ApJ, 328, 797

  35. [43]

    & Morris, M

    Mastrodemos, N. & Morris, M. 1998, ApJ, 497, 303

  36. [44]

    & Morris, M

    Mastrodemos, N. & Morris, M. 1999, ApJ, 523, 357

  37. [45]

    2023, A&A, 678, A175

    Matsumoto, K., Camps, P., Baes, M., et al. 2023, A&A, 678, A175

  38. [46]

    & Huggins, P

    Mauron, N. & Huggins, P. J. 2006, A&A, 452, 257

  39. [47]

    Merrill, P. W. 1952, ApJ, 116, 21

  40. [48]

    & Di Stefano, R

    Moe, M. & Di Stefano, R. 2017, ApJS, 230, 15

  41. [49]

    & Jura, M

    Morris, M. & Jura, M. 1983, ApJ, 264, 546

  42. [50]

    1985, A&A, 142, 107

    Morris, M., Lucas, R., & Omont, A. 1985, A&A, 142, 107

  43. [51]

    2006, in IAU Symposium, V ol

    Morris, M., Sahai, R., Matthews, K., et al. 2006, in IAU Symposium, V ol. 234, Planetary Nebulae in our Galaxy and Beyond, ed. M. J. Barlow & R. H. Mén- dez, 469–470

  44. [52]

    Olofsson, H., González Delgado, D., Kerschbaum, F., & Schöier, F. L. 2002, A&A, 391, 1053

  45. [53]

    J., Wurster, J., Tricco, T

    Price, D. J., Wurster, J., Tricco, T. S., et al. 2018, PASA, 35, e031

  46. [54]

    2006, A&A, 454

    Ramstedt, S., Schoeier, F., Olofsson, H., & Lundgren, A. 2006, A&A, 454

  47. [55]

    L., Olofsson, H., & Lundgren, A

    Ramstedt, S., Schöier, F. L., Olofsson, H., & Lundgren, A. A. 2008, A&A, 487, 645

  48. [56]

    Ramstedt, S., Vlemmings, W. H. T., Doan, L., et al. 2020, A&A, 640, A133

  49. [57]

    Saberi, M., Vlemmings, W. H. T., & De Beck, E. 2019, A&A, 625, A81

  50. [58]

    Saladino, M. I. & Pols, O. R. 2019, A&A, 629, A103

  51. [59]

    I., Pols, O

    Saladino, M. I., Pols, O. R., & Abate, C. 2019, A&A, 626, A68

  52. [60]

    I., Pols, O

    Saladino, M. I., Pols, O. R., van der Helm, E., Pelupessy, I., & Portegies Zwart, S. 2018, A&A, 618, A50 Schöier, F., van der Tak, F., van Dishoeck, E., & Black, J. 2010, LAMDA: Leiden Atomic and Molecular Database, Astrophysics Source Code Library, record ascl:1010.077 Schöie...

  53. [61]

    & Härm, R

    Schwarzschild, M. & Härm, R. 1965, ApJ, 142, 855

  54. [62]

    Siess, L., Homan, W., Toupin, S., & Price, D. J. 2022, A&A, 667, A75

  55. [63]

    1978, Physical processes in the interstellar medium (Wiley)

    Spitzer, L. 1978, Physical processes in the interstellar medium (Wiley)

  56. [64]

    Teyssier, D., Hernandez, R., Bujarrabal, V ., Yoshida, H., & Phillips, T. G. 2006, A&A, 450, 167

  57. [65]

    & Jorissen, A

    Theuns, T. & Jorissen, A. 1993, MNRAS, 265, 946

  58. [66]

    Tielens, A. G. G. M. 2005, The Physics and Chemistry of the Interstellar Medium (Cambridge University Press) van Winckel, H. 2003, ARA&A, 41, 391

  59. [67]

    F., & Black, J

    Visser, R., van Dishoeck, E. F., & Black, J. H. 2009, A&A, 503, 323

  60. [68]

    2003, Dust in the Galactic Environment, 2nd edn., Series in astron- omy and astrophysics (IOP)

    Whittet, D. 2003, Dust in the Galactic Environment, 2nd edn., Series in astron- omy and astrophysics (IOP)

  61. [69]

    C., Balakrishnan, N., & Forrey, R

    Yang, B., Stancil, P. C., Balakrishnan, N., & Forrey, R. C. 2010, ApJ, 718, 1062 Article number, page 17 of 22 A&A proofs: manuscript no. aa54905-25corr Appendix A: Comparison with SKIRT To further validate and check the accuracy of our results, we compare them by generating s...

Pith tools

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