With upcoming missions such as the James Webb Space Telescope (JWST), the European Extremely Large Telescope (ELT), and the Atmospheric Remote-sensing Infrared Exoplanet Large-survey (ARIEL), we soon will be on the verge of detecting and characterizing Earth-like exoplanetary atmospheres for the first time. These planets are most likely to be found around smaller and cooler K- and M-type stars. However, recent observations showed that their radiation environment might be much harsher than that of the Sun. Thus, the exoplanets are most likely exposed to an enhanced stellar radiation environment, which could affect their habitability, for example, in the form of a hazardous flux of energetic particles. Knowing the stellar radiation field, and being able to model the radiation exposure on the surface of a planet is crucial to assess its habitability. In this study, we present 3D magnetohydrodynamic (MHD)-based model efforts investigating M-stars, focusing on V374 Peg, Proxima Centauri, and LHS 1140, chosen because of their diverse astrospheric quantities. We show that V374 Peg has a much larger astrosphere (ASP) than our Sun, while Proxima Centauri and LHS 1140 most likely have ASPs comparable or even much smaller than the heliosphere, respectively. Based on a 1D transport model, for the first time, we provide numerical estimates of the modulation of Galactic cosmic rays (GCRs) within the three ASPs. We show that the impact of GCRs on the Earth-like exoplanets Proxima Centauri b and LHS 1140 b cannot be neglected in the context of exoplanetary habitability.

The AGB star $\pi^{1}$ Gruis has a known companion (at a separation of ~400 AU) which cannot explain the strong deviations from the spherical symmetry of the CSE. Recently, hydrodynamic simulations of mass transfer in closer binary systems have successfully reproduced the spiral-shaped CSEs found around a handful of sources. There is growing evidence for an even closer, undetected companion complicating the case of $\pi^{1}$ Gruis further. The improved spatial resolution allows for the investigation of the complex circumstellar morphology and the search for imprints on the CSE of the third component. We have observed the 12CO J=3-2 line emission from $\pi^{1}$ Gruis using both the compact and extended array of Atacama Large Millimeter/submillimeter Array (ALMA). The interferometric data has furthermore been combined with data from the ALMA total power (TP) array. The imaged brightness distribution has been used to constrain a non-local, non-LTE 3D radiative transfer model of the CSE. The high-angular resolution ALMA data have revealed the first example of a source on the AGB where both a faster bipolar outflow and a spiral pattern along the orbital plane can be seen in the gas envelope. The spiral can be traced in the low- to intermediate velocity, from 13 to 25 km s$^{-1}$, equatorial torus. The largest spiral-arm separation is $\approx$5".5 and consistent with a companion with an orbital period of $\approx$330 yrs and a separation of less than 70 AU. The kinematics of the bipolar outflow is consistent with it being created during a mass-loss eruption where the mass-loss rate from the system increased by at least a factor of 5 during 10-15 yrs. The spiral pattern is the result of an undetected companion. The bipolar outflow is the result of a rather recent mass-loss eruption event.

We present Atacama Large Millimeter-Submillimeter Array (ALMA) observations of CK Vulpeculae which is identified with "Nova Vulpeculae 1670". They trace obscuring dust in the inner regions of the associated nebulosity. The dust forms two cocoons, each extending ~5 arcsec north and south of the presumed location of the central star. Brighter emission is in a more compact east-west structure (2 arcsec by 1 arcsec) where the cocoons intersect. We detect line emission in NH$_2$CHO, CN, four organic molecules and C$^{17}$O. CN lines trace bubbles within the dusty cocoons; CH$_3$OH a north-south S-shaped jet; and other molecules a central cloud with a structure aligned with the innermost dust structure. The major axis of the overall dust and gas bubble structure has a projected inclination of ~24 degrees with respect to a 71 arcsec extended "hourglass" nebulosity, previously seen in H alpha. Three cocoon limbs align with dark lanes in the inner regions of the same H alpha images. The central 2 arcsec by 1 arcsec dust is resolved into a structure consistent with a warped dusty disc. The velocity structure of the jets indicates an origin at the centre of this disc and precession with an unknown period. Deceleration regions at both the northern and southern tips of the jets are roughly coincident with additional diffuse dust emission over regions approximately 2 arcsec across. These structures are consistent with a bipolar outflow expanding into surrounding high density material. We suggest that a white dwarf and brown dwarf merged between 1670 and 1672, with the observed structures and extraordinary isotopic abundances generated as a result.

We present 1D-MESA2HYDRO-3D, an open source, Python-based software tool that provides an accessible means of generating physically motivated initial conditions (ICs) for hydrodynamical simulations from 1-D stellar structure models. We test 1D-MESA2HYDRO-3D on five stellar models generated with the MESA stellar evolution code and verify its capacity as an IC generator with the Phantom smoothed-particle hydrodynamics code \citep{MESAIV, Phantom}. Consistency between the input density profiles, the 1D-MESA2HYDRO-3D-rendered particle distributions, and the state of the distributions after evolution over $10$ dynamical timescales is found for model stars ranging in structure and density from a radially extended supergiant to a white dwarf.

