Source author record

Christopher M. P. Russell

Christopher M. P. Russell appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2024arXiv

X-ray plasma flow and turbulence in the colliding winds of WR140

We analyse $\textit{XMM-Newton}$ RGS spectra of Wolf-Rayet (WR) 140, an archetype long-period eccentric WR+O colliding wind binary. We evaluate the spectra of O and Fe emission lines and find that the plasmas emitting these lines have the largest approaching velocities with the largest velocity dispersions between phases 0.935 and 0.968 where the inferior conjunction of the O star occurs. This behaviour is the same as that of the Ne line-emission plasma presented in our previous paper. We perform diagnosis of electron number density $n_{\rm e}$ using He-like triplet lines of O and Ne-like Fe-L lines. The former results in a conservative upper limit of $n_{\rm e} \lesssim 10^{10}$-10$^{12}$ cm$^{-3}$ on the O line-emission site, while the latter can not impose any constraint on the Fe line-emission site because of statistical limitations. We calculate the line-of-sight velocity and its dispersion separately along the shock cone. By comparing the observed and calculated line-of-sight velocities, we update the distance of the Ne line-emission site from the stagnation point. By assuming radiative cooling of the Ne line-emission plasma using the observed temperature and the local stellar wind density, we estimate the line-emission site extends along the shock cone by at most $\pm$58 per cent (phase 0.816) of the distance from the stagnation point. In this framework, excess of the observed velocity dispersion over the calculated one is ascribed to turbulence in the hot-shocked plasma at earlier orbital phases of 0.816, 0.912, and 0.935, with the largest velocity dispersion of 340-630 km s$^{-1}$ at phase 0.912.

preprint2022arXiv

Eta Carinae: an evolving view of the central binary, its interacting winds and its foreground ejecta

FUV spectra of Eta Car, recorded across two decades with HST/STIS, document multiple changes in resonant lines caused by dissipating extinction in our line of sight. The FUV flux has increased nearly ten-fold which has led to increased ionization of the multiple shells within the Homunculus and photo-destruction of molecular hydrogen. Comparison of observed resonant line profiles with CMFGEN model profiles allows separation of wind-wind collision and shell absorptions from the primary wind, P Cygni profiles.The dissipating occulter preferentially obscured the central binary and interacting winds relative to the very extended primary wind. We are now able to monitor changes in the colliding winds with orbital phase. High velocity transient absorptions occurred across the most recent periastron passage, indicating acceleration of the primary wind by the secondary wind which leads to a downstream, high velocity bowshock that is newly generated every orbital period. There is no evidence of changes in the properties of the binary winds.

preprint2022arXiv

NICER X-ray Observations of Eta Carinae During its Most Recent Periastron Passage

We report high-precision X-ray monitoring observations in the 0.4-10 keV band of the luminous, long-period colliding-wind binary Eta Carinae up to and through its most recent X-ray minimum/periastron passage in February 2020. Eta Carinae reached its observed maximum X-ray flux on 7 January 2020, at a flux level of $3.30 \times 10^{-10}$ ergs s$^{-1}$ cm$^{-2}$, followed by a rapid plunge to its observed minimum flux, $0.03 \times 10^{-10}$ ergs s$^{-1}$ cm$^{-2}$ near 17 February 2020. The NICER observations show an X-ray recovery from minimum of only $\sim$16 days, the shortest X-ray minimum observed so far. We provide new constraints of the "deep" and "shallow" minimum intervals. Variations in the characteristic X-ray temperature of the hottest observed X-ray emission indicate that the apex of the wind-wind "bow shock" enters the companion's wind acceleration zone about 81 days before the start of the X-ray minimum. There is a step-like increase in column density just before the X-ray minimum, probably associated with the presence of dense clumps near the shock apex. During recovery and after, the column density shows a smooth decline, which agrees with previous $N_{H}$ measurements made by SWIFT at the same orbital phase, indicating that changes in mass-loss rate are only a few percent over the two cycles. Finally, we use the variations in the X-ray flux of the outer ejecta seen by NICER to derive a kinetic X-ray luminosity of the ejecta of $\sim 10^{41}$ ergs s$^{-1}$ near the time of the "Great Eruption'.

