Context.
Red supergiants are observed to undergo vigorous mass loss. However, to date, no theoretical model has succeeded in explaining the origins of these objects’ winds. This strongly limits our ...understanding of red supergiant evolution and Type II-P and II-L supernova progenitor properties.
Aims.
We examine the role that vigorous atmospheric turbulence may play in initiating and determining the mass-loss rates of red supergiant stars.
Methods.
We analytically and numerically solve the equations of conservation of mass and momentum, which we later couple to an atmospheric temperature structure, to obtain theoretically motivated mass-loss rates. We then compare these to state-of-the-art empirical mass-loss rate scaling formulae as well as observationally inferred mass-loss rates of red supergiants.
Results.
We find that the pressure due to the characteristic turbulent velocities inferred for red supergiants is sufficient to explain the mass-loss rates of these objects in the absence of the normally employed opacity from circumstellar dust. Motivated by this initial success, we provide a first theoretical and fully analytic mass-loss rate prescription for red supergiants. We conclude by highlighting some intriguing possible implications of these rates for future studies of stellar evolution, especially in light of the lack of a direct dependence on metallicity.
Context. The behaviour of mass loss across the so-called bi-stability jump, where iron recombines from Fe IV to Fe III, is a key uncertainty in models of massive stars. Specifically, while an ...increase in mass loss is theoretically predicted, this has not yet been observationally confirmed. However, radiation-driven winds of hot massive stars are known to exhibit clumpy structures triggered by the line-deshadowing instability (LDI). This wind clumping severely affects empirical mass-loss rates inferred from ρ2-dependent spectral diagnostics. Thus, if clumping properties differ significantly for O and B supergiants across the bi-stability jump, this may help alleviate current discrepancies between theory and observations. Aims. We investigated with analyt ical and numerical tools how the onset of clumpy structures behave in the winds of O supergiants (OSG) and B supergiants (BSG) across the bi-stability jump. Methods. We derived a scaling relation for the linear growth rate of the LDI for a single optically thick line and applied it in the OSG and BSG regime. We ran 1D time-dependent line-driven instability simulations to study the non-linear evolution of the LDI in clumpy OSG and BSG winds. Results. Linear perturbation analysis for a single line shows that the LDI linear growth rate Ω scales strongly with stellar effective temperature and terminal wind speed: Ω∝v∞2Teff4 $\Omega \propto \varv_{\infty}^2 T_{\mathrm{eff}}^4$ Ω∝v∞2Teff4 . This implies significantly lower growth rates for (the cooler and slower) BSG winds than for OSG winds. This is confirmed by the non-linear simulations, which show significant differences in OSG and BSG wind structure formation, with the latter characterized by significantly weaker clumping factors and lower velocity dispersions. This suggests that lower correction factors due to clumping should be employed when deriving empirical mass-loss rates for BSGs on the cool side of the bi-stability jump. Moreover, the non-linear simulations provide a theoretical background towards explaining the general lack of observed intrinsic X-ray emission in single B-star winds.
Context. State of the art quantitative spectroscopy utilizes synthetic spectra to extract information from observations. For hot, massive stars, these synthetic spectra are calculated by means of 1D, ...spherically symmetric, NLTE atmosphere and spectrum-synthesis codes. Certain stellar atmospheres, however, show strong deviations from spherical symmetry, and need to be treated in 3D. Aims. We present and test a newly developed 3D radiative transfer code, tailored to the solution of the radiation field in rapidly expanding stellar atmospheres. We apply our code to the continuum transfer in wind-ablation models, and to the UV resonance line formation in magnetic winds. Methods. We have used a 3D finite-volume method for the solution of the time-independent equation of radiative transfer, to study continuum- and line-scattering problems, currently approximated by a two-level-atom. Convergence has been accelerated by coupling the formal solver to a non-local approximate Λ-iteration scheme. Particular emphasis has been put on careful tests, by comparing with alternative solutions for 1D, spherically symmetric model atmospheres. These tests allowed us to understand certain shortcomings of the methods, and to estimate limiting cases that can actually be calculated. Results. The typical errors of the converged source functions, when compared to 1D solutions, are of the order of 10–20%, and rapidly increase for optically thick (τ ≳ 10) continua, mostly due to the order of accuracy of our solution scheme. In circumstellar discs, the radiation temperatures in the (optically thin) transition region from wind to disc are quite similar to corresponding values in the wind. For MHD simulations of dynamical magnetospheres, the line profiles, calculated with our new 3D code, agree well with previous solutions using a 3D-SEI method. When compared with profiles resulting from the so-called analytic dynamical magnetosphere (ADM) model, however, significant differences become apparent. Conclusions. Due to similar radiation temperatures in the wind and the transition region to the disc, the same line-strength distribution can be applied within radiation hydrodynamic calculations for optically thick circumstellar discs in “accreting high-mass stars”. To properly describe the UV line formation in dynamical magnetospheres, the ADM model needs to be further developed, at least in a large part of the outer wind.
Context.
Knowledge about hot, massive stars is usually inferred from quantitative spectroscopy. To analyse non-spherical phenomena, the existing 1D codes must be extended to higher dimensions, and ...corresponding tools need to be developed.
Aims.
We present a 3D radiative transfer code that is capable of calculating continuum and line scattering problems in the winds of hot stars. By considering spherically symmetric test models, we discuss potential error sources, and indicate advantages and disadvantages by comparing different solution methods. Further, we analyse the ultra-violet (UV) resonance line formation in the winds of rapidly rotating O stars.
Methods.
We consider both a (simplified) continuum model including scattering and thermal sources, and a UV resonance line transition approximated by a two-level-atom. We applied the short-characteristics (SC) method, using linear or monotonic Bézier interpolations, for which monotonicity is of prime importance, to solve the equation of radiative transfer on a non-uniform Cartesian grid. To calculate scattering dominated problems, our solution method is supplemented by an accelerated Λ-iteration scheme using newly developed non-local operators.
Results.
For the spherical test models, the mean relative error of the source function is on the 5 − 20% level, depending on the applied interpolation technique and the complexity of the considered model. All calculated line profiles are in excellent agreement with corresponding 1D solutions. Close to the stellar surface, the SC methods generally perform better than a 3D finite-volume-method; however, they display specific problems in searchlight-beam tests at larger distances from the star. The predicted line profiles from fast rotating stars show a distinct behaviour as a function of rotational speed and inclination. This behaviour is tightly coupled to the wind structure and the description of gravity darkening and stellar surface distortion.
Conclusions.
Our SC methods are ready to be used for quantitative analyses of UV resonance line profiles. When calculating optically thick continua, both SC methods give reliable results, in contrast to the alternative finite-volume method.
Context.
Line-driven winds of hot, luminous stars are intrinsically unstable due to the line-deshadowing instability (LDI). In non-magnetic hot stars, the LDI leads to the formation of an ...inhomogeneous wind consisting of small-scale, spatially separated clumps that can have great effects on observational diagnostics. However, for magnetic hot stars the LDI generated structures, wind dynamics, and effects on observational diagnostics have not been directly investigated so far.
Aims.
We investigated the non-linear development of LDI generated structures and dynamics in a magnetic line-driven wind of a typical O-supergiant.
Methods.
We employed two-dimensional axisymmetric magnetohydrodynamic simulations of the LDI using the Smooth Source Function approximation for evaluating the assumed one-dimensional line force. To facilitate the interpretation of these magnetic models, they were compared with a corresponding non-magnetic LDI simulation as well as a magnetic simulation neglecting the LDI.
Results.
A central result obtained is that the wind morphology and wind clumping properties change strongly with increasing wind-magnetic confinement. Most notably, in magnetically confined flows, the LDI leads to large-scale, shellular sheets (‘pancakes’) that are quite distinct from the spatially separate, small-scale clumps in non-magnetic line-driven winds. We discuss the impact of these findings for observational diagnostic studies and stellar evolution models of magnetic hot stars.
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
Line-driven stellar winds from massive (OB) stars are subject to a strong line-deshadowing instability. Recently, spectropolarimetric surveys have collected ample evidence that a subset of ...Galactic massive stars hosts strong surface magnetic fields. We investigate here the propagation and stability of magnetoradiative waves in such a magnetized, line-driven wind. Our analytic, linear stability analysis includes line-scattering from the stellar radiation, and accounts for both radial and non-radial perturbations. We establish a bridging law for arbitrary perturbation wavelength after which we analyse separately the long- and short-wavelength limits. While long-wavelength radiative and magnetic waves are found to be completely decoupled, a key result is that short-wavelength, radially propagating Alfvén waves couple to the scattered radiation field and are strongly damped due to the line-drag effect. This damping of magnetic waves in a scattering-line-driven flow could have important effects on regulating the non-linear wind dynamics, and so might also have strong influence on observational diagnostics of the wind structure and clumping of magnetic line-driven winds.
The winds of hot, massive stars are variable from processes happening on both large and small spatial scales. A particular case of such wind variability is ‘discrete-absorption components’ (DACs) ...that manifest themselves as outward moving density features in UV resonance line spectra. Such DACs are believed to be caused by large-scale spiral-shaped density structures in the stellar wind. We consider novel 3-D radiation-hydrodynamic models of rotating hot star winds and study the emergence of co-rotating spiral structures due to a local (pseudo-)magnetic spot on the stellar surface. Subsequently, the hydrodynamic models are used to retrieve DAC spectral signatures in synthetic UV spectra created from a 3-D short-characteristics radiative transfer code.
Smoking has been associated with colorectal cancer (CRC) incidence and mortality in previous studies and might also be associated with prognosis after CRC diagnosis. However, current evidence on ...smoking in association with CRC prognosis is limited.
For this individual patient data meta-analysis, sociodemographic and smoking behavior information of 12 414 incident CRC patients (median age at diagnosis: 64.3 years), recruited within 14 prospective cohort studies among previously cancer-free adults, was collected at baseline and harmonized across studies. Vital status and causes of death were collected for a mean follow-up time of 5.1 years following cancer diagnosis. Associations of smoking behavior with overall and CRC-specific survival were evaluated using Cox regression and standard meta-analysis methodology.
A total of 5229 participants died, 3194 from CRC. Cox regression revealed significant associations between former hazard ratio (HR) = 1.12; 95 % confidence interval (CI) = 1.04–1.20 and current smoking (HR = 1.29; 95% CI = 1.04–1.60) and poorer overall survival compared with never smoking. Compared with current smoking, smoking cessation was associated with improved overall (HR<10 years = 0.78; 95% CI = 0.69–0.88; HR≥10 years = 0.78; 95% CI = 0.63–0.97) and CRC-specific survival (HR≥10 years = 0.76; 95% CI = 0.67–0.85).
In this large meta-analysis including primary data of incident CRC patients from 14 prospective cohort studies on the association between smoking and CRC prognosis, former and current smoking were associated with poorer CRC prognosis compared with never smoking. Smoking cessation was associated with improved survival when compared with current smokers. Future studies should further quantify the benefits of nonsmoking, both for cancer prevention and for improving survival among CRC patients, in particular also in terms of treatment response.
Context.
Classical Wolf-Rayet (WR) stars are direct supernova progenitors undergoing vigorous mass loss. Understanding the dense and fast outflows of such WR stars is thus crucial for understanding ...advanced stages of stellar evolution and the dynamical feedback of massive stars on their environments, and for characterizing the distribution of black hole masses.
Aims.
In this paper, we develop the first time-dependent, multidimensional, radiation-hydrodynamical models of the extended optically thick atmospheres and wind outflows of hydrogen-free classical WR stars.
Methods.
A flux-limiting radiation hydrodynamics approach is used on a finite volume mesh to model WR outflows. The opacities are described using a combination of tabulated Rosseland mean opacities and the enhanced line opacities expected within a supersonic flow.
Results.
For high-luminosity models, a radiation-driven, dense, supersonic wind is launched from deep subsurface regions associated with peaks in the Rosseland mean opacity. For a model with lower luminosity, on the other hand, the Rosseland mean opacity is not sufficient to sustain a net-radial outflow in the subsurface regions. Instead, what develops in this case, is a "standard" line-driven wind launched from the optically thin regions above an extended, moderately inflated, and highly turbulent atmosphere. We thus find here a natural transition from optically thick outflows of classical WR stars to optically thin winds of hot, compact subdwarfs; in our simulations, this transition occurs approximately at a luminosity that is ~40% of the Eddington luminosity. Because of the changing character of the wind-launching mechanism, this transition is also accompanied by a large drop (on the low-luminosity end) in the average mass-loss rate. Since the subsurface opacity peaks are further associated with convective instabilities, the flows are highly structured and turbulent, consisting of coexisting regions of outflowing, stagnated, and even pockets of infalling gas. Typical velocity dispersions in our 3D models are high, 100–300 km s
−1
, but the clumping factors are rather modest,
f
c1
≡ 〈
ρ
2
〉/〈
ρ
〉
2
~
2. We further find that, while the low-density gas in our simulations is strongly radiation-driven, the overdense structures are, after their initial launch, primarily advected outward by ram-pressure gradients. This inefficient radiative acceleration of dense "clumps" reflects the inverse dependence of line driving on mass density and leads to a general picture wherein high-density gas parcels move significantly slower than the mean and low-density wind material.