We present a detailed study of the 2016 eruption of nova V407 Lupi (ASASSN-16kt), including optical, near-infrared, X-ray, and ultraviolet data from SALT, SMARTS, SOAR, Chandra, Swift, and XMM-Newton. Timing analysis of the multiwavelength light-curves shows that, from 168 days post-eruption and for the duration of the X-ray supersoft source phase, two periods at 565 s and 3.57 h are detected. We suggest that these are the rotational period of the white dwarf and the orbital period of the binary, respectively, and that the system is likely to be an intermediate polar. The optical light-curve decline was very fast ($t_2 \leq$ 2.9 d), suggesting that the white dwarf is likely massive ($\gtrsim 1.25$ M$_{\odot}$). The optical spectra obtained during the X-ray supersoft source phase exhibit narrow, complex, and moving emission lines of He II, also characteristics of magnetic cataclysmic variables. The optical and X-ray data show evidence for accretion resumption while the X-ray supersoft source is still on, possibly extending its duration.

The symbiotic star R Aqr is part of a small sample of binary AGB stars observed with ALMA. The sample stars are: R Aqr, Mira, W Aql, and pi1 Gru. The sample covers a range in binary separation and wind properties, where R Aqr is the source with the smallest separation. The R Aqr binary pair consists of an M-type AGB star and a white dwarf at a separation of 45 mas, equivalent to about 10 AU at 218 pc. The aim of the ALMA study is to investigate the dependence of the wind shaping on the binary separation and to provide constraints for hydrodynamical binary interaction models. R Aqr is particularly interesting as the source with the smallest separation and a complex circumstellar environment that is strongly affected by the interaction between the two stars and by the high-energy radiation resulting from this interaction and from the hot white dwarf companion. The CO(J=3-2) line emission has been observed at 0.5" spatial resolution. The CO envelope around the binary pair is marginally resolved, showing what appears to be a rather complex distribution. The outer radius of the CO emitting region is estimated from the data and found to be about a factor of 10 larger than previously thought. This implies an average mass-loss rate during the past 100 yr of 2x10-7 Msun/yr, a factor of 45 less than previous estimates. The channel maps are presented and the molecular gas distribution is discussed and set into the context of what was previously known about the system. Additional molecular line emission detected within the bandwidth covered by the ALMA observations is also presented. Because of the limited extent of the emission, firm conclusions about the dynamical evolution of the system will have to wait for higher spatial resolution observations. However, the data presented here support the assumption that the mass-loss rate from the Mira star strongly varies and is focused on the orbital plane.

Context. Classical novae are thermonuclear explosions hosted by accreting white dwarfs in stellar binary systems. Material piles up on top of the white dwarf star under mildly degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope (mostly proton-capture reactions and beta-decays) heats the material up to peak temperatures ranging from 100 to 400 MK. In these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al), are ejected into the interstellar medium. Aims. To date, most of the efforts undertaken in the modeling of classical nova outbursts have focused on the early stages of the explosion and ejection, ignoring the interaction of the ejecta, first with the accretion disk orbiting the white dwarf, and ultimately with the secondary star. Methods. A suite of three-dimensional, SPH simulations of the interaction between the nova ejecta, the accretion disk, and the stellar companion have been performed to fill this gap, aimed at testing the influence of the different parameters (i.e., mass and velocity of the ejecta, mass and geometry of the accretion disk) on the dynamical and chemical properties of the system. Results. We discuss the conditions that lead to the disruption of the accretion disk and to mass loss from the binary system. In addition, we discuss the likelihood of chemical contamination of the stellar secondary induced by the impact with the nova ejecta and its potential effect on the next nova cycle.

We report on multiwavelength observations of nova SMCN 2016-10a. The present observational set is one of the most comprehensive for any nova in the Small Magellanic Cloud, including: low, medium, and high resolution optical spectroscopy and spectropolarimetry from SALT, FLOYDS, and SOAR; long-term OGLE $V$- and $I$- bands photometry dating back to six years before eruption; SMARTS optical and near-IR photometry from $\sim$ 11 days until over 280 days post-eruption; $Swift$ satellite X-ray and ultraviolet observations from $\sim$ 6 days until 319 days post-eruption. The progenitor system contains a bright disk and a main sequence or a sub-giant secondary. The nova is very fast with $t_2 \simeq$ 4.0 $\pm$ 1.0 d and $t_3 \simeq$ 7.8 $\pm$ 2.0 d in the $V$-band. If the nova is in the SMC, at a distance of $\sim$ 61 $\pm$ 10 kpc, we derive $M_{V,\mathrm{max}} \simeq - 10.5$ $\pm$ 0.5, making it the brightest nova ever discovered in the SMC and one of the brightest on record. At day 5 post-eruption the spectral lines show a He/N spectroscopic class and a FWHM of $\sim$ 3500 kms$^{-1}$ indicating moderately high ejection velocities. The nova entered the nebular phase $\sim$ 20 days post-eruption, predicting the imminent super-soft source turn-on in the X-rays, which started $\sim$ 28 days post-eruption. The super-soft source properties indicate a white dwarf mass between 1.2 M$_{\odot}$ and 1.3 M$_{\odot}$ in good agreement with the optical conclusions.