preprint2021arXiv

Conditions in the WR 140 wind-collision region revealed by the 1.083-micron He I line profile

We present spectroscopy of the P~Cygni profile of the 1.083-micron He I line in the WC7 + O5 colliding-wind binary (CWB) WR 140 (HD 193793), observed in 2008, before its periastron passage in 2009, and in 2016-17, spanning the subsequent periastron passage. Both absorption and emission components showed strong variations. The variation of the absorption component as the O5 star was occulted by the wind-collision region (WCR) sets a tight constraint on its geometry. While the sightline to the O5 star traversed the WCR, the strength and breadth of the absorption component varied significantly on time-scales of days. An emission sub-peak was observed on all our profiles. The variation of its radial velocity with orbital phase was shown to be consistent with formation in the WCR as it swung round the stars in their orbit. Modelling the profile gave a measure of the extent of the sub-peak forming region. In the phase range 0.93-0.99, the flux in the sub-peak increased steadily, approximately inversely proportionally to the stellar separation, indicating that the shocked gas in the WCR where the line was formed was adiabatic. After periastron, the sub-peak flux was anomalously strong and varied rapidly, suggesting formation in clumps down-stream in the WCR. For most of the time, its flux exceeded the 2-10-keV X-ray emission, showing it to be a significant coolant of the shocked wind.

preprint2020arXiv

3D simulations of clump formation in stellar wind collisions

The inner parsec of our Galaxy contains tens of Wolf-Rayet stars whose powerful outflows are constantly interacting while filling the region with hot, diffuse plasma. Theoretical models have shown that, in some cases, the collision of stellar winds can generate cold, dense material in the form of clumps. However, their formation process and properties are not well understood yet. In this work we present, for the first time, a statistical study of the clump formation process in unstable wind collisions. We study systems with dense outflows (${\sim}10^{-5}\rm\ M_{\odot}\ yr^{-1}$), wind speeds of $500$-$1500\rm\ km\ s^{-1}$, and stellar separations of ${\sim}20$-$200\rm\ au$. We develop 3D high resolution hydrodynamical simulations of stellar wind collisions with the adaptive-mesh refinement grid-based code Ramses. We aim to characterise the initial properties of clumps that form through hydrodynamic instabilities, mostly via the non-linear thin shell instability (NTSI). Our results confirm that more massive clumps are formed in systems whose winds are close to the transition between the radiative and adiabatic regimes. Increasing either the wind speed or the degree of asymmetry increases the dispersion of the clump mass and ejection speed distributions. Nevertheless, the most massive clumps are very light (${\sim}10^{-3}$-$10^{-2}\rm\ M_{\oplus}$), about three orders of magnitude less massive than theoretical upper limits. Applying these results to the Galactic Centre we find that clumps formed through the NTSI should not be heavy enough either to affect the thermodynamic state of the region or to survive for long enough to fall onto the central super-massive black hole.

preprint2019arXiv

Colliding Winds in and around the Stellar Group IRS 13E at the Galactic Center

IRS~13E is an enigmatic compact group of massive stars located in projection only 3.6 arcseconds away from Sgr A*. This group has been suggested to be bounded by an intermediate-mass black hole (IMBH). We present a multi-wavelength study of the group and its interplay with the environment. Based on Chandra observations, we find the X-ray spectrum of IRS~13E can be well characterized by an optically thin thermal plasma. The emission peaks between two strongly mass-losing Wolf-Rayet stars of the group. These properties can be reasonably well reproduced by simulated colliding winds of these two stars. However, this scenario under-predicts the X-ray intensity in outer regions. The residual emission likely results from the ram-pressure confinement of the IRS~13E group wind by the ambient medium and is apparently associated with a shell-like warm gas structure seen in Pa-alpha and in ALMA observations. These latter observations also show strongly peaked thermal emission with unusually large velocity spread between the two stars. These results indicate that the group is colliding with the bar of the dense cool gas mini-spiral around Sgr A*. The extended X-ray morphology of IRS~13E and its association with the bar further suggest that the group is physically much farther away than the projected distance from Sgr A*. The presence of an IMBH, while favorable to keep the stars bound together, is not necessary to explain the observed stellar and gas properties of IRS~13E.

