Refine
Year of publication
Document Type
- Article (134)
- Conference Proceeding (6)
- Other (6)
- Postprint (6)
- Doctoral Thesis (1)
- Review (1)
Keywords
- stars: massive (29)
- stars: Wolf-Rayet (25)
- stars: mass-loss (24)
- stars: atmospheres (21)
- stars: winds, outflows (21)
- stars: early-type (20)
- stars: evolution (13)
- Magellanic Clouds (11)
- X-rays: stars (9)
- binaries: close (9)
Institute
The long-term X-ray light curves of classical supergiant X-ray binaries and supergiant fast X-ray transients show relatively similar super-orbital modulations, which are still lacking a sound interpretation. We propose that these modulations are related to the presence of corotating interaction regions (CIRs) known to thread the winds of OB supergiants. To test this hypothesis, we couple the outcomes of three-dimensional (3D) hydrodynamic models for the formation of CIRs in stellar winds with a simplified recipe for the accretion onto a neutron star. The results show that the synthetic X-ray light curves are indeed modulated by the presence of the CIRs. The exact period and amplitude of these modulations depend on a number of parameters governing the hydrodynamic wind models and on the binary orbital configuration. To compare our model predictions with the observations, we apply the 3D wind structure previously shown to well explain the appearance of discrete absorption components in the UV time series of a prototypical B0.5I-type supergiant. Using the orbital parameters of IGRJ 16493-4348, which has the same B0.5I donor spectral type, the period and modulations in the simulated X-ray light curve are similar to the observed ones, thus providing support to our scenario. We propose that the presence of CIRs in donor star winds should be considered in future theoretical and simulation efforts of wind-fed X-ray binaries.
The information content of data on rotationally periodic recurrent discrete absorption components (DACs) in hot star wind emission lines is discussed. The data comprise optical depths tau(w,phi) as a function of dimensionless Doppler velocity w=(Deltalambda/lambda(0))(c/v(infinity)) and of time expressed in terms of stellar rotation angle phi. This is used to study the spatial distributions of density, radial and rotational velocities, and ionisation structures of the corotating wind streams to which recurrent DACs are conventionally attributed. The simplifying assumptions made to reduce the degrees of freedom in such structure distribution functions to match those in the DAC data are discussed and the problem then posed in terms of a bivariate relationship between tau(w, phi) and the radial velocity v(r)(r), transverse rotation rate Omega(r) and density rho(r, phi) structures of the streams. The discussion applies to cases where: the streams are equatorial; the system is seen edge on; the ionisation structure is approximated as uniform; the radial and transverse velocities are taken to be functions only of radial distance but the stream density is allowed to vary with azimuth. The last kinematic assumption essentially ignores the dynamical feedback of density on velocity and the relationship of this to fully dynamical models is discussed. The case of narrow streams is first considered, noting the result of Hamann et al. (2001) that the apparent acceleration of a narrow stream DAC is higher than the acceleration of the matter itself, so that the apparent slow acceleration of DACs cannot be attributed to the slowness of stellar rotation. Thus DACs either involve matter which accelerates slower than the general wind flow, or they are formed by structures which are not advected with the matter flow but propagate upstream (such as Abbott waves). It is then shown how, in the kinematic model approximation, the radial speed of the absorbing matter can be found by inversion of the apparent acceleration of the narrow DAC, for a given rotation law. The case of broad streams is more complex but also more informative. The observed tau(w,phi) is governed not only by v(r)(r) and Omega(r) of the absorbing stream matter but also by the density profile across the stream, determined by the azimuthal (phi(0)) distribution function F- 0(phi(0)) of mass loss rate around the stellar equator. When F-0(phi(0)) is fairly wide in phi(0), the acceleration of the DAC peak tau(w, phi) in w is generally slow compared with that of a narrow stream DAC and the information on v(r)(r), Omega(r) and F-0(phi(0)) is convoluted in the data tau(w, phi). We show that it is possible, in this kinematic model, to recover by inversion, complete information on all three distribution functions v(r)(r), Omega(r) and F- 0(phi(0)) from data on tau(w, phi) of sufficiently high precision and resolution since v(r)(r) and Omega(r) occur in combination rather than independently in the equations. This is demonstrated for simulated data, including noise effects, and is discussed in relation to real data and to fully hydrodynamic models
Two optically obscured Wolf-Rayet (WR) stars have been recently discovered by means of their infrared (IR) circumstellar shells, which show signatures of interaction with each other. Following the systematics of the WR star catalogues, these stars obtain the names WR 120bb and WR 120bc. In this paper, we present and analyse new near-IR, J-, H- and K-band spectra using the Potsdam Wolf-Rayet model atmosphere code. For that purpose, the atomic data base of the code has been extended in order to include all significant lines in the near-IR bands.
The spectra of both stars are classified as WN9h. As their spectra are very similar the parameters that we obtained by the spectral analyses hardly differ. Despite their late spectral subtype, we found relatively high stellar temperatures of 63 kK. The wind composition is dominated by helium, while hydrogen is depleted to 25 per cent by mass.
Because of their location in the Scutum-Centaurus Arm, WR 120bb and WR 120bc appear highly reddened, A(Ks) approximate to 2 mag. We adopt a common distance of 5.8 kpc to both stars, which complies with the typical absolute K-band magnitude for the WN9h subtype of -6.5 mag, is consistent with their observed extinction based on comparison with other massive stars in the region, and allows for the possibility that their shells are interacting with each other. This leads to luminosities of log(L/L-circle dot) = 5.66 and 5.54 for WR 120bb and WR 120bc, with large uncertainties due to the adopted distance.
The values of the luminosities of WR 120bb and WR 120bc imply that the immediate precursors of both stars were red supergiants (RSG). This implies in turn that the circumstellar shells associated with WR 120bb and WR 120bc were formed by interaction between the WR wind and the dense material shed during the preceding RSG phase.
We present an overview of four deep phase-constrained Chandra HETGS X-ray observations of delta Ori A. Delta Ori A is actually a triple system that includes the nearest massive eclipsing spectroscopic binary, delta Ori Aa, the only such object that can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, delta Ori Aa2, has a much lower X-ray luminosity than the brighter primary (delta Ori Aa1), delta Ori Aa provides a unique system with which to test the spatial distribution of the X-ray emitting gas around delta Ori Aa1 via occultation by the photosphere of, and wind cavity around, the X-ray dark secondary. Here we discuss the X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time of nearly 500 ks and covering nearly the entire binary orbit. The companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities obtained simultaneously with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectra. We find that the X-ray emission is dominated by embedded wind shock emission from star Aa1, with little contribution from the tertiary star Ab or the shocked gas produced by the collision of the wind of Aa1 against the surface of Aa2. We find a similar temperature distribution to previous X-ray spectrum analyses. We also show that the line half-widths are about 0.3-0.5 times the terminal velocity of the wind of star Aa1. We find a strong anti-correlation between line widths and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the strong and weak lines of Fe XVII and Ne X are inconsistent with model predictions, which may be an effect of resonance scattering.
Gamma 2 Velorum revisited
(1999)
Sk 183 is the visually brightest star in the N90 nebula, a young star-forming region in the Wing of the Small Magellanic Cloud (SMC). We present new optical spectroscopy from the Very Large Telescope which reveals Sk 183 to be one of the most massive O-type stars in the SMC. Classified as an O3-type dwarf on the basis of its nitrogen spectrum, the star also displays broadened He I absorption, which suggests a later type. We propose that Sk 183 has a composite spectrum and that it is similar to another star in the SMC, MPG 324. This brings the number of rare O2- and O3-type stars known in the whole of the SMC to a mere four. We estimate physical parameters for Sk 183 from analysis of its spectrum. For a single-star model, we estimate an effective temperature of 46 +/- 2 kK, a low mass-loss rate of similar to 10(-7) M-circle dot yr(-1), and a spectroscopic mass of 46(-8)(+ 9) M-circle dot (for an adopted distance modulus of 18.7 mag to the young population in the SMC Wing). An illustrative binary model requires a slightly hotter temperature (similar to 47.5 kK) for the primary component. In either scenario, Sk 183 is the earliest-type star known in N90 and will therefore be the dominant source of hydrogen-ionizing photons. This suggests Sk 183 is the primary influence on the star formation along the inner edge of the nebula.
We introduce a new survey of massive stars in the Galaxy and the Magellanic Clouds using the Fibre Large Array Multi- Element Spectrograph ( FLAMES) instrument at the Very Large Telescope ( VLT). Here we present observations of 269 Galactic stars with the FLAMES- Giraffe Spectrograph ( R similar or equal to 25 000), in fields centered on the open clusters NGC3293, NGC4755 and NGC6611. These data are supplemented by a further 50 targets observed with the Fibre- Fed Extended Range Optical Spectrograph ( FEROS, R = 48 000). Following a description of our scientific motivations and target selection criteria, the data reduction methods are described; of critical importance the FLAMES reduction pipeline is found to yield spectra that are in excellent agreement with less automated methods. Spectral classifications and radial velocity measurements are presented for each star, with particular attention paid to morphological peculiarities and evidence of binarity. These observations represent a significant increase in the known spectral content of NGC3293 and NGC4755, and will serve as standards against which our subsequent FLAMES observations in the Magellanic Clouds will be compared
We analyze the expansion of hydrogen-poor knots and filaments in the born-again planetary nebulae A30 and A78 based on Hubble Space Telescope ( HST) images obtained almost 20 yr apart. The proper motion of these features generally increases with distance to the central star, but the fractional expansion decreases, i.e., the expansion is not homologous. As a result, there is not a unique expansion age, which is estimated to be 610-950 yr for A30 and 600-1140 yr for A78. The knots and filaments have experienced complex dynamical processes: the current fast stellar wind is mass loaded by the material ablated from the inner knots; the ablated material is then swept up until it shocks the inner edges of the outer, hydrogen-rich nebula. The angular expansion of the outer filaments shows a clear dependence on position angle, indicating that the interaction of the stellar wind with the innermost knots channels the wind along preferred directions. The apparent angular expansion of the innermost knots seems to be dominated by the rocket effect of evaporating gas and by the propagation of the ionization front inside them. Radiation-hydrodynamical simulations show that a single ejection of material followed by a rapid onset of the stellar wind and ionizing flux can reproduce the variety of clumps and filaments at different distances from the central star found in A30 and A78.
Clumps in hot star winds can originate from shock compression due to the line driven instability. One-dimensional hydrodynamic simulations reveal a radial wind structure consisting of highly compressed shells separated by voids, and colliding with fast clouds. Two-dimensional simulations are still largely missing, despite first attempts. Clumpiness dramatically affects the radiative transfer and thus all wind diagnostics in the UV, optical, and in X-rays. The microturbulence approximation applied hitherto is currently superseded by a more sophisticated radiative transfer in stochastic media. Besides clumps, i.e. jumps in the density stratification, so-called kinks in the velocity law, i.e. jumps in dv/dr, play an eminent role in hot star winds. Kinks are a new type of radiative-acoustic shock, and propagate at super-Abbottic speed.
We discuss X-ray line formation in dense O star winds. A random distribution of wind shocks is assumed to emit X-rays that are partially absorbed by cooler wind gas. The cool gas resides in highly compressed fragments oriented perpendicular to the radial flow direction. For fully opaque fragments, we find that the blueshifted part of X-ray line profiles remains flat-topped even after severe wind attenuation, whereas the red part shows a steep decline. These box- type, blueshifted profiles resemble recent Chandra observations of the O3 star zeta Pup. For partially transparent fragments, the emission lines become similar to those from a homogeneous wind.
We discuss the connection between wind overloading and discrete absorption components in P Cygni line profiles from O stars. Overloading can create horizontal plateaus in the radial wind speed that cause the extra absorption in the line profile. The upstream propagation speed of these velocity plateaus is analyzed. The second part of the paper deals with X-ray emission from O stars. X-ray line profiles observed with Chandra and XMM are often symmetric, contrary to what is expected for lines from a homogeneous wind. We discuss the influence on line symmetry of photon escape channels in a strongly clumped wind.
Observations and theory suggest that line driven winds from hot stars and luminous accretion disks adopt a unique, critical solution which corresponds to maximum mass loss rate. We analyze the numerical stability of the infinite family of shallow wind solutions, which resemble solar wind breezes, and their transition to the critical wind. Shallow solutions are sub-critical with respect to radiative (or Abbott) waves. These waves can propagate upstream through shallow winds at high speeds. If the waves are not accounted for in the Courant time step, numerical runaway results. The outer boundary condition is equally important for wind stability. Assuming pure outflow conditions, as is done in the literature, triggers runaway of shallow winds to the critical solution or to accretion flow.
Context.
The supergiant ionized shell SMC-SGS 1 (DEM 167), which is located in the outer Wing of the Small Magellanic Cloud (SMC), resembles structures that originate from an energetic star-formation event and later stimulate star formation as they expand into the ambient medium. However, stellar populations within and surrounding SMC-SGS 1 tell a different story. Aims. We present a photometric study of the stellar population encompassed by SMC-SGS 1 in order to trace the history of such a large structure and its potential influence on star formation within the low-density, low-metallicity environment of the SMC.
Methods.
For a stellar population that is physically associated with SMC-SGS 1, we combined near-ultraviolet (NUV) photometry from the Galaxy Evolution Explorer with archival optical (V-band) photometry from the ESO Danish 1.54 m Telescope. Given their colors and luminosities, we estimated stellar ages and masses by matching observed photometry to theoretical stellar isochrone models. Results. We find that the investigated region supports an active, extended star-formation event spanning similar to 25-40 Myr ago, as well as continued star formation into the present. Using a standard initial mass function, we infer a lower bound on the stellar mass from this period of similar to 3 x 10(4) M-circle dot, corresponding to a star-formation intensity of similar to 6 x 10(-3) M-circle dot kpc(-2) yr(-1).
Conclusions.
The spatial and temporal distributions of young stars encompassed by SMC-SGS 1 imply a slow, consistent progression of star formation over millions of years. Ongoing star formation, both along the edge and interior to SMC-SGS 1, suggests a combined stimulated and stochastic mode of star formation within the SMC Wing. We note that a slow expansion of the shell within this low-density environment may preserve molecular clouds within the volume of the shell, leaving them to form stars even after nearby stellar feedback expels local gas and dust.
The NGC 346 young stellar system and associated N66 giant H ii region in the Small Magellanic Cloud are the nearest example of a massive star-forming event in a low metallicity (Z approximate to 0.2Z (circle dot)) galaxy. With an age of less than or similar to 3 Myr this system provides a unique opportunity to study relationships between massive stars and their associated H ii region. Using archival data, we derive a total H alpha luminosity of L(H alpha) = 4.1 x 10(38) erg s(-1) corresponding to an H-photoionization rate of 3 x 10(50) s(-1). A comparison with a predicted stellar ionization rate derived from the more than 50 known O-stars in NGC 346, including massive stars recently classified from Hubble Space Telescope far-ultraviolet (FUV) spectra, indicates an approximate ionization balance. Spectra obtained with SALT suggest the ionization structure of N66 could be consistent with some leakage of ionizing photons. Due to the low metallicity, the FUV luminosity from NGC 346 is not confined to the interstellar cloud associated with N66. Ionization extends through much of the spatial extent of the N66 cloud complex, and most of the cloud mass is not ionized. The stellar mass estimated from nebular L(H alpha) appears to be lower than masses derived from the census of resolved stars which may indicate a disconnect between the formation of high and low mass stars in this region. We briefly discuss implications of the properties of N66 for studies of star formation and stellar feedback in low metallicity environments.
Aims. To close this gap, we perform a comparative analysis of the optical companion in two important systems: IGR J175442619 (SFXT) and Vela X-1 (SGXB). We analyze the spectra of each star in detail and derive their stellar and wind properties. As a next step, we compare the wind parameters, giving us an excellent chance of recognizing key differences between donor winds in SFXTs and SGXBs. Methods. We use archival infrared, optical and ultraviolet observations, and analyze them with the non-local thermodynamic equilibrium (NLTE) Potsdam Wolf-Rayet model atmosphere code. We derive the physical properties of the stars and their stellar winds, accounting for the influence of X-rays on the stellar winds. Results. We find that the stellar parameters derived from the analysis generally agree well with the spectral types of the two donors: O9I (IGR J17544-2619) and B0.5Iae (Vela X-1). The distance to the sources have been revised and also agree well with the estimations already available in the literature. In IGR J17544-2619 we are able to narrow the uncertainty to d = 3.0 +/- 0.2 kpc. From the stellar radius of the donor and its X-ray behavior, the eccentricity of IGR J17544-2619 is constrained to e < 0.25. The derived chemical abundances point to certain mixing during the lifetime of the donors. An important difference between the stellar winds of the two stars is their terminal velocities (v(infinity) = 1500 km s(-1) in IGR J17544-2619 and v(infinity) = 700 km s(-1) in Vela X-1), which have important consequences on the X-ray luminosity of these sources. Conclusions. The donors of IGR J17544-2619 and Vela X-1 have similar spectral types as well as similar parameters that physically characterize them and their spectra. In addition, the orbital parameters of the systems are similar too, with a nearly circular orbit and short orbital period. However, they show moderate differences in their stellar wind velocity and the spin period of their neutron star which has a strong impact on the X-ray luminosity of the sources. This specific combination of wind speed and pulsar spin favors an accretion regime with a persistently high luminosity in Vela X-1, while it favors an inhibiting accretion mechanism in IGR J17544-2619. Our study demonstrates that the relative wind velocity is critical in class determination for the HMXBs hosting a supergiant donor, given that it may shift the accretion mechanism from direct accretion to propeller regimes when combined with other parameters.
Context. HD 93129A was classified as the earliest O-type star in the Galaxy (O2 If*) and is considered as the prototype of its spectral class. However, interferometry shows that this object is a binary system, while recent observations even suggest a triple configuration. None of the previous spectral analyses of this object accounted for its multiplicity. With new high-resolution UV and optical spectra, we have the possibility to reanalyze this key object, taking its binary nature into account for the first time. Aims. We aim to derive the fundamental parameters and the evolutionary status of HD 93129A, identifying the contributions of both components to the composite spectrum Results. Despite the similar spectral types of the two components, we are able to find signatures from each of the components in the combined spectrum, which allows us to estimate the parameters of both stars. We derive log(L/L-circle dot) = 6.15, T-eff = 52 kK, and log (M)over dot = -4.7[M-circle dot yr(-1)] for the primary Aa, and log(L/L-circle dot) = 5.58, T-eff = 45 kK, and log (M)over dot = -5.8 [M(circle dot)yr(-1)] for the secondary Ab. Conclusions. Even when accounting for the binary nature, the primary of HD 93129A is found to be one of the hottest and most luminous O stars in our Galaxy. Based on the theoretical decomposition of the spectra, we assign spectral types O2 If* and O3 III(f*) to components Aa and Ab, respectively. While we achieve a good fit for a wide spectral range, specific spectral features are not fully reproduced. The data are not sufficient to identify contributions from a hypothetical third component in the system.
We present the first non-LTE atmosphere models for WRstars that incorporate a self-consistent solution of the hydrodynamic equations. The models take iron-group line-blanketing and clumping into account, and compute the hydrodynamic structure of a radiatively driven wind consistently with the non-LTE radiation transport in the co-moving frame. We construct a self-consistent wind model that reproduces all observed properties of an early-type WCstar (WC5). We find that the WR-type mass-loss is initiated at high optical depth by the so-called "Hot Iron Bump" opacities (Fe IX- XVI). The acceleration of the outer wind regions is due to iron-group ions of lower excitation in combination with C and O. Consequently, the wind structure shows two acceleration regions, one close to the hydrostatic wind base in the optically thick part of the atmosphere, and another farther out in the wind. In addition to the radiative acceleration, the "Iron Bump" opacities are responsible for an intense heating of deep atmospheric layers. We find that the observed narrow O VI emission lines in the optical spectra of WC stars originate from this region. From their dependence on the clumping factor we gain important information about the location where the density inhomogeneities in WR-winds start to develop
Recent non-LTE models for expanding atmospheres, accounting for iron group line-blanketing and clumping, show a radiative acceleration which supplies a large part of the driving force of WR winds. Aiming at the calculation of fully consistent wind models, we developed a method to include the solution of the hydrodynamic equations into our code, taking into account the radiation pressure from the comoving-frame radiation transport. In the present work we discuss the resulting wind acceleration for WR- and O star models, and demonstrate the effects of clumping. In addition, we present a consistent hydrodynamic non-LTE model for the O-star zeta Puppis, which is calculated under consideration of complex model atoms of H, He, C, N, O, Si and the iron group elements. In its present state this model fails to reproduce the observed mass loss rate - probably due to still incomplete atomic data.
HST UV and optical spectra of the early-type [WC] star SMP 61 in the LMC are analyzed by means of line blanketed non-LTE models for expanding atmospheres. The known distance to the LMC allows a reliable determination of the stellar parameters. The low iron surface abundance of the object possibly indicates a preceding evolution through a very late thermal pulse (VLTP).
The blue compact H II galaxy CTS 1026 shows very strong WR emission features around 4686 AA and 5800 AA. We present high S/N optical spectra of the nucleus of this object. Byanalysis of the WR profile shapes, we determine the dominant spectral types and the WN/WC ratio in the starforming region. The ratio WR/O is determined via standard nebular diagnostics.
We describe the treatment of iron group line-blanketing in non-LTE model atmospheres for WR stars. As an example, a blanketed model for the early-type WC star WR 111 is compared to its un-blanketed counterpart. Blanketing affects the ionization structure and the emergent flux distribution of our models. The radiation pressure, as computed within our models, falls short by only a factor of two to provide the mechanical power of the WR wind.
We present XMM-Newton and Chandra observations of the born-again planetary nebula A 30. These X-ray observations reveal a bright unresolved source at the position of the central star whose X-ray luminosity exceeds by far the model expectations for photospheric emission and for shocks within the stellar wind. We suggest that a “born-again hot bubble” may be responsible for this X-ray emission. Diffuse X-ray emission associated with the petal-like features and one of the H-poor knots seen in the optical is also found. The weakened emission of carbon lines in the spectrum of the diffuse emission can be interpreted as the dilution of stellar wind by mass-loading or as the detection of material ejected during a very late thermal pulse.
The planetary nebula A30 is believed to have undergone a very late thermal pulse resulting in the ejection of knots of hydrogen-poor material. Using multi-epoch Hubble Space Telescope images, we have detected the angular expansion of these knots and derived an age of 850(-150)(+280) yr. To investigate the spectral and spatial properties of the soft X-ray emission detected by ROSAT, we have obtained Chandra and XMM-Newton deep observations of A30. The X-ray emission from A30 can be separated into two components: a point source at the central star and diffuse emission associated with the hydrogen-poor knots and the cloverleaf structure inside the nebular shell. To help us assess the role of the current stellar wind in powering this X-ray emission, we have determined the stellar parameters and wind properties of the central star of A30 using a non-LTE model fit to its optical and UV spectra. The spatial distribution and spectral properties of the diffuse X-ray emission are highly suggestive that it is generated by the post-born-again and present fast stellar winds interacting with the hydrogen-poor ejecta of the born-again event. This emission can be attributed to shock-heated plasma, as the hydrogen-poor knots are ablated by the stellar winds, under which circumstances the efficient mass loading of the present fast stellar wind raises its density and damps its velocity to produce the observed diffuse soft X-rays. Charge transfer reactions between the ions of the stellar winds and material of the born-again ejecta have also been considered as a possible mechanism for the production of diffuse X-ray emission, and upper limits on the expected X-ray production by this mechanism have been derived. The origin of the X-ray emission from the central star of A30 is puzzling: shocks in the present fast stellar wind and photospheric emission can be ruled out, while the development of a new, compact hot bubble confining the fast stellar wind seems implausible.
We report the first-ever discovery of a Wolf-Rayet (WR) star in the Large Magellanic Cloud via detection of a circular shell with the Spitzer Space Telescope. Follow-up observations with Gemini-South resolved the central star of the shell into two components separated from each other by a parts per thousand 2 arcsec (or a parts per thousand 0.5 pc in projection). One of these components turns out to be a WN3 star with H and He lines both in emission and absorption (we named it BAT99 3a using the numbering system based on extending the Breysacher et al. catalogue). Spectroscopy of the second component showed that it is a B0 V star. Subsequent spectroscopic observations of BAT99 3a with the du Pont 2.5-m telescope and the Southern African Large Telescope revealed that it is a close, eccentric binary system, and that the absorption lines are associated with an O companion star. We analysed the spectrum of the binary system using the non-LTE Potsdam WR (powr) code, confirming that the WR component is a very hot (a parts per thousand 90 kK) WN star. For this star, we derived a luminosity of log L/ L-aS (TM) = 5.45 and a mass-loss rate of 10(- 5.8) M-aS (TM) yr(- 1), and found that the stellar wind composition is dominated by helium with 20 per cent of hydrogen. Spectroscopy of the shell revealed an He iii region centred on BAT99 3a and having the same angular radius (a parts per thousand 15 arcsec) as the shell. We thereby add a new example to a rare class of high-excitation nebulae photoionized by WR stars. Analysis of the nebular spectrum showed that the shell is composed of unprocessed material, implying that the shell was swept-up from the local interstellar medium. We discuss the physical relationship between the newly identified massive stars and their possible membership of a previously unrecognized star cluster.
We report the discovery of two new Galactic candidate luminous blue variable (LBV) stars via detection of circular shells (typical of confirmed and candidate LBVs) and follow-up spectroscopy of their central stars. The shells were detected at 22 mu m in the archival data of the Mid-Infrared All Sky Survey carried out with the Wide-field Infrared Survey Explorer (WISE). Follow-up optical spectroscopy of the central stars of the shells conducted with the renewed Southern African Large Telescope (SALT) showed that their spectra are very similar to those of the well-known LBVs P Cygni and AG Car, and the recently discovered candidate LBV MN112, which implies the LBV classification for these stars as well. The LBV classification of both stars is supported by detection of their significant photometric variability: one of them brightened in the R and I bands by 0.68 +/- 0.10 and 0.61 +/- 0.04 mag, respectively, during the last 1318 years, while the second one (known as Hen 3-1383) varies its B, V, R, I and Ks brightnesses by similar or equal to 0.50.9 mag on time-scales from 10 d to decades. We also found significant changes in the spectrum of Hen 3-1383 on a time-scale of similar or equal to 3 months, which provides additional support for the LBV classification of this star. Further spectrophotometric monitoring of both stars is required to firmly prove their LBV status. We discuss a connection between the location of massive stars in the field and their fast rotation, and suggest that the LBV activity of the newly discovered candidate LBVs might be directly related to their possible runaway status.
We report the serendipitous discovery of a ring nebula around a candidate Wolf-Rayet (WR) star, HBHA 4202-22, in Cygnus using the Spitzer Space Telescope archival data. Our spectroscopic follow-up observations confirmed the WR nature of this star (we named it WR 138a) and showed that it belongs to the WN8-9h subtype. We thereby add a new example to the known sample of late WN stars with circumstellar nebulae. We analysed the spectrum of WR 138a by using the Potsdam Wolf-Rayet (PoWR) model atmospheres, obtaining a stellar temperature of 40 kK. The stellar wind composition is dominated by helium with 20 per cent of hydrogen. The stellar spectrum is highly reddened and absorbed (EB- V = 2.4 mag, A(V) = 7.4 mag). Adopting a stellar luminosity of log L/L-circle dot = 5.3, the star has a mass-loss rate of 10-4.7 M- circle dot yr-1, and resides in a distance of 4.2 kpc. We measured the proper motion for WR 138a and found that it is a runaway star with a peculiar velocity of similar or equal to 50 km s-1. Implications of the runaway nature of WR 138a for constraining the mass of its progenitor star and understanding the origin of its ring nebula are discussed.
We report the discovery of a new Wolf-Rayet star in Aquila via detection of its circumstellar nebula (reminiscent of ring nebulae associated with late WN stars) using the Spitzer Space Telescope archival data. Our spectroscopic follow-up of the central point source associated with the nebula showed that it is a WN7h star (we named it WR121b). We analysed the spectrum of WR 121b by using the Potsdam Wolf-Rayet model atmospheres, obtaining a stellar temperature of similar or equal to 50 kK. The stellar wind composition is dominated by helium with similar to 20 per cent of hydrogen. The stellar spectrum is highly reddened [E(B - V) = 2.85 mag]. Adopting an absolute magnitude of M-v = 5.7, the star has a luminosity of log L/L-circle dot = 5.75 and a mass-loss rate of 10(-4.7)M(circle dot)yr(-1), and resides at a distance of 6.3 kpc. We searched for a possible parent cluster of WR 121b and found that this star is located at similar or equal to 1 degrees from the young star cluster embedded in the giant HII region W43 (containing a WN7+a/OB? star - WR121a). We also discovered a bow shock around the O9.5III star ALS 9956, located at similar or equal to 0 degrees.5 from the cluster. We discuss the possibility that WR121b and ALS 9956 are runaway stars ejected from the cluster in W43.
Context:
The first directly detected gravitational waves (GW 150914) were emitted by two coalescing black holes (BHs) with masses of ≈ 36 M⊙ and ≈ 29 M⊙. Several scenarios have been proposed to put this detection into an astrophysical context. The evolution of an isolated massive binary system is among commonly considered models.
Aims:
Various groups have performed detailed binary-evolution calculations that lead to BH merger events. However, the question remains open as to whether binary systems with the predicted properties really exist. The aim of this paper is to help observers to close this gap by providing spectral characteristics of massive binary BH progenitors during a phase where at least one of the companions is still non-degenerate.
Methods:
Stellar evolution models predict fundamental stellar parameters. Using these as input for our stellar atmosphere code (Potsdam Wolf-Rayet), we compute a set of models for selected evolutionary stages of massive merging BH progenitors at different metallicities.
Results:
The synthetic spectra obtained from our atmosphere calculations reveal that progenitors of massive BH merger events start their lives as O2-3V stars that evolve to early-type blue supergiants before they undergo core-collapse during the Wolf-Rayet phase. When the primary has collapsed, the remaining system will appear as a wind-fed high-mass X-ray binary. Based on our atmosphere models, we provide feedback parameters, broad band magnitudes, and spectral templates that should help to identify such binaries in the future.
Conclusions:
While the predicted parameter space for massive BH binary progenitors is partly realized in nature, none of the known massive binaries match our synthetic spectra of massive BH binary progenitors exactly. Comparisons of empirically determined mass-loss rates with those assumed by evolution calculations reveal significant differences. The consideration of the empirical mass-loss rates in evolution calculations will possibly entail a shift of the maximum in the predicted binary-BH merger rate to higher metallicities, that is, more candidates should be expected in our cosmic neighborhood than previously assumed.
Context. Wolf-Rayet (WR) stars have a severe impact on their environments owing to their strong ionizing radiation fields and powerful stellar winds. Since these winds are considered to be driven by radiation pressure, it is theoretically expected that the degree of the wind mass-loss depends on the initial metallicity of WR stars.
Aims. Following our comprehensive studies of WR stars in the Milky Way, M31, and the LMC, we derive stellar parameters and mass-loss rates for all seven putatively single WN stars known in the SMC. Based on these data, we discuss the impact of a low-metallicity environment on the mass loss and evolution of WR stars.
Methods. The quantitative analysis of the WN stars is performed with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The physical properties of our program stars are obtained from fitting synthetic spectra to multi-band observations.
Results. In all SMC WN stars, a considerable surface hydrogen abundance is detectable. The majority of these objects have stellar temperatures exceeding 75 kK, while their luminosities range from 10(5.5) to 10(6.1) L-circle dot. The WN stars in the SMC exhibit on average lower mass-loss rates and weaker winds than their counterparts in the Milky Way, M31, and the LMC.
Conclusions. By comparing the mass-loss rates derived for WN stars in different Local Group galaxies, we conclude that a clear dependence of the wind mass-loss on the initial metallicity is evident, supporting the current paradigm that WR winds are driven by radiation. A metallicity effect on the evolution of massive stars is obvious from the HRD positions of the SMC WN stars at high temperatures and high luminosities. Standard evolution tracks are not able to reproduce these parameters and the observed surface hydrogen abundances. Homogeneous evolution might provide a better explanation for their evolutionary past.
The study of massive stars in different metallicity environments is a central topic of current stellar research. The spectral analysis of massive stars requires adequate model atmospheres. The computation of such models is difficult and time-consuming. Therefore, spectral analyses are greatly facilitated if they can refer to existing grids of models. Here we provide grids of model atmospheres for OB-type stars at metallicities corresponding to the Small and Large Magellanic Clouds, as well as to solar metallicity. In total, the grids comprise 785 individual models. The models were calculated using the state-of-the-art Potsdam Wolf-Rayet (PoWR) model atmosphere code. The parameter domain of the grids was set up using stellar evolution tracks. For all these models, we provide normalized and flux-calibrated spectra, spectral energy distributions, feedback parameters such as ionizing photons, Zanstra temperatures, and photometric magnitudes. The atmospheric structures (the density and temperature stratification) are available as well. All these data are publicly accessible through the PoWR website.
The Wolf-Rayet stars in the Large Magellanic Cloud - A comprehensive analysis of the WN class
(2014)
Context. Massive stars, although being important building blocks of galaxies, are still not fully understood. This especially holds true for Wolf-Rayet (WR) stars with their strong mass loss, whose spectral analysis requires adequate model atmospheres. Aims. Following our comprehensive studies of the WR stars in the Milky Way, we now present spectroscopic analyses of almost all known WN stars in the LMC.
Methods. For the quantitative analysis of the wind-dominated emission-line spectra, we employ the Potsdam Wolf-Rayet (PoWR) model atmosphere code. By fitting synthetic spectra to the observed spectral energy distribution and the available spectra (ultraviolet and optical), we obtain the physical properties of 107 stars.
Results. We present the fundamental stellar and wind parameters for an almost complete sample of WN stars in the LMC. Among those stars that are putatively single, two different groups can be clearly distinguished. While 12% of our sample are more luminous than 10(6) L-circle dot and contain a significant amount of hydrogen, 88% of the WN stars, with little or no hydrogen, populate the luminosity range between log (L/L-circle dot) = 5.3 ... 5.8.
Conclusions. While the few extremely luminous stars (log (L/L-circle dot) > 6), if indeed single stars, descended directly from the main sequence at very high initial masses, the bulk of WN stars have gone through the red-supergiant phase. According to their luminosities in the range of log (L/L-circle dot) = 5.3 ... 5.8, these stars originate from initial masses between 20 and 40 M-circle dot. This mass range is similar to the one found in the Galaxy, i.e. the expected metallicity dependence of the evolution is not seen. Current stellar evolution tracks, even when accounting for rotationally induced mixing, still partly fail to reproduce the observed ranges of luminosities and initial masses. Moreover, stellar radii are generally larger and effective temperatures correspondingly lower than predicted from stellar evolution models, probably due to subphotospheric inflation.
A detailed and comprehensive study of the Wolf-Rayet stars of the nitrogen sequence (WN
stars) in the Small Magellanic Cloud (SMC) and the Large Magellanic Cloud (LMC) is presented.
We derived the fundamental stellar and wind parameters for more than 100 massive stars, encompassing almost the whole WN population in the Magellanic Clouds (MCs). The observations are fitted with synthetic spectra, using the PotsdamWolf-Rayet model atmosphere
code (PoWR). For this purpose, large grids of line-blanket models for different metallicities have been calculated, covering a wide range of stellar temperatures, mass-loss rates, and hydrogen abundances. Our comprehensive sample facilitates statistical studies of the WN properties in the MCs without selection bias. To investigate the impact of the low LMC metallicity and the even lower SMC metallicity, we compare our new results to previous analyses of the Galactic WN population and the late type WN stars from M31. Based on these studies we derived an empirical relation between the WN mass-loss rates and the metallicity. Current stellar evolution tracks, even when accounting for rotationally induced mixing, partly fail to reproduce the observed ranges of luminosities and initial masses.
SwSt 1 (PN G001.5-06.7) is a bright and compact planetary nebula containing a late [WC]-type central star. Previous studies suggested that the nebular and stellar lines are slowly changing with time. We studied new and archival optical and ultraviolet spectra of the object. The [O III] 4959 and 5007 angstrom to H beta line flux ratios decreased between about 1976 and 1997/2015. The stellar spectrum also shows changes between these epochs. We modelled the stellar and nebular spectra observed at different epochs. The analyses indicate a drop of the stellar temperature from about 42 kK to 40.5 kK between 1976 and 1993. We do not detect significant changes between 1993 and 2015. The observations show that the star performed a loop in the H-R diagram. This is possible when a shell source is activated during its post-AGB evolution. We infer that a late thermal pulse (LTP) experienced by a massive post-AGB star can explain the evolution of the central star. Such a star does not expand significantly as the result of the LTP and does not became a born-again red giant. However, the released energy can remove the tiny H envelope of the star.
The non-LTE radiative transfer problem requires the consistent solution of two sets of equations: the radiative transfer equations, which couple the spatial points, and the equations of the statistical equilibrium, which couple the frequencies. The "Accelerated Lambda Iteration" (ALI) method allows for an iterative scheme, in which both sets of equations are solved in turn. For moving atmospheres the radiative transfer is preferably formulated in the co-moving frame-of-reference, which leads to a partial differential equation. "Classical" numerical solution methods are based on differencing schemes. For better numerical stability, we prefer "short characteristics" integration methods. Iron line blanketing is accounted for by means of the "superlevel" concept. In contrast to static atmospheres, the frequencies can not be re-ordered in the moving case because of the frequency coupling from Doppler shifts. One of our future aims is the coupling of elaborated radiative transfer calculations with the hydrodynamical equations in order to understand the driving of strong stellar winds, especially from Wolf-Rayet stars.
Stellar winds appear as a persistent feature of hot stars, irrespective of their wide range of different luminosities, masses, and chemical composition. Among the massive stars, the Wolf-Rayet types show considerably stronger mass loss than the O stars. Among hot low-mass stars, stellar winds are seen at central stars of planetary nebulae, where again the hydrogen-deficient stars show much stronger winds than those central stars with "normal" composition. We also studied mass-loss from a few extreme helium stars and sdOs. Their mass-loss rate roughly follows the same proportionality with luminosity to the power 1.5 as the massive O stars. This relation roughly marks a lower limit for the mass loss from hot stars of all kinds, and provides evidence that radiation pressure on spectral lines is the basic mechanism at work. For certain classes of stars the mass-loss rates lie significantly above this relation, for reasons that are not yet fully understood. Mass loss from low-mass stars may affect their evolution, by reducing the envelope mass, and can easily prevent diffusion from establishing atmospheric abundance patterns. In close binary systems, their winds can feed the accretion onto a companion.
The emission-line dominated spectra of Wolf-Rayet stars are formed in expanding layers of their atmosphere, i.e. in their strong stellar wind. Adequate modeling of such spectra has to face a couple of difficulties. Because of the supersonic motion, the radiative transfer is preferably formulated in the co-moving frame. The strong deviations from local thermodynamical equilibrium (LTE) require to solve the equations of statistical equilibrium for the population numbers, accounting for many hundred atomic energy levels and thousands of line transitions. Moreover, millions of lines from iron-group elements must be taken into account for their blanketing effect. Model atmospheres of the described kind can reproduce the observed WR spectra satisfyingly, and have been widely applied for corresponding spectral analyses.
Spectral lines formed in stellar winds from OB stars are observed to exhibit profile variations. Discrete Absorption Components (DACs) show a remarkably slow wavelength drift with time. In a straightforward interpretation, this is in sharp contradiction to the steep velocity law predicted by the radiation-driven wind theory, and by semi- empirical profile fitting. In the present paper we re-discuss the interpretation of the drift rate. We show that the Co- rotating Interaction Region (CIR) model for the formation of DACs does not explain their slow drift rate as a consequence of rotation. On the contrary, the apparent acceleration of a spectral CIR feature is even higher than for the corresponding kinematical model without rotation. However, the observations can be understood by distinguishing between the velocity field of the matter flow, and the velocity law for the motion of the patterns in which the DAC features are formed. If the latter propagate upstream against the matter flow, the resulting wavelength drift mimics a much slower acceleration although the matter is moving fast. Additional to the DACs, a second type of recurrent structures is present in observed OB star spectra, the so-called modulations. In contrast to the DACs, these structures show a steep acceleration compatible with the theoretically predicted velocity law. We see only two possible consistent scenarios. Either, the wind is accelerated fast, and the modulations are formed in advected structures, while the DACs come from structures which are propagating upstream. Or, alternatively, steep and shallow velocity laws may co-exist at the same time in different spatial regions or directions of the wind.
Context. Very massive stars pass through the Wolf-Rayet (WR) stage before they finally explode. Details of their evolution have not yet been safely established, and their physics are not well understood. Their spectral analysis requires adequate model atmospheres, which have been developed step by step during the past decades and account in their recent version for line blanketing by the millions of lines from iron and iron-group elements. However, only very few WN stars have been re-analyzed by means of line-blanketed models yet. Aims. The quantitative spectral analysis of a large sample of Galactic WN stars with the most advanced generation of model atmospheres should provide an empirical basis for various studies about the origin, evolution, and physics of the Wolf-Rayet stars and their powerful winds. Methods. We analyze a large sample of Galactic WN stars by means of the Potsdam Wolf-Rayet (PoWR) model atmospheres, which account for iron line blanketing and clumping. The results are compared with a synthetic population, generated from the Geneva tracks for massive star evolution. Results. We obtain a homogeneous set of stellar and atmospheric parameters for the GalacticWN stars, partly revising earlier results. Conclusions. Comparing the results of our spectral analyses of the Galactic WN stars with the predictions of the Geneva evolutionary calculations, we conclude that there is rough qualitative agreement. However, the quantitative discrepancies are still severe, and there is no preference for the tracks that account for the effects of rotation. It seems that the evolution of massive stars is still not satisfactorily understood.
Comprehensive spectral analyses of the Galactic Wolf-Rayet stars of the nitrogen sequence (i.e. the WN subclass) have been performed in a previous paper. However, the distances of these objects were poorly known. Distances have a direct impact on the "absolute" parameters, such as luminosities and mass-loss rates. The recent Gaia Data Release (DR2) of trigonometric parallaxes includes nearly all WN stars of our Galactic sample. In the present paper, we apply the new distances to the previously analyzed Galactic WN stars and rescale the results accordingly. On this basis, we present a revised catalog of 55 Galactic WN stars with their stellar and wind parameters. The correlations between mass-loss rate and luminosity show a large scatter, for the hydrogen-free WN stars as well as for those with detectable hydrogen. The slopes of the log L - log M correlations are shallower than found previously. The empirical Hertzsprung-Russell diagram (HRD) still shows the previously established dichotomy between the hydrogen-free early WN subtypes that are located on the hot side of the zero-age main sequence (ZAMS), and the late WN subtypes, which show hydrogen and reside mostly at cooler temperatures than the ZAMS (with few exceptions). However, with the new distances, the distribution of stellar luminosities became more continuous than obtained previously. The hydrogen-showing stars of late WN subtype are still found to be typically more luminous than the hydrogen-free early subtypes, but there is a range of luminosities where both subclasses overlap. The empirical HRD of the Galactic single WN stars is compared with recent evolutionary tracks. Neither these single-star evolutionary models nor binary scenarios can provide a fully satisfactory explanation for the parameters of these objects and their location in the HRD.
Grids of model atmospheres for Wolf-Rayet stars of the nitrogen sequence (WN subclass) are presented. The calculations account for the expansion of the atmosphere, non-LTE, clumping, and line blanketing from iron-group elements. Observed spectra of single Galactic WN stars can in general be reproduced consistently by this generation of models. The parameters of the presented model grids cover the whole relevant range of stellar temperatures and mass-loss rates. We point out that there is a degeneracy of parameters for very thick winds; their spectra tend to depend only on the ratio $L/{dot M}^{4/3}$. Abundances of the calculated grids are for Galactic WN stars without hydrogen and with 20% hydrogen (by mass), respectively. Model spectra and fluxes are available via internet (http://www.astro.physik.uni- potsdam.de/PoWR.html).
Most Central Stars of Planetary Nebulae exhibit a spectrum of a hydrogen-rich hot star with little or no stellar wind. About 20 % of the CSPN, however, show entirely different spectra dominated by bright and broad emission lines of carbon, oxygen and helium, resembling the so-called Wolf-Rayet (WR) spectral class originally established for massive, Pop. I stars. These spectra indicate a hydrogen-deficient surface composition and, at the same time, strong mass-loss. As the WR spectra are formed entirely in a dense stellar wind, their spectral analysis requires adequate modelling. Corresponding Non-LTE model atmospheres have been developed in the last decade and became more and more sophisticated. They have been applied yet for analyzing almost all available WR-type CSPN spectra, establishing the stellar parameters. The obtained surface abundances are not understandable in terms of "classical" evolutionary calculations, but agree in principle with the advanced models for AGB evolution which account consistently for diffusive mixing and nuclear burning. The underabundance of iron, which we established in a recent study of a WC-type central star (LMC-SMP 61), gives indirect evidence that neutron-capture synthesis has converted Fe into s-process elements.
The Potsdam Non-LTE code for expanding atmospheres, which accounts for clumping and iron-line blanketing, has been used to establish a grid of model atmospheres for WC stars. A parameter degeneracy is discovered for early-type WC models which do not depend on the "stellar temperature". 15 galactic WC4-7 stars are analyzed, showing a very uniform carbon abundance (He:C=55:40) with only few exceptions.
WR Central Stars
(2003)
Spectrum formation in clumped stellar winds : consequences for the analyses of Wolf-Rayet spectra
(1998)
The spectra of 18 WN stars in the Large Magellanic Cloud (LMC) are quantitatively analyzed by means of "standard" Wolf-Rayet model atmospheres, using the helium and nitrogen lines as well as the spectral energy distribution. The hydrogen abundance is also determined. Carbon is included for a subset of 4 stars. The studied sample covers all spectral subtypes (WN2 ... WN9) and also includes one WN/WC transition object. The luminosities of the program stars span a wide range ( L/Lsun = 5.0 ... 6.5). Due to the given LMC membership, these results are free from uncertainties inferred from the distance. 50 % of the studied stars (both, late and early WN subtypes) have rather low luminosity (L/Lsun < 5.5). This puts tough constraints on their evolutionary formation. If coming from single stars, it provides evidence for strong internal mixing processes. The empirical mass-loss rates are scaled down by a factor of about two due to the impact of clumping, compared to previous studies adopting homogeneous winds. There is no obvious strong correlation between the mass-loss rates and other parameters like luminosity, temperature and composition. The stellar parameters for the present LMC sample are not systematically different from those of the Galactic WN stars studied previously with the same techniques, in contrast to the expected metallicity effects.
The analysis of Wolf-Rayet spectra requires adequate model atmospheres which treat the non-LTE radiation transfer in a spherically expanding medium. Present state-of-the-art calculations account for complex model atoms with, typically, a few hundred energy levels and a few thousand spectral lines of He and CNO elements. In the most recent version of our model code, blanketing by millions of lines from iron-group elements is also included. These models have been widely applied for the spectral analysis of WN stars in the Galaxy and LMC. WN spectra can be well reproduced in most cases. WC stars have not yet been analyzed comprehensively, because the agreement with observations becomes satisfactory only when line-blanketed models are applied. The introduction of inhomogeneities (clumping), although treated in a rough approximation, has significantly improved the fit between synthetic and observed spectra with respect to the electron-scattering wings of strong lines. The mass-loss rates obtained from spectral analyses become smaller by a factor 2-3 if clumping is accounted for. A pre-specified velocity law is adopted for our models, but the radiation pressure can be evaluated from our detailed calculation and can be compared a posteriori with the required wind acceleration. Surprisingly we find that the line-blanketed models are not far from being hydrodynamically consistent, thus indicating that radiation pressure is probably the main driving force for the mass-loss from WR stars.
Non-LTE models of WR winds
(2000)
Modeling expanding atmospheres is a difficult task because of the extreme non-LTE situation, the need to account for complex model atoms, especially for the iron-group elements with their millions of lines, and because of the supersonic expansion. Adequate codes have been developed e.g. by Hillier (CMFGEN), the Munich group (Puls, Pauldrach), and in Potsdam (PoWR code, Hamann et al.). While early work was based on the assumption of a smooth and homogeneous spherical stellar wind, the need to account for clumping became obvious about ten years ago. A relatively simple first-order clumping correction was readily implemented into the model codes. However, its simplifying assumptions are severe. Most importantly, the clumps are taken to be optically thin at all frequencies (”microclumping”). We discuss the consequences of this approximation and describe an approach to account for optically thick clumps (“macroclumping”). First results demonstrate that macroclumping can generally reduce the strength of spectral features, depending on their optical thickness. The recently reported discrepancy between the Hα diagnostic and the Pv resonance lines in O star spectra can be resolved without decreasing the mass-loss rates, when macroclumping is taken into account.
To investigate if magnetic fields are present in Wolf-Rayet stars, we selected a few stars in the Galaxy and one in the Large Magellanic Cloud (LMC). We acquired low-resolution spectropolarimetric observations with the European Southern Observatory FORS 2 (FOcal Reducer low dispersion Spectrograph) instrument during two different observing runs. During the first run in visitor mode, we observed the LMC Wolf-Rayet star BAT99 7 and the stars WR 6, WR 7, WR 18, and WR 23 in our Galaxy. The second run in service mode was focused on monitoring the star WR 6. Linear polarization was recorded immediately after the observations of circular polarization. During our visitor observing run, the magnetic field for the cyclically variable star WR 6 was measured at a significance level of 3.3 sigma (< B-z > = 258 +/- 78 G). Among the other targets, the highest value for the longitudinal magnetic field, < B-z > = 327 +/- 141 G, was measured in the LMC star BAT99 7. Spectropolarimetric monitoring of the star WR 6 revealed a sinusoidal nature of the < B-z > variations with the known rotation period of 3.77 d, significantly adding to the confidence in the detection. The presence of the rotation-modulated magnetic variability is also indicated in our frequency periodogram. The reported field magnitude suffers from significant systematic uncertainties at the factor of 2 level, in addition to the quoted statistical uncertainties, owing to the theoretical approach used to characterize it. Linear polarization measurements showed no line effect in the stars, apart from WR 6. BAT99 7, WR 7, and WR 23 do not show variability of the linear polarization over two nights.
As WR 6 is a putatively single WN4 star, and is relatively bright (V = 6.9), it is an ideal case for studying the wind mechanisms in these extremely luminous stars. To obtain higher resolution spectra at higher energy (above 1 keV) than previously obtained with the XMM/Newton RGS, we have observed WR 6 with the Chandra High Energy Transmission Grating Spectrometer for 450 ks. We have resolved emission lines of S, Si, Mg, Ne, and Fe, which all show a “fin"-shaped prole, characteristic of a self-absorbed uniformly expanding shell. Steep blue edges gives robust maximal expansion velocities of about 2000 km/s, somewhat larger than the 1700km/s derived from UV lines. The He-like lines all indicate that X-ray emitting plasmas are far from the photosphere – even at the higher energies where opacity is lowest { as was also the case for the longer wavelength lines observed with XMM-Newton/RGS. Abundances determined from X-ray spectral modeling indicate enhancements consistent with nucleosynthesis. The star was also variable in X-rays and in simultaneous optical photometry obtained with Chandra aspect camera, but not coherently with the optically known period of 3.765 days.
With a deep Chandra/HETGS exposure of WR 6, we have resolved emission lines whose profiles show that the X-rays originate from a uniformly expanding spherical wind of high X-ray-continuum optical depth. The presence of strong helium-like forbidden lines places the source of X-ray emission at tens to hundreds of stellar radii from the photosphere. Variability was present in X-rays and simultaneous optical photometry, but neither were correlated with the known period of the system or with each other. An enhanced abundance of sodium revealed nuclear-processed material, a quantity related to the evolutionary state of the star. The characterization of the extent and nature of the hot plasma in WR 6 will help to pave the way to a more fundamental theoretical understanding of the winds and evolution of massive stars.
We obtained four pointings of over 100 ks each of the well-studied Wolf-Rayet star WR 6 with the XMM-Newton satellite. With a first paper emphasizing the results of spectral analysis, this follow-up highlights the X-ray variability clearly detected in all four pointings. However, phased light curves fail to confirm obvious cyclic behavior on the well-established 3.766 day period widely found at longer wavelengths. The data are of such quality that we were able to conduct a search for event clustering in the arrival times of X-ray photons. However, we fail to detect any such clustering. One possibility is that X-rays are generated in a stationary shock structure. In this context we favor a corotating interaction region (CIR) and present a phenomenological model for X-rays from a CIR structure. We show that a CIR has the potential to account simultaneously for the X-ray variability and constraints provided by the spectral analysis. Ultimately, the viability of the CIR model will require both intermittent long-term X-ray monitoring of WR 6 and better physical models of CIR X-ray production at large radii in stellar winds.
We have obtained spectropolarimetric observations of two Wolf-Rayet stars, WR 135 (WC8) and WR 136 (WN6), with the 6-m Russian telescope in July 2009 and July 2010. We have studied the He II 5412 angstrom line region, which contains also the C IV 5469 angstrom line (for WR 135 only). Our goals were to investigate the rapid line-profile variability (LPV) in WR star spectra and to search for magnetic fields. We find small amplitude emission peaks moving from the center of He II line to its wings during the night in spectra of both stars. These emission peaks are likely a signature of accelerating clumps in the stellar wind. We obtained upper limits of the magnetic field strength: approximate to 200G for WR 135 and approximate to 50G for WR 136. (C) 2011 WILEY-VCH Verlag GmbH&Co. KGaA, Weinheim
Many Central Stars of Planetary Nebulae are very similar to massive Wolf-Rayet stars of the carbon sequence with respect to their spectra, chemical composition and wind properties. Therefore their study opens an additional way towards the understanding of the Wolf-Rayet phenomenon. While the study of Line Profile Variation will be difficult, espescially for the very compact early types, the comparision with other hydrogen-deficient Central Stars illuminates the driving mechanism of their winds. We speculate that at least two ingredients are needed. The ionization of their atmpospheres has to be stratified to enable multi-scattering processes and the amount of carbon and oxygen has to be high (more than a few percent by mass).
Spectral analyses of central of planetary nebulae of early WC-type / NGC 6751 and Sanduleak 3
(1997)
We present our technique for solving the equations of radiation transfer in spherically expanding atmospheres. To ensure an efficient treatment of the Thomson scattering, the mean intensity J is derived by solving the moment equations in turn with the angle-dependent transfer equation. The latter provide the Eddington factors. Two different methods for the solution of the angle dependent equation are compared. Thereby the integration along short characteristics turned out to be superior in our context over the classical differencing scheme. The method is the basis of a non-LTE code suitable for the atmospheres of hot stars with high mass-loss.
We have obtained time-resolved observations of line-profile variations of the two Wolf-Rayet stars WR 135 and WR 111. The spectra, taken during two consecutive nights, cover a broad range from 4470 to 6590 Ang. The profile variability of the C iii emission line at 5696 Ang in WR 135 is shown in detail. The principal difficulties to constrain the velocity law from the frequency drift of discrete spectral features is discussed, emphasizing the crucial dependence on the adopted location of the line-emission region, and the possible necessity to distinguish between the motion of structures and the flow of the matter. - Full access to the observational data is provided via anonymous file transfer.
It has been observationally established that winds of hot massive stars have highly variable characteristics. The variability evident in the winds is believed to be caused by structures on a broad range of spatial scales. Small-scale structures (clumping) in stellar winds of hot stars are possible consequence of an instability appearing in their radiation hydrodynamics. To understand how clumping may influence calculation of theoretical spectra, different clumping properties and their 3D nature have to be taken into account. Properties of clumping have been examined using our 3D radiative transfer calculations. Effects of clumping for the case of the B[e] phenomenon are discussed.
Already for decades it has been known that the winds of massive stars are inhomogeneous (i.e. clumped). To properly model observed spectra of massive star winds it is necessary to incorporate the 3-D nature of clumping into radiative transfer calculations. In this paper we present our full 3-D Monte Carlo radiative transfer code for inhomogeneous expanding stellar winds. We use a set of parameters to describe dense as well as the rarefied wind components. At the same time, we account for non-monotonic velocity fields. We show how the 3-D density and velocity wind inhomogeneities strongly affect the resonance line formation. We also show how wind clumping can solve the discrepancy between P v and H alpha mass-loss rate diagnostics.
Context. Metal-poor massive stars are assumed to be progenitors of certain supernovae, gamma-ray bursts, and compact object mergers that might contribute to the early epochs of the Universe with their strong ionizing radiation. However, this assumption remains mainly theoretical because individual spectroscopic observations of such objects have rarely been carried out below the metallicity of the Small Magellanic Cloud. Aims. Here we explore the predictions of the state-of-the-art theories of stellar evolution combined with those of stellar atmospheres about a certain type of metal-poor (0.02 Z(circle dot)) hot massive stars, the chemically homogeneously evolving stars that we call Transparent Wind Ultraviolet INtense (TWUIN) stars. Methods. We computed synthetic spectra corresponding to a broad range in masses (20 130 M-circle dot) and covering several evolutionary phases from the zero-age main-sequence up to the core helium-burning stage. We investigated the influence of mass loss and wind clumping on spectral appearance and classified the spectra according to the Morgan-Keenan (MK) system. Results. We find that TWUIN stars show almost no emission lines during most of their core hydrogen-burning lifetimes. Most metal lines are completely absent, including nitrogen. During their core helium-burning stage, lines switch to emission, and even some metal lines (oxygen and carbon, but still almost no nitrogen) are detected. Mass loss and clumping play a significant role in line formation in later evolutionary phases, particularly during core helium-burning. Most of our spectra are classified as an early-O type giant or supergiant, and we find Wolf-Rayet stars of type WO in the core helium-burning phase. Conclusions. An extremely hot, early-O type star observed in a low-metallicity galaxy could be the result of chemically homogeneous evolution and might therefore be the progenitor of a long-duration gamma-ray burst or a type Ic supernova. TWUIN stars may play an important role in reionizing the Universe because they are hot without showing prominent emission lines during most of their lifetime.
Macroclumping in WR 136
(2015)
Macroclumping proved to resolve the discordance between different mass-loss rate diagnostics for O-type stars, in particular between Hα and the P v resonance lines. In this paper, we report first results from a corresponding investigation for WR stars. We apply our detailed 3-D Monte Carlo (MC) line formation code to the P v resonance doublet and show, for the Galactic WNL star WR136, that macroclumping is require to bring this line in accordance with the mass-loss rate derived from the emission-line spectrum.
Clumping in Galactic WN stars : a comparison of mass loss rates from UV/optical & radio diagnostics
(2007)
The mass loss rates and other parameters for a large sample of Galactic WN stars have been revised by Hamann et al. (2006), using the most up-to date Potsdam Wolf-Rayet (PoWR) model atmospheres. For a sub-sample of these stars exist measurements of their radio free-free emission. After harmonizing the adopted distance and terminal wind velocities, we compare the mass loss rates obtained from the two diagnostics. The differences are discussed as a possible consequence of different clumping contrast in the line-forming and radio-emitting regions.
The Quintuplet, one of three massive stellar clusters in the Galactic center (GC), is located about 30 pc in projection from Sagittarius A*. We aim at the construction of the Hertzsprung-Russell diagram (HRD) of the cluster to study its evolution and to constrain its star-formation history. For this purpose we use the most complete spectral catalog of the Quintuplet stars. Based on the K-band spectra we determine stellar temperatures and luminosities for all stars in the catalog under the assumption of a uniform reddening towards the cluster. We find two groups in the resulting HRD: early-type OB stars and late-type KM stars, well separated from each other. By comparison with Geneva stellar evolution models we derive initial masses exceeding 8 M-circle dot for the OB stars. In the HRD these stars are located along an isochrone corresponding to an age of about 4 Myr. This confirms previous considerations, where a similar age estimate was based on the presence of evolved Wolf-Rayet stars in the cluster. We derive number ratios for the various spectral subtype groups (e.g. N-WR/N-O, N-WC/N-WN) and compare them with predictions of population synthesis models. We find that an instantaneous burst of star formation at about 3.3 to 3.6 Myr ago is the most likely scenario to form the Quintuplet cluster. Furthermore, we apply a mass-luminosity relation to construct the initial mass function (IMF) of the cluster. We find indications for a slightly top-heavy IMF. The late-type stars in the LHO catalog are red giant branch (RGB) stars or red supergiants (RSGs) according to their spectral signatures. Under the assumption that they are located at about the distance of the Galactic center we can derive their luminosities. The comparison with stellar evolution models reveals that the initial masses of these stars are lower than 15 M-circle dot implying that they needed about 15 Myr (RSG) or even more than 30 Myr (RGB) to evolve into their present stage. It might be suspected that these late-type stars do not physically belong to the Quintuplet cluster. Indeed, most of them disqualify as cluster members because their radial velocities differ too much from the cluster average. Nevertheless, five of the brightest RGB/RSG stars from the LHO catalog share the mean radial velocity of the Quintuplet, and thus remain highly suspect for being gravitationally bound members. If so, this would challenge the cluster formation and evolution scenario.
Context: Three very massive clusters are known to reside in the Galactic center region, the Arches cluster, the Quintuplet cluster, and the central parsec cluster, each of them rich in young hot stars. With new infrared instruments, this region is no longer obscured for the observer. Aims: For understanding these very massive clusters, it is essential to know their stellar inventory. We provide comprehensive spectroscopic data for the stellar population of the Quintuplet cluster that will form the basis of subsequent spectral analyses. Methods. Spectroscopic observations of the Quintuplet cluster were obtained with the Integral Field Spectrograph SINFONI-SPIFFI at the ESO-VLT. The inner part of the Quintuplet cluster was covered by 22 slightly overlapping fields, each of them of 8 '' x 8 '' in size. The spectral range comprises the near-IR K-band from 1.94 to 2.45 mu m. The 3D data cubes of the individual fields were flux-calibrated and combined to one contiguous cube, from which the spectra of all detectable point sources were extracted. Results: We present a catalog of 160 stellar sources in the inner part of the Quintuplet cluster. The flux-calibrated K-band spectra of 98 early-type stars and 62 late-type stars are provided as Online Material. Based on these spectra, we assign spectral types to all detected sources and derive synthetic Ks-band magnitudes. Our sample is complete to about the 13th K-magnitude. We report the detection of two hitherto unknown Wolf- Rayet stars of late WC type (WC9 or later). Radial velocities are measured and employed to assess the cluster membership. The quantitative analysis of the early-type spectra will be the subject of a subsequent paper.
The quintuplet cluster III. Hertzsprung-Russell diagram and cluster age (vol 540, pg A14, 2012)
(2014)
High-mass X-ray binaries consist of a massive donor star and a compact object. While several of those systems have been well studied in X-rays, little is known for most of the donor stars as they are often heavily obscured in the optical and ultraviolet regime. There is an opportunity to observe them at infrared wavelengths, however. The goal of this study is to obtain the stellar and wind parameters of the donor star in the X1908+075 high-mass X-ray binary system with a stellar atmosphere model to check whether previous studies from X-ray observations and spectral morphology lead to a sufficient description of the donor star. We obtained H-and K-band spectra of X1908+075 and analysed them with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. For the first time, we calculated a stellar atmosphere model for the donor star, whose main parameters are: M-spec = 15 +/- 6 M-circle dot, T-* = 23(-3)(+6) kK, log g(eff) = 3.0 +/- 0.2 and log L/L-circle dot = 4.81 +/- 0.25. The obtained parameters point towards an early B-type (B0-B3) star, probably in a supergiant phase. Moreover we determined a more accurate distance to the system of 4.85 +/- 0.50 kpc than the previously reported value.
The development of infrared observational facilities has revealed a number of massive stars in obscured environments throughout the Milky Way and beyond. The determination of their stellar and wind properties from infrared diagnostics is thus required to take full advantage of the wealth of observations available in the near and mid infrared. However, the task is challenging. This session addressed some of the problems encountered and showed the limitations and successes of infrared studies of massive stars.
General Discussion
(2007)
The B fields in OB stars (BOB) survey is an ESO large programme collecting spectropolarimetric observations for a large number of early-type stars in order to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. As of July 2014, a total of 98 objects were observed over 20 nights with FORS2 and HARPSpol. Our preliminary results indicate that the fraction of magnetic OB stars with an organised, detectable field is low. This conclusion, now independently reached by two different surveys, has profound implications for any theoretical model attempting to explain the field formation in these objects. We discuss in this contribution some important issues addressed by our observations (e.g., the lower bound of the field strength) and the discovery of some remarkable objects.
We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution spectral data set of the delta Ori Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of approximate to 479 ks and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range of 5-25 is is confirmed, with a maximum amplitude of about +/- 15% within a single approximate to 125 ks observation. Periods of 4.76 and 2.04 days are found in the total X-ray flux, as well as an apparent overall increase in the flux level throughout the nine-day observational campaign. Using 40 ks contiguous spectra derived from the original observations, we investigate the variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at phi = 0.0 when the secondary delta Ori Aa2 is at the inferior conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission line width variability to the presence of a wind cavity created by a wind-wind collision, which is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary wind, thus producing phase-locked X-ray variability.
It is commonly adopted that X-rays from O stars are produced deep inside the stellar wind, and transported outwards through the bulk of the expanding matter which attenuates the radiation and affects the shape of emission line profiles. The ability of the X-ray observatories Chandra and XMM-Newton to resolve these lines spectroscopically provided a stringent test for the theory of the X-ray production. It turned out that none of the existing models was able to fit the observations consistently. The possible caveat of these models was the underlying assumption of a smooth stellar wind. Motivated by the evidence that the stellar winds are in fact structured, we present a 2-D numerical model of a stochastic, inhomogeneous wind. Small parcels of hot, X-ray emitting gas are permeated by cool, absorbing wind material which is compressed into thin shell fragments. Wind fragmentation alters the radiative transfer drastically, compared to homogeneous models of the same mass-loss rate. X-rays produced deep inside the wind, which would be totally absorbed in a homogeneous flow, can effectively escape from a fragmented wind. The wind absorption becomes wavelength independent if the individual fragments are optically thick. The X-ray line profiles are flat-topped in the blue part and decline steeply in the red part for the winds with a short acceleration zone. For the winds where the acceleration extends over significant distances, the lines can appear nearly symmetric and only slightly blueshifted, in contrast to the skewed, triangular line profiles typically obtained from homogeneous wind models of high optical depth. We show that profiles from a fragmented wind model can reproduce the observed line profiles from zeta Orionis. The present numerical modeling confirms the results from a previous study, where we derived analytical formulae from a statistical treatment
Archival X-ray spectra of the four prominent single, non-magnetic O stars zeta Pup, zeta Ori, xi Per and zeta Oph, obtained in high resolution with Chandra HETGS/MEG have been studied. The resolved X-ray emission line profiles provide information about the shocked, hot gas which emits the X-radiation, and about the bulk of comparably cool stellar wind material which partly absorbs this radiation. In this paper, we synthesize X-ray line profiles with a model of a clumpy stellar wind. We find that the geometrical shape of the wind inhomogeneities is important: better agreement with the observations can be achieved with radially compressed clumps than with spherical clumps. The parameters of the model, i.e. chemical abundances, stellar radius, mass-loss rate and terminal wind velocity, are taken from existing analyses of UV and optical spectra of the programme stars. On this basis, we also calculate the continuum-absorption coefficient of the cool-wind material, using the Potsdam Wolf-Rayet (POWR) model atmosphere code. The radial location of X-ray emitting gas is restricted from analysing the FIR line ratios of helium-like ions. The only remaining free parameter of our model is the typical distance between the clumps; here, we assume that at any point in the wind there is one clump passing by per one dynamical time-scale of the wind. The total emission in a model line is scaled to the observation. There is a good agreement between synthetic and observed line profiles. We conclude that the X-ray emission line profiles in O stars can be explained by hot plasma embedded in a cool wind which is highly clumped in the form of radially compressed shell fragments.
We present the first high-resolutionX-ray spectrum of a putatively singleWolf-Rayet (WR) star. 400 ks observations of WR 6 by the XMM-Newton telescope resulted in a superb quality high-resolution X-ray spectrum. Spectral analysis reveals that the X-rays originate far out in the stellar wind, more than 30 stellar radii from the photosphere, and thus outside the wind acceleration zone where the line-driving instability (LDI) could create shocks. The X-ray emitting plasma reaches temperatures up to 50 MK and is embedded within the unshocked, "cool" stellar wind as revealed by characteristic spectral signatures. We detect a fluorescent Fe line at approximate to 6.4 keV. The presence of fluorescence is consistent with a two-component medium, where the cool wind is permeated with the hot X-ray emitting plasma. The wind must have a very porous structure to allow the observed amount of X-rays to escape. We find that neither the LDI nor any alternative binary scenario can explain the data. We suggest a scenario where X-rays are produced when the fast wind rams into slow "sticky clumps" that resist acceleration. Our new data show that the X-rays in single WR star are generated by some special mechanism different from the one operating in the O-star winds.
We investigate the connections between the magnetic fields and the X-ray emission from massive stars. Our study shows that the X-ray properties of known strongly magnetic stars are diverse: while some comply to the predictions of the magnetically confined wind model, others do not. We conclude that strong, hard, and variable X-ray emission may be a sufficient attribute of magnetic massive stars, but it is not a necessary one. We address the general properties of X-ray emission from "normal" massive stars, especially the long standing mystery about the correlations between the parameters of X-ray emission and fundamental stellar properties. The recent development in stellar structure modeling shows that small-scale surface magnetic fields may be common. We suggest a "hybrid" scenario which could explain the X-ray emission from massive stars by a combination of magnetic mechanisms on the surface and shocks in the stellar wind. The magnetic mechanisms and the wind shocks are triggered by convective motions in sub-photospheric layers. This scenario opens the door for a natural explanation of the well established correlation between bolometric and X-ray luminosities.