Showing 1–43 of 43 results for author: Mohamed, S

Particle acceleration at the bow shock of runaway star LS 2355: non-thermal radio emission but no $γ$-ray counterpart

Authors: J. van den Eijnden, S. Mohamed, F. Carotenuto, S. Motta, P. Saikia, D. R. A. Williams-Baldwin

Abstract: Massive stars that travel at supersonic speeds can create bow shocks as their stellar winds interact with the surrounding interstellar medium. These bow shocks - prominent sites for mechanical feedback of individual massive stars - are predominantly observed in the infrared band. Confirmed high-energy emission from stellar bow shocks has remained elusive and confirmed radio counterparts, while ris… ▽ More Massive stars that travel at supersonic speeds can create bow shocks as their stellar winds interact with the surrounding interstellar medium. These bow shocks - prominent sites for mechanical feedback of individual massive stars - are predominantly observed in the infrared band. Confirmed high-energy emission from stellar bow shocks has remained elusive and confirmed radio counterparts, while rising in recent years, remain rare. Here, we present an in-depth multi-wavelength exploration of the bow shock driven by LS 2355, focusing on its non-thermal properties. Using the most-recent Fermi source catalogue, we rule out its previously-proposed association with an unidentified $γ$-ray source. Furthermore, we use deep ASKAP observations from the Rapid ASKAP Continuum Survey and the Evolutionary Map of the Universe survey to identify a non-thermal radio counterpart: the third spectrally confirmed non-thermal bow shock counterpart after BD +43$^{\rm o}$ 3654 and BD +60$^{\rm o}$ 2522. We finally use WISE IR data and Gaia to study the surrounding ISM and update the motion of LS 2355. Specifically, we derive a substantially reduced stellar velocity, $v_* = 7.0\pm2.5$ km/s, compared to previous estimates. The observed non-thermal properties of the bow shock can be explained by an interaction between the wind of LS 2355 and a dense HII region, at a magnetic field close to the maximum magnetic field strength allowed by the compressibility of the ISM. Similar to earlier works, we find that the thermal radio emission of the shocked ISM is likely to be substantially suppressed for it to be consistent with the observed radio spectrum. △ Less

Submitted 29 June, 2024; originally announced July 2024.

Comments: Accepted for publication in MNRAS

The messy death of a multiple star system and the resulting planetary nebula as observed by JWST

Authors: Orsola De Marco, Muhammad Akashi, Stavros Akras, Javier Alcolea, Isabel Aleman, Philippe Amram, Bruce Balick, Elvire De Beck, Eric G. Blackman, Henri M. J. Boffin, Panos Boumis, Jesse Bublitz, Beatrice Bucciarelli, Valentin Bujarrabal, Jan Cami, Nicholas Chornay, You-Hua Chu, Romano L. M. Corradi, Adam Frank, Guillermo Garcia-Segura, D. A. Garcia-Hernandez, Jorge Garcia-Rojas, Veronica Gomez-Llanos, Denise R. Goncalves, Martin A. Guerrero , et al. (44 additional authors not shown)

Abstract: Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures,… ▽ More Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures, likely shaped by a low-mass companion orbiting the central star at 40-60 AU. The images also reveal a mid-IR excess at the central star interpreted as a dusty disk, indicative of an interaction with another, closer companion. Including the previously known, A-type visual companion, the progenitor of the NGC3132 PN must have been at least a stellar quartet. The JWST images allow us to generate a model of the illumination, ionisation and hydrodynamics of the molecular halo, demonstrating the power of JWST to investigate complex stellar outflows. Further, new measurements of the A-type visual companion allow us to derive the value for the mass of the progenitor of a central star to date with excellent precision: 2.86+/-0.06 Mo. These results serve as path finders for future JWST observations of PNe providing unique insight into fundamental astrophysical processes including colliding winds, and binary star interactions, with implications for supernovae and gravitational wave systems. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 32 pages, 5 figures for the main article. 12 pages 8 figures for the supplementary material

Journal ref: Nature Astronomy, 2022, Vol. 6, p. 1421

Radio detections of IR-selected runaway stellar bow shocks

Authors: J. Van den Eijnden, P. Saikia, S. Mohamed

