Papers by christopher russell

arXiv (Cornell University), 2015
Colliding Wolf-Rayet (WR) winds produce thermal X-ray emission widely observed by X-ray telescope... more 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.

Proceedings of the International Astronomical Union, Jul 1, 2010
The colliding wind binary (CWB) systems η Carinae and WR140 provide unique laboratories for X-ray... more 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 apply 3D hydrodynamic simulations of the wind-wind collision using smoothed particle hydrodynamics (SPH). We find adiabatic simulations that account for the absorption of X-rays from an assumed point source of X-ray emission at the apex of the wind-collision shock cone can closely match the 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 account for the extended emission and absorption along the full wind-wind interaction shock front. For WR140, the computed X-ray light curves again match 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 radiative forces via an anti-gravity approach into the SPH code are also discussed.

Proceedings of the International Astronomical Union, Nov 1, 2016
Eta Carinae is one of the most massive observable binaries. Yet determination of its orbital and ... more Eta Carinae is one of the most massive observable binaries. Yet determination of its orbital and physical parameters is hampered by obscuring winds. However the effects of the strong, colliding winds changes with phase due to the high orbital eccentricity. We wanted to improve measures of the orbital parameters and to determine the mechanisms that produce the relatively brief, phase-locked minimum as detected throughout the electromagnetic spectrum. We conducted intense monitoring of the He ii λ4686 line in η Carinae for 10 months in the year 2014, gathering ∼300 high S/N spectra with ground-and space-based telescopes. We also used published spectra at the FOS4 SE polar region of the Homunculus, which views the minimum from a different direction. We used a model in which the He ii λ4686 emission is produced by two mechanisms: a) one linked to the intensity of the wind-wind collision which occurs along the whole orbit and is proportional to the inverse square of the separation between the companion stars; and b) the other produced by the 'bore hole' effect which occurs at phases across the periastron passage. The opacity (computed from 3D SPH simulations) as convolved with the emission reproduces the behavior of equivalent widths both for direct and reflected light. Our main results are: a) a demonstration that the He ii λ4686 light curve is exquisitely repeatable from cycle to cycle, contrary to previous claims for large changes; b) an accurate determination of the longitude of periastron, indicating that the secondary star is 'behind' the primary at periastron, a dispute extended over the past decade; c) a determination of the time of periastron passage, at ∼4 days after the onset of the deep light curve minimum; and d) show that the minimum is simultaneous for observers at different lines of sight, indicating that it is not caused by an eclipse of the secondary star, but rather by the immersion of the wind-wind collision interior to the inner wind of the primary.

The extremely massive (> 90 M O ) and luminous (= 5 x 106 Lo) star Eta Carinae, with its spectacu... more The extremely massive (> 90 M O ) and luminous (= 5 x 106 Lo) star Eta Carinae, with its spectacular bipolar "Homunculus" nebula, comprises one of the most remarkable and intensely observed stellar systems in the galaxy. However, many of its underlying physical parameters remain a mystery. Multiwavelength variations observed to occur every 5.54 years are interpreted as being due to the collision of a massive wind from the primary star with the fast, less dense wind of a hot companion star in a highly elliptical (e -0.9) orbit. Using three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the binary wind-wind collision in Eta Car, together with radiative transfer codes, we compute synthetic spectral images of [Fe 111] emission line structures and compare them to existing Hubble Space TelescopelSpace Telescope Imaging Spectrograph (HST/STIS) observations. We are thus able, for the first time, to constrain the absolute orientation of the binary orbit on the sky. An orbit with an inclination of i -40°, an argument of periapsis w -255°, and a projected orbital axis with a position angle of -312° east of north provides the best fit to the observations, implying that the orbital axis is closely aligned in 3-1) space with the Homunculus symmetry axis, and that the companion star orbits clockwise on the sky relative to the primary.

Monthly Notices of the Royal Astronomical Society: Letters, Jul 21, 2008
The very massive star system η Carinae exhibits regular 5.54 yr (2024 d) period disruptive events... more 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.

