3D modelling of the colliding winds in Eta Carinae - evidence for radiative inhibition

The X-ray emission from the super-massive star Eta Carinae 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 degrees, consistent with the polar axis of the Homunculus nebula (Smith 2006), and the projection of the observer's line-of-sight onto the orbital plane has an angle of ~ 0 - 30 degrees in the prograde direction on the apastron side of the semi-major axis. However, there are significant discrepancies between the observed and model lightcurves and spectra through the X-ray minimum. In particular, the hard flux in our synthetic spectra is an order of magnitude greater than observed. Further calculations reveal that radiative inhibition significantly reduces the preshock velocity of the companion wind. As a consequence the hard X-ray emission is quenched, but it is unclear whether the long duration of the minimum is due solely to this mechanism alone. Models incorporating a collapse/disruption of the WCR and/or reduced preshock companion wind velocities bring the predicted emission and the observations into much better agreement (abridged).Comment: 20 pages, 24 Figures, accepted to MNRA

The 2.35 year itch of Cyg OB2 #9. II. Radio monitoring

Cyg OB2 #9 is one of a small set of non-thermal radio emitting massive O-star binaries. The non-thermal radiation is due to synchrotron emission in the colliding-wind region. Cyg OB2 #9 was only recently discovered to be a binary system and a multi-wavelength campaign was organized to study its 2011 periastron passage. We report here on the results of the radio observations obtained in this monitoring campaign. We used the Expanded Very Large Array (EVLA) radio interferometer to obtain 6 and 20 cm continuum fluxes. The observed radio light curve shows a steep drop in flux sometime before periastron. The fluxes drop to a level that is comparable to the expected free-free emission from the stellar winds, suggesting that the non-thermal emitting region is completely hidden at that time. After periastron passage, the fluxes slowly increase. We introduce a simple model to solve the radiative transfer in the stellar winds and the colliding-wind region, and thus determine the expected behaviour of the radio light curve. From the asymmetry of the light curve, we show that the primary has the stronger wind. This is somewhat unexpected if we use the astrophysical parameters based on theoretical calibrations. But it becomes entirely feasible if we take into account that a given spectral type - luminosity class combination covers a range of astrophysical parameters. The colliding-wind region also contributes to the free-free emission, which can help to explain the high values of the spectral index seen after periastron passage. Combining our data with older Very Large Array (VLA) data allows us to derive a period P = 860.0 +- 3.7 days for this system. With this period, we update the orbital parameters that were derived in the first paper of this series.Comment: 10 pages, 4 figures, accepted for publication in A&

X-ray emission from the double-binary OB-star system QZ Car (HD 93206)

X-ray observations of the double-binary OB-star system QZ Car (HD 93206) obtained with the Chandra X-ray Observatory over a period of roughly 2 years are presented. The orbit of systems A (O9.7 I+b2 v, PA = 21 d) and B (O8 III+o9 v, PB = 6 d) are reasonably well sampled by the observations, allowing the origin of the X-ray emission to be examined in detail. The X-ray spectra can be well fitted by an attenuated three temperature thermal plasma model, characterised by cool, moderate, and hot plasma components at kT ~ 0.2, 0.7, and 2 keV, respectively, and a circumstellar absorption of ~ 0.2 x 10^22 cm-2. Although the hot plasma component could be indicating the presence of wind-wind collision shocks in the system, the model fluxes calculated from spectral fits, with an average value of ~ 7 x 10^-13 erg s-1 cm-2, do not show a clear correlation with the orbits of the two constituent binaries. A semi-analytical model of QZ Car reveals that a stable momentum balance may not be established in either system A or B. Yet, despite this, system B is expected to produce an observed X-ray flux well in excess of the observations. If one considers the wind of the O8 III star to be disrupted by mass transfer the model and observations are in far better agreement, which lends support to the previous suggestion of mass-transfer in the O8 III + o9 v binary. We conclude that the X-ray emission from QZ Car can be reasonably well accounted for by a combination of contributions mainly from the single stars and the mutual wind-wind collision between systems A and B.Comment: 11 pages, 7 figures. Accepted for the ApJS Special Issue on the Chandra Carina Complex Project (CCCP), scheduled for publication in May 2011. All 16 CCCP Special Issue papers are available at http://cochise.astro.psu.edu/Carina_public/special_issue.html through 2011 at leas

Future sea level change from Antarctica's Lambert-Amery glacial system

Future global mean sea level (GMSL) change is dependent on the complex response of the Antarctic ice sheet to ongoing changes and feedbacks in the climate system. The Lambert-Amery glacial system has been observed to be stable over the recent period yet is potentially at risk of rapid grounding line retreat and ice discharge given that a significant volume of its ice is grounded below sea level, making its future contribution to GMSL uncertain. Using a regional ice sheet model of the Lambert-Amery system, we find that under a range of future warming and extreme scenarios, the simulated grounding line remains stable and does not trigger rapid mass loss from grounding line retreat. This allows for increased future accumulation to exceed the mass loss from ice dynamical changes. We suggest that the Lambert-Amery glacial system will remain stable or gain ice mass and mitigate a portion of potential future sea level rise over the next 500 years, with a range of +3.6 to −117.5 mm GMSL equivalent

Variable millimetre radiation from the colliding-wind binary Cygnus OB2 #8A

