Context. The evolution of intermediate and low-mass stars on the asymptotic giant branch is dominated by their strong dust-driven winds. More massive stars evolve into red supergiants with a similar ...envelope structure and strong wind. These stellar winds are a prime source for the chemical enrichment of the interstellar medium. Aims. We aim to (1) set up simple and general analytical expressions to estimate mass-loss rates of evolved stars, and (2) from those calculate estimates for the mass-loss rates of the asymptotic giant branch, red supergiant, and yellow hypergiant stars in our galactic sample. Methods. The rotationally excited lines of carbon monoxide (CO) are a classic and very robust diagnostic in the study of circumstellar envelopes. When sampling different layers of the circumstellar envelope, observations of these molecular lines lead to detailed profiles of kinetic temperature, expansion velocity, and density. A state-of-the-art, nonlocal thermal equilibrium, and co-moving frame radiative transfer code that predicts CO line intensities in the circumstellar envelopes of late-type stars is used in deriving relations between stellar and molecular-line parameters, on the one hand, and mass-loss rate, on the other. These expressions are applied to our extensive CO data set to estimate the mass-loss rates of 47 sample stars. Results. We present analytical expressions for estimating the mass-loss rates of evolved stellar objects for 8 rotational transitions of the CO molecule and thencompare our results to those of previous studies. Our expressions account for line saturation and resolving of the envelope, thereby allowing accurate determination of very high mass-loss rates. We argue that, for estimates based on a single rotational line, the CO(2–1) transition provides the most reliable mass-loss rate. The mass-loss rates calculated for the asympotic giant branch stars range from 4 × 10-8 M⊙ yr-1 up to 8 × 10-5 M⊙ yr-1. For red supergiants they reach values between 2 × 10-7 M⊙ yr-1 and 3 × 10-4 M⊙ yr-1. The estimates for the set of CO transitions allow time variability to be identified in the mass-loss rate. Possible mass-loss-rate variability is traced for 7 of the sample stars. We find a clear relation between the pulsation periods of the asympotic giant branch stars and their derived mass-loss rates, with a levelling off at ~3 × 10-5 M⊙ yr-1 for periods exceeding 850 days. Conclusions.
CONFIRMATION OF CIRCUMSTELLAR PHOSPHINE AGUNDEZ, M; Cernicharo, J; Decin, L ...
Astrophysical journal. Letters,
08/2014, Letnik:
790, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Phosphine (PH sub(3)) was tentatively identified a few years ago in the carbon star envelopes IRC +10216 and CRL 2688 from observations of an emission line at 266.9 GHz attributable to the J = 1-0 ...rotational transition. We report the detection of the J = 2-1 rotational transition of PH sub(3) in IRC +10216 using the HIFI instrument on board Herschel, which definitively confirms the identification of PH sub(3). Radiative transfer calculations indicate that infrared pumping in excited vibrational states plays an important role in the excitation of PH sub(3) in the envelope of IRC +10216, and that the observed lines are consistent with phosphine being formed anywhere between the star and 100 Rlow * from the star, with an abundance of 10 super(-8) relative to H sub(2). The detection of PH sub(3) challenges chemical models, none of which offer a satisfactory formation scenario. Although PH sub(3) holds just 2% of the total available phosphorus in IRC +10216, it is, together with HCP, one of the major gas phase carriers of phosphorus in the inner circumstellar layers, suggesting that it could also be an important phosphorus species in other astronomical environments. This is the first unambiguous detection of PH sub(3) outside the solar system, and is a further step toward a better understanding of the chemistry of phosphorus in space.
Context.
Vigorous mass loss in the classical Wolf-Rayet (WR) phase is important for the late evolution and final fate of massive stars.
Aims.
We develop spherically symmetric time-dependent and ...steady-state hydrodynamical models of the radiation-driven wind outflows and associated mass loss from classical WR stars.
Methods.
The simulations are based on combining the opacities typically used in static stellar structure and evolution models with a simple parametrised form for the enhanced line opacity expected within a supersonic outflow.
Results.
Our simulations reveal high mass-loss rates initiated in deep and hot, optically thick layers around
T
≈ 200 kK. The resulting velocity structure is non-monotonic and can be separated into three phases: (i) an initial acceleration to supersonic speeds (caused by the static opacity), (ii) stagnation and even deceleration, and (iii) an outer region of rapid re-acceleration (by line opacity). The characteristic structures seen in converged steady-state simulations agree well with the outflow properties of our time-dependent models.
Conclusions.
By directly comparing our dynamic simulations to corresponding hydrostatic models, we explicitly demonstrate that the need to invoke extra energy transport in convectively inefficient regions of stellar structure and evolution models, in order to prevent drastic inflation of static WR envelopes, is merely an artefact of enforcing a hydrostatic outer boundary. Moreover, the dynamically inflated inner regions of our simulations provide a natural explanation for the often-found mismatch between predicted hydrostatic WR radii and those inferred from spectroscopy; by extrapolating a monotonic
β
-type velocity law from the observable supersonic regions to the invisible hydrostatic core, spectroscopic models likely overestimate the core radius by a factor of a few. Finally, we contrast our simulations with alternative recent WR wind models based on co-moving frame (CMF) radiative transfer to compute the radiation force. Since CMF transfer currently cannot handle non-monotonic velocity fields, the characteristic deceleration regions found here are avoided in such simulations by invoking an ad hoc very high degree of clumping.
ABSTRACT
Electromagnetic radiation plays a crucial role in various physical and chemical processes. Hence, almost all astrophysical simulations require some form of radiative transfer model. Despite ...many innovations in radiative transfer algorithms and their implementation, realistic radiative transfer models remain very computationally expensive, such that one often has to resort to approximate descriptions. The complexity of these models makes it difficult to assess the validity of any approximation and to quantify uncertainties on the model results. This impedes scientific rigour, in particular, when comparing models to observations, or when using their results as input for other models. We present a probabilistic numerical approach to address these issues by treating radiative transfer as a Bayesian linear regression problem. This allows us to model uncertainties on the input and output of the model with the variances of the associated probability distributions. Furthermore, this approach naturally allows us to create reduced-order radiative transfer models with a quantifiable accuracy. These are approximate solutions to exact radiative transfer models, in contrast to the exact solutions to approximate models that are often used. As a first demonstration, we derive a probabilistic version of the method of characteristics, a commonly-used technique to solve radiative transfer problems.
ABSTRACT
Chemical modelling of asymptotic giant branch (AGB) outflows is typically focused on either non-thermodynamic equilibrium chemistry in the inner region or photon-driven chemistry in the ...outer region. We include, for the first time, a comprehensive dust–gas chemistry in our AGB outflow chemical kinetics model, including both dust–gas interactions and grain-surface chemistry. The dust is assumed to have formed in the inner region, and follows an interstellar-like dust-size distribution. Using radiative transfer modelling, we obtain dust temperature profiles for different dust types in an O-rich and a C-rich outflow. We calculate a grid of models, sampling different outflow densities, drift velocities between the dust and gas, and dust types. Dust–gas chemistry can significantly affect the gas-phase composition, depleting parent and daughter species and increasing the abundance of certain daughter species via grain-surface formation followed by desorption/sputtering. Its influence depends on four factors: outflow density, dust temperature, initial composition, and drift velocity. The largest effects are for higher density outflows with cold dust and O-rich parent species, as these species generally have a larger binding energy. At drift velocities larger than ∼10 km s−1, ice mantles undergo sputtering; however, they are not fully destroyed. Models with dust–gas chemistry can better reproduce the observed depletion of species in O-rich outflows. When including colder dust in the C-rich outflows and adjusting the binding energy of CS, the depletion in C-rich outflows is also better reproduced. To best interpret high-resolution molecular line observations from AGB outflows, dust–gas interactions are needed in chemical kinetics models.
We report on the successful search forCO(2-1) and (3-2) emission associated with OH/IR stars in the Galactic bulge. We observed a sample of eight extremely red asymptotic giant branch stars with the ...Atacama Pathfinder EXperiment telescope and detected seven. The sources were selected at a sufficiently high galactic latitude to avoid interference by interstellar CO, which hampered previous studies of inner galaxy stars. To study the nature of our sample and the mass-loss, we constructed the spectral energy distribution (SEDs) from photometric data and Spitzer IRS spectroscopy. In a first step, we apply radiative transfer modelling to fit the SEDs and obtain luminosities and dust mass-loss rates (MLRs). Through dynamical modelling, we then retrieve the total MLR and the gas-to-dust ratios. We derived variability periods of our stars. The luminosities range between approximately 4000 and 5500 L⊙and periods are shorter than 700 d. The total MLR ranges between 10-5and 10-4M⊙yr-1. Comparison with evolutionary models shows that the progenitor mass ≈ 1.5 M⊙, similar to the bulge Miras of intermediate age (3 Gyr). The gas-to-dust ratios are between 100 and 400 and are similar to what is found for OH/IR stars in the galactic disc. One star, IRAS 17347-2319, has a very short period of approximately 300 d that may be decreasing further. It may belong to a class of Mira variables with a sudden change in period as observed in some Galactic objects. It would be the first example of an OH/IR star in this class and deserves further follow-up observations.
We provide an update for the astrometric solution for the Type II supernova progenitor Betelgeuse using the revised Hipparcos Intermediate Astrometric Data (HIAD) of van Leeuwen, combined with ...existing VLA and new e-MERLIN and ALMA positions. The 2007 Hipparcos refined abscissa measurements required the addition of so-called Cosmic Noise of 2.4 mas to find an acceptable 5-parameter stochastic solution. We find that a measure of radio Cosmic Noise should also be included for the radio positions because surface inhomogeneities exist at a level significant enough to introduce additional intensity centroid uncertainty. Combining the 2007 HIAD with the proper motions based solely on the radio positions leads to a parallax of mas ( pc), smaller than the Hipparcos 2007 value of 6.56 0.83 mas ( pc). Furthermore, combining the VLA and new e-MERLIN and ALMA radio positions with the 2007 HIAD, and including radio Cosmic Noise of 2.4 mas, leads to a nominal parallax solution of 4.51 0.80 mas ( pc), which, while only different from the 2008 solution of Harper et al., is 2.6 different from the solution of van Leeuwen. An accurate and precise parallax for Betelgeuse is always going to be difficult to obtain because it is small compared to the stellar angular diameter ( mas). We outline an observing strategy utilizing future mm and sub-mm high-spatial resolution interferometry that must be used if substantial improvements in the precision and accuracy of the parallax and distance are to be achieved.
Context. On the asymptotic giant branch (AGB) low- and intermediate-mass stars eject a large fraction of their envelope, but the mechanism driving these outflows is still poorly understood. For ...oxygen-rich AGB stars, the wind is thought to be driven by radiation pressure caused by scattering of radiation off dust grains. Aims. We study the photosphere, the warm molecular layer, and the inner wind of the close-by oxygen-rich AGB star R Doradus. We focus on investigating the spatial distribution of the dust grains that scatter light and whether these grains can be responsible for driving the outflow of this star. Methods. We use high-angular-resolution images obtained with SPHERE/ZIMPOL to study R Dor and its inner envelope in a novel way. We present observations in filters V, cntHα, and cnt820 and investigate the surface brightness distribution of the star and of the polarised light produced in the inner envelope. Thanks to second-epoch observations in cntHα, we are able to see variability on the stellar photosphere. We study the polarised-light data using a continuum-radiative-transfer code that accounts for direction-dependent scattering of photons off dust grains. Results. We find that in the first epoch the surface brightness of R Dor is asymmetric in V and cntHα, the filters where molecular opacity is stronger, while in cnt820 the surface brightness is closer to being axisymmetric. The second-epoch observations in cntHα show that the morphology of R Dor has changed completely in a timespan of 48 days to a more axisymmetric and compact configuration. This variable morphology is probably linked to changes in the opacity provided by TiO molecules in the extended atmosphere. The observations show polarised light coming from a region around the central star. The inner radius of the region from where polarised light is seen varies only by a small amount with azimuth. The value of the polarised intensity, however, varies by between a factor of 2.3 and 3.7 with azimuth for the different images. We fit the radial profile of the polarised intensity using a spherically symmetric model and a parametric description of the dust density profile, ρ(r) = ρ°r− n. On average, we find exponents of − 4.5 ± 0.5 that correspond to a much steeper density profile than that of a wind expanding at constant velocity. The dust densities we derive imply an upper limit for the dust-to-gas ratio of ~2 × 10-4 at 5.0 R⋆. Considering all the uncertainties in observations and models, this value is consistent with the minimum values required by wind-driving models for the onset of a wind, of ~3.3 × 10-4. However, if the steep density profile we find extends to larger distances from the star, the dust-to-gas ratio will quickly become too small for the wind of R Dor to be driven by the grains that produce the scattered light.
Context. In the inner regions of asymptotic giant branch (AGB) outflows, several molecules have been detected with abundances much higher than those predicted from thermodynamic equilibrium chemical ...models. The presence of the majority of these species can be explained by shock-induced non-equilibrium chemical models, where shocks caused by the pulsating star take the chemistry out of equilibrium in the inner region. Moreover, a non-uniform density structure has been detected in several AGB outflows. Both large-scale structures, such as spirals and disks, and small-scale density inhomogeneities or clumps have been observed. These structures may also have a considerable impact on the circumstellar chemistry. A detailed parameter study on the quantitative effects of a non-homogeneous outflow has so far not been performed. Aims. We examine the effects of a non-uniform density distribution within an AGB outflow on its chemistry by considering a stochastic, clumpy density structure. Methods. We implement a porosity formalism for treating the increased leakage of light associated with radiation transport through a clumpy, porous medium. We then use this method to examine the effects from the altered UV radiation field penetration on the chemistry, accounting also for the increased reaction rates of two-body processes in the overdense clumps. The specific clumpiness is determined by three parameters: the characteristic length scale of the clumps at the stellar surface, the clump volume filling factor, and the inter-clump density contrast. In this paper, the clumps are assumed to have a spatially constant volume filling factor, which implies that they expand as they move outward in the wind. Results. We present a parameter study of the effect of clumping and porosity on the chemistry throughout the outflow. Both the higher density within the clumps and the increased UV radiation field penetration have an important impact on the chemistry, as they both alter the chemical pathways throughout the outflow. The increased amount of UV radiation in the inner region leads to photodissociation of parent species, releasing the otherwise deficient elements. We find an increased abundance in the inner region of all species not expected to be present assuming thermodynamic equilibrium chemistry, such as HCN in O-rich outflows, H2O in C-rich outflows, and NH3 in both. Conclusions. A non-uniform density distribution directly influences the chemistry throughout the AGB outflow, both through the density structure itself and through its effect on the UV radiation field. Species not expected to be present in the inner region of the outflow assuming thermodynamic equilibrium chemistry are now formed in this region, including species that are not formed in greater abundance by shock-induced non-equilibrium chemistry models. Outflows whose clumps have a large overdensity and that are very porous to the interstellar UV radiation field yield abundances comparable to those observed in O-rich and C-rich outflows for most of the unexpected species investigated. The inner wind abundances of H2O in C-rich outflows and of NH3 in O-rich and C-rich outflows are however underpredicted.