preprint2019arXiv

Stellar winds pump the heart of the Milky Way

The central super-massive black hole of the Milky Way, Sgr A*, accretes at a very low rate making it a very underluminous galactic nucleus. Despite the tens of Wolf-Rayet stars present within the inner parsec supplying ${\sim}10^{-3}\rm\ M_{\odot}\ yr^{-1}$ in stellar winds, only a negligible fraction of this material ($<10^{-4}$) ends up being accreted onto Sgr A*. The recent discovery of cold gas (${\sim}10^4\rm\ K$) in its vicinity raised questions about how such material could settle in the hostile (${\sim}10^7\rm\ K$) environment near Sgr A*. In this work we show that the system of mass-losing stars blowing winds can naturally account for both the hot, inefficient accretion flow, as well as the formation of a cold disk-like structure. We run hydrodynamical simulations using the grid-based code Ramses starting as early in the past as possible to observe the state of the system at the present time. Our results show that the system reaches a quasi-steady state in about ${\sim}500\rm\ yr$ with material being captured at a rate of ${\sim}10^{-6}\rm\ M_{\odot}\ yr^{-1}$ at scales of ${\sim}10^{-4}\rm\ pc$, consistent with the observations and previous models. However, on longer timescales ($\gtrsim3000\rm\ yr$) the material accumulates close to the black hole in the form of a disk. Considering the duration of the Wolf-Rayet phase (${\sim}10^5\rm\ yr$), we conclude that this scenario likely has already happened, and could be responsible for the more active past of Sgr A*, and/or its current outflow. We argue that the hypothesis of the mass-losing stars being the main regulator of the activity of the black hole deserves further consideration.

preprint2016arXiv

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

The evolved, massive highly eccentric binary system, eta Carinae, 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 eta 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 Delta 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 star's 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 eta Car observations. The column density to the hardest emission component, NH ~1e24 cm-2, marked one of the highest values ever observed for eta Car, strongly suggesting 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. If the non-detection by NuSTAR is caused by absorption, the power-law source must be small and located very near the WWC apex. Alternatively, it may be that the power-law source is not related to either eta Car or the GeV gamma-ray source.

preprint2016arXiv

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

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 $\sim$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 $\dot{M}_A\approx8.5\times10^{-4}$ M$_\odot$ yr$^{-1}$, the model spectra well reproduce the observation as the strong- and no-coupling spectra bound the CCE observation for longitude of periastron $ω\approx252^\circ$, and bound/converge on the observation for $ω\approx90^\circ$. 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.

preprint2016arXiv

To $v_\infty$ and beyond! The He I absorption variability across the 2014.6 periastron passage of $η$ Carinae

We have monitored the massive binary star $η$ Carinae with the CTIO/SMARTS 1.5~m telescope and CHIRON spectrograph from the previous apastron passage of the system through the recent 2014.6 periastron passage. Our monitoring has resulted in a large, homogeneous data set with an unprecedented time-sampling, spectral resolving power, and signal-to-noise. This allowed us to investigate temporal variability previously unexplored in the system and discover a kinematic structure in the P Cygni absorption troughs of neutral helium wind lines. The features observed occurred prior to the periastron passage and are seen as we look through the trailing arm of the wind-wind collision shock cone. We show that the bulk of the variability is repeatable across the last five periastron passages, and that the absorption occurs in the inner 230 AU of the system. In addition, we found an additional, high-velocity absorption component super-imposed on the P Cygni absorption troughs that has been previously un-observed in these lines, but which bears resemblance to the observations of the He~I $λ$10830 Å feature across previous cycles. Through a comparison of the current smoothed particle hydrodynamical simulations, we show that the observed variations are likely caused by instabilities in the wind-wind collision region in our line-of-sight, coupled with stochastic variability related to clumping in the winds.