The star V766 Cen (=HR 5171A) was originally classified as a yellow hypergiant but lately found to more likely be a 27-36 Msun red supergiant (RSG). Recent observations indicated a close eclipsing companion in the contact or common-envelope phase. Here, we aim at imaging observations of V766 Cen to confirm the presence of the close companion. We used near-infrared H -band aperture synthesis imaging at three epochs in 2014, 2016, and 2017, employing the PIONIER instrument at the Very Large Telescope Interferometer (VLTI). The visibility data indicate a mean Rosseland angular diameter of 4.1+/-0.8 mas, corresponding to a radius of 1575+/-400 Rsun. The data show an extended shell (MOLsphere) of about 2.5 times the Rosseland diameter, which contributes about 30% of the H-band flux. The reconstructed images at the 2014 epoch show a complex elongated structure within the photospheric disk with a contrast of about 10%. The second and third epochs show qualitatively and quantitatively different structures with a single very bright and narrow feature and high contrasts of 20-30%. This feature is located toward the south-western limb of the photospheric stellar disk. We estimate an angular size of the feature of 1.7+/-0.3 mas, corresponding to a radius of 650+/-150 Rsun, and giving a radius ratio of 0.42+0.35/-0.10} compared to the primary stellar disk. We interpret the images at the 2016 and 2017 epochs as showing the close companion, or a common envelope toward the companion, in front of the primary. At the 2014 epoch, the close companion is behind the primary and not visible. Instead, the structure and contrast at the 2014 epoch are typical of a single RSG harboring giant photospheric convection cells. The companion is most likely a cool giant or supergiant star with a mass of 5+15/-3 Msun.

The S-type asymptotic giant branch (AGB) star $\pi^{1}$ Gruis has a known companion at a separation of $\approx$400 AU. The envelope structure, including an equatorial torus and a fast bipolar outflow, is rarely seen in the AGB phase and is particularly unexpected in such a wide binary system. Therefore a second, closer companion has been suggested, but the evidence is not conclusive. The new ALMA $^{12}$CO and $^{13}$CO $J$=3-2 data, together with previously published $^{12}$CO $J$=2-1 data from the Submillimeter Array (SMA), and the $^{12}$CO $J$=5-4 and $J$=9-8 lines observed with Herschel/Heterodyne Instrument for the Far-Infrared (HIFI), is modeled with the 3D non-LTE radiative transfer code SHAPEMOL. The data analysis clearly confirms the torus-bipolar structure. The 3D model of the CSE that satisfactorily reproduces the data consists of three kinematic components: a radially expanding torus with velocity slowly increasing from 8 to 13 km s$^{-1}$ along the equator plane; a radially expanding component at the center with a constant velocity of 14 km s$^{-1}$; and a fast, bipolar outflow with velocity proportionally increasing from 14 km s$^{-1}$ at the base up to 100 km s$^{-1}$ at the tip, following a linear radial dependence. The results are used to estimate an average mass-loss rate during the creation of the torus of 7.7$\times$10$^{-7}$ M$_{\odot}$ yr$^{-1}$. The total mass and linear momentum of the fast outflow are estimated at 7.3$\times$10$^{-4}$ M$_{\odot}$ and 9.6$\times$10$^{37}$ g cm s$^{-1}$, respectively. The momentum of the outflow is in excess (by a factor of about 20) of what could be generated by radiation pressure alone, in agreement with recent findings for more evolved sources. The best-fit model also suggests a $^{12}$CO/$^{13}$CO abundance ratio of 50. Possible shaping scenarios for the gas envelope are discussed

The CO(J=3-2) emission from the CSE of the binary S-type AGB star W Aql has been observed at subarcsecond resolution using ALMA. The aim of this paper is to investigate the wind properties of the AGB star and to analyse how the known companion has shaped the CSE. The average mass-loss rate during the creation of the detected CSE is estimated through modelling, using the ALMA brightness distribution and previously published single-dish measurements as observational constraints. The ALMA observations are presented and compared to the results from a 3D smoothed particle hydrodynamics (SPH) binary interaction model with the same properties as the W Aql system and with two different orbital eccentricities. Three-dimensional radiative transfer modelling is performed and the response of the interferometer is modelled and discussed. The estimated average mass-loss rate of W~Aql agrees with previous results. The size of the emitting region is consistent with photodissociation models. The CO(J=3-2) emission is dominated by a smooth component overlayed with two weak arc patterns with different separations. The larger pattern is predicted by the binary interaction model with separations of 10" and therefore likely due to the known companion. It is consistent with a binary orbit with low eccentricity. The smaller separation pattern is asymmetric and coincides with the dust distribution, but the separation timescale (200 yrs) is not consistent with any known process of the system. The separation of the known companions of the system is large enough to not have a very strong effect on the circumstellar morphology. The density contrast across the envelope of a binary with an even larger separation will not be easily detectable, even with ALMA, unless the orbit is strongly asymmetric or the AGB star has a much larger mass-loss rate.

Gamma-ray binaries are a small subclass of the high mass X-ray binary population which exhibit emission across the whole electromagnetic spectrum. We present radial velocities of 1FGL J1018.6$-$5856 based on observations obtained with the Southern African Large Telescope (SALT). We combine our measurements with those published in the literature to get a broad phase coverage. The mass function obtained supports a neutron star compact object, although a black hole mass is possible for very low inclination angles. The improved phase coverage allows constraints to be placed on the orbital eccentricity (e = 0.31 $\pm$ 0.16), which agrees with estimates from high-energy data.

