Issue |
A&A
Volume 521, October 2010
|
|
---|---|---|
Article Number | L55 | |
Number of page(s) | 4 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201015252 | |
Published online | 14 October 2010 |
LETTER TO THE EDITOR
Signature of wide-spread clumping in B supergiant winds
R. K. Prinja1 - D. L. Massa2
1 - Department of Physics & Astronomy, UCL, Gower Street, London
WC1E 6BT, UK
2 - Space Telescope Science Institute, 3700 San Marino Drive, Baltimore, MD 21218, USA
Received 22 June 2010 / Accepted 15 July 2010
Abstract
Aims. We seek to establish additional observational
signatures of the effects of clumping in OB star winds. The action
of clumping on strategic wind-formed spectral lines is tested to steer
the development of models for clumped winds and thus improve the
reliability of mass-loss determinations for massive stars.
Methods. The Si IV
1400
resonance line doublets of B0 to B5 supergiants are analysed
using empirical line-synthesis models. The focus is on decoding
information on wind clumping from measurements of ratios of the radial
optical depths (
)
of the red and blue components of the Si IV doublet.
We exploit in particular the fact that the two doublet components are
decoupled and formed independently for targets with relatively low wind
terminal velocities.
Results. Line-synthesis analyses reveal that the mean ratio of
of the blue to red Si IV components are rarely close to the canonical value of
2 (expected from atomic constants), and spread instead over a range of values between
1
and 2. These results are interpreted in terms of a photosphere
that is partially obscured by optically thick structures in the
outflowing gas.
Conclusions. The spectroscopic signatures established in this
study demonstrate the wide-spread existence of wind clumping in
B supergiants. The additional information in unsaturated doublet
profiles provides a means to quantify the porosity of the winds.
Key words: stars: early-type - stars: mass-loss - ultraviolet: stars
1 Introduction
There are current serious uncertainties in the mass-loss and thus
energy feedback processes of massive OB stars. Over the past
decade evidence has accumulated from a broad collection of
results to strongly challenge the current model of hot star mass-loss
via stellar winds. Of particular impact on mass-loss
determinations is the presence of (small-scale) clumping
and (large-scale) structure, resulting in porous winds. The
characteristics of the X-ray emission (e.g. Cohen 2008), excess flux at IR and mm wavelengths (Runacres & Blomme 1996), discordance between FUV (sensitive to density, )
and H
(sensitive to
)
measurements of mass-loss rates (e.g. Massa et al. 2003; Bouret et al. 2005; Fullerton et al. 2006; Puls et al. 2006)
collectively suggest the winds are highly clumped, with a consequent
reduction in actual mass-loss rates. Given that some of these studies
question mass-loss rates of massive stars at the order-of-magnitude
level, there is the highest urgency to robustly investigate this
discordance since it has potentially far reaching consequences for the
evolution and fate of massive stars (which is largely dependent on
mass-loss) and galactic chemical evolution (where mass-loss drives
feedback into the interstellar medium).
Furthermore, recent theoretical predictions of the effects of structure in porous winds (Oskinova et al. 2007; Sundqvist et al. 2010), ionization by shock produced XUV line radiation (Waldron & Cassinelli 2010) and X-rays (Krticka et al. 2009) have amplified the need to verify via additional methods the compelling FUV studies cited above, which were based on the P V
1118, 1128 doublet.
2 The Si IV
1400 resonance line doublet ratios
In this letter we report on additional, untapped information on the
nature of wind clumping that can be extracted from an analysis of the
doublet ratios of Si IV
1393.76,
1402.77 in B supergiants. The AGN community have previously
recognised the potential of ratios of doublet components as a signature
of clumping (e.g. Ganguly et al. 1999). The reason is that for a source obscured by a uniform distribution
of gas, the ratio of optical depths determined from the observed
absorption lines would equal the ratio of the oscillator strengths of
the two doublet components. However, in the case of a source that
is partially obscured by totally opaque clumps, the ratios of
observed optical depth will be 1, since they will depend only on
the fraction
of the source which is uncovered.
We exploit this diagnostic here (in an extension to Massa et al. 2008): in a homogeneous massive star wind, the ratio of the optical depths of the widely spaced Si IV
1400 doublet would be the usual ratio of the oscillator strengths of the two components (i.e.
2.01).
However, in an OB star wind composed of very optically thick
clumps embedded in a transparent medium, the radial optical depth
of an absorption line depends only on the geometric covering factor of
the clumps. Therefore in this case reducing the optical depths
of the clumps by a factor of 2 (as implied in the Si IV doublets)
makes minimal difference in the absorption. Thus in this clumped case
the ratio of doublet optical depths measured from the absorption lines
would lie between
1 and 2.
Table 1: Si IV SEI model profile fitting parameters for the core sample.
3 The B supergiant sample
The dataset of Si IV resonance line doublets
analysed in this study is derived from the B0 to
B5 supergiants discussed by Prinja et al. (2005), using high resolution International Ultraviolet Explorer (IUE) satellite spectra. In particular our core sample of targets are B supergiants that (i) show well-developed but unsaturated Si IV absorption troughs; and (ii) have wind terminal velocities that are less than 0.5 of the Si IV doublet separation; i.e. stars with
970 km s-1.
This is an important selection criterion since it means that the two
components of the doublet can be treated as radiatively decoupled and
can be analysed separately i.e. each component as a singlet. The
core sample of 25 stars fitting our primary criteria above are
listed in Table 1.
To extract physical information from the line profiles we employed the
``Sobolev with Exact Integration'' (SEI) code of Lamers et al. (1987), but with the modifications described by Massa et al. (2003), which permit the radial optical depth (
)
to be treated as 21 discrete, independent and variable bins
in velocity space. The bins are then adjusted to match the observed
profile by a non-linear least squares procedure. The often complex
photospheric spectrum in the region of Si IV
1400 was included in the SEI fitting by using the IUE spectra of dwarf B stars listed in Prinja et al. (2002). (The B dwarf spectra
do not exhibit Si IV wind lines and were rotationally broadened to match the target B supergiants.)
![]() |
Figure 1:
Examples of SEI line synthesis fits to Si IV |
Open with DEXTER |
By model fitting the blue and red Si IV components as singlets, we are then able to extract two independent sets of
,
which can be inspected for signatures of wind clumping
(as justified in Sect. 2). The SEI model fitting
parameters, including
,
velocity law index (
)
and turbulent velocity parameter (
)
are also listed in Table 1, for the 25 B supergiants for which the Si IV doublet components can be considered as decoupled.
4 The evidence for clumping
Three representative examples of the Si IV blue and red components independently matched by SEI models are shown in Fig. 1,
covering the spectral type range B0.7 to B4. The
corresponding radial optical depths are also plotted for each case. The
ratio of the blue-to-red optical depths, averaged between 0.2 to
0.9 ,
is shown for the full sample of 25 stars in Fig. 2 (filled circles; listed in Table 1). For reliability, we only consider optical depths in the interval 0.3
5.0, to avoid weaker and saturated lines.
We assign an error of
5% in the
values
adopted in each velocity bin. It is clear that in almost all cases
the best fits for the decoupled doublet components yield optical depth
ratios less than 2.01 (the expected value from atomic
constants). The majority of
ratios lie between
1.0 and 1.5, and the overall mean is
1.46 (s.d.
0.31). As argued in Sect. 2, blue-to-red Si IV doublet
ratios significantly lower than
2 imply that the B supergiant winds are porous and contain very optically thick clumps. Note in Fig. 2
that the ratios do not depend on the effective temperature of the
relatively wide B0 to B5 spectral type range, and we
therefore rule out any systematic
effects that would cause all of the derived ratios to differ in the
same manner. Our results clearly indicate that the weaker (red)
component of the Si IV doublet is stronger than expected from smooth wind models.
![]() |
Figure 2:
Filled circles - the mean optical depth ratio between 0.2 and 0.9 |
Open with DEXTER |
An alternative measure of the departure from the smooth wind predictions is listed in Table 1, where the product of mass-loss rate Si3+ ion fraction (
qi) is provided for each doublet component, as determined from the SEI fits (see formulation in Massa et al. 2003). (The values in Table 1 are integrated over 0.2
0.9.) In almost all cases the ratios of the blue to red
measurements lie between 0.5 and 1.0. These results differ from the optical depth ratios shown in Fig. 2 because
(where
contains stellar and atomic data which are the same for both
components). This relation has two consequences. First, the ratios
of the
scale as
times
the ratio of the mean optical depths, which accounts for difference
in range. Second, the determinations of the mean values in
Table 1 are weighted differently, which accounts for the fact that mean
ratio entry is not exactly 2 times
.
4.1 Additional data sample
The evidence for widespread wind clumping in B supergiants presented here (Fig. 2) is based on the robust case of well separated, decoupled Si IV components, i.e. for stars that have
970 km s-1. There are of course many additional B supergiants with available UV data that have
in excess of
1000 km s-1. As the terminal velocity in these cases exceeds more than half the Si IV doublet
separation, their doublet components cannot be matched independently as
``singlets'' (as above). We can nevertheless extract some
complementary information from these higher
stars.
From the original B supergiant IUE sample of Prinja et al. (2005), we selected an additional 41 B0 to B2 stars and fitted the (now coupled) Si IV doublets with the SEI models (see Massa et al. 2008).
However, in these cases the ratio of the doublet oscillator
strengths was allowed to be a free parameter in the least squares fits.
The best fits to the sample of overlapping Si IV doublets yielded f-value ratios that also lie between 1 and 2 in almost all cases. The results are also shown in Fig. 2 (open circles) as a function of
.
For a fixed mass-loss rate, velocity law and ion fraction, the
ratio of blue-to-red oscillator strengths is equivalent to the ratio of
optical depths. Therefore the results shown in Fig. 2 are consistent with the original sample (Table 1 and Fig. 2 filled circles) of low
,
decoupled doublets. Taken together, we conclude that the two datasets
make a compelling case for wide-spread clumping and porosity in
B supergiant winds.
5 Discussion
We have demonstrated that the Si IV
1400 doublets
of B supergiants contain additional signatures of clumping in the
stellar winds. Line synthesis models of 25 B0 to B5 stars for
which the Si IV blue and red components can be treated as decoupled (``singlets''), yield
blue-to-red optical depth ratios between 1 and 2, i.e. mostly well below the value
2
expected from atomic constants and a smooth wind. These results may be
interpreted as evidence for optically thick clumps in the wind,
covering only a fraction of the stellar source. Additional results
derived from B stars with overlapping Si IV doublets are consistent with this conclusion.
The optical depth ratios (Fig. 2)
are not dependent on the effective temperature. The large scatter
observed in the ratio for a given spectral type bin is also not
correlated with the projected stellar rotation velocity or mass-loss
rate. Aside from inclination effects connected to large-scale
structures such as co-rotating interaction regions, the scatter in
Fig. 2 may
reflect transient and episodic events in the wind. It may be
expected that for highly porous and
extremely clumped cases, the temporal effect due to the
time-dependent arrangement of clumps would be more pronounced.
As an illustration, we examined the IUE time-series dataset of HD 47240 (B1 Ib; Table 1) described by Prinja et al. (2002), to test for temporal variations in the mean blue-to-red optical depth ratios. Thirteen Si IV spectra (SWP48835 to 48953), spanning 16 days,
were matched by SEI models according to the methods for decoupled
doublets outlined in Sect. 3. The mean blue-to-red
ratio measured over 0.2 to 0.9
varies considerably during the 16 day monitoring of HD 47240, with values between
1.2 to 2.0, i.e. essentially the full range of scatter seen in Fig. 2 at a given spectral type. Consequently, we suspect that much of the scatter seen in Fig. 2
results from time-variable covering factors. The spectroscopic changes
modelled here in HD 47240 are empirically due to time-variable
velocity dependent structure in the UV absorption line profiles.
Several previous time-series studies of variations in the
UV resonance lines of OB stars have interpreted this temporal
behaviour as due to the evolution and velocity migration of large-scale
structure in the winds, perhaps due to density inhomogeneities or
velocity plateaus. The substantial UV line profile changes evident
in OB stars are not diagnostics of the effects of small-scale mirco-clumping. The inference here then may be that the variable blue-to-red
ratios
derived from the UV resonance lines of an individual OB star
reflect the porosity of a macro-clumped wind. We now plan additional
studies to derive information about the spatial extent and temporal
evolution of the obscuring structure, which may be decoded from the
suspected wavelength dependence of observed changes in the mean
ratios.
![]() |
Figure 3:
Plot of the observed P V
|
Open with DEXTER |
To a large degree the current debate concerning clumping in porous
winds, and the consequent over-estimation of clumping independent
mass-loss rates, is founded on observational and theoretical
studies of O-type stars (see e.g. references
in Sect. 1). The line doublets diagnostics that we have
demonstrated here for B supergiants can to some extent also be
tested on O-type stars. Since the terminal velocities of O-type stars
are mostly higher than for B supergiants, only the method
for overlapping doublets employed in Sect. 4.1 is feasible for
O stars, i.e. allowing the f-ratio
of the doublet to vary from the value expected from atomic constants.
A demonstration of this application is shown in Fig. 3 for the O4 I star, Pup from Massa et al. (2008).
We compare here matches to the FUV P V
1117, 1128 doublet P Cygni profile for (i) the f-ratio fixed to the value for the atomic constants (=2.02); and (ii) allowing the f-ratio
to be a free parameter; a value of 1.84 is adopted for the
final least squares fit. This 10% change in the f-ratio
undoubtedly improves the match, and is consistent with the evidence for
clumping in B supergiants reported in this letter. Finally, we
anticipate that the clumping signature contained in the optical depth
ratios of doublet lines may also be exploited
in other low
settings, such as the fast winds of PN central stars and OB stars in the Magellanic Clouds.
D.M. acknowledges support from NASA ADP contract NNH08CD06C to SGT, Inc.
References
- Bouret, J.-C., Lanz, T., & Hillier, D. J. 2005, A&A, 438, 301 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cohen, D. J. 2008, in Massive Stars as Cosmic Engines, ed. F. Bresolin, P. A. Crowther, & J. Puls (Cambridge University Press), Proc. IAU Symp., 250, 17 [Google Scholar]
- Fullerton, A. W., Massa, D. L., & Prinja, R. K. 2006, ApJ, 637, 1025 [NASA ADS] [CrossRef] [Google Scholar]
- Ganguly, R., Eracleous, M., Charlton, J. C., & Chruchill, C. W. 1999, AJ, 117, 2594 [NASA ADS] [CrossRef] [Google Scholar]
- Krticka, J., Feldmeier, A., Oskinova, L. M., Kubat, J., & Hamann, W.-R. 2009, A&A, 508, 841 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lamers, H. J. G. L. M., Cerruti-Sola, M., & Perinotto, M. 1987, ApJ, 314, 726 [NASA ADS] [CrossRef] [Google Scholar]
- Massa, D., Fullerton, A. W., Sonneborn, G., & Hutchings, J. B. 2003, ApJ, 586, 996 [NASA ADS] [CrossRef] [Google Scholar]
- Massa, D. L., Prinja, R. K., & Fullerton, A. W. 2008, in Clumping in Hot-Star Winds, ed. W.-R. Hamann, A. Feldmeier, & L. M. Oskinova, Univ.-Verl., 147 [Google Scholar]
- Oskinova, L. M., Hamann, W.-R., & Feldmeier, A. 2007, A&A, 476, 1331 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Prinja, R. K., Massa, D., & Fullerton, A. W. 2002, A&A, 388, 587 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Prinja, R. K., Massa, D., & Searle, S. C. 2005, A&A, 430, L41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Puls, J., Markova, N., Scuderi, S., et al. 2006, A&A, 454, 625 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Runacres, M. C., & Blomme, R. 1996, A&A, 309, 544 [NASA ADS] [Google Scholar]
- Searle, S. C., Prinja, R. K., Massa, D., & Ryans, R. 2008, A&A, 481, 777 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sundqvist, J., Puls, J., & Feldmeier, A. 2010, A&A, 510, 11 [Google Scholar]
- Waldron, W. L., & Cassinelli, J. P. 2010, ApJ, 711, L30 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Table 1: Si IV SEI model profile fitting parameters for the core sample.
All Figures
![]() |
Figure 1:
Examples of SEI line synthesis fits to Si IV |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Filled circles - the mean optical depth ratio between 0.2 and 0.9 |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Plot of the observed P V
|
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.