preprint2015arXiv

Hydrodynamic and radiative transfer modeling of X-ray emission from colliding WR winds: WR 140 & the Galactic center

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.

preprint2014arXiv

3D Dynamical Modeling of Wind Accretion in Cyg X-3

Cyg X-3 is a high mass X-ray binary consisting of a Wolf-Rayet star and a compact object in a very short orbital period of 4.8h. The only confirmed microquasar with high energy gamma-ray emission, Cyg X-3 provides a unique opportunity to study the relationship between the accretion power and the power in high energy emission. Because of a compact orbit and a slow Wolf-Rayet wind, the flow structure around the compact object is thought to be strongly affected by the orbital motion, details of which can be obtained only by numerical simulations. In this paper, we report on the results from 3D hydrodynamic simulations of the wind accretion in Cyg X-3. For simplicity we adopt an anti-gravity-like force that emulates the radiative acceleration consistent with the beta-velocity wind. Due to the rapid orbital motion, the flow around the compact object has large density gradients. As a result, the accretion rate onto the compact object is significantly lower than that of the Bondi-Hoyle-Lyttleton rate. We also calculate the model X-ray light curve. Although it roughly agrees with the observed light curve, more detailed modeling is needed for detailed comparison.

preprint2014arXiv

3D Hydrodynamic & Radiative Transfer Models of X-ray Emission from Colliding Wind Binaries

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 $\sim$3000 km s$^{-1}$, and their collision leads to hot plasma (up to $\sim10^8$K) that emit thermal X-rays (up to $\sim$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.

preprint2014arXiv

Suzaku Monitoring of Hard X-ray Emission from Eta 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.

preprint2014arXiv

V444 Cyg 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 XMM-Newton observations. In addition, we have characterized the intrinsic, 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 opening 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.

preprint2013arXiv

Constraints on decreases in Eta 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 results from large- (r=1545au) and small- (r=155au) domain, 3D smoothed particle hydrodynamic (SPH) simulations of Eta Car's colliding winds for 3 Eta A mass-loss rates (2.4, 4.8, and 8.5 x 10^-4 M_sun/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 Eta A's mass-loss rate 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. The SPH simulations together with 1D radiative transfer models of Eta A's spectra reveal that a factor of 2 or more drop in Eta A's mass-loss rate should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in Eta A's mass-loss rate was likely by a factor < 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 Eta A's mass-loss rate may be responsible, but changes in the wind/stellar parameters of Eta B cannot yet be fully ruled out. We suggest observations during Eta Car's next periastron in 2014 to further test for decreases in Eta A's mass-loss rate. If Eta A's mass-loss rate is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.

preprint2011arXiv

Constraining the Absolute Orientation of Eta Carinae's Binary Orbit: A 3-D 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 Car's 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 θ 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 = +38 degrees, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i \approx 130 to 145 degrees, θ \approx -15 to +30 degrees, and an orbital axis projected on the sky at a PA \approx 302 to 327 degrees 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 B, thus orbits clockwise on the sky and is on the observer's 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.

preprint2011arXiv

X-ray Modeling of η Carinae and WR140 from SPH Simulations

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 X-ray light curves and spectra, we apply 3D hydrodynamic simulations of the wind-wind collision using smoothed particle hydrodynamics (SPH), with the recent improvements of radiative cooling and the acceleration of the stellar winds according to a β law. For both systems, the 2-10 keV RXTE light curves are well-reproduced in absolute units for most phases, but the light curve dips associated with the periastron passages are not well matched. In WR140, the dip is too weak, and in η Carinae, the large difference in wind speeds of the two stars leads to a hot, post-periastron bubble that produces excess emission toward the end of the X-ray minimum.

preprint2011arXiv

X-ray Modeling of η Carinae and WR140 from SPH Simulations

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.