Aims. Our goal is to characterize the intermediate age, detached shell carbon star U Antliae morphologically and physically in order to study the mass-loss evolution after a possible thermal pulse. Methods. High spatial resolution ALMA observations of unprecedented quality in thermal CO lines allow us to derive first critical spatial and temporal scales and constrain modeling efforts to estimate mass-loss rates for both the present day as well as the ejection period of the detached shell. Results. The detached shell is remarkably thin, overall spherically symmetric, and shows a barely resolved filamentary substructure possibly caused by instabilities in the interaction zone of winds with different outflow velocities. The expansion age of the detached shell is of the order of 2700 years and its overall width indicates a high expansion-velocity and high mass-loss period of only a few hundred years at an average mass-loss rate of $\approx$ 10$^{-5}$ $M_\odot$ yr$^{-1}$. The post-high-mass-loss-rate-epoch evolution of U Ant shows a significant decline to a substantially lower gas expansion velocity and a mass-loss rate amounting to 4$\times$10$^{-8}$ $M_\odot$ yr$^{-1}$, at present being consistent with evolutionary changes as predicted for the period between thermal pulses.

We present near-infrared interferometry of the carbon-rich asymptotic giant branch (AGB) star R Sculptoris. The visibility data indicate a broadly circular resolved stellar disk with a complex substructure. The observed AMBER squared visibility values show drops at the positions of CO and CN bands, indicating that these lines form in extended layers above the photosphere. The AMBER visibility values are best fit by a model without a wind. The PIONIER data are consistent with the same model. We obtain a Rosseland angular diameter of 8.9+-0.3 mas, corresponding to a Rosseland radius of 355+-55 Rsun, an effective temperature of 2640+-80 K, and a luminosity of log L/Lsun=3.74+-0.18. These parameters match evolutionary tracks of initial mass 1.5+-0.5 Msun and current mass 1.3+-0.7 Msun. The reconstructed PIONIER images exhibit a complex structure within the stellar disk including a dominant bright spot located at the western part of the stellar disk. The spot has an H-band peak intensity of 40% to 60% above the average intensity of the limb-darkening-corrected stellar disk. The contrast between the minimum and maximum intensity on the stellar disk is about 1:2.5. Our observations are broadly consistent with predictions by dynamic atmosphere and wind models, although models with wind appear to have a circumstellar envelope that is too extended compared to our observations. The detected complex structure within the stellar disk is most likely caused by giant convection cells, resulting in large-scale shock fronts, and their effects on clumpy molecule and dust formation seen against the photosphere at distances of 2-3 stellar radii.

We present an ongoing, systematic search for extragalactic infrared transients, dubbed SPIRITS --- SPitzer InfraRed Intensive Transients Survey. In the first year, using Spitzer/IRAC, we searched 190 nearby galaxies with cadence baselines of one month and six months. We discovered over 1958 variables and 43 transients. Here, we describe the survey design and highlight 14 unusual infrared transients with no optical counterparts to deep limits, which we refer to as SPRITEs (eSPecially Red Intermediate Luminosity Transient Events). SPRITEs are in the infrared luminosity gap between novae and supernovae, with [4.5] absolute magnitudes between -11 and -14 (Vega-mag) and [3.6]-[4.5] colors between 0.3 mag and 1.6 mag. The photometric evolution of SPRITEs is diverse, ranging from < 0.1 mag/yr to > 7 mag/yr. SPRITEs occur in star-forming galaxies. We present an in-depth study of one of them, SPIRITS 14ajc in Messier 83, which shows shock-excited molecular hydrogen emission. This shock may have been triggered by the dynamic decay of a non-hierarchical system of massive stars that led to either the formation of a binary or a proto-stellar merger.

We report spectroscopic and photometric follow-up of the peculiar nova V5852~Sgr (discovered as OGLE-2015-NOVA-01), which exhibits a combination of features from different nova classes. The photometry shows a flat-topped light curve with quasi-periodic oscillations, then a smooth decline followed by two fainter recoveries in brightness. Spectroscopy with the Southern African Large Telescope shows first a classical nova with an Fe II or Fe IIb spectral type. In the later spectrum, broad emissions from helium, nitrogen and oxygen are prominent and the iron has faded which could be an indication to the start of the nebular phase. The line widths suggest ejection velocities around $1000\,{\rm km\,s^{-1}}$. The nova is in the direction of the Galactic bulge and is heavily reddened by an uncertain amount. The $V$ magnitude 16 days after maximum enables a distance to be estimated and this suggests that the nova may be in the extreme trailing stream of the Sagittarius dwarf spheroidal galaxy. If so it is the first nova to be detected from that, or from any dwarf spheroidal galaxy. Given the uncertainty of the method and the unusual light curve we cannot rule out the possibility that it is in the bulge or even the Galactic disk behind the bulge.

