X-ray Modeling of \eta\ Carinae and WR140 from SPH Simulations
christopher russell2011, arXiv (Cornell University)
Sign up for access to the world's latest research
checkGet notified about relevant papers
checkSave papers to use in your research
checkJoin the discussion with peers
checkTrack your impact
Abstract
The colliding wind binary (CWB) systems η Carinae and WR140 provide unique laboratories for X-ray astrophysics. Their wind-wind collisions produce hard X-rays that have been monitored extensively by several X-ray telescopes, including RXTE. To interpret these RXTE X-ray light curves, we model the wind-wind collision using 3D smoothed particle hydrodynamics (SPH) simulations. Adiabatic simulations that account for the absorption of X-rays from an assumed point source at the apex of the wind-collision shock cone by the distorted winds can closely match the observed 2-10keV RXTE light curves of both η Car and WR140. This point-source model can also explain the early recovery of η Car's X-ray light curve from the 2009.0 minimum by a factor of 2-4 reduction in the mass loss rate of η Car. Our more recent models relax the point-source approximation and account for the spatially extended emission along the wind-wind interaction shock front. For WR140, the computed X-ray light curve again matches the RXTE observations quite well. But for η Car, a hot, post-periastron bubble leads to an emission level that does not match the extended X-ray minimum observed by RXTE. Initial results from incorporating radiative cooling and radiatively-driven wind acceleration via a new anti-gravity approach into the SPH code are also discussed.
Related papers
X-Ray modeling of η Carinae & WR 140 from SPH simulations
Proceedings of the International Astronomical Union, 2010
The colliding wind binary (CWB) systems η Carinae and WR140 provide unique laboratories for X-ray astrophysics. Their wind-wind collisions produce hard X-rays that have been monitored extensively by several X-ray telescopes, including RXTE. To interpret these RXTE X-ray light curves, we model the wind-wind collision using 3D smoothed particle hydrodynamics (SPH) simulations. Adiabatic simulations that account for the absorption of X-rays from an assumed point source at the apex of the wind-collision shock cone by the distorted winds can closely match the observed 2-10keV RXTE light curves of both η Car and WR140. This point-source model can also explain the early recovery of η Car's X-ray light curve from the 2009.0 minimum by a factor of 2-4 reduction in the mass loss rate of η Car. Our more recent models relax the point-source approximation and account for the spatially extended emission along the wind-wind interaction shock front. For WR140, the computed X-ray light curve again matches the RXTE observations quite well. But for η Car, a hot, post-periastron bubble leads to an emission level that does not match the extended X-ray minimum observed by RXTE. Initial results from incorporating radiative cooling and radiatively-driven wind acceleration via a new anti-gravity approach into the SPH code are also discussed.
Radio, X-ray, and gamma-ray Emission Models of the Colliding Winds Binary WR 140
2006
We use hydrodynamical models of the wind-collision region (WCR) in the archetype colliding-wind system WR140 to determine the spatial and spectral distribution of the radio, X-ray and gamma-ray emission from shock accelerated electrons. Our calculations are for orbital phase 0.837 when the observed radio emission is close to maximum. Using the observed thermal X-ray emission together with the radio emission to constrain the mass-loss rates, we find that the O-star mass-loss rate is consistent with recent reductions for O4-5 supergiants. We demonstrate that radio VLBI observations of the WCR fail to constrain the opening angle. The observed low frequency turnover at ~3 GHz in the radio emission is due to free-free absorption, since models based on the Razin effect have an unacceptably large fraction of energy in non-thermal particles. The index of the non-thermal electron energy distribution is flatter than the canonical value for diffusive shock acceleration, namely p<2. Several mechanisms are discussed that could lead to such an index. Tighter constraints on p and the nature of the shocks in WR140 will be obtained from future observations at MeV and GeV energies, for which we generally predict lower fluxes than previous work. Since the high stellar photon fluxes prevent the acceleration of electrons beyond gamma > 1e5-1e6, TeV emission from CWB systems will provide unambiguous evidence of pion-decay emission from accelerated ions. We finish by commenting on the emission and physics of the multiple wind collisions in dense stellar clusters, paying particular attention to the Galactic Centre (abridged).
X-ray emission from hydrodynamical simulations in non-LTE wind models
Astronomy and Astrophysics, 2009
Hot stars are sources of X-ray emission originating in their winds. Although hydrodynamical simulations that are able to predict this X-ray emission are available, the inclusion of X-rays in stationary wind models is usually based on simplifying approximations. To improve this, we use results from time-dependent hydrodynamical simulations of the line-driven wind instability (seeded by the base perturbation) to derive the analytical approximation of X-ray emission in the stellar wind. We use this approximation in our non-LTE wind models and find that an improved inclusion of X-rays leads to a better agreement between model ionization fractions and those derived from observations. Furthermore, the slope of the LX − L relation is in better agreement with observations, however the X-ray luminosity is underestimated by a factor of three. We propose a possible solution for this discrepancy.
Hydrodynamic and radiative transfer modeling of X-ray emission from colliding WR winds: WR 140 & the Galactic center
arXiv (Cornell University), 2015
Colliding Wolf-Rayet (WR) winds produce thermal X-ray emission widely observed by X-ray telescopes. In wide WR+O binaries, such as WR 140, the X-ray flux is tied to the orbital phase, and is a direct probe of the winds' properties. In the Galactic center, $\sim$30 WRs orbit the super massive black hole (SMBH) within $\sim$10", leading to a smorgasbord of wind-wind collisions. To model the X-ray emission of WR 140 and the Galactic center, we perform 3D hydrodynamic simulations to trace the complex gaseous flows, and then carry out 3D radiative transfer calculations to compute the variable X-ray spectra. The model WR 140 RXTE light curve matches the data well for all phases except the X-ray minimum associated with periastron, while the model spectra agree with the RXTE hardness ratio and the shape of the Suzaku observations throughout the orbit. The Galactic center model of the Chandra flux and spectral shape match well in the region r$<$3", but the model flux falls off too rapidly beyond this radius.
Radio, X-ray, and gamma-ray emission models of the colliding-wind binary WR140
Mon Notic Roy Astron Soc, 2006
We use hydrodynamical models of the wind-collision region (WCR) in the archetype colliding-wind system WR 140 to determine the spatial and spectral distribution of the radio, X-ray and γ-ray emission from shock accelerated electrons. Our calculations are for orbital phase 0.837 when the observed radio emission is close to maximum. Using the observed thermal X-ray emission at this phase in conjunction with the radio emission to constrain the mass-loss rates, we find that the O-star mass-loss rate is consistent with the reduced estimates for O4-5 supergiants by Fullerton, Massa & Prinja (2006), and the wind momentum ratio, η = 0.02. This is independent of the opening angle deduced from radio VLBI observations of the WCR that we demonstrate fail to constrain the opening angle. We show that the turnover at ∼ 3 GHz in the radio emission is due to freefree absorption, since models based on the Razin effect have an unacceptably large fraction of energy in non-thermal electrons. We find the spectral index of the nonthermal electron energy distribution is flatter than the canonical value for diffusive shock acceleration (DSA), namely p < 2. Several mechanisms are discussed that could lead to such an index. Our inability to obtain fits to the radio data with p > 2 does not exclude the possibility of shock modification, but stronger evidence than currently exists is necessary for its support. Tighter constraints on p and the nature of the shocks in WR 140 will be obtained from future observations at MeV and GeV energies, for which we generally predict lower fluxes than previous work. Since the high stellar photon fluxes prevent the acceleration of electrons beyond γ ∼ > 10 5 − 10 6 , TeV emission from colliding-wind binary (CWB) systems will provide unambiguous evidence of pion-decay emission from accelerated ions. We finish by commenting on the emission and physics of the multiple wind collisions in dense stellar clusters, paying particular attention to the Galactic Center.
Probing the wind-wind collision in 2 Velorum with high-resolution Chandra X-ray spectroscopy: evidence for sudden radiative braking and non-equilibrium ionization
Mon Notic Roy Astron Soc, 2004
We present a new analysis of an archived Chandra HETGS X-ray spectrum of the WR+O colliding wind binary Gamma Velorum. The spectrum is dominated by emission lines from astrophysically abundant elements: Ne, Mg, Si, S and Fe. From a combination of broad-band spectral analysis and an analysis of line flux ratios we infer a wide range of temperatures in the X-ray emitting plasma (~4-40 MK). As in the previously published analysis, we find the X-ray emission lines are essentially unshifted, with a mean FWHM of 1240 +/- 30 km/s. Calculations of line profiles based on hydrodynamical simulations of the wind-wind collision predict lines that are blueshifted by a few hundred km/s. The lack of any observed shift in the lines may be evidence of a large shock-cone opening half-angle (> 85 degrees), and we suggest this may be evidence of sudden radiative braking. From the R and G ratios measured from He-like forbidden-intercombination-resonance triplets we find evidence that the Mg XI emission originates from hotter gas closer to the O star than the Si XIII emission, which suggests that non-equilibrium ionization may be present.
Modelling the thermal X-ray emission around the Galactic centre from colliding Wolf-Rayet winds
Proceedings of the International Astronomical Union, 2016
The Galactic centre is a hotbed of astrophysical activity, with the injection of wind material from ~30 massive Wolf-Rayet (WR) stars orbiting within 12″ of the super-massive black hole (SMBH) playing an important role. Hydrodynamic simulations of such colliding and accreting winds produce a complex density and temperature structure of cold wind material shocking with the ambient medium, creating a large reservoir of hot, X-ray-emitting gas. This work aims to confront the 3Ms of Chandra X-ray Visionary Program (XVP) observations of this diffuse emission by computing the X-ray emission from these hydrodynamic simulations of the colliding WR winds, amid exploring a variety of SMBH feedback mechanisms. The major success of the model is that it reproduces the spectral shape from the 2″–5″ ring around the SMBH, where most of the stellar wind material that is ultimately captured by Sgr A* is shock-heated and thermalised. This naturally explains that the hot gas comes from colliding WR win...
Theoretical X-ray line profiles from colliding wind binaries
Monthly Notices of the Royal Astronomical Society, 2003
We present theoretical X-ray line profiles from a range of model colliding wind systems. In particular, we investigate the effects of varying the stellar mass-loss rates, the wind speeds, and the viewing orientation. We find that a wide range of theoretical line profile shapes is possible, varying with orbital inclination and phase. At or near conjunction, the lines have approximately Gaussian profiles, with small widths (HWHM ∼ 0.1v ∞ ) and definite blue-or redshifts (depending on whether the star with the weaker wind is in front or behind). When the system is viewed at quadrature, the lines are generally much broader (HWHM ∼ v ∞ ), flat-topped and unshifted. Local absorption can have a major effect on the observed profiles -in systems with mass-loss rates of a few times 10 −6 M ⊙ yr −1 the lower energy lines (E < ∼ 1 keV) are particularly affected. This generally results in blueward-skewed profiles, especially when the system is viewed through the dense wind of the primary. The orbital variation of the line widths and shifts is reduced in a low inclination binary. The extreme case is a binary with i = 0 • , for which we would expect no line profile variation.
An X-ray Survey of Colliding Wind Binaries
2012
We have compiled a list of 36 O+O and 89 Wolf-Rayet binary candidates in the Milky Way and Magellanic clouds detected with the Chandra, XMM-Newton and ROSAT satellites to probe the connection between their X-ray properties and their system characteristics. Of the WR binaries with published parameters, all but two have kT > 0.9 keV. The most X-ray luminous WR binaries are typically very long period systems. The WR binaries show a nearly four-order of magnitude spread in X-ray luminosity, even among among systems with very similar WR primaries. Among the O+O binaries, short-period systems generally have soft X-ray spectra and longer period systems show harder X-ray spectra again with a large spread in L X /L bol .
Probing the wind-wind collision in γ 2 Velorum with high-resolution Chandra X-ray spectroscopy: evidence for sudden radiative braking and non-equilibrium ionization
Monthly Notices of the Royal Astronomical Society, 2005
We present a new analysis of an archived Chandra HETGS X-ray spectrum of the WR+O colliding wind binary γ 2 Velorum. The spectrum is dominated by emission lines from astrophysically abundant elements: Ne, Mg, Si, S and Fe. From a combination of broad-band spectral analysis and an analysis of line flux ratios we infer a wide range of temperatures in the X-ray-emitting plasma (∼4-40 MK). As in the previously published analysis, we find the X-ray emission lines are essentially unshifted, with a mean FWHM of 1240 ± 30 km s −1 .
Related topics
Related papers
Probing Colliding Wind Binaries with High-Resolution X-ray Spectra
X-ray line profiles represent a new way of studying the winds of massive stars. In particular, they enable us to probe in detail the wind-wind collision in colliding wind binaries, providing new insights into the structure and dynamics of the X-ray-emitting regions. We present the key results of new analyses of high-resolution Chandra X-ray spectra of two important colliding wind systems, Gamma Velorum and WR140. The lines of Gamma Vel are essentially unshifted from their rest wavelengths, which we suggest is evidence of a wide shock opening angle, indicative of sudden radiative braking. The widths of the lines of WR140 are correlated with ionization potential, implying non-equilibrium ionization. The implications of these results for the radio emission from these systems are discussed, as are some of the future directions for X-ray line profile modelling of colliding wind binaries.
3D modelling of the colliding winds in η Carinae - evidence for radiative inhibition
Monthly Notices of the Royal Astronomical Society, 2009
The X-ray emission from the super-massive star η Car is simulated using a three dimensional model of the wind-wind collision. In the model the intrinsic X-ray emission is spatially extended and energy dependent. Absorption due to the unshocked stellar winds and the cooled postshock material from the primary LBV star is calculated as the intrinsic emission is ray-traced along multiple sightlines through the 3D spiral structure of the circumstellar environment. The observable emission is then compared to available X-ray data, including the lightcurve observed by the Rossi X-ray Timing Explorer (RXTE) and spectra observed by XMM-Newton. The orientation and eccentricity of the orbit are explored, as are the wind parameters of the stars and the nature and physics of their close approach. Our modelling supports a viewing angle with an inclination of ≃ 42 • , consistent with the polar axis of the Homunculus nebula , and the projection of the observer's line-of-sight onto the orbital plane has an angle of ≃ 0 − 30 • in the prograde direction on the apastron side of the semi-major axis.
Particle Re-Acceleration in Colliding Winds Systems? Radio, X-ray, and gamma-ray Emission Models of WR 140
We present calculations of the spatial and spectral distribution of the radio, X-ray and γ-ray emission from shock accelerated electrons in the wind-collision region (WCR) of WR140. Our calculations are for orbital phase 0.837 when the observed radio emission is close to maximum. Using the observed thermal X-ray emission at this phase in conjunction with the radio emission to constrain the mass-loss rates, we find that the O-star mass-loss rate is consistent with the reduced estimates for O4-5 supergiants by Fullerton, , and the wind momentum ratio, η = 0.02.
Radio, X-ray, and γ-ray emission models of the colliding-wind binary WR 140
Monthly Notices of the Royal Astronomical Society, 2006
We use hydrodynamical models of the wind-collision region in the archetype colliding-wind system WR 140 to determine the spatial and spectral distributions of the radio, X-ray, and γ -ray emission from shock-accelerated electrons. Our calculations are for orbital phase 0.837 when the observed radio emission is close to maximum. Using the observed thermal X-ray emission at this phase in conjunction with the radio emission to constrain the mass-loss rates, we find that the O star mass-loss rate is consistent with the reduced estimates for O4-5 supergiants by Fullerton, Massa & Prinja, and the wind-momentum ratio, η = 0.02. This is independent of the opening angle deduced from radio very long baseline interferometry observations of the WCR that we demonstrate fail to constrain the opening angle.
Colliding stellar winds structure and X-ray emission
Monthly Notices of the Royal Astronomical Society
We investigate the structure and X-ray emission from the colliding stellar winds in massive star binaries. We find that the opening angle of the contact discontinuity (CD) is overestimated by several formulae in the literature at very small values of the wind momentum ratio, η. We find also that the shocks in the primary (dominant) and secondary winds flare by ≈20 • compared to the CD, and that the entire secondary wind is shocked when η 0.02. Analytical expressions for the opening angles of the shocks, and the fraction of each wind that is shocked, are provided. We find that the X-ray luminosity L x ∝ η, and that the spectrum softens slightly as η decreases.
Non-thermal X-ray and gamma-ray emission from the colliding wind binary WR 140
WR 140 is the archetype long-period colliding wind binary (CWB) system, and is well known for dramatic variations in its synchrotron emission during its 7.9-yr, highly eccentric orbit. This emission is thought to arise from relativistic electrons accelerated at the global shocks bounding the wind-collision region (WCR). The presence of non-thermal electrons and ions should also give rise to X-ray and gamma-ray emission from several separate mechanisms, including inverse-Compton cooling, relativistic bremsstrahlung, and pion decay. We describe new calculations of this emission and make some preliminary predictions for the new generation of gamma-ray observatories. We determine that WR 140 will likely require several Megaseconds of observation before detection with INTEGRAL, but should be a reasonably strong source for GLAST.
Modelling the<i>RXTE</i>light curve of Carinae from a 3D SPH simulation of its binary wind collision
Monthly Notices of the Royal Astronomical Society: Letters, 2008
The very massive star system η Carinae exhibits regular 5.54 yr (2024 d) period disruptive events in wavebands ranging from the radio to X-ray. There is a growing consensus that these events likely stem from periastron passage of an (as yet) unseen companion in a highly eccentric (e ∼ 0.9) orbit. This Letter presents 3D smoothed particle hydrodynamics (SPH) simulations of the orbital variation of the binary wind-wind collision, and applies these to modelling the X-ray light curve observed by the Rossi X-ray Timing Explorer (RXTE). By providing a global 3D model of the phase variation of the density of the interacting winds, the simulations allow computation of the associated variation in X-ray absorption, presumed here to originate from near the apex of the wind-wind interaction cone. We find that the observed RXTE light curve can be readily fitted if the observer's line of sight is within this cone along the general direction of apastron. Specifically, the data are well fitted by an assumed inclination i = 45 • for the orbit's polar axis, which is thus consistent with orbital angular momentum being along the inferred polar axis of the Homunculus nebula. The fits also constrain the position angle φ that an orbital-plane projection makes with the apastron side of the semimajor axis, strongly excluding positions φ < 9 • along or to the retrograde side of the axis, with the best-fitting position given by φ = 27 • . Overall the results demonstrate the utility of a fully 3D dynamical model for constraining the geometric and physical properties of this complex colliding wind binary system.
Modelling the RXTE light curve of Carinae from a 3D SPH simulation of its binary wind collision
Monthly Notices of the Royal Astronomical Society: Letters, 2008
The very massive star system η Carinae exhibits regular 5.54-year (2024-day) period disruptive events in wavebands ranging from the radio to X-ray. There is a growing consensus that these events likely stem from periastron passage of an (as yet) unseen companion in a highly eccentric (ǫ ∼ 0.9) orbit. This paper presents three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the orbital variation of the binary wind-wind collision, and applies these to modeling the X-ray light curve observed by the Rossi X-ray Timing Explorer (RXTE). By providing a global 3-D model of the phase variation of the density of the interacting winds, the simulations allow computation of the associated variation in X-ray absorption, presumed here to originate from near the apex of the wind-wind interaction cone. We find that the observed RXTE light curve can be readily fit if the observer's line of sight is within this cone along the general direction of apastron. Specifically, the data are well fit by an assumed inclination i = 45 • for the orbit's polar axis, which is thus consistent with orbital angular momentum being along the inferred polar axis of the Homunculus nebula. The fits also constrain the position angle φ that an orbital-plane projection makes with the apastron side of the semi-major axis, strongly excluding positions φ < 9 • along or to the retrograde side of the axis, with the best fit position given by φ = 27 • . Overall the results demonstrate the utility of a fully 3-D dynamical model for constraining the geometric and physical properties of this complex colliding-wind binary system.
X-ray and radio emission from colliding stellar winds
The collision of the hypersonic winds in early-type binaries produces shock heated gas, which radiates thermal X-ray emission, and relativistic electrons, which emit nonthermal radio emission. We review our current understanding of the emission in these spectral regions and discuss models which have been developed for the interpretation of this emission. Physical processes which affect the resulting emission are reviewed and ideas for the future are noted.
Preprint typeset using LATEX style emulateapj v. 11/10/09 MODELING BROADBAND X-RAY ABSORPTION OF MASSIVE STAR WINDS
2010
We present a method for computing the net transmission of X-rays emitted by shock-heated plasma distributed throughout a partially optically thick stellar wind from a massive star. We find the transmission by an exact integration of the formal solution, assuming that the emitting plasma and absorbing plasma are mixed at a constant mass ratio above some minimum radius, below which there is assumed to be no emission. This model is more realistic than either the slab absorption associated with a corona at the base of the wind or the exospheric approximation that assumes that all observed X-rays are emitted without attenuation from above the radius of optical depth unity. Our model is implemented in XSPEC as a pre-calculated table that can be coupled to a user-defined table of the wavelength dependent wind opacity. We provide a default wind opacity model that is more representative of real wind opacities than the commonly used neutral interstellar medium (ISM) tabulation. Preliminary modeling of Chandra grating data indicates that the X-ray hardness trend of OB stars with spectral subtype can largely be understood as a wind absorption effect.
Loading Preview
Sorry, preview is currently unavailable. You can download the paper by clicking the button above.