The J = 5.5 → 4.5 and J = 5 → 4 transitions of PO and PN, respectively, have been imaged in the envelope of hypergiant star VY Canis Majoris (VY CMa) using the Atacama Large Millimeter/submillimeter ...Array with angular resolutions of 0.″2 and 1.″5 and data from the Submillimeter Telescope of the Arizona Radio Observatory. These maps are the first high-fidelity images of PO and PN in a circumstellar envelope. Both molecules are primarily present in a spherical, star-centered region with a radius ∼60 R* (0.″5), indicating formation by LTE chemistry and then condensation into grains. PN, however, shows additional, fan-shaped emission 2″ southwest of the star, coincident with dust features resolved by Hubble Space Telescope (HST), as well as four newly identified distinct structures 1″–2″ toward the north, east, and west (Cloudlets I–IV), not visible in HST images. The “SW Fan” and the cloudlets are also prominent in the J = 5.5 → 4.5 transition of NS. The correlation of PN with NS, SiO, and dust knots in the SW Fan suggests a formation in shocked gas enhanced with nitrogen. Excess nitrogen is predicted to favor PN synthesis over PO. Abundances for PN and PO in the spherical source are f ∼ 4.4 × 10−8 and 1.4 × 10−7, respectively, relative to H2. Given a cosmic abundance of phosphorus, an unusually high fraction (∼35%) is contained in PO and PN. Alternatively, the stellar winds may be enriched in P (and N) by dredge-up from the interior of VY CMa.
Context. Evolved low-mass stars lose a significant fraction of their mass through stellar winds. While the overall morphology of the stellar wind structure during the asymptotic giant branch (AGB) ...phase is thought to be roughly spherically symmetric, the morphology changes dramatically during the post-AGB and planetary nebula phase, during which bipolar and multi-polar structures are often observed. Aims. We aim to study the inner wind structure of the closest well-known AGB star CW Leo. Different diagnostics probing different geometrical scales have implied a non-homogeneous mass-loss process for this star: dust clumps are observed at milli-arcsec scale, a bipolar structure is seen at arcsecond-scale, and multi-concentric shells are detected beyond 1′′. Methods. We present the first ALMA Cycle 0 band 9 data around 650 GHz (450 μm) tracing the inner wind of CW Leo. The full-resolution data have a spatial resolution of 0.̋42 × 0.̋24, allowing us to study the morpho-kinematical structure of CW Leo within ~6′′. Results. We have detected 25 molecular emission lines in four spectral windows. The emission of all but one line is spatially resolved. The dust and molecular lines are centered around the continuum peak position, which is assumed to be dominated by stellar emission. The dust emission has an asymmetric distribution with a central peak flux density of ~2 Jy. The molecular emission lines trace different regions in the wind acceleration region and imply that the wind velocity increases rapidly from about 5 R⋆, almost reaching the terminal velocity at ~11 R⋆. The images prove that vibrational lines are excited close to the stellar surface and that SiO is a parent molecule. The channel maps for the brighter lines show a complex structure; specifically, for the 13CO J = 6–5 line, different arcs are detected within the first few arcseconds. The curved structure in the position-velocity (PV) map of the 13CO J = 6–5 line can be explained by a spiral structure in the inner wind of CW Leo, probably induced by a binary companion. From modelling the ALMA data, we deduce that the potential orbital axis for the binary system lies at a position angle of ~10–20° to the north-east and that the spiral structure is seen almost edge-on. We infer an orbital period of 55 yr and a binary separation of 25 au (or ~8.2 R⋆). We tentatively estimate that the companion is an unevolved low-mass main-sequence star. Conclusions. A scenario of a binary-induced spiral shell can explain the correlated structure seen in the ALMA PV images of CW Leo. Moreover, this scenario can also explain many other observational signatures seen at different spatial scales and in different wavelength regions, such as the bipolar structure and the almost concentric shells. ALMA data hence for the first time provide the crucial kinematical link between the dust clumps seen at milli-arcsecond scale and the almost concentric arcs seen at arcsecond scale.
Complex stellar winds from evolved stars
Stars less than eight times the mass of the Sun end their lives as planetary nebulae, structures of ionized gas thrown off by the star and heated by the ...exposed stellar core. Planetary nebulae are often bipolar in shape or contain complex morphological features such as rings or spirals. Decin
et al.
observed the stellar winds of 14 stars during their asymptotic giant branch (AGB) phase of stellar evolution, which immediately precedes the planetary nebula phase. They found morphologies in the AGB winds similar to planetary nebulae and demonstrated that they are produced by the influence of a binary companion on the AGB wind.
Science
, this issue p.
1497
Complex morphologies in the winds of asymptotic giant branch stars and planetary nebulae are produced by binary interactions.
Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their winds exhibit distinct nonspherical geometries with morphological similarities to planetary nebulae (PNe). We infer that the same physics shapes both AGB winds and PNe; additionally, the morphology and AGB mass-loss rate are correlated. These characteristics can be explained by binary interaction. We propose an evolutionary scenario for AGB morphologies that is consistent with observed phenomena in AGB stars and PNe.
Red supergiants are the most common final evolutionary stage of stars that have initial masses between 8 and 35 times that of the Sun
. During this stage, which lasts roughly 100,000 years
, red ...supergiants experience substantial mass loss. However, the mechanism for this mass loss is unknown
. Mass loss may affect the evolutionary path, collapse and future supernova light curve
of a red supergiant, and its ultimate fate as either a neutron star or a black hole
. From November 2019 to March 2020, Betelgeuse-the second-closest red supergiant to Earth (roughly 220 parsecs, or 724 light years, away)
-experienced a historic dimming of its visible brightness. Usually having an apparent magnitude between 0.1 and 1.0, its visual brightness decreased to 1.614 ± 0.008 magnitudes around 7-13 February 2020
-an event referred to as Betelgeuse's Great Dimming. Here we report high-angular-resolution observations showing that the southern hemisphere of Betelgeuse was ten times darker than usual in the visible spectrum during its Great Dimming. Observations and modelling support a scenario in which a dust clump formed recently in the vicinity of the star, owing to a local temperature decrease in a cool patch that appeared on the photosphere. The directly imaged brightness variations of Betelgeuse evolved on a timescale of weeks. Our findings suggest that a component of mass loss from red supergiants
is inhomogeneous, linked to a very contrasted and rapidly changing photosphere.
Context.
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.
Aims.
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.
Methods.
We ran a small grid of high-resolution polytropic models with the smoothed particle hydrodynamics (SPH) numerical code P
HANTOM
to examine the 3D density structure of the AGB outflow in the orbital and meridional plane and around the poles. By constructing a basic toy model of the gravitational acceleration due to the companion, we analysed the terminal velocity reached by the outflow in the simulations.
Results.
We find that models with a stellar companion, large binary separation and high wind speed obtain a wind morphology in the orbital plane consisting of a single spiral structure, of which the two edges diverge due to a velocity dispersion caused by the gravitational slingshot mechanism. In the meridional plane the spiral manifests itself as concentric arcs, reaching all latitudes. When lowering the wind velocity and/or the binary separation, the morphology becomes more complex: in the orbital plane a double spiral arises, which is irregular for the closest systems, and the wind material gets focussed towards the orbital plane, with the formation of an equatorial density enhancement (EDE) as a consequence. Lowering the companion mass from a stellar to a planetary mass, reduces the formation of density perturbations significantly.
Conclusions.
With this 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.
Abstract
The
J
= 2 → 1 transition of CO near 230 GHz and the
J
= 3 → 2 line of HCN at 265 GHz have been imaged in the envelope of the red hypergiant star, VY Canis Majoris (VY CMa), using the Atacama ...Large Millimeter Array (ALMA) with angular resolutions 0.″2–1.″5; single-dish data were added to provide sensitivity up to 30″. These images reveal a far more complex envelope, with previously unseen outflows extending 4″–9″ from the star. These new structures include an arc-like outflow with an angular separation of ∼9″ northeast from the stellar position (“NE Arc”), twin fingerlike features approximately 4″ to the north/northeast (“NE Extension”), and a roughly spherical region observed ∼7″ E of the star (“E Bubble”). The NE Arc appears to be decelerating from base (
V
LSR
∼ 7 km s
−1
) to tip (
V
LSR
∼ 18 km s
−1
), while the NE Extension is blueshifted with
V
LSR
∼ −7 km s
−1
. Among the new features, HCN is only detected in the NE Arc. In addition, known structures Arc 1, Arc 2, and NW Arc, as well as other features closer to the star, are closely replicated in CO, suggesting that the gas and dust are well mixed. The CO spectra are consistent with the kinematic picture of VY CMa derived from HST data. Arc 2, however, has added complexity. Preliminary results from CO suggest
12
C/
13
C ∼ 22–38 across the envelope. The additional presence of at least three major episodic mass ejection events significantly broadens the current perspective of the envelope structure and mass-loss history of VY CMa.
Context. Phosphorus-bearing compounds have only been studied in the circumstellar environments of the asymptotic giant branch star IRC +10 216 and the protoplanetary nebula CRL 2688, both carbon-rich ...objects, and the oxygen-rich red supergiant VY CMa. The current chemical models cannot reproduce the high abundances of PO and PN derived from observations of VY CMa. No observations have been reported of phosphorus in the circumstellar envelopes of oxygen-rich asymptotic giant branch stars. Aims. We aim to set observational constraints on the phosphorous chemistry in the circumstellar envelopes of oxygen-rich asymptotic giant branch stars, by focussing on the Mira-type variable star IK Tau. Methods. Using the IRAM 30 m telescope and the Submillimeter Array, we observed four rotational transitions of PN (J = 2−1,3−2,6−5,7−6) and four of PO (J = 5/2−3/2,7/2−5/2,13/2−11/2,15/2−13/2). The IRAM 30 m observations were dedicated line observations, while the Submillimeter Array data come from an unbiased spectral survey in the frequency range 279−355 GHz. Results. We present the first detections of PN and PO in an oxygen-rich asymptotic giant branch star and estimate abundances X(PN/H2) ≈ 3 × 10-7 and X(PO/H2) in the range 0.5−6.0 × 10-7. This is several orders of magnitude higher than what is found for the carbon-rich asymptotic giant branch star IRC +10 216. The diameter (≲0.′′7) of the PN and PO emission distributions measured in the interferometric data corresponds to a maximum radial extent of about 40 stellar radii. The abundances and the spatial occurrence of the molecules are in very good agreement with the results reported for VY CMa. We did not detect PS or PH3 in the survey. Conclusions. We suggest that PN and PO are the main carriers of phosphorus in the gas phase, with abundances possibly up to several 10-7. The current chemical models cannot account for this, underlining the strong need for updated chemical models that include phosphorous compounds.
Context.
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, which are believed to be formed primarily via an 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.
Aims.
We aim to improve our understanding of the wind morphology of eccentric asymptotic giant branch (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.
Methods.
Using the smoothed particle hydrodynamics (SPH) code
PHANTOM
, we generated 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.
Results.
Compared to the zero-eccentricity situation, we find that for low eccentricities, for the case of a high wind velocity, and hence limited interaction between the wind and the companion, the standard two-edged spiral structure that dominates the shape of the wind in the orbital plane is only minimally affected. When the wind speed is lower, strong compression of the wind material by the companion occurs, causing a high-pressure region around the companion which shapes the wind into an irregular spiral. In extreme cases, with low wind velocity and high eccentricity, these instabilities grow to such proportion that they cause high-speed ejections of matter along the orbital plane, shaping the wind into a highly irregular morphology. In more eccentric orbits, the amplitude of the phase-dependent wind-companion interaction increases significantly, introducing additional complexities that make the outbursts even more energetic, leading in some cases to high-speed polar flows of matter. Further, the orbital motion of the stars tends to flatten the global density distribution of the models with no instabilities. We distinguish global flattening from an equatorial density enhancement, the latter being formed by a strong gravitational interaction of the companion with the wind particles. We classify the resulting morphologies according to these new definitions, and find that (i) all low-velocity models have an equatorial density enhancement and (ii), in general, the flattening increases for decreasing wind velocity, until the low wind velocity results in high-energy outflows that clear away the flattening.
Conclusions.
We conclude that for models with a high wind velocity, the short interaction with the companion results in a regular spiral morphology, which 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.
Context
. Recent observations suggest the presence of clouds in exoplanet atmospheres, but they have also shown that certain chemical species in the upper atmosphere might not be in chemical ...equilibrium. Present and future interpretation of data from, for example, CHEOPS, JWST, PLATO, and Ariel require a combined understanding of the gas-phase and the cloud chemistry.
Aims
. The goal of this work is to calculate the two main cloud formation processes, nucleation, and bulk growth consistently from a non-equilibrium gas phase. The aim is also to explore the interaction between a kinetic gas-phase and cloud microphysics.
Methods
. The cloud formation is modelled using the moment method and kinetic nucleation, which are coupled to a gas-phase kinetic rate network. Specifically, the formation of cloud condensation nuclei is derived from cluster rates that include the thermochemical data of (TiO
2
)
N
from
N
= 1 to 15. The surface growth of nine bulk Al, Fe, Mg, O, Si, S, and Ti binding materials considers the respective gas-phase species through condensation and surface reactions as derived from kinetic disequilibrium. The effect of the completeness of rate networks and the time evolution of the cloud particle formation is studied for an example exoplanet, HD 209458 b.
Results
. A consistent, fully time-dependent cloud formation model in chemical disequilibrium with respect to nucleation, bulk growth, and the gas-phase is presented and first test cases are studied. This model shows that cloud formation in exoplanet atmospheres is a fast process. This confirms previous findings that the formation of cloud particles is a local process. Tests on selected locations within the atmosphere of the gas-giant HD 209458 b show that the cloud particle number density and volume reach constant values within 1 s. The complex kinetic polymer nucleation of TiO
2
confirms results from classical nucleation models. The surface reactions of SiOs and SiO
2
s can create a catalytic cycle that dissociates H
2
to 2 H, resulting in a reduction of the CH
4
number densities.
ABSTRACT
Asymptotic giant branch (AGB) stars shed a significant amount of their mass in the form of a stellar wind, creating a vast circumstellar envelope (CSE). Owing to the ideal combination of ...relatively high densities and cool temperatures, CSEs serve as rich astrochemical laboratories. While the chemical structure of AGB outflows has been modelled and analysed in detail for specific physical setups, there is a lack of understanding regarding the impact of changes in the physical environment on chemical abundances. A systematic sensitivity study is necessary to comprehend the nuances in the physical parameter space, given the complexity of the chemistry. This is crucial for estimating uncertainties associated with simulations and observations. In this work, we present the first sensitivity study of the impact of varying outflow densities and temperature profiles on the chemistry. With the use of a chemical kinetics model, we report on the uncertainty in abundances, given a specific uncertainty on the physical parameters. In addition, we analyse the molecular envelope extent of parent species and compare our findings to observational studies. Mapping the impact of differences in physical parameters throughout the CSE on the chemistry is a strong aid to observational studies.