Free Access
Issue
A&A
Volume 550, February 2013
Article Number A55
Number of page(s) 9
Section Stellar atmospheres
DOI https://doi.org/10.1051/0004-6361/201118544
Published online 24 January 2013

© ESO, 2013

1. Introduction

Since the launch of the first X-ray observatories, it has been known that O-type stars and early B-type stars are X-ray sources. Harnden et al. (1979), Seward et al. (1979), and Cassinelli et al. (1981) described the first observations with Einstein. The investigations were continued with ASCA (Corcoran et al. 1994) and ROSAT (Haberl & White 1993; Berghöfer et al. 1996). Although the origin of the emission was initially a matter of debate, it is now generally believed that instabilities and clumping in the strong line-driven stellar wind is the mechanism that heats small parts of the wind (e.g., Lucy & Solomon 1970; Lucy & White 1980; Lucy 1982; Owocki et al. 1988; Dessart & Owocki 2005). The launches of Chandra and XMM-Newton in 1999 with high-resolution instruments aboard have offered the possibility to study the X-ray spectra in more detail, resulting in, e.g., estimates of line ratios, line profiles, temperatures, emission measures, and abundances (Waldron & Cassinelli 2001, 2007; Cohen et al. 2006; Kahn et al. 2001; Raassen et al. 2008). Generally, the overall spectra have been described by means of collisional ionization equilibrium (CIE) models with several temperature components, often complemented by analysis of individual line profiles, line fluxes, and line flux ratios of He-like ions. Based on RGS observations of ζ Ori, Pollock (2007) suggests an alternative mechanism due to the non-equilibrium ion-ion interactions far out in the wind. Since extending the high-resolution coverage to the longer wavelengths available with the Chandra LETGS should provide new constraints on theory, it has been proposed to record the spectrum of δ Ori, the softest of the X-ray brightest O-stars observed with ROSAT. This paper discusses the LETGS data obtained in the context of individual line diagnostics to obtain plasma characteristics, as well as overall thermal models.

The O-star δ Ori (HD 36486) is the western star in the belt of Orion, named Mintaka, with a parallax of 4.71(0.58) mas, i.e., at a distance of 212 pc (SIMBAD: Van Leeuwen 2007)1. The position in J2000 coordinates is RA = 05 32 00.400and Dec = −00 17 56.74. It is a bright X-ray source and suffers lower interstellar absorption than ζ Ori. It is part of a triple star system, containing of δ OriA, δ OriB, and δ OriC. Here we focus on the “single star” δ OriA, which in itself is also a triple star system (Aa1, Aa2, and Ab). The tertiary star (Ab) is too far away to play a role in our observation. The primary (Aa1), a supergiant, is of spectral type O9.5II, with a temperature of 33 000 K. The secondary (Aa2) is of spectral type B0.5III, with a temperature of 27 000 K (Voels et al. 1989; Tarasov et al. 1995). They form a binary system with a period of 5.7325 days (Harvey et al. 1987). Different stellar mass values are given in the literature. Harvin et al. (2002) give Mp = 10.3 M as average with Rp = 11 R and for the secondary star Ms = 5.2 M and Rs = 4 R, while Mayer et al. (2010) suggest values of Mp = 25 M with Rs = 16−17 R. The primary and secondary stars are separated by 33 R, equivalent to 3 Rp. The properties of δ Ori, which we have applied, are given in Table 1.

Table 1

Adopted properties of δ Ori.

The spectrum of δ Ori was studied in X-rays with earlier satellites by Cassinelli & Swank (1983), Haberl & White (1993), and Corcoran et al. (1994) and most recently at high resolution by Waldron & Cassinelli (2007) and Miller et al. (2002) with the Chandra HETGS. They focused on line profiles and line widths, together with ratios of forbidden and intercombination lines in He-like ions.

The set up of this paper is as follows. The emphasis of the paper is on determining the location of the X-ray emitting plasma and its temperature in relation to the stellar surface and the UV-radiation field. Therefore we need individual line fluxes. In Sect. 2 the observation and the lightcurve are discussed. We show individual line fluxes, line wavelengths, and line broadenings in Sect. 3. Based on the individual line fluxes of He-like ions, the location relative to the stellar surface, as well as the temperature of the plasma, is derived. They are given in Sect. 4. We discuss a multitemperature model fit to the total spectrum with a CIE-model in Sect. 5. As suggested by the light curve, we also make a time split, dividing the total observation of 98 ks into the first part of 49 ks and the second part of 49 ks.

2. Observation and lightcurve

2.1. Data and pipelining

The X-ray spectrum of δ Ori was obtained by means of HRC-S in combination with the Low-Energy-Transmission-Grating (LETG) on board Chandra on November 9, 2007 during 97 ks (JD = 54 413.92705–54 415.05066). The log of the observation is shown in Table 2.