Recent interpretations of narrow, variable absorption lines detected in some Type Ia supernovae suggest that their progenitors are surrounded by dense, circumstellar material. Similar variations detected in the symbiotic recurrent nova system RS Oph, which undergoes thermonuclear outbursts every ~20 years, making it an ideal candidate to investigate the origin of these lines. To this end, we present simulations of multiple mass transfer-nova cycles in RS Oph. We find that the quiescent mass transfer produces a dense, equatorial outflow, i.e., concentrated towards the binary orbital plane, and an accretion disc forms around the white dwarf. The interaction of a spherical nova outburst with these aspherical circumstellar structures produces a bipolar outflow, similar to that seen in HST imaging of the 2006 outburst. In order to produce an ionization structure that is consistent with observations, a mass-loss rate of $5 \times 10^{-7}\,\mathrm{M}_{\odot}\,\mathrm{yr}^{-1}$ from the red giant is required. The simulations also produce a polar accretion flow, which may explain the broad wings of the quiescent H {\alpha} line and hard X-rays. By comparing simulated absorption line profiles to observations of the 2006 outburst, we are able to determine which components arise in the wind and which are due to the novae. We explore the possible behaviour of absorption line profiles as they may appear should a supernova occur in a system like RS Oph. Our models show similarities to supernovae like SN 2006X, but require a high mass-loss rate, $\dot{M} \sim 10^{-6}$ to $10^{-5}\,\mathrm{M_\odot}\,\mathrm{yr}^{-1}$, to explain the variability in SN 2006X.

Mid-infrared arcs of dust emission are often seen near ionizing stars within HII regions. A possible explanations for these arcs is that they could show the outer edges of asymmetric stellar wind bubbles. We use two-dimensional, radiation-hydrodynamics simulations of wind bubbles within HII regions around individual stars to predict the infrared emission properties of the dust within the HII region. We assume that dust and gas are dynamically well-coupled and that dust properties (composition, size distribution) are the same in the HII region as outside it, and that the wind bubble contains no dust. We post-process the simulations to make synthetic intensity maps at infrared wavebands using the TORUS code. We find that the outer edge of a wind bubble emits brightly at 24um through starlight absorbed by dust grains and re-radiated thermally in the infrared. This produces a bright arc of emission for slowly moving stars that have asymmetric wind bubbles, even for cases where there is no bow shock or any corresponding feature in tracers of gas emission. The 24um intensity decreases exponentially from the arc with increasing distance from the star because the dust temperature decreases with distance. The size distribution and composition of the dust grains has quantitative but not qualitative effects on our results. Despite the simplifications of our model, we find good qualitative agreement with observations of the HII region RCW120, and can provide physical explanations for any quantitative differences. Our model produces an infrared arc with the same shape and size as the arc around CD -38 11636 in RCW120, and with comparable brightness. This suggests that infrared arcs around O stars in HII regions may be revealing the extent of stellar wind bubbles, although we have not excluded other explanations.

We present a study of Hen 2-155 and Hen 2-161, two planetary nebulae which bear striking morphological similarities to other planetary nebulae known to host close-binary central stars. Both central stars are revealed to be photometric variables while spectroscopic observations confirm that Hen 2-155 is host to a double-eclipsing, post-common-envelope system with an orbital period of 3h33m making it one of the shortest period binary central stars known. The observations of Hen 2-161 are found to be consistent with a post-common-envelope binary of period ~1 day. A detailed model of central star of Hen 2-155, is produced, showing the nebular progenitor to be a hot, post-AGB remnant of approximately 0.62 Msol, consistent with the age of the nebula, and the secondary star to be an M dwarf whose radius is almost twice the expected ZAMS radius for its mass. In spite of the small numbers, all main-sequence companions, of planetary nebulae central stars, to have had their masses and radii constrained by both photometric and spectroscopic observations have also been found to display this "inflation". The cause of the "inflation" is uncertain but is probably related to rapid accretion, immediately before the recent common-envelope phase, to which the star has not yet thermally adjusted. The chemical composition of both nebulae is also analysed, showing both to display elevated abundance discrepancy factors. This strengthens the link between elevated abundance discrepancy factors and close binarity in the nebular progenitor.

Interstellar bubbles around O stars are driven by a combination of the star's wind and ionizing radiation output. The wind contribution is uncertain because the boundary between the wind and interstellar medium is difficult to observe. Mid-infrared observations (e.g., of the H II region RCW 120) show arcs of dust emission around O stars, contained well within the H II region bubble. These arcs could indicate the edge of an asymmetric stellar wind bubble, distorted by density gradients and/or stellar motion. We present two-dimensional, radiation-hydrodynamics simulations investigating the evolution of wind bubbles and H II regions around massive stars moving through a dense (n=3000 cm^{-3}), uniform medium with velocities ranging from 4 to 16 km/s. The H II region morphology is strongly affected by stellar motion, as expected, but the wind bubble is also very aspherical from birth, even for the lowest space velocity considered. Wind bubbles do not fill their H II regions (we find filling factors of 10-20%), at least for a main sequence star with mass M~30 Msun. Furthermore, even for supersonic velocities the wind bow shock does not significantly trap the ionization front. X-ray emission from the wind bubble is soft, faint, and comes mainly from the turbulent mixing layer between the wind bubble and the H II region. The wind bubble radiates <1 per cent of its energy in X-rays; it loses most of its energy by turbulent mixing with cooler photoionized gas. Comparison of the simulations with the H II region RCW 120 shows that its dynamical age is <=0.4 Myr and that stellar motion <=4 km/s is allowed, implying that the ionizing source is unlikely to be a runaway star but more likely formed in situ. The region's youth, and apparent isolation from other O or B stars, makes it very interesting for studies of massive star formation and of initial mass functions.

