Christopher M. P. Russell

RECENT X-RAY VARIABILITY OF η CARINAE: THE QUICK ROAD TO RECOVERY

We report continued monitoring of the superluminous binary system η Car by the Proportional Counter Array on the Rossi X-ray Timing Observatory (RXTE) through the 2009 X-ray minimum. The RXTE campaign shows that the minimum began on 2009 January 16, consistent with the phasings of the two previous minima, and overall, the temporal behavior of the X-ray emission was similar to that observed by RXTE in the previous two cycles. However, important differences did occur. The 2-10 keV X-ray flux and X-ray hardness decreased in the 2.5 year interval leading up to the 2009 minimum compared to the previous cycle. Most intriguingly, the 2009 X-ray minimum was about 1 month shorter than either of the previous two minima. During the egress from the 2009 minimum the X-ray hardness increased markedly as it had during egress from the previous two minima, although the maximum X-ray hardness achieved was less than the maximum observed after the two previous recoveries. We suggest that the cycle-to-cycle variations, especially the unexpectedly early recovery from the 2009 X-ray minimum, might have been the result of a decline in η Cars wind momentum flux produced by a drop in η Cars mass loss rate or wind terminal velocity (or some combination), though if so the change in wind momentum flux required to match the X-ray variation is surprisingly large.

Constraining the absolute orientation of η Carinae’s binary orbit: a 3D dynamical model for the broad [Fe iii] emission

We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission observed in Eta Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) simulations of Eta Cars binary colliding winds. Radiative transfer codes are used to generate synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA(theta) that the orbital plane projection of the line-of-sight makes with the apastron side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously reproduce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38deg, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i approx. = 130deg to 145deg, Theta approx. = -15deg to +30deg, and an orbital axis projected on the sky at a P A approx. = 302deg to 327deg east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, Eta(sub B), thus orbits clockwise on the sky and is on the observers side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modeling to determine the stellar masses.

Constraints on decreases in η Carinae's mass-loss from 3D hydrodynamic simulations of its binary colliding winds