The data were pipelined with CIAO3.4 with CALDB version 3.4.1 and prepared for SPEX by the programs extract and crele developed at SRON. The spectrum derived from this observation is the sum of the +1 and −1 order. The background, which is quite dominant in the higher wavelength area of LETG spectra, has been subtracted. The spectrum contains 31 500 counts in the spectral range from 5 to 175 Å.

thumbnail Fig. 1

The zeroth-order lightcurve of δ Ori, observed with LETGS, and rebinned to time bins of 15 700 s. The curve indicates a minimum at the beginning of the observation and a flattening of the increase at the end of the observation.

Table 2

Observation log of the data of δ Ori.

thumbnail Fig. 2

Line shapes of the O viii line at 18.97 Å (top panel) and the C vi line at 33.74 Å in the spectrum of δ Ori. Data are in red/gray, the fitted Gaussian model line folded through the response matrix is given in black. The rest wavelength is also shown. Keep in mind that the width is also influenced by the instrumental line spread function.

2.2. Zeroth-order lightcurve

In other wavelength ranges (optical and IUE) fluctuations due to eclipsing have been observed in the lightcurve of δ Ori. These changes are about .12 m for the primary eclipse and .07m for the secondary eclipse (Harvin et al. 2002). This corresponds to a decrease of about 11% during the primary eclipse and of about 6% during the secondary eclipse. Based on the zeroth-order image a lightcurve has been constructed. Figure 1 shows the zeroth-order lightcurve of this observation of δ Ori. The total number of counts in this zeroth-order image is 26 350 counts. No periodicity could be established in the lightcurve, as this observation is too short to cover several stellar revolutions. However, after rebinning to time bins of about 16 ks the lightcurve shows some variation in count rate from 0.264(0.004) to 0.278(0.004) counts/s during this observation. This is a change of about 5%. During the observation the phase of the binary system was 0.83−0.02, derived from Harvin et al. (2002), or 0.87–0.07, as derived from Mayer et al. (2010). This implies that the secondary is in front of the primary during the observation. In later sections we investigate a possible relation between the behavior of the lightcurve and the observed spectrum by cutting the observation into a first part from 0 to 49 ks and a second part from 49 to 98 ks.

Table 3

Line positions and line fluxes in δ Ori.

Table 4

Higher-order lines and their line fluxes in the spectrum of δ Ori.

3. Individual line fluxes

The individual line fluxes of the LETGS spectrum of δ Ori were measured. For each observed line a Gaussian profile was folded through the response matrix and fit to the observed line profile, establishing the line positions, line fluxes, as well as the linewidths. A powerlaw was added to describe the continuum near the line features. The line profiles are symmetric. This is shown in Fig. 2 for the O VIII and C VI lines. This agrees with observations of δ Ori by Miller et al. (2002), when applying HETGS. Seeming asymmetries in the data are due to asymmetries in the instrumental line spread function (LSF). It also influences the observed widths of the lines. The situation in ζ Ori is different. In that star the lines are asymmetric and broadened by themselves as shown by Pollock (2007).

The measured wavelengths (Å) and line fluxes (in units 10-4 photons cm-2 s-1 and in units 10-13 erg cm-2 s-1) are collected in Table 3, together with the theoretical wavelengths from laboratory measurements collected by Kelly (1987) and the NIST4 database2. The stronger lines have a low value for the (1σ) uncertainty. Their uncertainties are a few mÅ.

There are very few line features above 50 Å, and they are often due to higher-order lines, especially third-order lines, which have their strong first-order component between 10 and 35 Å. The second-order lines of Fe xvii at 16.777, 17.054, and 17.097 Å contribute a bit to the wings of the C vi line at 33.736 Å. The third-order Ne x resonance line at 12.134 Å is blended by the S xii transition at 36.500 Å, and the second-order O viii resonance line contributes to the wing of the S xi and Si xiii feature at 37.787 Å. Some free unblended higher-order lines are given in Table 4. The higher order lines are implemented in the model fits of Sects. 5.1 and 5.2. First-order lines in the wavelength range above 50 Å belong to lower ionization stages of silicon and iron, and their optimal ionization temperature is around 0.1 keV.