We have mapped the CO(3-2) line emission around the Mira AB system at 0.5 resolution using the Atacama Large Millimeter/submillimeter Array (ALMA). The CO map shows amazing complexity. The circumstellar gas has been shaped by different dynamical actors during the evolution of the system and several morphological components can be identified. The companion is marginally resolved in continuum emission and is currently at 0.487$\pm$0.006 separation. In the main line component, centered on the stellar velocity, spiral arcs around Mira A are found. The spiral appears to be relatively flat and oriented in the orbital plane. An accretion wake behind the companion is clearly visible and the projected arc separation is of order 5''. In the blue wing of the line emission, offset from the main line, several large ($\sim$5-10''), opposing arcs are found. We tentatively suggest that this structure is created by the wind of Mira B blowing a bubble in the expanding envelope of Mira A.

Betelgeuse, a nearby red supergiant, is a runaway star with a powerful stellar wind that drives a bow shock into its surroundings. This picture has been challenged by the discovery of a dense and almost static shell that is three times closer to the star than the bow shock and has been decelerated by some external force. The two physically distinct structures cannot both be formed by the hydrodynamic interaction of the wind with the interstellar medium. Here we report that a model in which Betelgeuse's wind is photoionized by radiation from external sources can explain the static shell without requiring a new understanding of the bow shock. Pressure from the photoionized wind generates a standing shock in the neutral part of the wind and forms an almost static, photoionization-confined shell. Other red supergiants should have significantly more massive shells than Betelgeuse, because the photoionization-confined shell traps up to 35 per cent of all mass lost during the red supergiant phase, confining this gas close to the star until it explodes. After the supernova explosion, massive shells dramatically affect the supernova lightcurve, providing a natural explanation for the many supernovae that have signatures of circumstellar interaction.

We discuss two important effects for the astrospheres of runaway stars: the propagation of ionizing photons far beyond the astropause, and the rapid evolution of massive stars (and their winds) near the end of their lives. Hot stars emit ionizing photons with associated photoheating that has a significant dynamical effect on their surroundings. 3D simulations show that HII regions around runaway O stars drive expanding conical shells and leave underdense wakes in the medium they pass through. For late O stars this feedback to the interstellar medium is more important than that from stellar winds. Late in life, O stars evolve to cool red supergiants more rapidly than their environment can react, producing transient circumstellar structures such as double bow shocks. This provides an explanation for the bow shock and linear bar-shaped structure observed around Betelgeuse.

The circumstellar medium around massive stars is strongly impacted by stellar winds, radiation, and explosions. We use numerical simulations of these interactions to constrain the current properties and evolutionary history of various stars by comparison with observed circumstellar structures. Two- and three-dimensional simulations of bow shocks around red supergiant stars have shown that Betelgeuse has probably only recently evolved from a blue supergiant to a red supergiant, and hence its bow shock is very young and has not yet reached a steady state. We have also for the first time investigated the magnetohydrodynamics of the photoionised H II region around the nearby runaway O star Zeta Oph. Finally, we have calculated a grid of models of bow shocks around main sequence and evolved massive stars that has general application to many observed bow shocks, and which forms the basis of future work to model the explosions of these stars into their pre-shaped circumstellar medium.

In this Letter, we explore the hypothesis that the smooth appearance of bow shocks around some red supergiants (RSGs) might be caused by the ionization of their winds by external sources of radiation. Our numerical simulations of the bow shock generated by IRC-10414 (the first-ever RSG with an optically detected bow shock) show that the ionization of the wind results in its acceleration by a factor of two, which reduces the difference between the wind and space velocities of the star and makes the contact discontinuity of the bow shock stable for a range of stellar space velocities and mass-loss rates. Our best fit model reproduces the overall shape and surface brightness of the observed bow shock and suggests that the space velocity and mass-loss rate of IRC-10414 are $\approx$50 ${\rm km} \, {\rm s}^{-1}$ and $\approx$$10^{-6}$ $M_\odot \, {\rm yr}^{-1}$, respectively, and that the number density of the local ISM is $\approx$3 ${\rm cm}^{-3}$. It also shows that the bow shock emission comes mainly from the shocked stellar wind. This naturally explains the enhanced nitrogen abundance in the line-emitting material, derived from the spectroscopy of the bow shock. We found that photoionized bow shocks are $\approx$15$-$50 times brighter in optical line emission than their neutral counterparts, from which we conclude that the bow shock of IRC-10414 must be photoionized.