Abstract: Massive stars moving at supersonic peculiar velocities through the interstellar medium (ISM) can create bow shocks, arc-like structures at the interface between the stellar wind and the ISM. Many such bow shocks have been detected and catalogued at IR wavelengths, but detections in other wavebands remain rare. Strikingly, while electrons are expected to be accelerated in the bow shock and their no… ▽ More Massive stars moving at supersonic peculiar velocities through the interstellar medium (ISM) can create bow shocks, arc-like structures at the interface between the stellar wind and the ISM. Many such bow shocks have been detected and catalogued at IR wavelengths, but detections in other wavebands remain rare. Strikingly, while electrons are expected to be accelerated in the bow shock and their non-thermal emission may include synchrotron emission at low frequencies, only two massive runaway stellar bow shocks have to date been detected in the radio band. Here, we examine a sample of fifty IR-detected bow shocks from the E-BOSS catalogues in recently released radio images from the Rapid ASKAP Continuum Survey (RACS). We identify three confident and three likely counterparts, as well as three inconclusive candidates requiring confirmation via follow-up observations. These detections significantly increase the number of known radio massive stellar bow shocks and highlight the advantage of dedicated searches with current and next-generation radio telescopes. We investigate the underlying radio emission mechanism for these radio sources, finding a mix of free-free-dominated and synchrotron-dominated systems. We also discuss the non-detected targets by putting constraints on their emission properties and investigating their detectability with future observations. Finally, we propose several future avenues of research to advance the study and understanding of bow shocks at radio frequencies. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Comments: 17 pages, 12 figures, 3 tables. Resubmitted to MNRAS after final round of very minor referee comments

3D models of the circumstellar environments of evolved stars: Formation of multiple spiral structures

Authors: Elias Aydi, Shazrene Mohamed

Abstract: We present 3D hydrodynamic models of the interaction between the outflows of evolved, pulsating, Asymptotic Giant Branch (AGB) stars and nearby ($< 3$ stellar radii) sub-stellar companions ($M_{\mathrm{comp}} \lesssim 40$ M$_J$). Our models show that due to resonances between the orbital period of the companion and the pulsation period of the AGB star, multiple spiral structures can form; the shoc… ▽ More We present 3D hydrodynamic models of the interaction between the outflows of evolved, pulsating, Asymptotic Giant Branch (AGB) stars and nearby ($< 3$ stellar radii) sub-stellar companions ($M_{\mathrm{comp}} \lesssim 40$ M$_J$). Our models show that due to resonances between the orbital period of the companion and the pulsation period of the AGB star, multiple spiral structures can form; the shocks driven by the pulsations are enhanced periodically in different regions as they encounter the denser material created by the sub-stellar companion's wake. We discuss the properties of these spiral structures and the effect of the companion parameters on them. We also demonstrate that the gravitational potential of the nearby companion enhances the mass loss from the AGB star. For more massive ($M_{\mathrm{comp}} > 40$ M$_J$) and more distant companions ($> 4$ stellar radii), a single spiral arm forms. We discuss the possibility of observing these structures with the new generations of high-resolution, high-sensitivity instruments, and using them to `find' sub-stellar companions around bright, evolved stars. Our results also highlight possible structures that could form in our solar system when the Sun turns into an AGB star. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: 31 pages, 30 figures, accepted for publication at MNRAS

L-GALAXIES 2020: The formation and chemical evolution of stellar haloes in Milky Way Analogues and galaxy clusters

Authors: Geoff G. Murphy, Robert M. Yates, Shazrene Mohamed