Some stronger lines are collected in Table 5 and shown in Fig. 3. This table shows the line shifts and the line broadenings. The line shifts show velocities between −160 and +160 km s-1. The earth’s orbital velocity in the direction of δ Ori during the observation and the velocity of δ Ori itself (Wilson 1963) are both about 16 km s-1, but they have opposite signs and are canceling. The obtained velocities are comparable with those given by Miller et al. (2002), based on observations by HETGS. The half widths at half maximum (HWHMs) of the strong lines result in a velocity of about 840 km s-1, below the terminal velocity of 2000 km s-1. Our value is higher than the values given by Miller et al. (2002), based on HETGS. This might be due to an inconsistency between the instrumental LSFs of LETGS and HETGS. However, our value agrees with Wojdowski & Schulz (2005), who also used HETGS. The ultimate velocity at the end of the wing, corrected for the instrumental width are also given. These values are about 1850 km s-1 with an uncertainty of 200 km s-1. The lines in the separated spectra obtained from the two time intervals (0–49 ks and 49–98 ks) were measured, and no significant differences between the line measurements of the two time intervals have been established.

4. He-like line ratios

The individual line fluxes of the line triplets in the He-like ions, consisting of a resonance line (r), an intercombination line (i), and a forbidden line (f), are tools for density and temperature diagnostics (Gabriel & Jordan 1969), because the upper level of the forbidden line (1s2s 3S1) is depopulated by collisions in favor of the upper level of the intercombination line (1s2p 3P1). On the other hand, the same effect of depopulation can be caused by a strong UV radiation field (Blumenthal et al. 1972). The latter gives the size of the UV field flux and therewith an indication of the distance to the stellar surface where the X-ray emitting ions are formed. While for late type stars the collisions and therewith the electron density are the dominant processes, for hot, early type stars the strong UV radiation field (the location relative to the stellar surface) is most important. This has become a well-proven diagnostic tool for establishing the radial location of the X-ray source in the stellar wind of hot stars (Miller et al. 2002; Leutenegger et al. 2006; Waldron & Cassinelli 2007). The appropriate features in the He-like ions of C v, N vi, O vii, Ne ix, Mg xi, and Si xiii lie in the range from 5 to 45 Å, which is covered by LETGS.

The line fluxes of these special cases have been measured a slightly different way from the other lines given in Table 3. After the final overall multitemperature CIE-model fit (see Sect. 5), limited wavelength ranges were selected just around the relevant line features of the He-like ions. Then the lines from the ions under investigation were ignored from the model. In this way the continuum is based on the multitemperature fit and possible blends are taken into account. From that point on the same procedure was used as in Sect. 4: a Gaussian was folded through the response matrix and fitted to the line profiles of each ion line individually. The line flux ratios (f/i) obtained in this way are collected in Table 6. They might differ from ratios derived from Table 3, as in that table total feature fluxes, polluted by blends and satellite lines are given with their dominant identifications. The line flux ratios are compared with those given by Miller et al. (2002) using HETGS. Their values and ours are in good agreement. The ratios have been applied to a Planck UV-radiation field of 33 000 K at the wavelengths of the 1s2s3S–1s2p3P transitions of the corresponding ions, in combination with scaled formulas given by Blumenthal et al. (1972). This UV radiation field decreases with distance from the stellar surface by a dilution factor: W ( r ) = 1 2   [ 1 ( 1 ( R r ) 2 ) 1 / 2 ] · \begin{equation} W(r)=\frac{1}{2}~\left[1-\left(1-\left(\frac{R_{*}}{r}\right)^2\right)^{1/2}\right]\cdot \end{equation}(1)As a result of the radial dependence of the radiation field, the observed f/i ratio can be used to derive the radial location of the He-like ions that are producing the observed fir lines.

thumbnail Fig. 3

Part of the spectrum of δ Ori. The most dominant lines are indicated. Data are in red (gray), the CIE-model in black.

Table 5

Wavelength shifts and half width at half maximum (HWHM) values for some strong isolated lines in the spectrum of δ Ori.

Table 6

He-like and H-like over He-like line flux ratios and corresponding distances in R and temperatures.

thumbnail Fig. 4

The f/i ratio versus the distance to the stellar surface in stellar radii (R ∗ ).

The influence of the UV radiation field on the population of the 1s2s 3S1 and 1s2p 3P1 levels is reflected in the f/i ratio, hence in the distance to the stellar surface as shown in Fig. 4. The shapes of the curves are the same as those given by Miller et al. (2002) for Mg xi, Ne ix, and O vii. The curves of N vi and C v are new and established as a result of our investigation. The bold intervals on the trajectory show the average distance from the stellar surface (including the 1σ statistical error) where the ion is formed. From this figure it is clear that the hotter ions (Mg xi and Ne ix) are formed closer to the stellar surface than the cooler ones (C v and N vi). In the short wavelengths our results are in good agreement with Fig. 4 in the paper by Miller et al. (2002) and with the values of Mg xi and O vii by Waldron & Cassinelli (2007), based on HETGS observations. At the longer wavelengths, N vi and C v (outside the HETGS range) show that the X-rays are formed much farther out in the wind than the wind acceleration zone, which stops at about 10 R. The N vi is located at 49(13)   R and C v at 75(14)   R with temperatures of 1.6 MK and 1.0 MK respectively (see Table 6).