The Astrophysical Journal, Aug 19, 2015
We present time-resolved and phase-resolved variability studies of an extensive X-ray highresolut... more We present time-resolved and phase-resolved variability studies of an extensive X-ray highresolution spectral dataset of the δ Ori Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of ≈479 ks and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range 5-25 Å is confirmed, with maximum amplitude of about ±15% within a single ≈125 ks observation. Periods of 4.76d and 2.04d are found in the total X-ray flux, as well as an apparent overall increase in flux level throughout the 9-day observational campaign. Using 40 ks contiguous spectra derived from the original observations, we investigate variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at φ=0.0 when the secondary δ Ori Aa2 is at inferior conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission line width variability to the presence of a wind cavity created by a wind-wind collision, which is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary wind, thus producing phase locked X-ray variability.

arXiv (Cornell University), Oct 8, 2011
The colliding wind binary (CWB) systems η Carinae and WR140 provide unique laboratories for X-ray... more 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.

arXiv (Cornell University), Sep 20, 2021
WR 140 is a long-period, highly eccentric Wolf-Rayet star binary system with exceptionally wellde... more WR 140 is a long-period, highly eccentric Wolf-Rayet star binary system with exceptionally welldetermined orbital and stellar parameters. Bright, variable X-ray emission is generated in shocks produced by the collision of the winds of the WC7pd+O5.5fc component stars. We discuss the variations in the context of the colliding-wind model using broad-band spectrometry from the RXTE, Swift, and NICER observatories obtained over 20 years and nearly 1000 observations through 3 consecutive 7.94-year orbits including 3 periastron passages. The X-ray luminosity varies as expected with the inverse of the stellar separation over most of the orbit: departures near periastron are produced when cooling shifts to excess optical emission in C III λ5696 in particular. We use X-ray absorption to estimate mass-loss rates for both stars and to constrain the system morphology. The absorption maximum coincides closely with inferior conjunction of the WC star and provides evidence of the ion-reflection mechanism that underlies the formation of collisionless shocks governed by magnetic fields probably generated by the Weibel instability. Comparisons with K-band emission and He I λ10830 absorption show that both are correlated after periastron with the asymmetric X-ray absorption. Dust appears within a few days of periastron suggesting formation within shocked gas near the stagnation point. X-ray flares seen in η Car have not occurred in WR 140, suggesting the absence of large-scale wind inhomogeneities. Relatively constant soft emission revealed during the X-ray minimum is probably not from recombining plasma entrained in outflowing shocked gas.