Context. Massive binaries have stellar winds that collide. In the colliding-wind region, various physically interesting processes occur, leading to enhanced X-ray emission, non-thermal radio emission, as well as non-thermal X-rays and gamma-rays. Non-thermal radio emission (due to synchrotron radiation) has so far been observed at centimetre wavelengths. At millimetre wavelengths, the stellar winds and the colliding-wind region emit more thermal free-free radiation, and it is expected that any non-thermal contribution will be difficult or impossible to detect. Aims. We aim to determine if the material in the colliding-wind region contributes substantially to the observed millimetre fluxes of a colliding-wind binary. We also try to distinguish the synchrotron emission from the free-free emission. Methods. We monitored the massive binary Cyg OB2 #8A at 3 mm with the NOrthern Extended Millimeter Array (NOEMA) interferometer of the Institut de Radioastronomie Millimétrique (IRAM). The data were collected in 14 separate observing runs (in 2014 and 2016), and provide good coverage of the orbital period. Results. The observed millimetre fluxes range between 1.1 and 2.3 mJy, and show phase-locked variability, clearly indicating that a large part of the emission is due to the colliding-wind region. A simple synchrotron model gives fluxes with the correct order of magnitude, but with a maximum that is phase-shifted with respect to the observations. Qualitatively this phase shift can be explained by our neglect of orbital motion on the shape of the colliding-wind region. A model using only free-free emission results in only a slightly worse explanation of the observations. Additionally, on the map of our observations we also detect the O6.5 III star Cyg OB2 #8B, for which we determine a 3 mm flux of 0.21 ± 0.033 mJy. Conclusions. The question of whether synchrotron radiation or free-free emission dominates the millimetre fluxes of Cyg OB2 #8A remains open. More detailed modelling of this system, based on solving the hydrodynamical equations, is required to give a definite answer

Waiting in the Wings: Reflected X-ray Emission from the Homunculus Nebula

We report the first detection of X-ray emission associated with the Homunculus Nebula which surrounds the supermassive star Eta Carinae. The emission is characterized by a temperature in excess of 100 MK, and is consistent with scattering of the time-delayed X-ray flux associated with the star. The nebular emission is bright in the northwestern lobe and near the central regions of the Homunculus, and fainter in the southeastern lobe. We also report the detection of an unusually broad Fe K fluorescent line, which may indicate fluorescent scattering off the wind of a companion star or some other high velocity outflow. The X-ray Homunculus is the nearest member of the small class of Galactic X-ray reflection nebulae, and the only one in which both the emitting and reflecting sources are distinguishable.Comment: 17 pages, 5 figures, accepted by Ap

Chandra X-ray Grating Spectrometry of Eta Carinae near X-ray Minimum: I. Variability of the Sulfur and Silicon Emission Lines

We report on variations in important X-ray emission lines in a series of Chandra grating spectra of the supermassive colliding wind binary star eta Car, including key phases around the X-ray minimum/periastron passage in 2003.5. The X-rays arise from the collision of the slow, dense wind of eta Car with the fast, low-density wind of an otherwise hidden companion star. The X-ray emission lines provide the only direct measure of the flow dynamics of the companion's wind along the wind-wind collision zone. We concentrate here on the silicon and sulfur lines, which are the strongest and best resolved lines in the X-ray spectra. Most of the line profiles can be adequately fit with symmetric Gaussians with little significant skewness. Both the silicon and sulfur lines show significant velocity shifts and correlated increases in line widths through the observations. The R = forbidden-to-intercombination ratio from the Si XIII and S XV triplets is near or above the low-density limit in all observations, suggesting that the line-forming region is > 1.6 stellar radii from the companion star, and that the emitting plasma may be in a non-equilibrium state. We show that simple geometrical models cannot simultaneously fit both the observed centroid variations and changes in line width as a function of phase. We show that the observed profiles can be fitted with synthetic profiles with a reasonable model of the emissivity along the wind-wind collision boundary. We use this analysis to help constrain the line formation region as a function of orbital phase, and the orbital geometry. Subject headings: X-rays: stars -stars: early-type-stars: individual (q Car

Spiralling out of control: 3D hydrodynamical modelling of the colliding winds in η \eta\thinspaceCarinae

Three dimensional (3D) adaptive-mesh refinement (AMR) hydrodynamical simulations of the wind-wind collision between the enigmatic super-massive star \etacar and its mysterious companion star are presented which include radiative driving of the stellar winds, gravity, optically-thin radiative cooling, and orbital motion. Simulations with static stars with a periastron passage separation reveal that the preshock companion star's wind speed is sufficiently reduced that radiative cooling in the postshock gas becomes important, permitting the runaway growth of non-linear thin shell (NTSI) instabilities which massively distort the WCR. However, large-scale simulations which include the orbital motion of the stars, show that orbital motion reduces the impact of radiative inhibition, and thus increases the acquired preshock velocities. As such, the postshock gas temperature and cooling time see a commensurate increase, and sufficient gas pressure is preserved to stabilize the WCR against catastrophic instability growth. We then compute synthetic X-ray spectra and lightcurves and find that, compared to previous models, the X-ray spectra agree much better with {\it XMM-Newton} observations just prior to periastron. The narrow width of the 2009 X-ray minimum can also be reproduced. However, the models fail to reproduce the extended X-ray mimimum from previous cycles. We conclude that the key to explaining the extended X-ray minimum is the rate of cooling of the companion star's postshock wind. If cooling is rapid then powerful NTSIs will heavily disrupt the WCR. Radiative inhibition of the companion star's preshock wind, albeit with a stronger radiation-wind coupling than explored in this work, could be an effective trigger.Comment: 25 pages, 20 figures, accepted for publication in Ap