Only an upper value has been determined for Si xiii. This feature is better imaged by HETGS, and therefore the distance of its plasma to the stellar surface is excluded from Table 6 and Fig. 4. Also included in Table 6 are the temperature indicators G = (f + i)/r of the He-like ions, as well as the ratios of the H-like resonance lines and the He-like line triplets. They result in temperatures for the X-ray emitting plasma, related to the location of the specific ions in the plasma.

5. Multitemperature fitting

Along with the individual line analysis we investigated the LETG-spectrum of δ Ori using a multithermal model for optically thin plasma in CIE as implemented in SPEX (Kaastra et al. 1996a) in combination with MEKAL (Mewe et al. 1985, 1995; Kaastra et al. 1996b)3. It provides a calculation with thousands of lines and a model electron continuum. The ionization equilibrium is based on calculations by Arnaud & Rothenflug (1985) and Arnaud & Raymond (1992), the latter especially for iron. Application of a non equilibrium ionization model did not improve the description of the spectrum.

Table 7

The 3-T fit to the LETGS spectrum of δ Ori, applying a CIE-model to the total (time) spectrum.

The fits were applied using the C-statistics approach. For the transition probabilities of the model line features, an uncertainty of 15% was assumed, a value that accords with the NIST database4. The first calculation was fitted to the data, taking the errors on the data count rate into account. In the follow up of the fitting procedure the errors on the model count rate were applied to avoid over-influences of incidental low count rate bins.

5.1. A 3-T model

We started the fitting procedure with more fixed temperature components, out of which only three components were significant. From that point on a three temperature (3-T) fit was applied for the CIE-model. The values of the fit are collected in Table 7. This table shows the interstellar column density NH, the obtained temperatures, the emission measures, the luminosities of the different temperature bins, absorption column densities for the three temperature components in the wind, abundances, line shifts (z), and line broadening (vmic).

Three statistically well-determined temperature components at 0.10, 0.23, and 0.65 keV have been established. The two coolest components correspond to the lines of the He-like ions, discussed in Sect. 4. According to conclusions from that section, the hottest component (0.65 keV), that is responsible for Si xiv and S xv is formed deeper in the wind.

The results of the total 3-T fit to the spectrum are shown in Fig. 5. The lines and continua from the three individual temperature components are shown in the three top panels of Fig. 6. The bottom panel shows the data with the calculated continuum as sum of the three continua of these components.

thumbnail Fig. 5

The spectrum of δ Ori (red/gray), together with the model (black) derived from the multitemperature fit applying the parameters given in Table 7.

thumbnail Fig. 6

The spectrum of δ Ori (red/gray), together with the three individual temperature contributions to the model (black), applying the parameters given in Table 7. They imply the line and continuum radiation from the individual components. The bottom panel shows the model electron continuum (black). This implies the sum of the three continua.

The tabulated emission measure is defined as EM = n e n H V, \begin{equation} EM=n_{\rm e}n_{\rm H}V, \end{equation}(2)in which nH = 0.85 × ne.

Knowledge of the electron density offers the possibility to calculate the X-ray emitting volume. The electron density for hot stars with a spherically symmetric stellar wind can be obtained via the following continuity equation: = 4 π R 2 ρ ( R ) V ( R ) . \begin{equation} \dot M = 4\pi R^2 \rho(R) V(R). \end{equation}(3)In this formula is the mass loss rate (1.07 × 10-6 M/year, Lamers & Leitherer 1993) for δ Ori, R is the distance from the stellar surface, ρ(R) the density at distance R, and V(R) the wind velocity at distance R, related to the terminal velocity (V) by V ( R ) = V × ( 1 R / R ) β , \begin{equation} V(R) = V_{\infty}\times (1 - R_*/R)^{\beta}, \end{equation}(4)in which β is a coefficient between 0.85 for subgiants and 1.1 for supergiants. We use a value of 1.0 for β.

The InterStellar Medium (ISM) NH was fitted and found to be 1.5 × 1020 cm-2. Three additional wind absorption parameters were applied using the hot module in SPEX. It describes a warm plasma, as the CIE-model, while calculating the continuum absorption, as well as the line absorption. Emission lines with the highest oscillator strength are mostly influenced. The temperatures of this hot plasma were coupled to the temperatures of the components in the CIE-model. The hottest component, which is formed deeper in the wind (see Sect. 4), results in the highest absorption column density. No absorption for the cool wind (33 000 K) could be determined. This implies that the overall cool stellar wind absorption is not taken into account and might have been taken over partly by the fitted ISM-parameter NH.