arXiv (Cornell University), May 19, 2014
Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars ... more Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars in these systems possess powerful stellar winds with speeds up to ∼3000 km s -1 , and their collision leads to hot plasma (up to ∼ 10 8 K) that emit thermal X-rays (up to ∼10 keV). Many X-ray telescopes have observed CWBs, including Suzaku, and our work aims to model these X-ray observations. We use 3D smoothed particle hydrodynamics (SPH) to model the wind-wind interaction, and then perform 3D radiative transfer to compute the emergent X-ray flux, which is folded through X-ray telescopes' response functions to compare directly with observations. In these proceedings, we present our models of Suzaku observations of the multi-year-period, highly eccentric systems η Carinae and WR 140. The models reproduce the observations well away from periastron passage, but only η Carinae's X-ray spectrum is reproduced at periastron; the WR 140 model produces too much flux during this more complicated phase.
Nature Astronomy, Jul 2, 2018
Cosmic-ray acceleration has been a long-standing mystery and despite more than a century of study... more Cosmic-ray acceleration has been a long-standing mystery and despite more than a century of study, we still do not have a complete census of acceleration mechanisms. The collision of strong stellar winds in massive binary systems creates powerful shocks, which have been expected to produce high-energy cosmic-rays through Fermi acceleration at the shock interface. The accelerated particles should collide with stellar photons or ambient material, producing non-thermal emission observable in X-rays and γ-rays [3,. The supermassive binary star η Car drives the strongest colliding wind shock in the solar neighborhood . Observations with non-focusing high-energy observatories indicate a high energy source near η Car, but have been unable to conclusively identify η Car as the source because of their rela-

arXiv (Cornell University), Oct 8, 2011
The colliding wind binary (CWB) systems η Carinae and WR140 provide unique laboratories for X-ray... more 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.

arXiv (Cornell University), Jul 28, 2015
Eta Carinae (η Car) is the nearest example of a supermassive, superluminous, unstable star. Mass ... more Eta Carinae (η Car) is the nearest example of a supermassive, superluminous, unstable star. Mass loss from the system is critical in shaping its circumstellar medium and in determining its ultimate fate. Eta Car currently loses mass via a dense, slow stellar wind and possesses one of the largest mass loss rates known. It is prone to episodes of extreme mass ejection via eruptions from some as-yet unspecified cause; the best examples of this are the large-scale eruptions which occurred in 19th century. Eta Car is a colliding wind binary in which strong variations in X-ray emission and in other wavebands are driven by the violent collision of the wind of η Car-A and the fast, less dense wind of an otherwise hidden companion star. X-ray variations are the simplest diagnostic we have to study the wind-wind collision and allow us to measure the state of the stellar mass loss from both stars. We present the X-ray lightcurve over the last 20 years from ROSAT observations and monitoring with the Rossi X-ray Timing Explorer and the X-ray Telescope on the Swift satellite. We compare and contrast the behavior of the X-ray emission from the system over that timespan, including surprising variations during the 2014 X-ray minimum.

arXiv (Cornell University), May 19, 2014
Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars ... more Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars in these systems possess powerful stellar winds with speeds up to ∼3000 km s -1 , and their collision leads to hot plasma (up to ∼ 10 8 K) that emit thermal X-rays (up to ∼10 keV). Many X-ray telescopes have observed CWBs, including Suzaku, and our work aims to model these X-ray observations. We use 3D smoothed particle hydrodynamics (SPH) to model the wind-wind interaction, and then perform 3D radiative transfer to compute the emergent X-ray flux, which is folded through X-ray telescopes' response functions to compare directly with observations. In these proceedings, we present our models of Suzaku observations of the multi-year-period, highly eccentric systems η Carinae and WR 140. The models reproduce the observations well away from periastron passage, but only η Carinae's X-ray spectrum is reproduced at periastron; the WR 140 model produces too much flux during this more complicated phase.

Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars ... more Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars in these systems possess powerful stellar winds with speeds up to∼3000 km s−1, and their collision leads to hot plasma (up to ∼ 108K) that emit thermal X-rays (up to ∼10 keV). Many X-ray telescopes have observed CWBs, including Suzaku, and our work aims to model these X-ray observations. We use 3D smoothed particle hydrodynamics (SPH) to model the wind-wind interaction, and then perform 3D radiative transfer to compute the emergent X-ray flux, which is folded through X-ray telescopes ’ response functions to compare directly with observations. In these proceedings, we present our models of Suzaku observations of the multi-year-period, highly eccentric systems η Carinae and WR 140. The models reproduce the observations well away from periastron passage, but only η Carinae’s X-ray spectrum is reproduced at periastron; the WR 140 model produces too much flux during this more complicated phase.

Monthly Notices of the Royal Astronomical Society, 2017
We report on the first multi-colour precision light curve of the bright Wolf-Rayet binary γ 2 Vel... more We report on the first multi-colour precision light curve of the bright Wolf-Rayet binary γ 2 Velorum, obtained over six months with the nanosatellites in the BRITE-Constellation fleet. In parallel, we obtained 488 high-resolution optical spectra of the system. In this first report on the data sets, we revise the spectroscopic orbit and report on the bulk properties of the colliding winds. We find a dependence of both the light curve and excess emission properties that scales with the inverse of the binary separation. When analysing the spectroscopic properties in combination with the photometry, we find that the phase dependence is caused only by excess emission in the lines, and not from a changing continuum. We also detect a narrow, high-velocity absorption component from the He I λ5876 transition, which appears twice in the orbit. We calculate smoothed-particle hydrodynamical simulations of the colliding winds and can accurately associate the absorption from He I to the leading and trailing arms of the wind shock cone passing tangentially through our line of sight. The simulations also explain the general strength and kinematics of the emission excess observed in wind lines such as C III λ5696 of the system. These results represent the first in a series of investigations into the winds and properties of γ 2 Velorum through multi-technique and multi-wavelength observational campaigns.

Proceedings of the International Astronomical Union, 2016
The Galactic centre is a hotbed of astrophysical activity, with the injection of wind material fr... more 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...

The Astrophysical Journal, 2016
η Carinae is an extremely massive binary system in which rapid spectrum variations occur near per... more η Carinae is an extremely massive binary system in which rapid spectrum variations occur near periastron. Most notably, near periastron the He ii λ4686 line increases rapidly in strength, drops to a minimum value, then increases briefly before fading away. To understand this behavior, we conducted an intense spectroscopic monitoring of the He ii λ4686 emission line across the 2014.6 periastron passage using ground-and space-based telescopes. Comparison with previous data confirmed the overall repeatability of EW (He ii λ4686), the line radial velocities, and the timing of the minimum, though the strongest peak was systematically larger in 2014 than in 2009 by 26%. The EW (He ii λ4686) variations, combined with other measurements, yield an orbital period 2022.7 ± 0.3 d. The observed variability of the EW (He ii λ4686) was reproduced by a model in which the line flux primarily arises at the apex of the windwind collision and scales inversely with the square of the stellar separation, if we account for the excess emission as the companion star plunges into the hot inner layers of the primary's atmosphere, and including absorption from the disturbed primary wind between the source and the observer. This model constrains the orbital inclination to 135 • -153 • , and the longitude of periastron to 234 • -252 • . It also suggests that periastron passage occurred on T 0 = 2456874.4 ± 1.3 d. Our model also reproduced EW (He ii λ4686) variations from a polar view of the primary star as determined from the observed He ii λ4686 emission scattered off the Homunculus nebula.

Monthly Notices of the Royal Astronomical Society, 2016
The X-ray emission of η Carinae shows multiple features at various spatial and temporal scales. T... more The X-ray emission of η Carinae shows multiple features at various spatial and temporal scales. The central constant emission (CCE) component is centred on the binary and arises from spatial scales much smaller than the bipolar Homunculus nebula, but likely larger than the central wind-wind collision region between the stars as it does not vary over the ∼2-3 month X-ray minimum when it can be observed. Using large-scale 3D smoothed particle hydrodynamics (SPH) simulations, we model both the colliding-wind region between the stars, and the region where the secondary wind collides with primary wind ejected from the previous periastron passage. The simulations extend out to one hundred semimajor axes and make two limiting assumptions (strong coupling and no coupling) about the influence of the primary radiation field on the secondary wind. We perform 3D radiative transfer calculations on the SPH output to synthesize the X-ray emission, with the aim of reproducing the CCE spectrum. For the preferred primary mass-loss rate ṀA ≈ 8.5 × 10 -4 M yr -1 , the model spectra well reproduce the observation as the strong-and no-coupling spectra bound the CCE observation for longitude of periastron ω ≈ 252 • , and bound/converge on the observation for ω ≈ 90 • . This suggests that η Carinae has moderate coupling between the primary radiation and secondary wind, that both the region between the stars and the comoving collision on the backside of the secondary generate the CCE, and that the CCE cannot place constraints on the binary's line of sight. We also discuss comparisons with common X-ray fitting parameters.

The Astrophysical Journal, 2015
Eclipsing systems of massive stars allow one to explore the properties of their components in gre... more Eclipsing systems of massive stars allow one to explore the properties of their components in great detail. We perform a multi-wavelength, non-LTE analysis of the three components of the massive multiple system δ Ori A, focusing on the fundamental stellar properties, stellar winds, and X-ray characteristics of the system. The primary's distance-independent parameters turn out to be characteristic for its spectral type (O9.5 II), but usage of the Hipparcos parallax yields surprisingly low values for the mass, radius, and luminosity. Consistent values follow only if δ Ori lies at about twice the Hipparcos distance, in the vicinity of the σ-Orionis cluster. The primary and tertiary dominate the spectrum and leave the secondary only marginally detectable. We estimate the V-band magnitude difference between primary and secondary to be V 2 . 8 m 1 - , which agrees with hydrodynamic predictions, and provides a consistent picture along the X-ray, UV, optical, and radio spectral domains.
Uploads
Papers by christopher russell