In circumstellar gas, the complex organic molecule methanol has been found almost exclusively around young stellar objects, and is thus regarded as a signpost of recent star formation. Here we report the first probable detection of methanol around an evolved high-mass star, in the complex circumstellar environment around the Luminous Blue Variable $\eta$ Carinae, while using ALMA to investigate molecular cloud conditions traced by CO (2-1) in an orbit phase of the massive binary preceding the 2020 periastron. Favoring methanol over a $^{13}$CS alternative, the emission originates from hot ($T_{\rm{gas}} \simeq$ 700 K) material, $\sim$2$''$ (0.02 pc) across, centered on the dust-obscured binary in contrast to the CO which traces inner layers of the extended massive equatorial torus, and is accompanied by prominent absorption in a cooler ($T_{\rm{gas}} \simeq$ 110 K) layer of gas. We also report detections of water in $Herschel$/HIFI observations at 557 GHz and 988 GHz. The methanol abundance is several to 50 times higher than observed towards several lower mass stars, while water abundances are similar to those observed in cool, dense molecular clouds. The very high methanol:water abundance ratio in the core of $\eta$ Carinae may suggest methanol formation processes similar to Fischer-Tropsch-type catalytic reactions on dust grains. These observations prove that complex molecule formation can occur in the chemically evolved environments around massive stars in the end stages of their evolution, given sufficient gas densities and shielding conditions as may occur in material around massive interacting companions and merger remnants.

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 gamma-rays. The supermassive binary star eta Carinae drives the strongest colliding wind shock in the solar neighborhood. Observations with non-focusing high-energy observatories indicate a high energy source near eta Carinae, but have been unable to conclusively identify eta Carinae as the source because of their relatively poor angular resolution. Here we present the first direct focussing observations of the non-thermal source in the extremely hard X-ray band, which is found to be spatially coincident with the star within several arc-seconds. These observations show that the source of non-thermal X-rays varies with the orbital phase of the binary, and that the photon index of the emission is similar to that derived through analysis of the gamma-ray spectrum. This is conclusive evidence that the high-energy emission indeed originates from non-thermal particles accelerated at colliding wind shocks.