Classical barium stars are binary systems which consist of a late-type giant enriched in carbon and slow neutron capture (s-process) elements and an evolved white dwarf (WD) that is invisible at optical wavelengths. The youngest observed barium stars are surrounded by planetary nebulae (PNe), ejected soon after the wind accretion of polluted material when the WD was in its preceeding asymptotic giant branch (AGB) phase. Such systems are rare but powerful laboratories for studying AGB nucleosynthesis as we can measure the chemical abundances of both the polluted star and the nebula ejected by the polluter. Here we present evidence for a barium star in the PN Hen 2-39. The polluted giant is very similar to that found in WeBo 1. It is a cool (Teff=4250 +/- 150 K) giant enhanced in carbon ([C/H]=0.42 +/- 0.02 dex) and barium ([Ba/Fe]=1.50 +/- 0.25 dex). A spectral type of C-R3 C_24 nominally places Hen 2-39 amongst the peculiar early R-type carbon stars, however the barium enhancement and likely binary status mean that it is more likely to be a barium star with similar properties, rather than a true member of this class. An AGB star model of initial mass 1.8 Msun and a relatively large carbon pocket size can reproduce the observed abundances well, provided mass is transferred in a highly conservative way from the AGB star to the polluted star (e.g. wind Roche-lobe overflow). The nebula exhibits an apparent ring morphology in keeping with the other PNe around barium stars (WeBo 1 and A 70) and shows a high degree of ionization implying the presence of an invisible hot pre-WD companion that will require confirmation with UV observations. In contrast to A 70, the nebular chemical abundance pattern is consistent with non-Type I PNe, in keeping with the trend found from nebular s-process studies that non-Type I PNe are more likely to be s-process enhanced. (abridged)

[abridged] The 12CO/13CO ratio is often used as a measure of the 12C/13C ratio in the circumstellar environment, carrying important information about the stellar nucleosynthesis. External processes can change the 12CO and 13CO abundances, and spatially resolved studies of the 12CO/13CO ratio are needed to quantify the effect of these processes on the globally determined values. Additionally, such studies provide important information on the conditions in the circumstellar environment. The detached-shell source R Scl, displaying CO emission from recent mass loss, in a binary-induced spiral structure as well as in a clumpy shell produced during a thermal pulse, provides a unique laboratory for studying the differences in CO isotope abundances throughout its recent evolution. We observed both the 12CO(J=3-2) and the 13CO(J=3-2) line using ALMA. We find significant variations in the 12CO/13CO intensity ratios and consequently in the abundance ratios. The average CO isotope abundance ratio is at least a factor three lower in the shell (~19) than that in the present-day (<300 years) mass loss (>60). Additionally, variations in the ratio of more than an order of magnitude are found in the shell itself. We attribute these variations to the competition between selective dissociation and isotope fractionation in the shell, of which large parts cannot be warmer than ~35 K. However, we also find that the 12CO/13CO ratio in the present-day mass loss is significantly higher than the 12C/13C ratio determined in the stellar photosphere from molecular tracers (~19). The origin of this discrepancy is still unclear, but we speculate that it is due to an embedded source of UV-radiation that is primarily photo-dissociating 13CO. This radiation source could be the hitherto hidden companion. Alternatively, the UV-radiation could originate from an active chromosphere of R Scl itself....

Carbon-enhanced metal-poor (CEMP) stars count for 9-25% of all the very metal-poor stars of the halo. In at least some CEMP stars the chemical enrichment is believed to be due to wind mass transfer in the past from an AGB donor star on to a low-mass companion. However, binary population synthesis models predict much lower CEMP fractions. As an alternative to the canonical Bondi-Hoyle-Lyttleton (BHL) wind accretion model, recent hydrodynamical simulations suggest an efficient mode of wind mass transfer, called wind Roche-lobe overflow (WRLOF), can reproduce observations of AGB winds in binary systems. We use our population synthesis model to test the consequences of WRLOF on a population of CEMP stars. Compared to previous predictions based on the BHL model we find a modest increase of the fraction of CEMP stars and substantial differences in the distributions of carbon and periods in the population of CEMP stars.

Betelgeuse, the bright red supergiant (RSG) in Orion, is a runaway star. Its supersonic motion through the interstellar medium has resulted in the formation of a bow shock, a cometary structure pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuse's bow shock. We show that the bow shock morphology depends substantially on the growth timescale for Rayleigh-Taylor versus Kelvin-Helmholtz instabilities. We discuss our models in light of the recent Herschel, GALEX and VLA observations. If the mass in the bow shock shell is low (~few x 0.001 Msun), as seems to be implied by the AKARI and Herschel observations, then Betelgeuse's bow shock is very young and is unlikely to have reached a steady state. The circular, smooth bow shock shell is consistent with this conclusion. We further discuss the implications of our results, in particular, the possibility that Betelgeuse may have only recently entered the RSG phase.

The nearby red supergiant (RSG) Betelgeuse has a complex circumstellar medium out to at least 0.5 parsecs from its surface, shaped by its mass-loss history within the past 0.1 Myr, its environment, and its motion through the interstellar medium (ISM). In principle its mass-loss history can be constrained by comparing hydrodynamic models with observations. Observations and numerical simulations indicate that Betelgeuse has a very young bow shock, hence the star may have only recently become a RSG. To test this possibility we calculated a stellar evolution model for a single star with properties consistent with Betelgeuse. We incorporated the resulting evolving stellar wind into 2D hydrodynamic simulations to model a runaway blue supergiant (BSG) undergoing the transition to a RSG near the end of its life. The collapsing BSG wind bubble induces a bow shock-shaped inner shell which at least superficially resembles Betelgeuse's bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG wind's interaction with the ISM. We investigate whether this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuse's bow shock.

