Inhomogeneous Stellar Wind Research Articles

Analysis of X-Ray Emission from OB Stars. IV. About X-Ray Emission from Inhomogeneous Winds of OB Stars

We study the origin of X-ray emission from OB stars due to collisions of stellar winds and/or inhomogeneities in the winds. The low-resolution X-ray spectra of a big sample of OB stars were fitted by both the stationary APEC/MEKAL models and by this model with an additional PSHOCK component describing the nonstationary X-ray emission. These spectra were also described by two-temperature PSHOCK models. More than ∼50% of considered spectra can be described by the above-mentioned model combinations including the PSHOCK model and the quality of the fits appears to be better for O stars. The plasma temperature of the PSHOCK component is about 1–5 keV with the ionization timescale τ u ∼ 108–1013 s cm−3. The temperature of the PSHOCK component increases with the momentum and kinetic energy of the stellar wind by a power law with an index ∼0.12–0.14. Such dependencies were not revealed through modeling by the stationary APEC/MEKAL models only. At the same time the X-ray luminosity of OB stars depends on momentum and kinetic energy of their winds similarly either for stationary or for nonstationary models. We conclude that many O stars and some B stars can be sources of the nonstationary X-rays formed in their inhomogeneous stellar wind.

Modeling the effects of clumpy winds in the high-energy light curves of γ-ray binaries

Context. High-mass gamma-ray binaries are powerful nonthermal galactic sources, some of them hosting a pulsar whose relativistic wind interacts with a likely inhomogeneous stellar wind. So far, modeling these sources including stellar wind inhomogeneities has been done using either simple analytical approaches or heavy numerical simulations, none of which allow for an exploration of the parameter space that is both reasonably realistic and general. Aims. Applying different semi-analytical tools together, we study the dynamics and high-energy radiation of a pulsar wind colliding with a stellar wind with different degrees of inhomogeneity to assess the related observable effects. Methods. We computed the arrival of clumps to the pulsar wind-stellar wind interaction structure using a Monte Carlo method and a phenomenological clumpy-wind model. The dynamics of the clumps that reach deep into the pulsar wind zone was computed using a semi-analytical approach. This approach allows for the characterization of the evolution of the shocked pulsar wind region in times much shorter than the orbital period. With this three-dimensional information about the emitter, we applied analytical adiabatic and radiative models to compute the variable high-energy emission produced on binary scales. Results. An inhomogeneous stellar wind induces stochastic hour-timescale variations in the geometry of the two-wind interaction structure on binary scales. Depending on the degree of stellar wind inhomogeneity, 10–100% level hour-scale variability in the X-rays and gamma rays is predicted, with the largest variations occurring roughly once per orbit. Conclusions. Our results, based on a comprehensive approach, show that present X-ray and future very-high-energy instrumentation can allow us to trace the impact of a clumpy stellar wind on the shocked pulsar wind emission in a gamma-ray binary.

Towards a Unified View of Inhomogeneous Stellar Winds in Isolated Supergiant Stars and Supergiant High Mass X-Ray Binaries

Massive stars, at least $\sim$ 10 times more massive than the Sun, have two key properties that make them the main drivers of evolution of star clusters, galaxies, and the Universe as a whole. On the one hand, the outer layers of massive stars are so hot that they produce most of the ionizing ultraviolet radiation of galaxies; in fact, the first massive stars helped to re-ionize the Universe after its Dark Ages. Another important property of massive stars are the strong stellar winds and outflows they produce. This mass loss, and finally the explosion of a massive star as a supernova or a gamma-ray burst, provide a significant input of mechanical and radiative energy into the interstellar space. These two properties together make massive stars one of the most important cosmic engines: they trigger the star formation and enrich the interstellar medium with heavy elements, that ultimately leads to formation of Earth-like rocky planets and the development of complex life. The study of massive star winds is thus a truly multidisciplinary field and has a wide impact on different areas of astronomy. [...] This detailed review summarises the current knowledge on the theory and observations of winds from massive stars, as well as on observations and accretion processes in wind-fed high mass X-ray binaries. The aim is to combine in the near future all available theoretical diagnostics and observational measurements to achieve a unified picture of massive star winds in isolated objects and in binary systems.

Non-thermal radiation from a pulsar wind interacting with an inhomogeneous stellar wind

Binaries hosting a massive star and a non-accreting pulsar are powerful non-thermal emitters due to the interaction of the pulsar and the stellar wind. The winds of massive stars are thought to be inhomogeneous, which could have an impact on the non-thermal emission. We study numerically the impact of the presence of inhomogeneities or clumps in the stellar wind on the high-energy non-thermal radiation of high-mass binaries hosting a non-accreting pulsar. We compute the trajectories and physical properties of the streamlines in the shocked pulsar wind without clumps, with a small clump, and with a large one. This information is used to compute the synchrotron and inverse Compton emission from the non-thermal populations, accounting also for the effect of gamma-ray absorption through pair creation. A specific study is done for PSR B1259-63/LS2883. When stellar wind clumps perturb the two-wind interaction region, the associated non-thermal radiation in the X-ray band,of synchrotron origin, and in the GeV-TeV band, of inverse Compton origin, is affected by several effects: (i) strong changes in the the plasma velocity direction that result in Doppler boosting factor variations, (ii) strengthening of the magnetic field that mainly enhances the synchrotron radiation, (iii) strengthening of the pulsar wind kinetic energy dissipation at the shock, potentially available for particle acceleration, and (iv) changes in the rate of adiabatic losses that affect the lower energy part of the non-thermal particle population. The radiation above 100 GeV detected, presumably, during the post-periastron crossing of the Be star disc in PSR B1259-63/LS2883, can be roughly reproduced assuming that the crossing of the disc is modeled as the encounter with a large inhomogeneity.

