A small fraction of the radiative flux emitted by hot stars is absorbed by their winds and redistributed towards longer wavelengths. This effect, which leads also to the heating of the stellar ...photosphere, is termed wind blanketing. For stars with variable winds, the effect of wind blanketing may lead to the photometric variability. We have studied the consequences of line driven wind instability and wind blanketing for the light variability of O stars. We combined the results of wind hydrodynamic simulations and of global wind models to predict the light variability of hot stars due to the wind blanketing and instability. The wind instability causes stochastic light variability with amplitude of the order of tens of millimagnitudes and a typical timescale of the order of hours for spatially coherent wind structure. The amplitude is of the order of millimagnitudes when assuming that the wind consists of large number of independent concentric cones. The variability with such amplitude is observable using present space borne photometers. We show that the simulated light curve is similar to the light curves of O stars obtained using BRITE and CoRoT satellites.
Massive stars lose a significant fraction of mass during their evolution. However, the corresponding mass-loss rates are rather uncertain, especially for evolved stars. To improve this, we calculated ...global line-driven wind models for Galactic B supergiants. Our models predict radial wind structure and particularly the mass-loss rates and terminal velocities directly from basic stellar parameters. The hydrodynamic structure of the flow is consistently determined from the photosphere in nearly hydrostatic equilibrium to supersonically expanding wind. The radiative force is derived from the solution of the radiative transfer equation in the comoving frame. We provide a simple formula that predicts theoretical mass-loss rates as a function of stellar luminosity and effective temperature. The mass-loss rate of B supergiants slightly decreases with temperature down to about 22.5 kK, where the region of recombination of Fe
IV
to Fe
III
starts to appear. In this region, which is about 5 kK wide, the mass-loss rate gradually increases by a factor of about 6. The increase of the mass-loss rate is associated with a gradual decrease of terminal velocities by a factor of about 2. We compared the predicted wind parameters with observations. While the observed wind terminal velocities are reasonably reproduced by the models, the situation with mass-loss rates is less clear. The mass-loss rates derived from observations that are uncorrected for clumping are by a factor of 3 to 9 higher than our predictions on cool and hot sides of the studied sample, respectively. These observations can be reconciled with theory assuming a temperature-dependent clumping factor that is decreasing toward lower effective temperatures. On the other hand, the mass-loss rate estimates that are not sensitive to clumping agree with our predictions much better. Our predictions are by a factor of about 10 lower than the values currently used in evolutionary models appealing for reconsideration of the role of winds in the stellar evolution.
Fraction of hot stars posses strong magnetic fields that channel their radiatively driven outflows. We study the influence of line splitting in the magnetic field (Zeeman effect) on the wind ...properties. We use our own global wind code with radiative transfer in the comoving frame to understand the influence of the Zeeman splitting on the line force. We show that the Zeeman splitting has a negligible influence on the line force for magnetic fields that are weaker than about 100 kG. This means that the wind mass-loss rates and terminal velocities are not affected by the magnetic line splitting for magnetic fields as are typically found on the surface of nondegenerate stars. Neither have we found any strong flux variability that would be due to the magnetically split line blanketing.
Some supernova (SN) explosions show evidence for an interaction with a pre-existing nonspherically symmetric circumstellar medium (CSM) in their light curves, spectral line profiles, and polarization ...signatures. The origin of this aspherical CSM is unknown, but binary interactions have often been implicated. To better understand the connection with binary stars and to aid in the interpretation of observations, we performed two-dimensional axisymmetric hydrodynamic simulations where an expanding spherical SN ejecta initialized with realistic density and velocity profiles collide with various aspherical CSM distributions. We consider CSM in the form of a circumstellar disk, colliding wind shells in binary stars with different orientations and distances from the SN progenitor, and bipolar lobes representing a scaled down version of the Homunculus nebula of
η
Car. We study how our simulations map onto observables, including approximate light curves, indicative spectral line profiles at late times, and estimates of a polarization signature. We find that the SN–CSM collision layer is composed of normal and oblique shocks, reflected waves, and other hydrodynamical phenomena that lead to acceleration and shear instabilities. As a result, the total shock heating power fluctuates in time, although the emerging light curve might be smooth if the shock interaction region is deeply embedded in the SN envelope. SNe with circumstellar disks or bipolar lobes exhibit late-time spectral line profiles that are symmetric with respect to the rest velocity and relatively high polarization. In contrast, SNe with colliding wind shells naturally lead to line profiles with asymmetric and time-evolving blue and red wings and low polarization. Given the high frequency of binaries among massive stars, the interaction of SN ejecta with a pre-existing colliding wind shell must occur and the observed signatures could be used to characterize the binary companion.
Context. In wind-powered X-ray binaries, the radiatively driven stellar wind from the primary may be inhibited by the X-ray irradiation. This creates the feedback that limits the X-ray luminosity of ...the compact secondary. Wind inhibition might be weakened by the effect of small-scale wind inhomogeneities (clumping) possibly affecting the limiting X-ray luminosity. Aims. We study the influence of X-ray irradiation on the stellar wind for different radial distributions of clumping. Methods. We calculate hot star wind models with external irradiation and clumping using our global wind code. The models are calculated for different parameters of the binary. We determine the parameters for which the X-ray wind ionization is so strong that it leads to a decrease of the radiative force. This causes a decrease of the wind velocity and even of the mass-loss rate in the case of extreme X-ray irradiation. Results. Clumping weakens the effect of X-ray irradiation because it favours recombination and leads to an increase of the wind mass-loss rate. The best match between the models and observed properties of high-mass X-ray binaries (HMXBs) is derived with radially variable clumping. We describe the influence of X-ray irradiation on the terminal velocity and on the mass-loss rate in a parametric way. The X-ray luminosities predicted within the Bondi-Hoyle-Lyttleton theory agree nicely with observations when accounting for X-ray irradiation. Conclusions. The ionizing feedback regulates the accretion onto the compact companion resulting in a relatively stable X-ray source. The wind-powered accretion model can account for large luminosities in HMXBs only when introducing the ionizing feedback. There are two possible states following from the dependence of X-ray luminosity on the wind terminal velocity and mass-loss rate. One state has low X-ray luminosity and a nearly undisturbed wind, and the second state has high X-ray luminosity and exhibits a strong influence of X-rays on the flow.
Context.
Fast rotating Wolf-Rayet stars are expected to be progenitors of long duration gamma-ray bursts. However, the observational test of this model is problematic. Spectral lines of Wolf-Rayet ...stars originate in expanding stellar wind, therefore a reliable spectroscopical determination of their rotational velocities is difficult. Intrinsic polarization of Wolf-Rayet stars due to the rotational modulation of the stellar wind may provide an indirect way to determine the rotational velocities of these stars. However, detailed wind models are required for this purpose.
Aims.
We determine the intrinsic polarization of Wolf-Rayet stars from hydrodynamical wind models as a function of rotational velocity.
Methods.
We used 2.5D hydrodynamical simulations to calculate the structure of rotating winds of Wolf-Rayet stars. The simulations account for the deformation of the stellar surface due to rotation, gravity darkening, and nonradial forces. From the derived models, we calculated the intrinsic stellar polarization. The mass loss rate was scaled to take realistic wind densities of Wolf-Rayet stars into account.
Results.
The hydrodynamical wind models predict a prolate wind structure, which leads to a relatively low level of polarization. Even relatively large rotational velocities are allowed by observational constrains. The obtained wind structure is similar to that obtained previously for rotating optically thin winds.
Conclusions.
Derived upper limits of rotational velocities of studied Wolf-Rayet stars are not in conflict with the model of long duration gamma-ray bursts.
Context. Massive stars are characterized by a significant loss of mass either via (nearly) spherically symmetric stellar winds or pre-explosion pulses, or by aspherical forms of circumstellar matter ...(CSM) such as bipolar lobes or outflowing circumstellar equatorial disks. Since a significant fraction of most massive stars end their lives by a core collapse, supernovae (SNe) are always located inside large circumstellar envelopes created by their progenitors. Aims. We study the dynamics and thermal effects of collision between expanding ejecta of SNe and CSM that may be formed during, for example, a sgBe star phase, a luminous blue variable phase, around PopIII stars, or by various forms of accretion. Methods. For time-dependent hydrodynamic modeling we used our own grid-based Eulerian multidimensional hydrodynamic code built with a finite volumes method. The code is based on a directionally unsplit Roe’s method that is highly efficient for calculations of shocks and physical flows with large discontinuities. Results. We simulate a SNe explosion as a spherically symmetric blast wave. The initial geometry of the disks corresponds to a density structure of a material that orbits in Keplerian trajectories. We examine the behavior of basic hydrodynamic characteristics, i.e., the density, pressure, velocity of expansion, and temperature structure in the interaction zone under various geometrical configurations and various initial densities of CSM. We calculate the evolution of the SN–CSM system and the rate of aspherical deceleration as well as the degree of anisotropy in density, pressure, and temperature distribution. Conclusions. Our simulations reveal significant asphericity of the expanding envelope above all in the case of dense equatorial disks. Our “low density” model however also shows significant asphericity in the case of the disk mass-loss rate Ṁcsd = 10−6 M⊙ yr−1. The models also show the zones of overdensity in the SN–disk contact region and indicate the development of Kelvin-Helmholtz instabilities within the zones of shear between the disk and the more freely expanding material outside the disk.
We provide mass-loss rate predictions for O stars from Large and Small Magellanic Clouds. We calculate global (unified, hydrodynamic) model atmospheres of main sequence, giant, and supergiant stars ...for chemical composition corresponding to Magellanic Clouds. The models solve radiative transfer equation in comoving frame, kinetic equilibrium equations (also known as NLTE equations), and hydrodynamical equations from (quasi-)hydrostatic atmosphere to expanding stellar wind. The models allow us to predict wind density, velocity, and temperature (consequently also the terminal wind velocity and the mass-loss rate) just from basic global stellar parameters. As a result of their lower metallicity, the line radiative driving is weaker leading to lower wind mass-loss rates with respect to the Galactic stars. We provide a formula that fits the mass-loss rate predicted by our models as a function of stellar luminosity and metallicity. On average, the mass-loss rate scales with metallicity as Ṁ ~ Z0.59. The predicted mass-loss rates are lower than mass-loss rates derived from Hα diagnostics and can be reconciled with observational results assuming clumping factor Cc = 9. On the other hand, the predicted mass-loss rates either agree or are slightly higher than the mass-loss rates derived from ultraviolet wind line profiles. The calculated P v ionization fractions also agree with values derived from observations for LMC stars with Teff ≤ 40 000 K. Taken together, our theoretical predictions provide reasonable models with consistent mass-loss rate determination, which can be used for quantitative study of stars from Magellanic Clouds.
Wind-fed high-mass X-ray binaries are powered by accretion of the radiatively driven wind of the luminous component on the compact star. Accretion-generated X-rays alter the ionization state of the ...wind. Because higher ionization states drive the wind less effectively, X-ray ionization may brake acceleration of the wind. This causes a decrease in the wind terminal velocity and mass flux in the direction toward the X-ray source. Here we study the effect of X-ray ionization on the stellar wind of B supergiants. We determine the binary parameters for which the X-ray irradiation significantly influences the stellar wind. This can be conveniently studied in diagrams that plot the optical depth parameter versus the X-ray luminosity. For low optical depths or for high X-ray luminosities, X-ray ionization leads to a disruption in the wind aimed toward the X-ray source. Observational parameters of high-mass X-ray binaries with B-supergiant components appear outside the wind disruption zone. The X-ray feedback determines the resulting X-ray luminosity. We recognize two states with a different level of feedback. For low X-ray luminosities, ionization is weak, and the wind is not disrupted by X-rays and flows at large velocities, consequently the accretion rate is relatively low. On the other hand, for high X-ray luminosities, the X-ray ionization disrupts the flow braking the acceleration, the wind velocity is low, and the accretion rate becomes high. These effects determine the X-ray luminosity of individual binaries. Accounting for the X-ray feedback, estimated X-ray luminosities reasonably agree with observational values. We study the effect of small-scale wind inhomogeneities (clumping), showing that clumping weakens the effect of X-ray ionization by increasing recombination and the mass-loss rate. This effect is particularly important in the region of the so-called bistability jump. We show that ultraluminous X-ray binaries with
L
X
≲ 10
40
erg s
−1
may be powered by accretion of a B-supergiant wind on a massive black hole.
Context.
Fast line-driven stellar winds play an important role in the evolution of planetary nebulae, even though they are relatively weak.
Aims.
We provide global (unified) hot star wind models of ...central stars of planetary nebulae. The models predict wind structure including the mass-loss rates, terminal velocities, and emergent fluxes from basic stellar parameters.
Methods.
We applied our wind code for parameters corresponding to evolutionary stages between the asymptotic giant branch and white dwarf phases for a star with a final mass of 0.569
M
⊙
. We study the influence of metallicity and wind inhomogeneities (clumping) on the wind properties.
Results.
Line-driven winds appear very early after the star leaves the asymptotic giant branch (at the latest for
T
eff
≈ 10 kK) and fade away at the white dwarf cooling track (below
T
eff
= 105 kK). Their mass-loss rate mostly scales with the stellar luminosity and, consequently, the mass-loss rate only varies slightly during the transition from the red to the blue part of the Hertzsprung–Russell diagram. There are the following two exceptions to the monotonic behavior: a bistability jump at around 20 kK, where the mass-loss rate decreases by a factor of a few (during evolution) due to a change in iron ionization, and an additional maximum at about
T
eff
= 40−50 kK. On the other hand, the terminal velocity increases from about a few hundreds of km s
−1
to a few thousands of km s
−1
during the transition as a result of stellar radius decrease. The wind terminal velocity also significantly increases at the bistability jump. Derived wind parameters reasonably agree with observations. The effect of clumping is stronger at the hot side of the bistability jump than at the cool side.
Conclusions.
Derived fits to wind parameters can be used in evolutionary models and in studies of planetary nebula formation. A predicted bistability jump in mass-loss rates can cause the appearance of an additional shell of planetary nebula.