The interacting binary Eta Carinae remains one of the most enigmatic massive stars in our Galaxy despite over four centuries of observations. In this work, its light curve from the ultraviolet to the near-infrared is analysed using spatially resolved HST observations and intense monitoring at the La Plata Observatory, combined with previously published photometry. We have developed a method to separate the central stellar object in the ground-based images using HST photometry and applying it to the more numerous ground-based data, which supports the hypothesis that the central source is brightening faster than the almost-constant Homunculus. After detrending from long-term brightening, the light curve shows periodic orbital modulation ($\Delta V$ $\sim$ 0.6 mag) attributed to the wind-wind collision cavity as it sweeps around the primary star and it shows variable projected area to our line-of-sight. Two quasi-periodic components with time scales of 2-3 and 8-10 yr and low amplitude, $\Delta V$ $<$ 0.2 mag, are superimposed on the brightening light-curve, being the only stellar component of variability found, which indicates minimal stellar instability. Moreover, the light curve analysis shows no evidence of `shell ejections' at periastron. We propose that the long-term brightening of the stellar core is due to the dissipation of a dusty clump in front of the central star, which works like a natural coronagraph. Thus, the central stars appear to be more stable than previously thought since the dominant variability originates from a changing circumstellar medium. We predict that the brightening phase, due mainly to dust dissipation, will be completed around 2032 $\pm$ 4 yr, when the star will be brighter than in the 1600's by up to $\Delta V$ $\sim$ 1 mag.

$\eta$ Car is a massive, eccentric binary with a rich observational history. We obtained the first high-cadence, high-precision light curves with the BRITE-Constellation nanosatellites over 6 months in 2016 and 6 months in 2017. The light curve is contaminated by several sources including the Homunculus nebula and neighboring stars, including the eclipsing binary CPD$-$59$^\circ$2628. However, we found two coherent oscillations in the light curve. These may represent pulsations that are not yet understood but we postulate that they are related to tidally excited oscillations of $\eta$ Car's primary star, and would be similar to those detected in lower-mass eccentric binaries. In particular, one frequency was previously detected by van Genderen et al. and Sterken et al. through the time period of 1974 to 1995 through timing measurements of photometric maxima. Thus, this frequency seems to have been detected for nearly four decades, indicating that it has been stable in frequency over this time span. These pulsations could help provide the first direct constraints on the fundamental parameters of the primary star if confirmed and refined with future observations.

3-D printing moves beyond interactive 3-D graphics and provides an excellent tool for both visual and tactile learners, since 3-D printing can now easily communicate complex geometries and full color information. Some limitations of interactive 3-D graphics are also alleviated by 3-D printable models, including issues of limited software support, portability, accessibility, and sustainability. We describe the motivations, methods, and results of our work on using 3-D printing (1) to visualize and understand the Eta Car Homunculus nebula and central binary system and (2) for astronomy outreach and education, specifically, with visually impaired students. One new result we present is the ability to 3-D print full-color models of Eta Car's colliding stellar winds. We also demonstrate how 3-D printing has helped us communicate our improved understanding of the detailed structure of Eta Car's Homunculus nebula and central binary colliding stellar winds, and their links to each other. Attached to this article are full-color 3-D printable files of both a red-blue Homunculus model and the Eta Car colliding stellar winds at orbital phase 1.045. 3-D printing could prove to be vital to how astronomer's reach out and share their work with each other, the public, and new audiences.

Context. The mass loss from massive stars is not understood well. Eta Car is a unique object for studying the massive stellar wind during the LBV phase. It is also an eccentric binary with a period of 5.54 yr. The nature of both stars is uncertain, although we know from X-ray studies that there is a wind-wind collision whose properties change with orbital phase. Methods. Observations of Eta Car were carried out with the ESO VLTI and the AMBER instrument between approximately five and seven months before the August 2014 periastron passage. Velocity-resolved aperture-synthesis images were reconstructed from the spectrally dispersed interferograms. Interferometric studies can provide information on the binary orbit, the primary wind, and the wind collision. Results. We present velocity-resolved aperture-synthesis images reconstructed in more than 100 different spectral channels distributed across the Br Gamma 2.166 micrometer emission line. The intensity distribution of the images strongly depends on wavelength. At wavelengths corresponding to radial velocities of approximately -140 to -376 km/s measured relative to line center, the intensity distribution has a fan-shaped structure. At the velocity of -277 km/s, the position angle of the symmetry axis of the fan is ~ 126 degree. The fan-shaped structure extends approximately 8.0 mas (~ 18.8 au) to the southeast and 5.8 mas (~ 13.6 au) to the northwest, measured along the symmetry axis at the 16% intensity contour. The shape of the intensity distributions suggests that the obtained images are the first direct images of the innermost wind-wind collision zone. Therefore, the observations provide velocity-dependent image structures that can be used to test three-dimensional hydrodynamical, radiative transfer models of the massive interacting winds of Eta Car.

Eta Carinae, the closest, active, massive binary containing a highly unstable LBV, exhibits expanding, compressed wind shells, seen in emission, that are spatially and spectrally resolved by \hst/\stis. Starting in June 2009, these structures were mapped across its 5.54-year, highly elliptical, binary orbit to follow temporal changes in the light of [Fe~III] 4659\AA\ and [Fe~II] 4815\AA. The emissions trace portions of fossil wind-shells, that were formed by wind-wind interactions across each cycle. Over the high-ionization state, dense arcs, photo-ionized by far ultraviolet radiation from the hot secondary, are seen in [Fe~III]. Other arcs, ionized by mid-ultraviolet radiation from the primary star, are seen in [Fe II]. The [Fe III] structures tend to be interior to [Fe II] structures that trace extensive, less disturbed primary wind. During the brief periastron passage when the secondary plunges deep into the primary's extremely dense wind, on the far side of primary star, high-ionization [Fe III] structures fade and reappear in [Fe II]. Multiple fossil wind-structures were traced across the 5.7-year monitoring interval. The strong similarity of the expanding [Fe II] shells suggests that the wind and photo-ionization properties of the massive binary have not changed substantially from one orbit to the next over the past several orbital cycles. These observations trace structures that can be used to test three-dimensional hydrodynamical and radiative-transfer models of massive, interacting winds. They also provide a baseline for following future changes in Eta Carinae, especially of its winds and photoionization properties.

We have monitored the massive binary star $\eta$ 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 $\lambda$10830 \AA\ 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.

The X-ray emission of $\eta$ 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 $\omega\approx252^\circ$, and bound/converge on the observation for $\omega\approx90^\circ$. This suggests that $\eta$ 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 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.

{\eta} Carinae is an extremely massive binary system in which rapid spectrum variations occur near periastron. Most notably, near periastron the He II $\lambda 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 $\lambda 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 $\lambda 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 $\lambda 4686$) variations, combined with other measurements, yield an orbital period $2022.7\pm0.3$ d. The observed variability of the EW(He II $\lambda 4686$) was reproduced by a model in which the line flux primarily arises at the apex of the wind-wind 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^\circ$-$153^\circ$, and the longitude of periastron to $234^\circ$-$252^\circ$. It also suggests that periastron passage occurred on $T_0 = 2456874.4\pm1.3$ d. Our model also reproduced EW(He II $\lambda 4686$) variations from a polar view of the primary star as determined from the observed He II $\lambda 4686$ emission scattered off the Homunculus nebula.

We present results of full 3D hydrodynamical and radiative transfer simulations of the colliding stellar winds in the massive binary system Eta Carinae. We accomplish this by applying the SimpleX algorithm for 3D radiative transfer on an unstructured Voronoi-Delaunay grid to recent 3D smoothed particle hydrodynamics (SPH) simulations of the binary colliding winds. We use SimpleX to obtain detailed ionization fractions of hydrogen and helium, in 3D, at the resolution of the original SPH simulations. We investigate several computational domain sizes and Luminous Blue Variable primary star mass-loss rates. We furthermore present new methods of visualizing and interacting with output from complex 3D numerical simulations, including 3D interactive graphics and 3D printing. While we initially focus on Eta Car, the methods employed can be applied to numerous other colliding wind (WR 140, WR 137, WR 19) and dusty 'pinwheel' (WR 104, WR 98a) binary systems. Coupled with 3D hydrodynamical simulations, SimpleX simulations have the potential to help determine the regions where various observed time-variable emission and absorption lines form in these unique objects.

Eta Carinae 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 eta 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.

Spectral observations of the massive colliding wind binary Eta Carinae show phase-dependent variations, in intensity and velocity, of numerous helium emission and absorption lines throughout the entire 5.54-year orbit. Approaching periastron, the 3D structure of the wind-wind interaction region (WWIR) gets highly distorted due to the eccentric ($e \sim 0.9$) binary orbit. The secondary star ($\eta_{\mathrm{B}}$) at these phases is located deep within the primary's dense wind photosphere. The combination of these effects is thought to be the cause of the particularly interesting features observed in the helium lines at periastron. We perform 3D radiative transfer simulations of $\eta$ Car's interacting winds at periastron. Using the SimpleX radiative transfer algorithm, we post-process output from 3D smoothed particle hydrodynamic simulations of the inner 150 au of the $\eta$ Car system for two different primary star mass-loss rates ($\dot{M}_{\eta_{\mathrm{A}}}$). Using previous results from simulations at apastron as a guide for the initial conditions, we compute 3D helium ionization maps. We find that, for higher $\dot{M}_{\eta_{\mathrm{A}}}$, $\eta_{\mathrm{B}}$ He$^{0+}$-ionizing photons are not able to penetrate into the pre-shock primary wind. He$^{+}$ due to $\eta_{\mathrm{B}}$ is only present in a thin layer along the leading arm of the WWIR and in a small region close to the stars. Lowering $\dot{M}_{\eta_{\mathrm{A}}}$ allows $\eta_{\mathrm{B}}$'s ionizing photons to reach the expanding unshocked secondary wind on the apastron side of the system, and create a low fraction of He$^{+}$ in the pre-shock primary wind. With apastron on our side of the system, our results are qualitatively consistent with the observed variations in strength and radial velocity of $\eta$ Car's helium emission and absorption lines, which helps better constrain the regions where these lines arise.

We present the first 3D prints of output from a supercomputer simulation of a complex astrophysical system, the colliding stellar winds in the massive (>120 M_Sun), highly eccentric (e ~ 0.9) binary star system Eta Carinae. We demonstrate the methodology used to incorporate 3D interactive figures into a PDF journal publication and the benefits of using 3D visualization and 3D printing as tools to analyze data from multidimensional numerical simulations. Using a consumer-grade 3D printer (MakerBot Replicator 2X), we successfully printed 3D smoothed particle hydrodynamics (SPH) simulations of Eta Carinae's inner (r ~ 110 au) wind-wind collision interface at multiple orbital phases. The 3D prints and visualizations reveal important, previously unknown 'finger-like' structures at orbital phases shortly after periastron (phi ~ 1.045) that protrude radially outward from the spiral wind-wind collision region. We speculate that these fingers are related to instabilities (e.g. thin-shell, Rayleigh-Taylor) that arise at the interface between the radiatively-cooled layer of dense post-shock primary-star wind and the fast (3000 km/s), adiabatic post-shock companion-star wind. The success of our work and easy identification of previously unrecognized physical features highlight the important role 3D printing and interactive graphics can play in the visualization and understanding of complex 3D time-dependent numerical simulations of astrophysical phenomena.

The highly eccentric binary system Eta Carinae shows numerous time-variable emission and absorption features. These observational signatures are the result of interactions between the complex three-dimensional (3D) wind-wind collision regions and photoionization by the luminous stars. Specifically, helium presents several interesting spectral features that provide important clues on the geometry and physical properties of the system and the individual stars. We use the SimpleX algorithm to post-process 3D smoothed particle hydrodynamics simulation output of the interacting winds in Eta Car in order to obtain the fractions of ionized helium assuming three different primary star mass-loss rates. The resultant ionization maps constrain the regions where helium is singly- and doubly-ionized. We find that reducing the primary's mass-loss rate increases the volume of He+. Lowering the primary mass-loss rate produces large variations in the volume of He+ in the pre-shock primary wind on the periastron side of the system. Our results show that binary orientations in which apastron is on our side of the system are more consistent with available observations. We suggest that small variations in the primary's mass-loss rate might explain the observed increase in HeI absorption in recent decades, although numerous questions regarding this scenario remain open. We also propose that the absence of broad HeI lines in the spectra of Eta Car between its 1890's eruption and 1944 might be explained by the companion star's He0+ ionizing photons not being able to penetrate the wind-wind interaction region, due to a higher primary mass-loss rate at that time (by a factor >2, compared to the present value).

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.

We investigate, using the modeling code SHAPE, the three-dimensional structure of the bipolar Homunculus nebula surrounding Eta Carinae, as mapped by new ESO VLT/X-Shooter observations of the H2 $\lambda=2.12125$ micron emission line. Our results reveal for the first time important deviations from the axisymmetric bipolar morphology: 1) circumpolar trenches in each lobe positioned point-symmetrically from the center and 2) off-planar protrusions in the equatorial region from each lobe at longitudinal (~55 degrees) and latitudinal (10-20 degrees) distances from the projected apastron direction of the binary orbit. The angular distance between the protrusions (~110 degrees) is similar to the angular extent of each polar trench (~130 degrees) and nearly equal to the opening angle of the wind-wind collision cavity (~110 degrees). As in previous studies, we confirm a hole near the centre of each polar lobe and no detectable near-IR H2 emission from the thin optical skirt seen prominently in visible imagery. We conclude that the interaction between the outflows and/or radiation from the central binary stars and their orientation in space has had, and possibly still has, a strong influence on the Homunculus. This implies that prevailing theoretical models of the Homunculus are incomplete as most assume a single star origin that produces an axisymmetric nebula. We discuss how the newly found features might be related to the Homunculus ejection, the central binary and the interacting stellar winds. We also include a 3D printable version of our Homunculus model.

Eta Carinae is an ideal astrophysical laboratory for studying massive binary interactions and evolution, and stellar wind-wind collisions. Recent three-dimensional (3D) simulations set the stage for understanding the highly complex 3D flows in $\eta$ Car. Observations of different broad high- and low-ionization forbidden emission lines provide an excellent tool to constrain the orientation of the system, the primary's mass-loss rate, and the ionizing flux of the hot secondary. In this work we present the first steps towards generating synthetic observations to compare with available and future HST/STIS data. We present initial results from full 3D radiative transfer simulations of the interacting winds in $\eta$ Car. We use the SimpleX algorithm to post-process the output from 3D SPH simulations and obtain the ionization fractions of hydrogen and helium assuming three different mass-loss rates for the primary star. The resultant ionization maps of both species constrain the regions where the observed forbidden emission lines can form. Including collisional ionization is necessary to achieve a better description of the ionization states, especially in the areas shielded from the secondary's radiation. We find that reducing the primary's mass-loss rate increases the volume of ionized gas, creating larger areas where the forbidden emission lines can form. We conclude that post processing 3D SPH data with SimpleX is a viable tool to create ionization maps for $\eta$ Car.

X-ray emission from the supermassive binary system Eta Carinae 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.9e-12 ergs 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 emissivity of the WWC shock.

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.

A series of three HST/STIS spectroscopic mappings, spaced approximately one year apart, reveal three partial arcs in [Fe II] and [Ni II] emissions moving outward from eta Carinae. We identify these arcs with the shell-like structures, seen in the 3D hydrodynamical simulations, formed by compression of the primary wind by the secondary wind during periastron passages.

We analyze spatially resolved spectroscopic observations of the Eta Carinae binary system obtained with HST/STIS. Eta Car is enshrouded by the dusty Homunculus nebula, which scatters light emitted by the central binary and provides a unique opportunity to study a massive binary system from different vantage points. We investigate the latitudinal and azimuthal dependence of H$\alpha$ line profiles caused by the presence of a wind-wind collision (WWC) cavity created by the companion star. Using two-dimensional radiative transfer models, we find that the wind cavity can qualitatively explain the observed line profiles around apastron. Regions of the Homunculus which scatter light that propagated through the WWC cavity show weaker or no H alpha absorption. Regions scattering light that propagated through a significant portion of the primary wind show stronger P Cygni absorption. Our models overestimate the H alpha absorption formed in the primary wind, which we attribute to photoionization by the companion, not presently included in the models. We can qualitatively explain the latitudinal changes that occur during periastron, shedding light on the nature of Eta Car's spectroscopic events. Our models support the idea that during the brief period of time around periastron when the primary wind flows unimpeded toward the observer, H alpha absorption occurs in directions toward the central object and Homunculus SE pole, but not toward equatorial regions close to the Weigelt blobs. We suggest that observed latitudinal and azimuthal variations are dominated by the companion star via the WWC cavity, rather than by rapid rotation of the primary star.

We present 2D radiative transfer modeling of the Eta Carinae binary system accounting for the presence of a wind-wind collision (WWC) cavity carved in the optically-thick wind of the primary star. By comparing synthetic line profiles with HST/STIS spectra obtained near apastron, we show that the WWC cavity has a strong influence on multi-wavelength diagnostics. This influence is regulated by the modification of the optical depth in the continuum and spectral lines. We find that H-alpha, H-beta, and Fe II lines are the most affected by the WWC cavity, since they form over a large volume of the primary wind. These spectral lines depend on latitude and azimuth since, according to the orientation of the cavity, different velocity regions of a spectral line are affected. For 2D models with orientation corresponding to orbital inclination angle 110deg < i < 140deg and longitude of periastron 210deg < omega < 330deg, the blueshifted and zero-velocity regions of the line profiles are the most affected. These orbital orientations are required to simultaneously fit the UV and optical spectrum of Eta Car, for a half-opening angle of the cavity in the range 50-70deg. We find that the excess P-Cygni absorption seen in H-alpha, H-beta and optical Fe II lines in spherical models becomes much weaker or absent in the 2D models, in agreement with the observations. The observed UV spectrum of Eta Car, dominated by Fe II absorption lines, is superbly reproduced by our 2D cavity models. Small discrepancies still remain, as H-gamma and H-delta absorptions are overestimated by our models. We suggest that photoionization of the wind of the primary by the hot companion star is responsible for the weak absorption seen in these lines. Our CMFGEN models indicate that the primary star has a mass-loss rate of 8.5x10e-4 Msun/yr and wind terminal velocity of 420 km/s around the 2000 apastron.

We present a new model for the behavior of scattered time-dependent, asymmetric near-UV emission from the nearby ejecta of {\eta} Car. Using a 3-D hydrodynamical simulation of {\eta} Car's binary colliding winds, we show that the 3-D binary orientation derived by Madura et al. (2012) is capable of explaining the asymmetric near-UV variability observed in the Hubble Space Telescope Advanced Camera for Surveys/High Resolution Camera (HST ACS/HRC) F220W images of Smith et al. (2004b). Models assuming a binary orientation with i ~ 130 to 145 degrees, {\omega} ~ 230 to 315 degrees, PAz ~ 302 to 327 degrees are consistent with the observed F220W near-UV images. We find that the hot binary companion does not significantly contribute to the near-UV excess observed in the F220W images. Rather, we suggest that a bore-hole effect and the reduction of Fe II optical depths inside the wind-wind collision cavity carved in the extended photosphere of the primary star lead to the time-dependent directional illumination of circum-binary material as the companion moves about in its highly elliptical orbit.

Eta Carinae is considered to be a massive colliding wind binary system with a highly eccentric (e \sim 0.9), 5.54-yr orbit. However, the companion star continues to evade direct detection as the primary dwarfs its emission at most wavelengths. Using three-dimensional (3-D) SPH simulations of Eta Car's colliding winds and radiative transfer codes, we are able to compute synthetic observables across multiple wavebands for comparison to the observations. The models show that the presence of a companion star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line profiles, as well as the ground-based photometric monitoring. Here, we focus on the Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in the dense primary wind, allowing increased escape of radiation from the hotter/deeper layers of the primary's extended wind photosphere. The results have important implications for interpretations of Eta Car's observables at multiple wavelengths.

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 {\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 = +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, {\theta} \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.

We report new HST/STIS observations that map the high-ionization forbidden line emission in the inner arcsecond of Eta Car, the first that fully image the extended wind-wind interaction region of the massive colliding wind binary. These observations were obtained after the 2009.0 periastron at orbital phases 0.084, 0.163, and 0.323 of the 5.54-year spectroscopic cycle. We analyze the variations in brightness and morphology of the emission, and find that blue-shifted emission (-400 to -200 km s-1) is symmetric and elongated along the northeast-southwest axis, while the red-shifted emission (+100 to +200 km s-1) is asymmetric and extends to the north- northwest. Comparison to synthetic images generated from a 3-D dynamical model strengthens the 3-D orbital orientation found by Madura et al. (2011), with an inclination i \approx 138\degree, argument of periapsis {\omega} \approx 270\degree, and an orbital axis that is aligned at the same PA on the sky as the symmetry axis of the Homunculus, 312\degree. We discuss the potential that these and future mappings have for constraining the stellar parameters of the companion star and the long-term variability of the system.

The colliding wind binary (CWB) systems \eta\ 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 \beta\ 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 \eta\ 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.

The colliding wind binary (CWB) systems \eta\ 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 \eta\ Car and WR140. This point-source model can also explain the early recovery of \eta\ Car's X-ray light curve from the 2009.0 minimum by a factor of 2-4 reduction in the mass loss rate of \eta\ 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 \eta\ 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.

We analyze interferometric measurements of the Luminous Blue Variable Eta Carinae with the goal of constraining the rotational velocity of the primary star and probing the influence of the companion. Using 2-D radiative transfer models of latitude-dependent stellar winds, we find that prolate wind models with a ratio of the rotational velocity (vrot) to the critical velocity (vcrit) of W=0.77-0.92, inclination angle of i=60-90 degrees, and position angle PA=108-142 degrees reproduce simultaneously K-band continuum visibilities from VLTI/VINCI and closure phase measurements from VLTI/AMBER. Interestingly, oblate models with W=0.73-0.90 and i=80-90 degrees produce similar fits to the interferometric data, but require PA=210-230 degrees. Therefore, both prolate and oblate models suggest that the rotation axis of the primary star is not aligned with the Homunculus polar axis. We also compute radiative transfer models of the primary star allowing for the presence of a cavity and dense wind-wind interaction region created by the companion star. We find that the wind-wind interaction has a significant effect on the K-band image mainly via free-free emission from the compressed walls and, for reasonable model parameters, can reproduce the VLTI/VINCI visibilities taken at phase 0.92-0.93. We conclude that the density structure of the primary wind can be sufficiently disturbed by the companion, thus mimicking the effects of fast rotation in the interferometric observables. Therefore, fast rotation may not be the only explanation for the interferometric observations. Intense temporal monitoring and 3-D modeling are needed to resolve these issues.

We report near-IR spectroscopic observations of the Eta Carinae massive binary system during 2008-2009 using VLT/CRIRES. We detect a strong, broad absorption wing in He I 10833 extending up to -1900 km/s across the 2009.0 spectroscopic event. Archival HST/STIS ultraviolet and optical data shows a similar high-velocity absorption (up to -2100 km/s) in the UV resonance lines of Si IV 1394, 1403 across the 2003.5 event. UV lines from low-ionization species, such as Si II 1527, 1533 and C II 1334, 1335, show absorption up to -1200 km/s, indicating that the absorption with v from -1200 to -2100 km/s originates in a region markedly faster and more ionized than the nominal wind of the primary star. Observations obtained at the OPD/LNA during the last 4 spectroscopic cycles (1989-2009) also display high-velocity absorption in He I 10833 during periastron. Based on the OPD/LNA dataset, we determine that material with v < -900 km/s is present in the phase range 0.976 < phi < 1.023 of the spectroscopic cycle, but absent in spectra taken at phi < 0.947 and phi > 1.049. Therefore, we constrain the duration of the high-velocity absorption to be 95 to 206 days (or 0.047 to 0.102 in phase). We suggest that the high-velocity absorption originates from shocked gas in the wind-wind collision zone, at distances of 15 to 45 AU in the line-of-sight to the primary star. Using 3-D hydrodynamical simulations of the wind-wind collision zone, we find that the dense high-velocity gas is in the line-of-sight to the primary star only if the binary system is oriented in the sky so that the companion is behind the primary star during periastron, corresponding to a longitude of periastron of omega ~ 240 to 270 degrees. We study a possible tilt of the orbital plane relative to the Homunculus equatorial plane and conclude that our data are broadly consistent with orbital inclinations in the range i=40 to 60 degrees.

We analyze the steady 1D flow equations for a rotating stellar wind based on a ``nozzle'' analogy for terms that constrain the local mass flux. For low rotation, we find the nozzle minimum occurs near the stellar surface, allowing a transition to a standard, CAK-type steep-acceleration solution; but for rotations > 75% of the critical rate, this inner nozzle minimum exceeds the global minimum, implying near-surface supercritical solutions would have an overloaded mass loss rate. In steady, analytic models in which the acceleration is assumed to be monotonically positive, this leads the solution to switch to a slow acceleration mode. However, time-dependent simulations using a numerical hydrodynamics code show that, for rotation rates 75 - 85% of critical, the flow can develop abrupt "kink" transitions from a steep acceleration to a decelerating solution. For rotations above 85% of critical, the hydrodynamic simulations confirm the slow acceleration, with the lower flow speed implying densities 5 - 30 times higher than the polar (or a nonrotating) wind. Still, when gravity darkening and 2D flow effects are accounted for, it seems unlikely that rotationally modified equatorial wind outflows could account for the very large densities inferred for the equatorial regions around B[e] supergiants.