Simulations of an inhomogeneous stellar wind interacting with a pulsar wind in a binary system

Binary systems containing a massive star and a non-accreting pulsar present strong interaction between the stellar and the pulsar winds. The properties of this interaction, which largely determine the non-thermal radiation in these systems, strongly depend on the structure of the stellar wind, which can be clumpy or strongly anisotropic, as in Be stars. We study numerically the influence of inhomogeneities in the stellar wind on the structure of the two-wind interaction region. We carried out for the first time axisymmetric, relativistic hydrodynamical simulations, with Lorentz factors of ~6 and accounting for the impact of instabilities, to study the impact in the two-wind interaction structure of an over-dense region of the stellar wind. We also followed the evolution of this over-dense region or clump as it faces the impact of the pulsar wind. For typical system parameters, and adopting a stellar wind inhomogeneity with a density contrast >~10, clumps with radii of a few percent of the binary size can significantly perturb the two-wind interaction region, pushing the two-wind interface to <~40% of the initial distance to the pulsar. After it is shocked, the inhomogeneity quickly expands and is disrupted when it reaches the smallest distance to the pulsar. It eventually fragments, being advected away from the binary system. The whole interaction region is quite unstable, and the shocked pulsar wind can strongly change under small perturbations. We confirm the sensitive nature of the two-wind interaction structure to perturbations, in particular when the stellar wind is inhomogeneous. For realistic over-dense regions of the stellar wind, the interaction region may shrink by a factor of a few, with the shocked flow presenting a complex spatial and temporal pattern. This can lead to strong variations in the non-thermal radiation.

Modulated gamma-ray emission from compact millisecond pulsar binary systems

A significant amount of the millisecond pulsars has been discovered within binary systems. In several such binary systems the masses of the companion stars have been derived allowing to distinguish two classes of objects, called the Black Widow and the Redback binaries. Pulsars in these binary systems are expected to produce winds which, colliding with stellar winds, create conditions for acceleration of electrons. These electrons should interact with the anisotropic radiation from the companion stars producing gamma-ray emission modulated with the orbital period of the binary system. We consider the interaction of a millisecond pulsar (MSP) wind with a very inhomogeneous stellar wind from the companion star within binary systems of the Black Widow and Redback types. It is expected that the pulsar wind should mix efficiently with the inhomogeneous stellar wind. Electrons accelerated in such mixed, turbulent winds can interact with the magnetic field and also strong radiation from the companion star producing not only synchrotron radiation but also gamma-rays in the the Inverse Compton process. Applying numerical methods, we calculated the GeV-TeV gamma-ray spectra and the light curves expected from some millisecond pulsar binary systems. It is concluded that energetic millisecond pulsar binary systems create a new class of TeV gamma-ray sources which could be detectable by the future Cherenkov arrays (e.g. CTA) and possibly also by the extensive campaigns with the present arrays (HESS, MAGIC, VERITAS). However, gamma-ray emission from the millisecond pulsar binary systems is predicted to have different features than those observed in the case of massive TeV gamma-ray binaries such as LS I 303 61 or LS 5039. The maximum in the TeV gamma-ray orbital light curve should appear when the MSP is behind the companion star.

Clumped stellar winds in supergiant high-mass X-ray binaries: X-ray variability and photoionization

The clumping of massive star winds is an established paradigm, which is confirmed by multiple lines of evidence and is supported by stellar wind theory. The purpose of this paper is to bridge the gap between detailed models of inhomogeneous stellar winds in single stars and the phenomenological description of donor winds in supergiant high‐mass X‐ray binaries (HMXBs). We use the results from time‐dependent hydrodynamical models of the instability in the line‐driven wind of a massive supergiant star to derive the time‐dependent accretion rate on to a compact object in the Bondi–Hoyle–Lyttleton approximation. The strong density and velocity fluctuations in the wind result in strong variability of the synthetic X‐ray light curves. The model predicts a large‐scale X‐ray variability, up to eight orders of magnitude, on relatively short time‐scales. The apparent lack of evidence for such strong variability in the observed HMXBs indicates that the details of the accretion process act to reduce the variability resulting from the stellar wind velocity and density jumps. We study the absorption of X‐rays in the clumped stellar wind by means of a two‐dimensional stochastic wind model. The monochromatic absorption in the cool stellar wind, depending on the orbital phase, is computed for realistic stellar wind opacity. We find that the absorption of X‐rays changes strongly at different orbital phases. The degree of the variability resulting from the absorption in the wind depends on the shape of the wind clumps, and this is stronger for oblate clumps. We address the photoionization in the clumped wind, and we show that the degree of ionization is affected by the wind clumping. We derive a correction factor for the photoionization parameter, and we show that the photoionization parameter is reduced by a factor compared to the smooth wind models with the same mass‐loss rate, where is the wind inhomogeneity parameter. We conclude that wind clumping must also be taken into account when comparing the observed and model spectra of the photoionized stellar wind.

X-ray emission lines from inhomogeneous stellar winds

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. None of the existing models was able to fit the observationed lines consistently. The possible caveat of these models was the assumption of a smooth stellar wind. Motivated by the various evidences 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 with 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 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.

Similarity Solution for a Supernova Remnant Evolving in an Inhomogeneous Stellar Wind

Similarity Solution for a Supernova Remnant Evolving in an Inhomogeneous Stellar Wind

Profile variability of Hα in the B3e Herbig star HD 200775 as indication of matter infall and inhomogeneous stellar wind

Significant profile variability of Hα on the timescale of months is investigated on the basis of high-dispersion CCD spectra of the B3e Herbig star HD 200775. Variations in position and shape of the central absorption component are considered to indicate complex structure of circumstellar environment, containing infalling streams and inhomogeneous stellar wind.