A slight shift in λ (≈ −2.3e-4) was applied to optimize the fitted wavelengths. A line broadening parameter vmic of 1039 km s-1 was added and fitted to the overall spectrum to take into account the broadening of the line features. The HWHM of a line is related to the vmic by the following formula (see also Sect. 3): HWHM / λ obs = ln 2 v mic c · \begin{equation} {\it HWHM}/\lambda_{\rm obs} =\frac{\sqrt{\ln 2}v_{\rm mic}}{c}\cdot \end{equation}(5)The X-ray luminosities of the three individual temperature components are given over the LETG energy range from 0.07 to 3.0 keV. For the luminosities no individual statistical errors have been determined. Their uncertainties are supposed to have the same relative values as those for the emission measures because the luminosity is mostly determined by the emission measures in this temperature range, and therefore LX ∝ EM. The total luminosities (summed over the three temperature components) in other energy bands (0.5–4 keV, 0.1–2.4 keV, and 0.4–4 keV) for Einstein, ROSAT, ASCA given in the papers by Cassinelli & Swank (1983), Haberl & White (1993), and Corcoran et al. (1994), respectively, do agree with our values when taking differences in interstellar absorption or distance into account.

The abundances are relative to solar photospheric values (Lodders et al. 2009). The absolute values of the abundances are strongly influenced by the obtained values of the emission measures. The product of abundances and emission measures (A × EM), however, is the robust quantity for the line spectrum, therefore, abundance-ratios are often given relative to an abundant element in the stellar atmosphere. Here the abundance of iron was fixed at its solar photospheric value. The abundance values agree with the value of  ≈ 0.9 that Simón-Díaz (2010) finds for oxygen and silicon abundances in the Orion star-forming region.

A fit to the spectrum taken by HETGS applying the same set of parameters needs a lower value for vmic (see also Sect. 3). This might be due to inconsistency between the LSFs of the used instruments. Also the description of the Si xiii and Si xiv line features is not satisfactory. This might be because these features are near the end of the calibrated area for both instruments.

To investigate the time dependence of the spectrum, as suggested by the light curve (see Fig. 1) we split the observation into two time intervals, 0–49 ks and 49–98 ks. The emission measure of the hottest component (formed closest to the stellar surface) is about 10% lower in the first time interval than in the last 49 ks. This might agree with some eclipsing by the companion star.

5.2. DEM modeling

To show the smooth connection between the separated temperature bins, we performed a continuous DEM-modeling for the spectrum of δ Ori by applying the regularization method (Kaastra et al. 1996b) in the SPEX-code. This method uses direct matrix inversion with the additional constraint that the second-order derivative of the solution with respect to the temperature is as smooth as possible. In this DEM modeling the abundances given in Table 7 have been used. The results are shown in Fig. 7.

thumbnail Fig. 7

DEM-modeling of δ Ori, based on the regularization method (Kaastra et al. 1996b) see text. The EM is given in units of 1054 cm-3.

The results of the fits, collected in Table 7 and Fig. 7, are in good agreement. The peak of our DEM agrees with the peaks found by Miller et al. (2002) and Wojdowski & Schulz (2005). However, our DEM extends more to lower temperatures and their values extend further to higher temperatures as a consequence of using two different instruments. The peak value given by Wojdowski & Schulz (2005) is higher.

As for the CIE-fit we also performed a DEM-modeling for the partial spectra of 0–49 ks and 49–98 ks, again applying the regularization method. These two DEM-modelings are shown in Fig. 8. Although the hottest component is again a bit stronger for the last 49 ks, no significant difference can be established between fits to the partial spectra.

6. Discussion and conclusions

The X-ray spectrum of the O-star δ Ori obtained with LETGS is rich in lines belonging to H-like and He-like ions of C, N, O, Ne, Mg, and Si. A weak He-like (S xv) triplet is also present. Other possible weak features present below 5 Å fall outside the calibrated bandwidth of the LETGS. Iron L-shell lines are dominant. The Ne-like transitions, between electron configurations of an inert noble gas, are favored for iron (Fe XVII) as well as for calcium (Ca XI), and nickel (Ni XIX).

thumbnail Fig. 8

DEM-modeling of δ Ori, based on the regularization method. The two DEMs show a modeling of the spectrum over the time interval 0–49 ks (red/gray) and 49–98 ks (black).

Above 50 Å a few new line features were identified, some of which are third-order lines associated with strong first-order features located between 15 and 35 Å. No lines of cool ions were detected at wavelengths higher than 70 Å. The main reasons for the absence of these line features is the absence of a noticeable volume of cool (0.05 keV) plasma and the interstellar absorption.