Abstract: We present an analysis of the formation and chemical evolution of stellar haloes around (a) Milky Way Analogue (MWA) galaxies and (b) galaxy clusters in the L-Galaxies 2020 semi-analytic model of galaxy evolution. Observed stellar halo properties are better reproduced when assuming a gradual stripping model for the removal of cold gas and stars from satellites, compared to an instantaneous strippi… ▽ More We present an analysis of the formation and chemical evolution of stellar haloes around (a) Milky Way Analogue (MWA) galaxies and (b) galaxy clusters in the L-Galaxies 2020 semi-analytic model of galaxy evolution. Observed stellar halo properties are better reproduced when assuming a gradual stripping model for the removal of cold gas and stars from satellites, compared to an instantaneous stripping model. The slope of the stellar mass -- metallicity relation for MWA stellar haloes is in good agreement with that observed in the local Universe. This extends the good agreement between L-Galaxies 2020 and metallicity observations from the gas and stars inside galaxies to those outside. Halo stars contribute on average only $\sim{}0.1$ per cent of the total circumgalactic medium (CGM) enrichment by $z=0$ in MWAs, ejecting predominantly carbon produced by AGB stars. Around a quarter of MWAs have a single `significant progenitor' with a mean mass of $\sim{}2.3\times{}10^{9}M_{\odot}$, similar to that measured for Gaia Enceladus. For galaxy clusters, L-Galaxies 2020 shows good correspondence with observations of stellar halo mass fractions, but slightly over-predicts stellar masses. Assuming a Navarro-Frenk-White profile for the stellar halo mass distribution provides the best agreement. On average, the intracluster stellar component (ICS) is responsible for 5.4 per cent of the total intracluster medium (ICM) iron enrichment, exceeding the contribution from the brightest cluster galaxy (BCG) by $z=0$. We find that considering gradual stripping of satellites and realistic radial profiles is crucial for accurately modelling stellar halo formation on all scales in semi-analytic models. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: 19 pages, 17 figures. Submitted to MNRAS, comments welcome

MeerKAT discovery of radio emission from the Vela X-1 bow shock

Authors: J. van den Eijnden, I. Heywood, R. Fender, S. Mohamed, G. R. Sivakoff, P. Saikia, T. D. Russell, S. Motta, J. C. A. Miller-Jones, P. A. Woudt

Abstract: Vela X-1 is a runaway X-ray binary system hosting a massive donor star, whose strong stellar wind creates a bow shock as it interacts with the interstellar medium. This bow shock has previously been detected in H$α$ and IR, but, similar to all but one bow shock from a massive runaway star (BD+43$^{\rm o}$3654), has escaped detection in other wavebands. We report on the discovery of $1.3$ GHz radio… ▽ More Vela X-1 is a runaway X-ray binary system hosting a massive donor star, whose strong stellar wind creates a bow shock as it interacts with the interstellar medium. This bow shock has previously been detected in H$α$ and IR, but, similar to all but one bow shock from a massive runaway star (BD+43$^{\rm o}$3654), has escaped detection in other wavebands. We report on the discovery of $1.3$ GHz radio emission from the Vela X-1 bow shock with the MeerKAT telescope. The MeerKAT observations reveal how the radio emission closely traces the H$α$ line emission, both in the bow shock and in the larger-scale diffuse structures known from existing H$α$ surveys. The Vela X-1 bow shock is the first stellar-wind-driven radio bow shock detected around an X-ray binary. In the absence of a radio spectral index measurement, we explore other avenues to constrain the radio emission mechanism. We find that thermal/free-free emission can account for the radio and H$α$ properties, for a combination of electron temperature and density consistent with earlier estimates of ISM density and the shock enhancement. In this explanation, the presence of a local ISM over-density is essential for the detection of radio emission. Alternatively, we consider a non-thermal/synchrotron scenario, evaluating the magnetic field and broad-band spectrum of the shock. However, we find that exceptionally high fractions ($\gtrsim 13$%) of the kinetic wind power would need to be injected into the relativistic electron population to explain the radio emission. Assuming lower fractions implies a hybrid scenario, dominated by free-free radio emission. Finally, we speculate about the detectability of radio bow shocks and whether it requires exceptional ISM or stellar wind properties. △ Less

Submitted 24 November, 2021; v1 submitted 19 November, 2021; originally announced November 2021.

Comments: 17 pages including appendices, 12 figures; accepted for publication in MNRAS. Updated version with minor corrections to reference list

The binary central star of the bipolar pre-planetary nebula IRAS 08005-2356 (V510 Pup)

Authors: Rajeev Manick, Brent Miszalski, Devika Kamath, Patricia A. Whitelock, Hans Van Winckel, Bruce J. Hrivnak, Brad N. Barlow, Shazrene Mohamed

Abstract: Current models predict that binary interactions are a major ingredient for the formation of bipolar planetary nebulae (PNe) and pre-planetary nebulae (PPNe). Despite years of radial velocity (RV) monitoring, the paucity of known binaries amongst the latter systems is insufficient to examine this relationship in detail. In this paper, we report on the discovery of a long period (P=2654$\pm$124 d) b… ▽ More Current models predict that binary interactions are a major ingredient for the formation of bipolar planetary nebulae (PNe) and pre-planetary nebulae (PPNe). Despite years of radial velocity (RV) monitoring, the paucity of known binaries amongst the latter systems is insufficient to examine this relationship in detail. In this paper, we report on the discovery of a long period (P=2654$\pm$124 d) binary at the centre of the Galactic bipolar PPN, IRAS 08005-2356 (V510 Pup) determined from long-term spectroscopic and near-infrared time series data. The spectroscopic orbit is fit with an eccentricity of 0.36$\pm$0.05 that is similar to other long period post-AGB binaries. Time resolved H$α$ profiles reveal high-velocity outflows (jets) with de-projected velocities up to 231$_{-27}^{+31}$ km s$^{-1}$ seen at phases when the luminous primary is behind the jet. The outflow traced by H$α$ is likely produced via accretion onto a main sequence companion for which we calculate a mass of 0.63$\pm$0.13 M$_\odot$. This discovery is one of the first cases of a confirmed binary PPN and demonstrates the importance of high-resolution spectroscopic monitoring surveys on large telescopes in revealing binarity among these systems. △ Less

Submitted 20 August, 2021; originally announced August 2021.

Comments: Accepted for publication in MNRAS, 10 Pages, 7 Figures

AT 2019qyl in NGC 300: Internal Collisions in the Early Outflow from a Very Fast Nova in a Symbiotic Binary

Authors: Jacob E. Jencson, Jennifer E. Andrews, Howard E. Bond, Viraj Karambelkar, David J. Sand, Schuyler D. van Dyk, Nadejda Blagorodnova, Martha L. Boyer, Mansi M. Kasliwal, Ryan M. Lau, Shazrene Mohamed, Robert Williams, Patricia A. Whitelock, Rachael C. Amaro, K. Azalee Bostroem, Yize Dong, Michael J. Lundquist, Stefano Valenti, Samuel D. Wyatt, Jamie Burke, Kishalay De, Saurabh W. Jha, Joel Johansson, César Rojas-Bravo, David A. Coulter , et al. (17 additional authors not shown)

Abstract: Nova eruptions, thermonuclear explosions on the surfaces of white dwarfs (WDs), are now recognized to be among the most common shock-powered astrophysical transients. We present the early discovery and rapid ultraviolet (UV), optical, and infrared (IR) temporal development of AT 2019qyl, a recent nova in the nearby Sculptor Group galaxy NGC 300. The light curve shows a rapid rise lasting… ▽ More Nova eruptions, thermonuclear explosions on the surfaces of white dwarfs (WDs), are now recognized to be among the most common shock-powered astrophysical transients. We present the early discovery and rapid ultraviolet (UV), optical, and infrared (IR) temporal development of AT 2019qyl, a recent nova in the nearby Sculptor Group galaxy NGC 300. The light curve shows a rapid rise lasting $\lesssim 1$ day, reaching a peak absolute magnitude of $M_V = -9.2$ mag, and a very fast decline, fading by 2 mag over 3.5 days. A steep dropoff in the light curves after 71 days and the rapid decline timescale suggest a low-mass ejection from a massive WD with $M_{\rm WD} \gtrsim 1.2~M_{\odot}$. We present an unprecedented view of the early spectroscopic evolution of such an event. Three spectra prior to the peak reveal a complex, multicomponent outflow giving rise to internal collisions and shocks in the ejecta of an He/N-class nova. We identify a coincident IR-variable counterpart in the extensive preeruption coverage of the transient location and infer the presence of a symbiotic progenitor system with an O-rich asymptotic-giant-branch donor star, as well as evidence for an earlier UV-bright outburst in 2014. We suggest that AT 2019qyl is analogous to the subset of Galactic recurrent novae with red-giant companions such as RS Oph and other embedded nova systems like V407 Cyg. Our observations provide new evidence that internal shocks between multiple, distinct outflow components likely contribute to the generation of the shock-powered emission from such systems. △ Less

Submitted 5 September, 2021; v1 submitted 22 February, 2021; originally announced February 2021.

Comments: 32 pages, 10 figures, accepted in ApJ

On the Diversity of M-Star Astrospheres and the Role of Galactic Cosmic Rays Within

Authors: Konstantin Herbst, Klaus Scherer, Stefan E. S. Ferreira, Lennart R. Baalmann, N. Eugene Engelbrecht, Horst Fichtner, Jens Kleimann, R. Du Toit Strauss, Daniel M. Moeketsi, Shazrene Mohamed

Abstract: 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. Ho… ▽ More 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. △ Less

Submitted 17 June, 2020; originally announced June 2020.

Comments: accepted at ApJ Letters

The extended molecular envelope of the asymptotic giant branch star $π^{1}$ Gruis as seen by ALMA II. The spiral-outflow observed at high-angular resolution

Authors: L. Doan, S. Ramstedt, W. H. T. Vlemmings, S. Mohamed, S. Höfner, E. De Beck, F. Kerschbaum, M. Lindqvist, M. Maercker, C. Paladini, M. Wittkowski

Abstract: The AGB star $π^{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 complica… ▽ More The AGB star $π^{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 $π^{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 $π^{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. △ Less

Submitted 27 November, 2019; v1 submitted 25 November, 2019; originally announced November 2019.

Comments: 12 pages, 11 figures

Density Conversion between 1-D and 3-D Stellar Models with 1D-MESA2HYDRO-3D

Authors: Meridith Joyce, Lianne Lairmore, Daniel J. Price, Thomas Reichardt, Shazrene Mohamed

Abstract: 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-part… ▽ More 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. △ Less

Submitted 21 July, 2019; originally announced July 2019.

Comments: Accepted to ApJ; 14 pages, 13 figures

ALMA reveals the aftermath of a white dwarf--brown dwarf merger in CK Vulpeculae

Authors: Stewart Eyres, Aneurin Evans, Albert Zijlstra, Adam Avison, Robert Gehrz, Marcin Hajduk, Sumner Starrfield, Shazrene Mohamed, Charles Woodward, R. Mark Wagner

Abstract: 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 arcs… ▽ More 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. △ Less

Submitted 11 October, 2019; v1 submitted 16 September, 2018; originally announced September 2018.

Comments: 10 pages, 4 figures, accepted for publication in MNRAS 14th September 2018

Multiwavelength observations of V407 Lupi (ASASSN-16kt) --- a very fast nova erupting in an intermediate polar

Authors: E. Aydi, M. Orio, A. P. Beardmore, J. -U. Ness, K. L. Page, N. P. M. Kuin, F. M. Walter, D. A. H. Buckley, S. Mohamed, P. Whitelock, J. P. Osborne, J. Strader, L. Chomiuk, M. J. Darnley, A. Dobrotka, A. Kniazev, B. Miszalski, G. Myers, N. Ospina, M. Henze, S. Starrfield, C. E. Woodward

Abstract: 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… ▽ More 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. △ Less

Submitted 2 July, 2018; originally announced July 2018.

Comments: 39 pages, 32 figures, 23 tables, accepted for publication in MNRAS

CO envelope of the symbiotic star R Aqr seen by ALMA

Authors: S. Ramstedt, S. Mohamed, T. Olander, W. H. T. Vlemmings, T. Khouri, S. Liljegren

Abstract: 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… ▽ More 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. △ Less

Submitted 13 June, 2018; originally announced June 2018.

Comments: 6 pages, 4 figures, accepted in A&A

Journal ref: A&A 616, A61 (2018)

Three-dimensional simulations of the interaction between the nova ejecta, the accretion disk, and the companion star

Authors: Joana Figueira, Jordi Jose, Enrique Garcia-Berro, Simon W. Campbell, Domingo Garcia-Senz, Shazrene Mohamed

Abstract: 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 temperat… ▽ More 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. △ Less

Submitted 22 December, 2017; originally announced December 2017.

Comments: Accepted for Publication in Astronomy and Astrophysics

Journal ref: A&A 613, A8 (2018)

Multiwavelength observations of nova SMCN 2016-10a --- one of the brightest novae ever observed

Authors: E. Aydi, K. L. Page, N. P. M. Kuin, M. J. Darnley, F. M. Walter, P. Mróz, D. Buckley, S. Mohamed, P. Whitelock, P. Woudt, S. C. Williams, M. Orio, R. E. Williams, A. P. Beardmore, J. P. Osborne, A. Kniazev, V. A. R. M. Ribeiro, A. Udalski, J. Strader, L. Chomiuk

Abstract: 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… ▽ More 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. △ Less

Submitted 10 October, 2017; originally announced October 2017.

Comments: Accepted for publication in MNRAS on 2017 October 10 (31 pages, 26 figures, 11 tables)

The extended molecular envelope of the asymptotic giant branch star $π^{1}$ Gruis as seen by ALMA I. Large-scale kinematic structure and CO excitation properties

Authors: L. Doan, S. Ramstedt, W. H. T. Vlemmings, S. Höfner, E. De Beck, F. Kerschbaum, M. Lindqvist, M. Maercker, S. Mohamed, C. Paladini, M. Wittkowski

Abstract: The S-type asymptotic giant branch (AGB) star $π^{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… ▽ More The S-type asymptotic giant branch (AGB) star $π^{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 △ Less

Submitted 27 September, 2017; originally announced September 2017.

Journal ref: A&A, Volume 605, September 2017

Multi-epoch VLTI-PIONIER imaging of the supergiant V766 Cen: Image of the close companion in front of the primary

Authors: M. Wittkowski, F. J. Abellan, B. Arroyo-Torres, A. Chiavassa, J. C. Guirado, J. M. Marcaide, A. Alberdi, W. J. de Wit, K. -H. Hofmann, A. Meilland, F. Millour, S. Mohamed, J. Sanchez-Bermudez

Abstract: 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 imagin… ▽ More 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. △ Less

Submitted 27 September, 2017; originally announced September 2017.

Comments: 4 pages, 2 figures plus 4 pages appendix. Accepted for publication in Astronomy and Astrophysics (A&A) as a letter

Journal ref: A&A 606, L1 (2017)

The circumstellar envelope around the S-type AGB star W Aql Effects of an eccentric binary orbit

Authors: S. Ramstedt, S. Mohamed, W. H. T. Vlemmings, T. Danilovich, M. Brunner, E. De Beck, E. M. L. Humphreys, M. Lindqvist, M. Maercker, H. Olofsson, F. Kerschbaum, G. Quintana-Lacaci

Abstract: 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 distributio… ▽ More 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. △ Less

Submitted 21 September, 2017; originally announced September 2017.

Comments: 10 pages, 8 figures

Journal ref: A&A 605, A126 (2017)

The orbit of the gamma-ray binary 1FGL J1018.6-5856

Authors: Itumeleng Monageng, Vanessa McBride, Lee Townsend, Alexei Kniazev, Shazrene Mohamed, Markus Böttcher

Abstract: 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 obta… ▽ More 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. △ Less

Submitted 21 August, 2017; originally announced August 2017.

Comments: 6 pages, 5 figures, accepted for publication in ApJ

Rings and filaments: The remarkable detached CO shell of U Antliae

Authors: F. Kerschbaum, M. Maercker, M. Brunner, M. Lindqvist, H. Olofsson, M. Mecina, E. De Beck, M. A. T. Groenewegen, E. Lagadec, S. Mohamed, C. Paladini, S. Ramstedt, W. H. T. Vlemmings, M. Wittkowski

Abstract: 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 m… ▽ More 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. △ Less

Submitted 9 August, 2017; originally announced August 2017.

Journal ref: A&A 605, A116 (2017)

Aperture synthesis imaging of the carbon AGB star R Sculptoris: Detection of a complex structure and a dominating spot on the stellar disk

Authors: M. Wittkowski, K. -H. Hofmann, S. Höfner, J. B. Le Bouquin, W. Nowotny, C. Paladini, J. Young, J. -P. Berger, M. Brunner, I. de Gregorio-Monsalvo, K. Eriksson, J. Hron, E. M. L. Humphreys, M. Lindqvist, M. Maercker, S. Mohamed, H. Olofsson, S. Ramstedt, G. Weigelt

Abstract: 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 va… ▽ More 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. △ Less

Submitted 8 February, 2017; originally announced February 2017.

Comments: 13 pages, 11 figures, 4 pages appendix, accepted for publication in Astronomy and Astrophysics (A&A)

Journal ref: A&A 601, A3 (2017)

SPIRITS: Uncovering Unusual Infrared Transients With Spitzer

Authors: Mansi M. Kasliwal, John Bally, Frank Masci, Ann Marie Cody, Howard E. Bond, Jacob E. Jencson, Samaporn Tinyanont, Yi Cao, Carlos Contreras, Devin A. Dykhoff, Samuel Amodeo, Lee Armus, Martha Boyer, Matteo Cantiello, Robert L. Carlon, Alexander C. Cass, David Cook, David T. Corgan, Joseph Faella, Ori D. Fox, Wayne Green, Robert Gehrz, George Helou, Eric Hsiao, Joel Johansson , et al. (22 additional authors not shown)

Abstract: 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 transient… ▽ More 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. △ Less

Submitted 4 January, 2017; originally announced January 2017.

V5852 Sgr: An Unusual Nova Possibly Associated with the Sagittarius Stream

Authors: E. Aydi, P. Mróz, P. A. Whitelock, S. Mohamed, Ł. Wyrzykowski, A. Udalski, P. Vaisanen, T. Nagayama, M. Dominik, A. Scholz, H. Onozato, R. E. Williams, S. T. Hodgkin, S. Nishiyama, M. Yamagishi, A. M. S. Smith, T. Ryu, A. Iwamatsu, I. Kawamata

Abstract: 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… ▽ More 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. △ Less

Submitted 8 June, 2016; originally announced June 2016.

Comments: Accepted for publication in MNRAS on June 8th, 2016 (11 pages, 8 figures, 5 tables)

Modelling the circumstellar medium in RS Ophiuchi and its link to Type Ia supernovae

Authors: Richard A. Booth, Shazrene Mohamed, Philipp Podsiadlowski

Abstract: 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 prese… ▽ More 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 α 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. △ Less

Submitted 11 January, 2016; originally announced January 2016.

Comments: 15 pages, 9 figures. MNRAS accepted

Detecting stellar-wind bubbles through infrared arcs in HII regions

Authors: Jonathan Mackey, Thomas J. Haworth, Vasilii V. Gvaramadze, Shazrene Mohamed, Norbert Langer, Tim J. Harries

Abstract: 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 regio… ▽ More 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. △ Less

Submitted 21 December, 2015; originally announced December 2015.

Comments: 16 pages, 11 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 586, A114 (2016)

The post-common envelope central stars of the planetary nebulae Henize 2-155 and Henize 2-161

Authors: David Jones, Henri M. J. Boffin, Pablo Rodríguez-Gil, Roger Wesson, Romano L. M. Corradi, Brent Miszalski, Shazrene Mohamed

Abstract: 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 i… ▽ More 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. △ Less

Submitted 10 June, 2015; v1 submitted 4 June, 2015; originally announced June 2015.

Comments: Accepted for publication in A&A. 14 main pages (10 pages of online material), 12 figures, 10 tables. Resubmission to correct table numbering

Journal ref: A&A 580, A19 (2015)

The wonderful complexity of the Mira AB system

Authors: S. Ramstedt, S. Mohamed, W. H. T. Vlemmings, M. Maercker, R. Montez, A. Baudry, E. De Beck, M. Lindqvist, H. Olofsson, E. M. L. Humphreys, A. Jorissen, F. Kerschbaum, A. Mayer, M. Wittkowski, N. L. J. Cox, E. Lagadec, M. L. Leal-Ferreira, C. Paladini, A. Pérez-Sánchez, S. Sacuto

Abstract: 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 emissi… ▽ More 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. △ Less

Submitted 6 October, 2014; originally announced October 2014.

Comments: Letter accepted in A&A

Journal ref: A&A 570, L14 (2014)

Wind bubbles within H II regions around slowly moving stars

Authors: Jonathan Mackey, Vasilii V. Gvaramadze, Shazrene Mohamed, Norbert Langer

Abstract: 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 cou… ▽ More 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. △ Less

Submitted 17 November, 2014; v1 submitted 30 September, 2014; originally announced October 2014.

Comments: 14 pages, 11 figures, accepted for publication in Astronomy and Astrophysics (new version corrects an error in the simulation postprocessing, figs 6,7,11 are modified slightly, conclusions unchanged)

Journal ref: A&A 573, A10 (2015)

Interacting supernovae from photoionization-confined shells around red supergiant stars

Authors: Jonathan Mackey, Shazrene Mohamed, Vasilii V. Gvaramadze, Rubina Kotak, Norbert Langer, Dominique M. -A. Meyer, Takashi J. Moriya, Hilding R. Neilson

Abstract: 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 hydrodyna… ▽ More 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. △ Less

Submitted 11 August, 2014; originally announced August 2014.

Comments: 22 pages, 8 figures, authors' version, accepted for publication in Nature; minor differences may exist compared with published version; final version available at http://dx.doi.org/10.1038/nature13522

Journal ref: Nature, 512, 282-285 (21 August 2014)

Effects of stellar evolution and ionizing radiation on the environments of massive stars

Authors: Jonathan Mackey, Norbert Langer, Shazrene Mohamed, Vasilii V. Gvaramadze, Hilding R. Neilson, Dominique M. -A. Meyer

Abstract: 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… ▽ More 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. △ Less

Submitted 31 July, 2014; originally announced July 2014.

Comments: 5 pages, 2 figures, conference proceedings in "Wind Bubbles, Astrospheres and the Heliosphere: Environments and Cosmic Rays", Bochum, Germany (Nov 2013). Published in ASTRA proceedings, an Open Access Journal for Refereed Proceedings in Extraterrestrial Research

Journal ref: ASTRA Proceedings, Volume 1, Issue 1, 2014, pp.61-63

The Circumstellar Medium of Massive Stars in Motion

Authors: Jonathan Mackey, Norbert Langer, Dominique M. -A. Meyer, Vasilii V. Gvaramadze, Shazrene Mohamed, Hilding R. Neilson, Andrea Mignone

Abstract: 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 Betelgeu… ▽ More 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. △ Less

Submitted 3 June, 2014; originally announced June 2014.

Comments: 7 pages, 7 figures, page format edited to save space, conference proceedings: NIC Symposium, 12-13 February 2014, Juelich, Germany (Eds. K. Binder, G. Muenster, M. Kremer)

Journal ref: Publication Series of the John von Neumann Institute for Computing (NIC), 2014, Volume 47 (ISBN 978-3-89336-933-1), pp. 77-85

On the stability of bow shocks generated by red supergiants: the case of IRC-10414

Authors: D. M. -A. Meyer, V. V Gvaramadze, N. Langer, J. Mackey, P. Boumis, S. Mohamed

Abstract: 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 facto… ▽ More 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. △ Less

Submitted 5 December, 2013; originally announced December 2013.

Comments: 5 pages, 5 figures. Accepted for publication in MNRAS Letters

SALT reveals the barium central star of the planetary nebula Hen 2-39

Authors: B. Miszalski, H. M. J Boffin, D. Jones, A. I. Karakas, J. Köppen, A. A. Tyndall, S. S. Mohamed, P. Rodríguez-Gil, M. Santander-García

Abstract: 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 gi… ▽ More 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) △ Less

Submitted 20 September, 2013; originally announced September 2013.

Comments: Accepted for publication in MNRAS. 15 pages, 8 figures and 8 tables

ALMA observations of the variable 12CO/13CO ratio around the asymptotic giant branch star R Sculptoris

Authors: W. H. T. Vlemmings, M. Maercker, M. Lindqvist, S. Mohamed, H. Olofsson, S. Ramstedt, M. Brunner, M. A. T. Groenewegen, F. Kerschbaum, M. Wittkowski

Abstract: [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… ▽ More [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.... △ Less

Submitted 16 August, 2013; originally announced August 2013.

Comments: 6 pages, 5 figures, online data available at http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=J/A+A/556/L1

Journal ref: 2013 Astronomy & Astrophysics, 556, 1

The elusive origin of Carbon-Enhanced Metal-Poor stars

Authors: C. Abate, O. R. Pols, R. G. Izzard, S. S. Mohamed, S. E. de Mink

Abstract: 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 (BH… ▽ More 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. △ Less

Submitted 9 April, 2013; originally announced April 2013.

Comments: 6 pages, 2 figures, accepted for publication in PoS, XII International Symposium on Nuclei in the Cosmos (5-10 August 2012, Cairns, Australia)

3D Smoothed Particle Hydrodynamics Models of Betelgeuse's Bow Shock

Authors: Shazrene Mohamed, Jonathan Mackey, Norbert Langer

Abstract: 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… ▽ More 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. △ Less

Submitted 18 March, 2013; originally announced March 2013.

Comments: 9 pages, 4 figures, Betelgeuse workshop, November 2012, Paris. To be published in the European Astronomical Society Publications Series, editors: Pierre Kervella, Thibaut Le Bertre & Guy Perrin

Journal ref: EAS Publications Series, Volume 60, 2013, pp.243-251

Numerical models for the circumstellar medium around Betelgeuse

Authors: Jonathan Mackey, Shazrene Mohamed, Hilding R. Neilson, Norbert Langer, Dominique M. -A. Meyer

Abstract: 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 ind… ▽ More 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. △ Less

Submitted 8 March, 2013; originally announced March 2013.

Comments: 5 pages, 3 figures; to appear in proceedings of the Betelgeuse 2012 Workshop, Paris, Nov. 2012

Journal ref: EAS Publications Series, Volume 60, 2013, pp.253-259

Wind Roche-lobe overflow: Application to carbon-enhanced metal-poor stars

Authors: C. Abate, O. R. Pols, R. G. Izzard, S. S. Mohamed, S. E. de Mink

Abstract: 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 popu… ▽ More 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. △ Less

Submitted 18 February, 2013; originally announced February 2013.

Comments: 12 pages, 9 figures, accepted for publication in A&A

Unexpectedly large mass loss during the thermal pulse cycle of the red giant R Sculptoris!

Authors: M. Maercker, S. Mohamed, W. H. T. Vlemmings, S. Ramstedt, M. A. T. Groenewegen, E. Humphreys, F. Kerschbaum, M. Lindqvist, H. Olofsson, C. Paladini, M. Wittkowski, I. de Gregorio-Monsalvo, L. -A. Nyman

Abstract: 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 asymptot… ▽ More 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. △ Less

Submitted 10 October, 2012; originally announced October 2012.

Comments: Accepted by Nature

Double bow shocks around young, runaway red supergiants: application to Betelgeuse

Authors: Jonathan Mackey, Shazrene Mohamed, Hilding R. Neilson, Norbert Langer, Dominique M. -A. Meyer

Abstract: 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… ▽ More 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. △ Less

Submitted 17 April, 2012; originally announced April 2012.

Comments: 6 pages, 5 figures, emulateapj format, accepted to ApJ letters, simulation animations at http://www.astro.uni-bonn.de/~jmackey/research.html#Betelgeuse

Journal ref: 2012, ApJ, 751, L10

3D Simulations of Betelgeuse's Bow Shock

Authors: S. Mohamed, J. Mackey, N. Langer

Abstract: 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… ▽ More 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. △ Less

Submitted 21 March, 2012; v1 submitted 7 September, 2011; originally announced September 2011.

Comments: Minor revisions, replaced Fig. 1, 15, and 16, added movies. For a pdf version with higher resolution, see A&A: Forthcoming

Herschel's view into Mira's head

Authors: A. Mayer, A. Jorissen, F. Kerschbaum, S. Mohamed, S. Van Eck, R. Ottensamer, J. A. D. L. Blommaert, L. Decin, M. A. T. Groenewegen, Th. Posch, B. Vandenbussche, C. Waelkens

Abstract: 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… ▽ More 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. △ Less

Submitted 18 June, 2011; originally announced June 2011.

Comments: 4 pages

Journal ref: Astronomy & Astrophysics, Volume 531, L4 (July 2011)