Carbon-enhanced metal-poor stars (CEMP) are observed as a substantial fraction of the very metal-poor stars in the Galactic halo. Most CEMP stars are also enriched in s-process elements and these are often found in binary systems. This suggests that the carbon enrichment is due to mass transfer in the past from an asymptotic giant branch (AGB) star on to a low-mass companion. Models of binary population synthesis are not able to reproduce the observed fraction of CEMP stars without invoking non-standard nucleosynthesis or a substantial change in the initial mass function. This is interpreted as evidence of missing physical ingredients in the models. Recent hydrodynamical simulations show that efficient wind mass transfer is possible in the case of the slow and dense winds typical of AGB stars through a mechanism called wind Roche-lobe overflow (WRLOF), which lies in between the canonical Bondi-Hoyle-Lyttleton (BHL) accretion and Roche-lobe overflow. WRLOF has an effect on the accretion efficiency of mass transfer and on the angular momentum lost by the binary system. The aim of this work is to understand the overall effect of WRLOF on the population of CEMP stars. To simulate populations of low-metallicity binaries we combined a synthetic nucleosynthesis model with a binary population synthesis code. In this code we implemented the WRLOF mechanism. We used the results of hydrodynamical simulations to model the effect of WRLOF on the accretion efficiency and we took the effect on the angular momentum loss into account by assuming a simple prescription. As a result the number of CEMP stars predicted by our model increases by a factor 1.2-1.8 compared to earlier results that consider the BHL prescription. Moreover, higher enrichments of carbon are produced and the final orbital period distribution is shifted towards shorter periods.

The asymptotic giant branch star R Sculptoris is surrounded by a detached shell of dust and gas. The shell originates from a thermal pulse during which the star undergoes a brief period of increased mass loss. It has hitherto been impossible to constrain observationally the timescales and mass-loss properties during and after a thermal pulse - parameters that determine the lifetime on the asymptotic giant branch and the amount of elements returned by the star. Here we report observations of CO emission from the circumstellar envelope and shell around R Sculptoris with an angular resolution of 1.3 arcsec. What was hitherto thought to be only a thin, spherical shell with a clumpy structure, is revealed to contain a spiral structure. Spiral structures associated with circumstellar envelopes have been seen previously, from which it was concluded that the systems must be binaries. Using the data, combined with hydrodynamic simulations, we conclude that R Sculptoris is a binary system that underwent a thermal pulse approximately 1800 years ago, lasting approximately 200 years. About 0.003 Msun of mass was ejected at a velocity of 14.3 km s-1 and at a rate approximately 30 times higher than the prepulse mass-loss rate. This shows that approximately 3 times more mass is returned to the interstellar medium during and immediately after a pulse than previously thought.

A significant fraction of massive stars are moving supersonically through the interstellar medium (ISM), either due to disruption of a binary system or ejection from their parent star cluster. The interaction of their wind with the ISM produces a bow shock. In late evolutionary stages these stars may undergo rapid transitions from red to blue and vice versa on the Hertzsprung-Russell diagram, with accompanying rapid changes to their stellar winds and bow shocks. Recent 3D simulations of the bow shock produced by the nearby runaway red supergiant (RSG) Betelgeuse, under the assumption of a constant wind, indicate that the bow shock is very young (<30000 years old), hence Betelgeuse may have only recently become a RSG. To test this possibility, we have calculated stellar evolution models for single stars which match the observed properties of Betelgeuse in the RSG phase. The resulting evolving stellar wind is incorporated into 2D hydrodynamic simulations in which we model a runaway blue supergiant (BSG) as it undergoes the transition to a RSG near the end of its life. We find that the collapsing BSG wind bubble induces a bow shock-shaped inner shell around the RSG wind that resembles Betelgeuse's bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG wind's interaction with the ISM. We suggest that this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuse's bow shock.

Betelgeuse, the bright, cool red supergiant in Orion, is moving supersonically relative to the local interstellar medium. The star emits a powerful stellar wind which collides with this medium, forming a cometary structure, a bow shock, pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuse's bow shock. The models include realistic low temperature cooling and cover a range of plausible interstellar medium densities and stellar velocities between 0.3 - 1.9 cm-3 and 28 - 73 km/s. We show that the flow dynamics and morphology of the bow shock differ substantially due to the preferential growth of Rayleigh-Taylor or Kelvin-Helmholtz instabilities in the models. The former dominate the models with slow stellar velocities resulting in a clumpy bow shock sub-structure, whereas the latter produce a smoother, more layered sub-structure in the fast models. If the mass in the bow shock shell is low, as seems to be implied by the AKARI luminosities (~0.003 Msun), then Betelgeuse's bow shock is very young and is unlikely to have reached a steady state. The circular nature of the bow shock shell is consistent with this conclusion. Thus, our results suggest that Betelgeuse entered the red supergiant phase only recently.

Herschel's PACS instrument observed the environment of the binary system Mira Ceti in the 70 and 160 micron bands. These images reveal bright structures shaped as five broken arcs and fainter filaments in the ejected material of Mira's primary star. The overall shape of the IR emission around Mira deviates significantly from the expected alignment with Mira's exceptionally high space velocity. The observed broken arcs are neither connected to each other nor are they of a circular shape; they stretch over angular ranges of 80 to 100 degrees. By comparing Herschel and GALEX data, we found evidence for the disruption of the IR arcs by the fast outflow visible in both Halpha and the far UV. Radial intensity profiles are derived, which place the arcs at distances of 6-85" (550 - 8000 AU) from the binary. Mira's IR environment appears to be shaped by the complex interaction of Mira's wind with its companion, the bipolar jet, and the ISM.