Recent work suggests that the mass-loss rate of the primary star Eta-A in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present result from large- (+/- 1545 au) and small- (+/- 155 au) domain, 3D smoothed particle hydrodynamics (SPH) simulations of Eta Cars colliding winds for three Eta-A mass-loss rates ( (dot-M(sub Eta-A) = 2.4, 4.8 and 8.5 10(exp 4) M(solar)/ yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling and radiative forces. We find that dot-M Eta-A greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star Eta-B switches from the adiabatic to the radiative-cooling regime during periastron passage (Phi approx.= 0.985-1.02). This switchover starts later and ends earlier the lower the value of dot-M Eta-A and is caused by the encroachment of the wind of Eta-A into the acceleration zone of Eta-Bs wind, plus radiative inhibition of Eta-Bs wind by Eta-A. The SPH simulations together with 1D radiative transfer models of Eta-As spectra reveal that a factor of 2 or more drop in dot-M EtaA should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in dot- M Eta-A was likely by a factor of approx. < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line of sight to the system lies close to the dense walls of the WWC zone. A modest decrease in dot-M Eta-A may be responsible, but changes in the wind/stellar parameter of Eta-B, while less likely, cannot yet be fully ruled out. We suggest observations during Eta-Cars next periastron in 2014 to further test for decreases in dot-M Eta-A. If dot-M Eta-A is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.

Modelling the RXTE light curve of η Carinae from a 3D SPH simulation of its binary wind collision

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 observe r’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 orb ital 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 si de 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 Emission from Eta Carinae near Periastron in 2009. I. A Two-state Solution

X-ray emission from the supermassive binary system η Car declines sharply around periastron. This X-ray minimum has two distinct phases—the lowest flux phase in the first ~3 weeks and a brighter phase thereafter. In 2009, the Chandra X-ray Observatory monitored the first phase five times and found the lowest observed flux at ~1.9 × 10–12 erg cm–2 s–1 (3-8 keV). The spectral shape changed such that the hard band above ~4 keV dropped quickly at the beginning and the soft band flux gradually decreased to its lowest observed value in ~2 weeks. The hard band spectrum had begun to recover by that time. This spectral variation suggests that the shocked gas producing the hottest X-ray gas near the apex of the wind-wind collision (WWC) is blocked behind the dense inner wind of the primary star, which later occults slightly cooler gas downstream. Shocked gas previously produced by the system at earlier orbital phases is suggested to produce the faint residual X-ray emission seen when the emission near the apex is completely blocked by the primary wind. The brighter phase is probably caused by the re-appearance of the WWC plasma, whose emissivity significantly declined during the occultation. We interpret this to mean that the X-ray minimum is produced by a hybrid mechanism of an occultation and a decline in the emissivity of the WWC shock. We constrain timings of superior conjunction and periastron based on these results.

A Coordinated X-Ray and Optical Campaign of the Nearest Massive Eclipsing Binary, δ Orionis Aa. III. Analysis of Optical Photometric (MOST) and Spectroscopic (Ground-based) Variations

We report on both high-precision photometry from the Microvariability and Oscillations of Stars (MOST) space telescope and ground-based spectroscopy of the triple system δ Ori A, consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7 days, and a more distant tertiary (O9 IV P > 400 years). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for three weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy, we have a well-constrained radial velocity (RV) curve of Aa1. While we are unable to recover RV variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries that show evidence for tidally induced pulsations.

Suzaku monitoring of hard X-ray emission from η Carinae over a single binary orbital cycle

The Suzaku X-ray observatory monitored the supermassive binary system Eta Carinae 10 times during the whole 5.5 year orbital cycle between 2005-2011. This series of observations presents the first long-term monitoring of this enigmatic system in the extremely hard X-ray band between 15-40 keV. During most of the orbit, the 15-25 keV emission varied similarly to the 2-10 keV emission, indicating an origin in the hard energy tail of the kT ~4 keV wind-wind collision (WWC) plasma. However, the 15-25 keV emission declined only by a factor of 3 around periastron when the 2-10 keV emission dropped by two orders of magnitude due probably to an eclipse of the WWC plasma. The observed minimum in the 15-25 keV emission occurred after the 2-10 keV flux had already recovered by a factor of ~3. This may mean that the WWC activity was strong, but hidden behind the thick primary stellar wind during the eclipse. The 25-40 keV flux was rather constant through the orbital cycle, at the level measured with INTEGRAL in 2004. This result may suggest a connection of this flux component to the gamma-ray source detected in this field. The Helium-like Fe Kalpha line complex at ~6.7 keV became strongly distorted toward periastron as seen in the previous cycle. The 5-9 keV spectra can be reproduced well with a two-component spectral model, which includes plasma in collision equilibrium (CE) and a plasma in non-equilibrium ionization (NEI) with tau ~1e11 cm-3 s-1. The NEI plasma increases in importance toward periastron.

V444 Cygni X-ray and polarimetric variability: Radiative and Coriolis forces shape the wind collision region

We present results from a study of the eclipsing, colliding- wind binary V444 Cyg that uses a combination of X-ray and optical spectropolarimetric methods to describe the 3-D nature of the shock and wind structure within the system. We have created the most complete X-ray light curve of V444 Cyg to date using 40 ksec of new data from Swift, and 200 ksec of new and archived XMMNewton observations. In addition, we have characterized the intri nsic, polarimetric phase-dependent behavior of the strongest optical emission lines using data obtained with the University of Wisconsin’s Half-Wave Spectropolarimeter. We have detected evidence of the Coriolis distortion of the wind-wind collision in the X- ray regime, which manifests itself through asymmetric behavior around the eclipses in the system’s X-ray light curves. The large openi ng angle of the X-ray emitting region, as well as its location (i.e. the WN wind does not collide with the O star, but rather its wind) are evidence of radiative braking/inhibition occurring within the system. Additionally, the polarimetric results show evidence of the cavity the wind-wind collision region carves out of the Wolf-Rayet star’s wind.

Eta Carinae's Thermal X-Ray Tail Measured with XMM-Newton and NuSTAR

The evolved, massive highly eccentric binary system, η Car, underwent a periastron passage in the summer of 2014. We obtained two coordinated X-ray observations with XMM-Newton and NuSTAR during the elevated X-ray flux state and just before the X-ray minimum flux state around this passage. These NuSTAR observations clearly detected X-ray emission associated with η Car extending up to ~50 keV for the first time. The NuSTAR spectrum above 10 keV can be fit with the bremsstrahlung tail from a kT ~6 keV plasma. This temperature is ΔkT ~2 keV higher than those measured from the iron K emission line complex, if the shocked gas is in collisional ionization equilibrium. This result may suggest that the companion stars pre-shock wind velocity is underestimated. The NuSTAR observation near the X-ray minimum state showed a gradual decline in the X-ray emission by 40% at energies above 5 keV in a day, the largest rate of change of the X-ray flux yet observed in individual η Car observations. The column density to the hardest emission component, N H ~1024 H cm-2, marked one of the highest values ever observed for η Car, strongly suggesting the increased obscuration of the wind-wind colliding X-ray emission by the thick primary stellar wind prior to superior conjunction. Neither observation detected the power-law component in the extremely hard band that INTEGRAL and Suzaku observed prior to 2011. The power-law source might have faded before these observations.

Modelling the Central Constant Emission X-ray component of η Carinae

The X-ray emission of