Based on the forbidden line over intercombination line ratio in He-like ions the average distance of the X-ray emitting plasma to the stellar surface has been determined. Depending on the ion (temperature), the plasma emitting region is between 2 and 80 R for ions from Mg xi to C v. The cooler plasma (lower ionization stages) is located farther out in the wind. For the hot ions, Mg xi, Ne ix, and O vii, the average X-ray emitting plasma is from 2 to 10 R, respectively. This is in accordance with observations by Waldron & Cassinelli (2007) and Miller et al. (2002) using HETGS. Using LETGS that extends to longer wavelengths, the cooler ions N vi and C v have been observed as part of this investigation. The X-rays from these ions are formed on average at 50 R and 75 R, respectively, much farther out in the wind than the acceleration zone, where the final velocity has been reached. Here the wind is expected to be steady. However, small deviations in velocity (weak shocks) will be able to produce the cooler ions.

The spectrum is described well by a CIE model. Three well-distinguished temperature components have been determined: 0.10, 0.23, and 0.65 keV. The first two are in good agreement with the two temperatures found by Haberl & White (1993). They did not determine a third temperature component. The two lower temperature components (0.098 and 0.228 keV) describe the line-forming regions of the cooler ions, N vi and C v, and the moderately ionized species Mg xi, Ne ix, O vii, and O viii, respectively, while the hottest component is supposed to produce even hotter ions deeper in the wind and a hot continuum.

Three emission measures are established, related to the obtained temperatures. Knowledge of the electron density offers the possibility to calculate the X-ray emitting volume. Applying the distances (see Table 6 and Fig. 4) and the wind velocity and density formulas 4 and 5, given in Sect. 5.1, we have calculated the electron densities for the regions where the specific ions produce X-ray radiation.

Because the distances given in Table 6 and Fig. 4 are average values with the statistical error on this average and not the limits of the X-ray emitting distance, it is difficult to calculate the total wind volume in which the conditions are optimal for the corresponding He-like ion emission. By assuming a wind volume from half to twice the obtained average distance we get the percentage of the wind that contributes to that ion. The emitting volumes are about 1% of the wind volume for the cooler plasma and less for the hotter ions.

The NH of the ISM is mainly based on the long wavelength part of the spectrum. Its value of 1.51(.06) × 1020 cm-2 is in good agreement with the value of 1.56 × 1020 cm-2 found by Jenkins et al. (1999) based on 57 IUE observations, and slightly lower than the value found by Haberl & White (1993) for their T = 0.1 keV component. Additional wind absorption column densities have been determined (see Table 7). They are related to the temperature components of 0.10, 0.23, and 0.65 keV, with distances to the stellar surface of about 60 R, 6 R, and 2 R, respectively, where the plasma is formed on average. In combination with formula 4 we calculated that about 7% of the wind contributes to the hot wind absorption. No absorption of the overall cool stellar wind could be determined. This means that this absorption component is not explicitly taken into account during the fitting procedure. However, the effect might have been taken over by the other absorption parameters.

The X-ray luminosities of the three temperature components in the total LETGS band (0.07–3 keV) are 7.10, 5.14, and 2.59 × 1031 erg/s. Two partial spectra were derived: one from the first 49 ks and another from the last 49 ks of the observation. This division was suggested by the light curve (see Fig. 1). No significant differences between the spectra, their temperature structures, their line shapes, or line fluxes have been established.

The obtained abundances are relative to solar photospheric values (Lodders et al. 2009). They were fit by fixing the iron abundance at its solar photospheric value and leaving the carbon, nitrogen, oxygen, neon, magnesium, and silicon abundances free to vary. Fixing the iron abundance was done because the abundances and emission measures are strongly anticorrelated, especially in spectra with a weak continuum to establish the emission measure. The obtained values for oxygen and silicon are in good agreement with values found by Simón-Díaz (2010) in the Orion star-forming region.

The lines are symmetric and their HWHMs are about 840 km s-1, far below the terminal velocity of 2000 km s-1, but higher than the values found by Miller et al. (2002) using HETGS. This might be due to inconsistencies between the instrumental line spread functions in the calibration files of HETGS and LETGS. The ultimate velocities, measured at the end of the wings, are about 1850 km s-1 with an uncertainty of 200 km s-1.

The measured line shifts given in Table 5, as well as the z-value given in Table 7, indicate that the velocity of the radiative plasma in the direction from or to the observer is low. This is well known for OB stars, especially among giants and supergiants (Waldron & Cassinelli 2007).

The conclusion should be that the radiation we observe is produced in plasma that has a very low velocity or moves perpendicular to the line of sight. There might be several explanations. The clumps, which collide with the wind and produce the X-ray radiation, have low velocities. In that case the relative velocity of the clumps to the wind is highest far out in the wind at terminal wind velocity. As a consequence the highest stages of ionization are produced far out in the wind, which contradicts our observation based on the forbidden-to-intercombination line ratios in He-like ions (see Fig. 4). Additionally, after the acceleration phase one does not expect velocity differences that are too high. The secondary star and its wind play important roles in the explanation of the δ Ori spectrum. However, this is discussed by Miller et al. (2002) and rejected. A third option is that the clumps screen their own X-ray radiation when they are in between the star and the observer. Only when the wind-clump collision is almost perpendicular to the line of sight is the radiation observable and spread over a narrow velocity range as derived from the HWHM. The effects of clumping and porosity on X-Ray emission-line profiles from hot-star winds have been described extensively by Owocki & Cohen (2006) and Oskinova et al. (2006).


