Red supergiant stars are surrounded by a gaseous and dusty circumstellar environment created by their mass loss which spreads heavy elements into the interstellar medium. The structure and the dynamics of this envelope are crucial to understand the processes driving the red supergiant mass loss and the shaping of the pre-supernova ejecta. We have observed the emission from the CO $J = 2-1$ line from the red supergiant star $\mu$~Cep with the NOEMA interferometer. In the line the synthesized beam was $0.92 \times 0.72$~arcsec ($590 \times 462$~au at 641~pc). The continuum map shows only the unresolved contribution of the free-free emission of the star chromosphere. The continuum-subtracted channel maps reveal a very inhomogeneous and clumpy circumstellar environment. In particular, we detected a bright CO clump, as bright as the central source in the line, at 1.80~arcsec south-west from the star, in the blue channel maps. After a deprojection of the radial velocity assuming two different constant wind velocities, the observations were modelled using the 3D radiative transfer code \textsc{lime} to derive the characteristics of the different structures. We determine that the gaseous clumps observed around $\mu$~Cep are responsible for a mass loss rate of $(4.9 \pm 1.0) \times 10^{-7}~{\rm M}_\odot\,{\rm yr}^{-1}$, in addition to a spatially unresolved wind component with an estimated mass-loss rate of $2.0 \times 10^{-6}~{\rm M}_\odot\,{\rm yr}^{-1}$. Therefore, the clumps have a significant role in $\mu$~Cep's mass loss ($\ge 25 \%$). We cannot exclude that the unresolved central outflow may be made of smaller unresolved clumps.

The material lost through stellar winds of Asymptotic Giant Branch (AGB) stars is one of the main contributors to the chemical enrichment of galaxies. The general hypothesis of the mass loss mechanism of AGB winds is a combination of stellar pulsations and radiative pressure on dust grains, yet current models still suffer from limitations. Among others, they assume chemical equilibrium of the gas, which may not be justified due to rapid local dynamical changes in the wind. This is important as it is the chemical composition that regulates the thermal structure of the wind, the creation of dust grains in the wind, and ultimately the mass loss by the wind. Using a self-consistent hydrochemical model, we investigated how non-equilibrium chemistry affects the dynamics of the wind. This paper compares a hydrodynamical and a hydrochemical dust-free wind, with focus on the chemical heating and cooling processes. No sustainable wind arises in a purely hydrodynamical model with physically reasonable pulsations. Moreover, temperatures are too high for dust formation to happen, rendering radiative pressure on grains impossible. A hydrochemical wind is even harder to initiate due to efficient chemical cooling. However, temperatures are sufficiently low in dense regions for dust formation to take place. These regions occur close to the star, which is needed for radiation pressure on dust to sufficiently aid in creating a wind. Extending this model self-consistently with dust formation and evolution, and including radiation pressure, will help to understand the mass loss by AGB winds.

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.

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.

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).

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.