No Arabic abstract
Bow-shaped mid-infrared emission regions have been discovered in satellite observations of numerous late-type O and early-type B stars with moderate velocities relative to the ambient interstellar medium. Previously, hydrodynamical bow shock models have been used to study this emission. It appears that such models are incomplete in that they neglect kinetic effects associated with long mean free paths of stellar wind particles, and the complexity of Weibel instability fronts. Wind ions are scattered in the Weibel instability and mix with the interstellar gas. However, they do not lose their momentum and most ultimately diffuse further into the ambient gas like cosmic rays, and share their energy and momentum. Lacking other coolants, the heated gas transfers energy to interstellar dust grains, which radiate it. This process, in addition to grain photo-heating, provides the energy for the emission. A weak R-type ionization front, formed well outside the infrared emission region, generally moderates the interstellar gas flow into the emission region. The theory suggests that the infrared emission process is limited to cases of moderate stellar peculiar velocities, evidently in accord with the observations.
Bow shocks and related density enhancements produced by the winds of massive stars moving through the interstellar medium provide important information regarding the motions of the stars, the properties of their stellar winds, and the characteristics of the local medium. Since bow shocks are aspherical structures, light scattering within them produces a net polarization signal even if the region is spatially unresolved. Scattering opacity arising from free electrons and dust leads to a distribution of polarized intensity across the bow shock structure. That polarization encodes information about the shape, composition, opacity, density, and ionisation state of the material within the structure. In this paper we use the Monte Carlo radiative transfer code SLIP to investigate the polarization created when photons scatter in a bow shock-shaped region of enhanced density surrounding a stellar source. We present results assuming electron scattering, and investigate the polarization behaviour as a function of optical depth, temperature, and source of photons for two different cases: pure scattering and scattering with absorption. In both regimes we consider resolved and unresolved cases. We discuss the implication of these results as well as their possible use along with observational data to constrain the properties of observed bow shock systems. In different situations and under certain assumptions, our simulations can constrain viewing angle, optical depth and temperature of the scattering region, and the relative luminosities of the star and shock.
We study the polarization produced by scattering from dust in a bow shock-shaped region of enhanced density surrounding a stellar source, using the Monte Carlo radiative transfer code SLIP. Bow shocks are structures formed by the interaction of the winds of fast-moving stars with the interstellar medium. Our previous study focused on the polarization produced in these structures by electron scattering; we showed that polarization is highly dependent on inclination angle and that multiple scattering changes the shape and degree of polarization. In contrast to electron scattering, dust scattering is wavelength-dependent, which changes the polarization behaviour. Here we explore different dust particle sizes and compositions and generate polarized spectral energy distributions for each case. We find that the polarization SED behaviour depends on the dust composition and grain size. Including dust emission leads to polarization changes with temperature at higher optical depth in ways that are sensitive to the orientation of the bow shock. In various scenarios and under certain assumptions, our simulations can constrain the optical depth and dust properties of resolved and unresolved bow shock-shaped scattering regions.Constraints on optical depth can provide estimates of local ISM density for observed bow shocks. We also study the impact of dust grains filling the region between the star and bow shock. We see that as the density of dust between the star and bow shock increases, the resulting polarization is suppressed for all the optical depth regimes.
Mass-loss rates and terminal wind velocities are key parameters that determine the kinetic wind energy and momenta of massive stars. Furthermore, accurate mass-loss rates determine the mass and rotational velocity evolution of mass stars, and their fates as neutron stars and black holes in function of metallicity (Z). Here we update our Monte Carlo mass-loss Recipe with new dynamically-consistent computations of the terminal wind velocity -- as a function of Z. These predictions are particularly timely as the HST ULLYSES project will observe ultraviolet spectra with blue-shifted P Cygni lines of hundreds of massive stars in the low-Z Large and Small Magellanic Clouds, as well as sub-SMC metallicity hosts. Around 35 000 K, we uncover a weak-wind dip and we present diagnostics to investigate its physics with ULLYSES and X-Shooter data. We discuss how the dip may provide important information on wind-driving physics, and how this is of key relevance towards finding a new gold-standard for OB star mass-loss rates. For B supergiants below the Fe IV to III bi-stability jump, the terminal velocity is found to be independent of Z and M, while the mass-loss rate still varies as $dot{M} propto Z^{0.85}$. For O-type stars above the bi-stability jump we find a terminal-velocity dependence of $v_{infty} propto Z^{0.19}$ and the Z-dependence of the mass-loss rate is found to be as shallow as $dot{M} propto Z^{0.42}$, implying that to reproduce the `heavy black holes from LIGO/VIRGO, the `low Z requirement becomes even more stringent than was previously anticipated.
By considering the advection and interaction of the vector momentum flux in highly supersonic spherically diverging winds, we derive a simple analytic description of the asymptotic opening angle of a wind-collision shock cone, in the approximation that the shocked gas is contained in a cone streaming out along a single characteristic opening angle. Both highly radiative and highly adiabatic limits are treated, and their comparison is the novel result. Analytic closed-form expressions are obtained for the inferred wind momentum ratios as a function of the observed shock opening angle, allowing the conspicuous shape of the asymptotic bow shock to be used as a preliminary constraint on more detailed modeling of the colliding winds. In the process, we explore from a general perspective the limitations in applying to the global shock geometry the so-called Dyson approximation, which asserts a local balance in the perpendicular ram pressure across the shock.
Bow-shock pulsar wind nebulae are a subset of pulsar wind nebulae that form when the pulsar has high velocity due to the natal kick during the supernova explosion. The interaction between the relativistic wind from the fast-moving pulsar and the interstellar medium produces a bow-shock and a trail, which are detectable in H$_{alpha}$ emission. Among such bow-shock pulsar wind nebulae, the Guitar Nebula stands out for its peculiar morphology, which consists of a prominent bow-shock head and a series of bubbles further behind. We present a scenario in which multiple bubbles can be produced when the pulsar encounters a series of density discontinuities in the ISM. We tested the scenario using 2-D and 3-D hydrodynamic simulations. The shape of the guitar nebula can be reproduced if the pulsar traversed a region of declining low density. We also show that if a pulsar encounters an inclined density discontinuity, it produces an asymmetric bow-shock head, consistent with observations of the bow-shock of the millisecond pulsar J2124-3358.