Acknowledgments

The SRON National Institute for Space Research is supported financially by NWO.

References

  1. Arnaud, M., & Raymond, R. 1992, ApJ, 398, 394 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arnaud, M., & Rothenflug, J. 1985, A&AS, 60, 425 [NASA ADS] [Google Scholar]
  3. Berghöfer, T.W., Schmitt, J. H. M. M., & Cassinelli, J. P. 1996, A&AS, 118, 481 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Blumenthal, G. R., Drake, G. W. F., & Tucker, W. H. 1972, ApJ, 172, 205 [NASA ADS] [CrossRef] [Google Scholar]
  5. Cassinelli, J. P., & Swank, J. H. 1983, ApJ, 271, 681 [NASA ADS] [CrossRef] [Google Scholar]
  6. Cassinelli, J. P., Waldron, W. L., Sanders, W. T., et al. 1981, ApJ, 250, 677 [NASA ADS] [CrossRef] [Google Scholar]
  7. Cohen, D. H., Leutenegger, M. A., Grizzard, K. T., et al. 2006, MNRAS, 368, 1905 [NASA ADS] [CrossRef] [Google Scholar]
  8. Corcoran, M. F., Waldron, W. L., MacFarlane, J. J., et al. 1994, ApJ, 436, L95 [NASA ADS] [CrossRef] [Google Scholar]
  9. Dessart, L., & Owocki, S. P. 2005, A&A, 437, 657 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Gabriel, A. H., & Jordan, C. 1969, MNRAS, 145, 241 [NASA ADS] [CrossRef] [Google Scholar]
  11. Haberl, F., & White, N. 1993, A&A, 280, 519 [NASA ADS] [Google Scholar]
  12. Harnden, F. R., Branduardi, G., Elvis, M., et al. 1979, ApJ, 234, L51 [NASA ADS] [CrossRef] [Google Scholar]
  13. Harvey, A. S., Stickland, D. J., Horwarth, I. D., & Zuiderwijk, E. L. 1987, Observatory, 107, 205 [NASA ADS] [Google Scholar]
  14. Harvin, J. A., Gies, D. R., Bagnuolo, W. G., et al. 2002, ApJ, 565, 1216 [NASA ADS] [CrossRef] [Google Scholar]
  15. Jenkins, E. B., Tripp, T. M., Woźniak, P. R., Ulysses, J. S., & Sonneborn, G. 1999, ApJ, 520, 182 [NASA ADS] [CrossRef] [Google Scholar]
  16. Kaastra, J. S., Mewe, R., & Nieuwenhuijzen, H. 1996a, in UV and X-ray Spectroscopy of Astrophysical and Laboratory Plasmas (Tokyo: Universal Academy Press, Inc.), eds. K. Yamashita, & T. Watanabe, 411 [Google Scholar]
  17. Kaastra, J. S., Mewe, R., Liedahl, D. A., et al. 1996b, A&A, 314, 547 [NASA ADS] [Google Scholar]
  18. Kahn, S. M., Leutenegger, M. A., Cottam, J., et al. 2001, A&A, 365, L312 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Kelly, R. L. 1987, J. Phys. Chem. Ref. Data, 16, 1 [CrossRef] [Google Scholar]
  20. Lamers, H. J. G. L. M., & Leitherer, C. 1993, ApJ, 412, 771 [NASA ADS] [CrossRef] [Google Scholar]
  21. Leutenegger, M. A., Paerels, F. B. S., Kahn, S. M., & Cohen, D. H. 2006, ApJ, 650, 1096 [NASA ADS] [CrossRef] [Google Scholar]
  22. Lodders, K., Palme, H., & Gail, H.-P. 2009, Landolt-Bornstein, New Series, Astronomy and Astrophysics (Berlin: Springer Verlag), VI/4B, 560 [Google Scholar]
  23. Lucy, L. B. 1982, ApJ, 255, 286 [NASA ADS] [CrossRef] [Google Scholar]
  24. Lucy, L. B., & Solomon, P. M. 1970, ApJ, 159, 879 [NASA ADS] [CrossRef] [Google Scholar]
  25. Lucy, L. B., & White, R. 1980, ApJ, 241, 300 [NASA ADS] [CrossRef] [Google Scholar]
  26. Mayer, P., Harmanec, P., Wolf, M., Božic, H., & Šlechta, M. 2010, A&A, 520, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Mewe, R., Gronenschild, E. H. B. M., & van den Oord, G. H. J. 1985, A&AS, 62, 1 [Google Scholar]
  28. Mewe, R., Kaastra, J. S., & Liedahl, D. A. 1995, Legacy, 6, 16 [Google Scholar]
  29. Miller, N. A., Cassinelli, J. P., Waldron, W. L., MacFarlane, J. J., & Cohen, D. H. 2002, ApJ, 577, 951 [NASA ADS] [CrossRef] [Google Scholar]
  30. Oskinova, L. M., Feldmeier, A., & Hamann, W.-R. 2006, MNRAS, 372, 313 [NASA ADS] [CrossRef] [Google Scholar]
  31. Owocki, S. P., & Cohen, D. H. 2006, ApJ, 648, 565 [NASA ADS] [CrossRef] [Google Scholar]
  32. Owocki, S. P., Castor, J. I., & Rybicki, G. B. 1988, ApJ, 335, 914 [NASA ADS] [CrossRef] [Google Scholar]
  33. Pollock, A. M. T. 2007, A&A, 463, 1111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Raassen, A. J. J., van der Hucht, K. A., Miller, N. A., & Cassinelli, J. P. 2008, A&A, 478, 513 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Seward, F. D., Forman, W. R., Giacconi, R., et al. 1979, ApJ, 234, L55 [NASA ADS] [CrossRef] [Google Scholar]
  36. Simón-Díaz, S. 2010, A&A, 510, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Tarasov, A. E., Harmanec, P., Horn, J., et al. 1995, A&AS, 110, 59 [NASA ADS] [Google Scholar]
  38. Van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Voels, S. A., Bohannan, B., Abbott, D. C., & Hummer, D. G. 1989, ApJ, 340, 1073 [NASA ADS] [CrossRef] [Google Scholar]
  40. Waldron, W. L., & Cassinelli, J. P. 2001, ApJ, 548, L45 [NASA ADS] [CrossRef] [Google Scholar]
  41. Waldron, W. L., & Cassinelli, J. P. 2007, ApJ, 668, 456 [NASA ADS] [CrossRef] [Google Scholar]
  42. Wilson, R. E. 1963, in General Catalogue Of Stellar Radial Velocities (Carnegie Institution of Washington publication), 601 [Google Scholar]
  43. Wojdowski, P. S., & Schulz, N. S. 2005, ApJ, 627, 953 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Adopted properties of δ Ori.

Table 2

Observation log of the data of δ Ori.

Table 3

Line positions and line fluxes in δ Ori.

Table 4

Higher-order lines and their line fluxes in the spectrum of δ Ori.

Table 5

Wavelength shifts and half width at half maximum (HWHM) values for some strong isolated lines in the spectrum of δ Ori.

Table 6

He-like and H-like over He-like line flux ratios and corresponding distances in R and temperatures.

Table 7

The 3-T fit to the LETGS spectrum of δ Ori, applying a CIE-model to the total (time) spectrum.

All Figures

thumbnail Fig. 1

The zeroth-order lightcurve of δ Ori, observed with LETGS, and rebinned to time bins of 15 700 s. The curve indicates a minimum at the beginning of the observation and a flattening of the increase at the end of the observation.

In the text
thumbnail Fig. 2

Line shapes of the O viii line at 18.97 Å (top panel) and the C vi line at 33.74 Å in the spectrum of δ Ori. Data are in red/gray, the fitted Gaussian model line folded through the response matrix is given in black. The rest wavelength is also shown. Keep in mind that the width is also influenced by the instrumental line spread function.

In the text
thumbnail Fig. 3

Part of the spectrum of δ Ori. The most dominant lines are indicated. Data are in red (gray), the CIE-model in black.

In the text
thumbnail Fig. 4

The f/i ratio versus the distance to the stellar surface in stellar radii (R ∗ ).

In the text
thumbnail Fig. 5

The spectrum of δ Ori (red/gray), together with the model (black) derived from the multitemperature fit applying the parameters given in Table 7.

In the text
thumbnail Fig. 6

The spectrum of δ Ori (red/gray), together with the three individual temperature contributions to the model (black), applying the parameters given in Table 7. They imply the line and continuum radiation from the individual components. The bottom panel shows the model electron continuum (black). This implies the sum of the three continua.

In the text
thumbnail Fig. 7

DEM-modeling of δ Ori, based on the regularization method (Kaastra et al. 1996b) see text. The EM is given in units of 1054 cm-3.

In the text
thumbnail Fig. 8

DEM-modeling of δ Ori, based on the regularization method. The two DEMs show a modeling of the spectrum over the time interval 0–49 ks (red/gray) and 49–98 ks (black).

In the text

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.