Search

Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars). Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution data. To help interpret the polarized signal, we produced synthetic maps of light scattering by dust, through 3D radiative transfer simulations with the RADMC3D code. Results. The degree of linear polarization (DoLP) observed by ZIMPOL spreads across several optical filters. We infer that it primarily probes dust located just outside of the point spread function, and in or near the plane of the sky, with a total optical depth close to unity in the line of sight, representing only a fraction of the total circumstellar dust. The maximum DoLP ranges from 0.03-0.38 depending on the source, fractions that can be reproduced by our 3D pilot models for grains composed of common dust species. The spatial structure of the DoLP shows a diverse set of shapes. Only for three sources do we note a correlation between the ALMA CO and SiO lines, which trace the gas density, and the DoLP, which traces the dust. Conclusion. The clumpiness of the DoLP and the lack of a consistent correlation between the gas and the dust location show that, in the inner circumstellar environment (CSE), dust formation occurs at very specific sites. This has potential consequences for the derived mass-loss rates and dust-to-gas ratio in the inner region of the CSE. Except for $\pi^1$~Gru and perhaps GY Aql, we do not detect interactions between the circumstellar wind and the hypothesized companions that shape the wind at larger scales. This suggests that the orbits of any other companions are tilted out of the plane of the sky., Comment: Accepted for publication in Astronomy & Astrophysics. 22 pages, 15 figures, 5 tables

In some semi-detached binary systems, the donor star may transfer mass to the companion at a very high rate. We propose that, at sufficiently high mass-transfer rates such that the accretion disk around the companion becomes geometrically thick (or advection-dominated) near the disk outer radius, a large fraction of the transferred mass will be lost through the outer Lagrangian (L2) point, as a result of the excessive energy generated by viscous heating that cannot be efficiently radiated away. A physical model is constructed where the L2 mass loss fraction is given by the requirement that the remaining material in the disk has Bernoulli number equal to the L2 potential energy. Our model predicts significant L2 mass loss at mass-transfer rates exceeding a few times 10^{-4} Msun/yr. An equatorial circum-binary outflow (CBO) is formed in these systems. Implications for the orbital evolution and the observational appearance are discussed. In particular, (1) rapid angular momentum loss from the system tends to shrink the orbit and hence may increase the formation rate of mergers and gravitational-wave sources; (2) photons from the hot disk wind are reprocessed by the CBO into longer wavelength emission in the infrared bands, consistent with Spitzer observations of some ultra-luminous X-ray sources., MNRAS in press, 14 pages, 7 figures, plus appendix. Codes available here: https://github.com/wenbinlu/L2massloss

The nearby M dwarf WX UMa has recently been detected at radio wavelengths with LOFAR. The combination of its observed brightness temperature and circular polarisation fraction suggests that the emission is generated via the electron-cyclotron maser instability. Two distinct mechanisms have been proposed to power such emission from low-mass stars: either a sub-Alfv\'enic interaction between the stellar magnetic field and an orbiting planet, or reconnection at the edge of the stellar magnetosphere. In this paper, we investigate the feasibility of both mechanisms, utilising the information about the star's surrounding plasma environment obtained from modelling its stellar wind. Using this information, we show that a Neptune-sized exoplanet with a magnetic field strength of 10-100 G orbiting at ~0.034 au can accurately reproduce the observed radio emission from the star, with corresponding orbital periods of 7.4 days. Due to the stellar inclination, a planet in an equatorial orbit is unlikely to transit the star. While such a planet could induce radial velocity semi-amplitudes from 7 to 396 m s$^{-1}$, it is unlikely that this signal could be detected with current techniques due to the activity of the host star. The application of our planet-induced radio emission model here illustrates its exciting potential as a new tool for identifying planet-hosting candidates from long-term radio monitoring. We also develop a model to investigate the reconnection-powered emission scenario. While this approach produces less favourable results than the planet-induced scenario, it nevertheless serves as a potential alternative emission mechanism which is worth exploring further., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS

We present results from new self-consistent 3D MHD simulations of the magnetospheres from massive stars with a dipole magnetic axis that has a non-zero obliquity angle ($\beta$) to the star's rotation axis. As an initial direct application, we compare the global structure of co-rotating disks for nearly aligned ($\beta=5^o$) versus half-oblique ($\beta=45^o$) models, both with moderately rapid rotation ($\sim$ 0.5 critical). We find that accumulation surfaces broadly resemble the forms predicted by the analytic Rigidly Rotating Magnetosphere (RRM) model, but the mass buildup to near the critical level for centrifugal breakout against magnetic confinement distorts the field from the imposed initial dipole. This leads to an associated warping of the accumulation surface toward the rotational equator, with the highest density concentrated in {\em wings} centered on the intersection between the magnetic and rotational equators. These MHD models can be used to synthesize rotational modulation of photometric absorption and H$\alpha$ emission for a direct comparison with observations., Comment: Accepted for publication in MNRAS

V906 Carinae was one of the best observed novae of recent times. It was a prolific dust producer and harboured shocks in the early evolving ejecta outflow. Here, we take a close look at the consequences of these early interactions through study of high-resolution Ultraviolet and Visual Echelle spectrograph spectroscopy of the nebular stage and extrapolate backwards to investigate how the final structure may have formed. A study of ejecta geometry and shaping history of the structure of the shell is undertaken following a spectral line $\rm {\small SHAPE}$ model fit. A search for spectral tracers of shocks in the nova ejecta is undertaken and an analysis of the ionized environment. Temperature, density, and abundance analyses of the evolving nova shell are presented.

During the third all-sky survey (eRASS3), eROSITA, the soft X-ray instrument aboard Spectrum-Roentgen-Gamma, detected a new hard X-ray transient, eRASSt J040515.6-745202, in the direction of the Magellanic Bridge. We arranged follow-up observations and searched for archival data to reveal the nature of the transient. Using X-ray observations with XMM-Newton, NICER, and Swift, we investigated the temporal and spectral behaviour of the source for over about 10 days. The X-ray light curve obtained from the XMM-Newton observation with an 28 ks exposure revealed a type-I X-ray burst with a peak bolometric luminosity of at least 1.4e37 erg/s. The burst energetics are consistent with a location of the burster at the distance of the Magellanic Bridge. The relatively long exponential decay time of the burst of 70 s indicates that it ignited in a H-rich environment. The non-detection of the source during the other eROSITA surveys, twelve and six months before and six months after eRASS3, suggests that the burst was discovered during a moderate outburst which reached 2.6e36 erg/s in persistent emission. During the NICER observations, the source showed alternating flux states with the high level at a similar brightness as during the XMM-Newton observation. This behaviour is likely caused by dips as also seen during the last hour of the XMM-Newton observation. Evidence for a recurrence of the dips with a period of 21.8 hr suggests eRASSt J040515.6-745202 is a low-mass X-ray binary (LMXB) system with an accretion disk seen nearly edge on. We identify a multi-wavelength counterpart to the X-ray source in UVW1 and g, r, i, and z images obtained by the optical/UV monitor on XMM-Newton and the Dark Energy Camera at the Cerro Tololo Inter-American Observatory. (abbreviated), Comment: Accepted for publication in Astronomy & Astrophysics. 10 pages, 13 figures

Luminous Blue Variables (LBVs) are massive stars that are likely to be a transitionary phase between O stars and hydrogen-free classical Wolf–Rayet stars. The variability of these stars has been an area of study for both professional and amateur astronomers for more than a century. In this paper, we present 5 yr of precision photometry of the classical LBV P Cygni taken with the BRITE-Constellation nanosatellites. We have analyzed these data with Fourier analysis to search for periodicities that could elucidate the drivers of variability for these stars. These data show some long-time-scale variability over the course of all six calendar years of observations, but the frequencies needed to reproduce the individual light curves are not consistent from 1 yr to the next. These results likely show that there is no periodic phenomenon present for P Cygni, meaning that the variability is largely stochastic. We interpret the data as being caused by internal gravity waves similar to those seen in other massive stars, with P Cygni exhibiting a larger amplitude and lower characteristic frequency than the main-sequence or blue supergiant stars previously studied. These results show evidence that LBVs may be an extrapolation of the blue supergiants, which have previously been shown to be an extension of main-sequence stars in the context of the stochastic low-frequency photometric variability.

Fil: Ghoreyshi, Maziar. University of Western Ontario. Department of Physics and Astronomy. Ontario, Canada. Fil: Jones, Carol E. University of Western Ontario. Department of Physics and Astronomy. Ontario, Canada. Fil: Granada, Anahi. Universidad Nacional de Rio Negro. Centro Interdisciplinario de Telecomunicaciones, Electrónica, Computación y Ciencia Aplicada. Rio Negro, Argentina. Fil: Granada, Anahi. Consejo Nacional de Investigaciones Científicas y Técnicas. Buenos Aires, Argentina. Circumstellar discs around Be stars are formed by the material ejected by the central star. This process removes excess angular momentum from the star as viscosity facilitates the mass and angular momentum transfer within the disc and its growth. The angular momentum loss rates (AMLR) of Be stars is a subject of debate in the literature. Through the modelling of the disc formation and dissipation phases observed from Be stars, their average AMLR can be determined and this is the goal of this work. We use the viscous decretion disc (VDD) model to provide a range of the average AMLR for Be stars and compare these rates with predicted values from the literature. We explore the reasons for discrepancies between the predicted values of average AMLR, using the VDD and Geneva stellar evolution models that were previously reported in literature and find that the largest differences occur when Be stars are rotating below their critical speeds. We show that the time over which the mass reservoir builds up is inversely proportional to the average AMLR. Also, we determine a revised value of the average AMLR for the Galactic Be star ω CMa of 4.7×1036gcm2s−2⁠, which is in better agreement with the values expected for a typical B2-type star. Finally, the effect of disc truncation due to the presence of a companion star is investigated and we find that this has a minimal effect on the average AMLR. Los discos circumestelares alrededor de las estrellas Be se forman por el material expulsado por la estrella central. Este proceso elimina el exceso de momento angular de la estrella, ya que la viscosidad facilita la transferencia de masa y momento angular dentro del disco y su crecimiento. Las tasas de pérdida de momento angular (AMLR) de las estrellas Be son objeto de debate en la literatura. A través de la modelización de las fases de formación y disipación del disco observadas en estrellas Be, se puede determinar su AMLR media y este es el objetivo de este trabajo. Utilizamos el modelo de disco de decrecimiento viscoso (VDD) para proporcionar un rango de la AMLR media de las estrellas Be y comparamos estas tasas con los valores predichos en la literatura. Exploramos las razones de las discrepancias entre los valores predichos de AMLR media, utilizando el modelo VDD y los modelos de evolución estelar de Ginebra de los que se había informado previamente en la literatura, y descubrimos que las mayores diferencias se producen cuando las estrellas Be están rotando por debajo de sus velocidades críticas. Demostramos que el tiempo durante el cual se acumula el depósito de masa es inversamente proporcional a la AMLR media. Además, determinamos un valor revisado de la AMLR media para la estrella galáctica Be ω CMa de 4,7×1036gcm2s-2, que concuerda mejor con los valores esperados para una estrella típica de tipo B2. Por último, se investiga el efecto del truncamiento del disco debido a la presencia de una estrella compañera y se observa que tiene un efecto mínimo sobre la AMLR media.

e present the first 3D radiation-hydrodynamic simulations on the formation and evolution of born-again planetary nebulae (PNe), with particular emphasis to the case of HuBi 1, the inside-out PN. We use the extensively tested GUACHO code to simulate the formation of HuBi 1 adopting mass-loss and stellar wind terminal velocity estimates obtained from observations presented by our group. We found that, if the inner shell of HuBi 1 was formed by an explosive very late thermal pulse (VLTP) ejecting material with velocities of ∼300 km s−1, the age of this structure is consistent with that of ≃200 yr derived from multi-epoch narrow-band imaging. Our simulations predict that, as a consequence of the dramatic reduction of the stellar wind velocity and photon ionizing flux during the VLTP, the velocity and pressure structure of the outer H-rich nebula are affected creating turbulent ionized structures surrounding the inner shell. These are indeed detected in Gran Telescopio Canarias MEGARA optical observations. Furthermore, we demonstrate that the current relatively low ionizing photon flux from the central star of HuBi 1 is not able to completely ionize the inner shell, which favours previous suggestions that its excitation is dominated by shocks. Our simulations suggest that the kinetic energy of the H-poor ejecta of HuBi 1 is at least 30 times that of the clumps and filaments in the evolved born-again PNe A 30 and A 78, making it a truly unique VLTP event. © 2021 The Author(s)., JAT thanks funding by Fundación Marcos Moshinsky (Mexico) and the Direcciónn General de Asuntos del Personal Académico (DGAPA) of the Universidad Nacional Autónoma de México (UNAM) project IA100720. VL gratefully acknowledges support from the CONACYT Research Fellowship program. BMM and MAG acknowledge support of the Spanish Ministerio de Ciencia, Innovación y Universidades (MCIU) grant PGC2018-102184-B-I00. AE acknowledges support from DGAPA-PAPIIT (UNAM) grant IN 109518. This work has made extensive use of NASA’s Astrophysics Data System., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

We focus on post-asymptotic giant branch (post-AGB) binaries and study the interaction between the different components of these complex systems. These components comprise the post-AGB primary, a main sequence secondary, a circumbinary disk, as well as a fast bipolar outflow (jet) launched by the companion. We obtained well-sampled time series of high resolution optical spectra over the last decade and these spectra provide the basis of our study. The jet is detected in absorption, at superior conjunction, when the line of sight towards the primary goes through the bipolar cone. Our spectral time series scan the jets during orbital motion. Our spatio-kinematic model is constrained by these dynamical spectra. We complement this with a radiative-transfer model in which the Balmer series are used to derive total mass-loss rates in the jets. The jets are found to be wide and display an angle-dependent density structure with a dense and slower outer region near the jet cone and a fast inner part along the jet symmetry axes. The deprojected outflow velocities confirm that the companions are main sequence companions. The total mass-loss rates are large (10^{-8} and 10^{-5}\,solar mass per year), from which we can infer that the mass-accretion rates onto the companion star must be high as well. The circumbinary disk is likely the main source for the accretion disk around the companion. All systems with full disks that start near the sublimation radius show jets, whereas for systems with evolved transition disks, this lowers to a detection rate of 50%. Objects without an infrared excess do not show jets. We conclude that jet creation in post-AGB binaries is a mainstream process. The interaction between the circumbinary disks and the central binary provide the needed accretion flow, but the presence of a circumbinary disk does not seem to be the only prerequisite to launch a jet., 17 pages, 12 figures, accepted by Astronomy and Astrophysics

Hot, massive (OB) stars experience strong line-driven stellar winds and mass loss. As the majority of efficient driving lines are metallic, the amount of wind driving and mass loss is dependent on the stellar metallicity Z. In addition, line-driven winds are intrinsically inhomogeneous and clumpy. However, to date, neither theoretical nor empirical studies of line-driven winds have investigated how such wind clumping may also depend on Z. We theoretically investigated the degree of wind clumping due to the line-deshadowing instability (LDI) as a function of Z. We performed two-dimensional hydrodynamic simulations of the LDI with an assumed one-dimensional radiation line force for a grid of O-star wind models with fixed luminosity, but with different metal contents by varying the accumulative line strength Qbar describing the total ensemble of driving lines. We find that, for this fixed luminosity, the amount of wind clumping decreases with metallicity. The decrease is clearly seen in the statistical properties of our simulations, but is nonetheless rather weak; a simple power-law fit for the dependence of the clumping factor f_cl = /^2 on metallicity yields f_cl ~ Z^(0.15 +/- 0.01). This implies that empirically derived power-law dependencies of mass-loss rate Mdot on metallicity -- which were previously inferred from spectral diagnostics effectively depending on Mdot*sqrt(f_cl) but without having any constraints on f_cl(Z) -- should be only modestly altered by clumping. We expect that this prediction can be directly tested using new data from the Hubble Space Telescope Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) project., Comment: accepted for publication in A&A

Mass loss remains a major uncertainty in stellar modelling. In low-mass stars, mass loss is most significant on the red giant branch (RGB), and will impact the star's evolutionary path and final stellar remnant. Directly measuring the mass difference of stars in various phases of evolution represents one of the best ways to quantify integrated mass loss. Globular clusters (GCs) are ideal objects for this. M4 is currently the only GC for which asteroseismic data exists for stars in multiple phases of evolution. Using K2 photometry, we report asteroseismic masses for 75 red giants in M4, the largest seismic sample in a GC to date. We find an integrated RGB mass loss of $\Delta\bar{M} = 0.17 \pm 0.01 ~\mathrm{M}_{\odot}$, equivalent to a Reimers' mass-loss coefficient of $\eta_R = 0.39$. Our results for initial mass, horizontal branch mass, $\eta_R$, and integrated RGB mass loss show remarkable agreement with previous studies, but with higher precision using asteroseismology. We also report the first detections of solar-like oscillations in early asymptotic giant branch (EAGB) stars in GCs. We find an average mass of $\bar{M}_{\text{EAGB}}=0.54 \pm 0.01 ~\mathrm{M}_{\odot}$, significantly lower than predicted by models. This suggests larger-than-expected mass loss on the horizontal branch. Alternatively, it could indicate unknown systematics in seismic scaling relations for the EAGB. We discover a tentative mass bi-modality in the RGB sample, possibly due to the multiple populations. In our red horizontal branch sample, we find a mass distribution consistent with a single value. We emphasise the importance of seismic studies of GCs since they could potentially resolve major uncertainties in stellar theory., Comment: 15 pages, 12 figures, accepted by MNRAS

We present a new asteroseismic analysis of the stars in the Globular Cluster (GC) M4 based on the data collected by the K2 mission. We report the detection of solar-like oscillation in 37 stars, 32 red giant branch (RGB) and 6 red horizontal branch (rHB) stars, the largest sample for this kind of study in GC up to date. Combining information from asteroseismology and multi-band photometry we estimate both the masses and the radii of our targets. Our estimates are in agreement with independent sources, serving as a crucial verification of asteroseismology in the low metallicity regime. As M4 is an old GC, it hosts multiple stellar populations differing in light-element abundances and in helium mass fraction. This generates a mass difference between the populations along the RGB, which in the case of M4 is estimated to be $0.017 M_\odot$. With this wealth of information we can assign population membership and estimate the average mass of the stellar populations, but the current uncertainties do not allow us to resolve this mass difference. The population membership and the seismic data of RGB and HB stars, allow us, however, to assess the integrated mass loss along the RGB of the first generation stars in the cluster. We obtain $\rm \Delta M=0.227 \pm0.028 M_\odot$, in good agreement with independent estimates. Finally, we observe the presence of a statistically significant mass-temperature gradient in the rHB stars. This represents the first direct, model-independent observation of the colour-temperature-mass correlation predicted by the theory., Comment: 14 pages, 11 figures, 3 tables. Accepted for publication in A&A

Mass loss by red giants is an important process to understand the final stages of stellar evolution and the chemical enrichment of the interstellar medium. Mass-loss rates are thought to be controlled by pulsation-enhanced dust-driven outflows. Here we investigate the relationships between mass loss, pulsations, and radiation, using 3213 luminous Kepler red giants and 135000 ASAS-SN semiregulars and Miras. Mass-loss rates are traced by infrared colours using 2MASS and WISE and by observed-to-model WISE fluxes, and are also estimated using dust mass-loss rates from literature assuming a typical gas-to-dust mass ratio of 400. To specify the pulsations, we extract the period and height of the highest peak in the power spectrum of oscillation. Absolute magnitudes are obtained from the 2MASS Ks band and the Gaia DR2 parallaxes. Our results follow. (i) Substantial mass loss sets in at pulsation periods above ~60 and ~100 days, corresponding to Asymptotic-Giant-Branch stars at the base of the period-luminosity sequences C' and C. (ii) The mass-loss rate starts to rapidly increase in semiregulars for which the luminosity is just above the Red-Giant-Branch tip and gradually plateaus to a level similar to that of Miras. (iii) The mass-loss rates in Miras do not depend on luminosity, consistent with pulsation-enhanced dust-driven winds. (iv) The accumulated mass loss on the Red Giant Branch consistent with asteroseismic predictions reduces the masses of red-clump stars by 6.3%, less than the typical uncertainty on their asteroseismic masses. Thus mass loss is currently not a limitation of stellar age estimates for galactic archaeology studies., 14 pages, 12 figures, accepted for publication in MNRAS

Observations of rapidly-rotating cool stars often show coronal slingshot prominences that remove mass and angular momentum when they are ejected. The derived masses of these prominences show a scatter of some two orders of magnitude. In order to investigate if this scatter could be intrinsic, we use a full magnetic cycle of solar magnetograms to model the coronal structure and prominence distribution in a young Sun, where we scale the field strength in the magnetograms with angular velocity according to $ B \propto \Omega^{-1.32} $. We reproduce both the observed prominence masses and their scatter. We show that both the field strength and the field geometry contribute to the prominence masses that can be supported and to the rate at which they are ejected. Predicted prominence masses follow the magnetic cycle, but with half the period, peaking both at cycle maximum and at cycle minimum. We show that mass loss rates in prominences are less than those predicted for the stellar wind. We also investigate the role of small-scale field that may be unresolved in typical stellar magnetograms. This provides only a small reduction in the predicted total prominence mass, principally by reducing the number of large magnetic loops that can support slingshot prominences. We conclude that the observed scatter in prominence masses can be explained by underlying magnetic cycles., Comment: 10 pages, 7 figures. Accepted for publication in MNRAS

Context. The fate of a massive star during the latest stages of its evolution is highly dependent on its mass-loss rate/geometry and therefore knowing the geometry of the circumstellar material close to the star and its surroundings is crucial. Aims. We aim to study the nature (i.e. geometry, rates) of mass-loss episodes. In this context, yellow hypergiants are great targets. Methods. We analyse a large set of optical/near-infrared data, in spectroscopic and photometric (X-shooter/VLT), spectropolarimetric (ISIS/WHT), and interferometric GRAVITY-AMBER/VLTI) modes, toward the yellow hypergiant IRAS 17163-3907. We present the first model-independent reconstructed images of IRAS 17163-3907 at these wavelengths at milli-arcsecond scales. Lastly, we apply a 2D radiative transfer model to fit the dereddened photometry and the radial profiles of published VISIR images at 8.59 {\mu}m, 11.85 {\mu}m and 12.81 {\mu}m simultaneously, adopting the revised Gaia distance (DR2). Results. The interferometric observables around 2 {\mu}m show that the Br{\gamma} emission is more extended and asymmetric than the Na i and the continuum emission. In addition to the two known shells surrounding IRAS 17163-3907 we report on the existence of a third hot inner shell with a maximum dynamical age of only 30 yr. Conclusions. The interpretation of the presence of Na i emission at closer distances to the star compared to Br{\gamma} has been a challenge in various studies. We argue that the presence of a pseudophotosphere is not needed, but it is rather an optical depth effect. The three observed distinct mass-loss episodes are characterised by different mass-loss rates and can inform the theories on mass-loss mechanisms, which is a topic still under debate. We discuss these in the context of photospheric pulsations and wind bi-stability mechanisms., Comment: 26 pages, 22 Figures, 6 tables, accepted for publication in A&A. Original abstract will be available after publication

The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $\mu$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust., Comment: 20 pages, 13 figures, 5 tables, accepted for publications in MNRAS

Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind properties of the massive stars in R136, in particular the parameters related to wind clumping. Methods: We simultaneously analyse optical and UV spectroscopy of 53 O-type and 3 WNh-stars using the FASTWIND model atmosphere code and a genetic algorithm. The models account for optically thick clumps and effects related to porosity and velocity-porosity, as well as a non-void interclump medium. Results: We obtain stellar parameters, surface abundances, mass-loss rates, terminal velocities and clumping characteristics and compare these to theoretical predictions and evolutionary models. The clumping properties include the density of the interclump medium and the velocity-porosity of the wind. For the first time, these characteristics are systematically measured for a wide range of effective temperatures and luminosities. Conclusions: We confirm a cluster age of 1.0-2.5 Myr and derive an initial stellar mass of $\geq 250 {\rm M}_\odot$ for the most massive star in our sample, R136a1. The winds of our sample stars are highly clumped, with an average clumping factor of $f_{\rm cl}=29\pm15$. We find tentative trends in the wind-structure parameters as a function of mass-loss rate, suggesting that the winds of stars with higher mass-loss rates are less clumped. We compare several theoretical predictions to the observed mass-loss rates and terminal velocities and find that none satisfactorily reproduces both quantities. The prescription of Krti\v{c}ka & Kub\'at (2018) matches best the observed mass-loss rates., Comment: Accepted for publication in A&A; Appendix I will not be included in the published version

12 pags., 10 figs. 2 tabs., Red supergiant stars (RSGs, Minit = 10-40Msun) are known to eject large amounts of material, as much as half of their initial mass during this evolutionary phase. However, the processes powering the mass ejection in low- and intermediate-mass stars do not work for RSGs and the mechanism that drives the ejection remains unknown. Different mechanisms have been proposed as responsible for this mass ejection but so far little is known about the actual processes taking place in these objects. Here we present high angular resolution interferometric ALMA maps of VY CMa continuum and molecular emission, which resolve the structure of the ejecta with unprecedented detail. The study of the molecular emission from the ejecta around evolved stars has been shown to be an essential tool in determining the characteristics of the mass loss ejections. Our aim is thus to use the information provided by these observations to understand the ejections undergone by VY CMa and to determine their possible origins. We inspected the kinematics of molecular emission observed. We obtained position-velocity diagrams and reconstructed the 3D structure of the gas traced by the different species. It allowed us to study the morphology and kinematics of the gas traced by the different species surrounding VY CMa. Two types of ejecta are clearly observed: extended, irregular, and vast ejecta surrounding the star that are carved by localized fast outflows. The structure of the outflows is found to be particularly flat. We present a 3D reconstruction of these outflows and proof of the carving. This indicates that two different mass loss processes take place in this massive star. We tentatively propose the physical cause for the formation of both types of structures. These results provide essential information on the mass loss processes of RSGs and thus of their further evolution., The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013) / ERC Grant Agreement n. 610256 (NANOCOSMOS). We would also like to thank the Spanish MINECO for funding support from grants CSD2009-00038, AYA2012-32032. M.A. also thanks for funding support from the Ramón y Cajal programme of Spanish MINECO (RyC2014-16277). This publication is part of the “I+D+i” (research, development, and innovation) projects PID2019-107115GB-C21, and PID2019-106110GBI00, and PID2020-117034RJ-I00, supported by the Spanish “Ministerio de Ciencia e Innovación” MCIN/AEI/10.13039/501100011033.

We present the first instalment of a deep imaging catalogue containing 58 True, Likely, and Possible extended PNe detected with the Isaac Newton Telescope Photometric H α Survey (IPHAS). The three narrow-band filters in the emission lines of H α, [N ii] λ6584 Å, and [O iii] λ5007 Å used for this purpose allowed us to improve our description of the morphology and dimensions of the nebulae. In some cases even the nature of the source has been reassessed. We were then able to unveil new macro- and micro-structures, which will without a doubt contribute to a more accurate analysis of these PNe. It has been also possible to perform a primary classification of the targets based on their ionization level. A Deep Learning classification tool has also been tested. We expect that all the PNe from the IPHAS catalogue of new extended planetary nebulae will ultimately be part of this deep H α, [N ii], and [O iii] imaging catalogue. © 2021 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., LS acknowledges support from PAPIIT grant IN-101819 (Mexico). MAG acknowledges support from grant AYA PGC2018-102184-B-I00 co-funded with FEDER funds. GR-L acknowledges support from CONACyT (grant 263373) and PRODEP (Mexico). JAT acknowledges funding by Dirección General de Asuntos del Personal Academico of the Universidad Nacional Autónoma de México (DGAPA, UNAM) project IA100720 and the Marcos Moshinsky Fundation (Mexico). AAZ acknowledges funding from STFC under grant number ST/T000414/1 and Newton grant ST/R006768/1, and from a Hung Hing Ying visiting professorship at the University of Hong Kong. DNFAI acknowledges funding under ‘Deep Learning for Classification of Astronomical Archives’ from the Newton-Ungku Omar Fund: F08/STFC/1792/2018., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

We investigate the impact of optically thick clumping on stellar wind diagnostics in O supergiants and constrain wind parameters associated with porosity in velocity space. This is the first time the effects of optically thick clumping have been investigated for a sample of massive hot stars, using models including a full optically thick clumping description. We re-analyse spectroscopic observations of a sample of eight O supergiants. Using a genetic algorithm wrapper around the NLTE atmosphere code FASTWIND we obtain simultaneous fits to optical and UV spectra and determine photospheric and wind properties and surface abundances. We provide empirical constraints on a number of wind parameters including the clumping factors, mass-loss rates and terminal wind velocities. Additionally, we establish the first systematic empirical constraints on velocity filling factors and interclump densities. These parameters describe clump distribution in velocity-space and density of the interclump medium in physical-space, respectively. We observe a mass-loss rate reduction of a factor of 3.6 compared to theoretical predictions from Vink et al. (2000), and mass-loss rates within a factor 1.4 of predictions from Bj\"orklund et al. (2021). We confirm that including optically thick clumping allows simultaneous fitting of recombination lines and resonance lines, including the unsaturated UV phosphorus lines (Pv 1118-1128), without reducing the phosphorus abundance. We find that, on average, half of the wind velocity field is covered by dense clumps. We also find that these clumps are 25 times denser than the average wind, and that the interclump medium is 3-10 times less dense than the mean wind. The former result agrees well with theoretical predictions, the latter suggests that lateral filling-in of radially compressed gas might be critical for setting the scale of the rarefied interclump matter., Comment: Accepted for publication in A&A; 35 pages, 20 figures

25 pags., 15 figs., 4 tabs., Aims. We aim to constrain the sizes of, and investigate deviations from spherical symmetry in, the CO circumstellar envelopes (CSEs) of 16 S-type stars, along with an additional 7 and 4 CSEs of C-type and M-type AGB stars, respectively. Methods. We map the emission from the CO J = 2-1 and 3-2 lines observed with the Atacama Compact Array (ACA) and its total power (TP) antennas, and fit with a Gaussian distribution in the uv- and image planes for ACA-only and TP observations, respectively. The major axis of the fitted Gaussian for the CO(2-1) line data gives a first estimate of the size of the CO-line-emitting CSE. We investigate possible signs of deviation from spherical symmetry by analysing the line profiles and the minor-to-major axis ratio obtained from visibility fitting, and by investigating the deconvolved images. Results. The sizes of the CO-line-emitting CSEs of low-mass-loss-rate (low-MLR) S-type stars fall between the sizes of the CSEs of C-stars, which are larger, and those of M-stars, which are smaller, as expected because of the differences in their respective CO abundances and the dependence of the photodissociation rate on this quantity. The sizes of the low-MLR S-type stars show no dependence on circumstellar density, as measured by the ratio of the MLR to terminal outflow velocity, irrespective of variability type. The density dependence steepens for S-stars with higher MLRs. While the CO(2-1) brightness distribution size of the low-density S-stars is in general smaller than the predicted photodissociation radius (assuming the standard interstellar radiation field), the measured size of a few of the high-density sources is of the same order as the expected photodissociation radius. Furthermore, our results show that the CO CSEs of most of the S-stars in our sample are consistent with a spherically symmetric and smooth outflow. For some of the sources, clear and prominent asymmetric features are observed which are indicative of intrinsic circumstellar anisotropy. Conclusions. As the majority of the S-type CSEs of the stars in our sample are consistent with a spherical geometry, the CO envelope sizes obtained in this paper will be used to constrain detailed radiative transfer modelling to directly determine more accurate MLR estimates for the stars in our sample. For several of our sources that present signs of deviation from spherical symmetry, further high-resolution observations would be necessary to investigate the nature of, and the physical processes behind, these asymmetrical structures. This will provide further insight into the mass-loss process and its related chemistry in S-type AGB stars., This paper makes use of the following ALMA data: ADS/JAO.ALMA#2018.1.01434.S; ADS/JAO.ALMA#2017.1.00595.S and ADS/JAO.ALMA#2012.1.00524.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. MA acknowledges support from the Nordic ALMA Regional Centre (ARC) node based at Onsala Space Observatory. The Nordic ARC node is funded through Swedish Research Council grant No 2017-00648. This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme under grant agreements No. 883867 [EXWINGS] and 730562 [RadioNet].

Cay, Y. Z., et al., We present the spectroscopic and photometric study of five intermediate-luminosity red transients (ILRTs), namely AT 2010dn, AT 2012jc, AT 2013la, AT 2013lb, and AT 2018aes. They share common observational properties and belong to a family of objects similar to the prototypical ILRT SN 2008S. These events have a rise time that is less than 15 days and absolute peak magnitudes of between-11.5 and-14.5 mag. Their pseudo-bolometric light curves peak in the range 0.5-9.0 × 1040 erg s-1 and their total radiated energies are on the order of (0.3-3) × 1047 erg. After maximum brightness, the light curves show a monotonic decline or a plateau, resembling those of faint supernovae IIL or IIP, respectively. At late phases, the light curves flatten, roughly following the slope of the 56Co decay. If the late-time power source is indeed radioactive decay, these transients produce 56Ni masses on the order of 10-4 to 10-3 M⊙. The spectral energy distribution of our ILRT sample, extending from the optical to the mid-infrared (MIR) domain, reveals a clear IR excess soon after explosion and non-negligible MIR emission at very late phases. The spectra show prominent H lines in emission with a typical velocity of a few hundred km s-1, along with Ca II features. In particular, the [Ca II] λ7291,7324 doublet is visible at all times, which is a characteristic feature for this family of transients. The identified progenitor of SN 2008S, which is luminous in archival Spitzer MIR images, suggests an intermediate-mass precursor star embedded in a dusty cocoon. We propose the explosion of a super-asymptotic giant branch star forming an electron-capture supernova as a plausible explanation for these events.

Context. Eta Carinae is a highly eccentric, massive binary system (semimajor axis ~15.5 au) with powerful stellar winds and a phase-dependent wind-wind collision (WWC) zone. The primary star, η Car A, is a luminous blue variable (LBV); the secondary, η Car B, is a Wolf-Rayet or O star with a faster but less dense wind. Aperture-synthesis imaging allows us to study the mass loss from the enigmatic LBV η Car. Understanding LBVs is a crucial step toward improving our knowledge about massive stars and their evolution. Aims. Our aim is to study the intensity distribution and kinematics of η Car’s WWC zone. Methods. Using the VLTI-MATISSE mid-infrared interferometry instrument, we perform Brα imaging of η Car’s distorted wind. Results. We present the first VLTI-MATISSE aperture-synthesis images of η Car A’s stellar windin several spectral channels distributed across the Brα 4.052 μm line (spectral resolving power R ~ 960). Our observations were performed close to periastron passage in February 2020 (orbital phase ~ 14.0022). The reconstructed iso-velocity images show the dependence of the primary stellar wind on wavelength or line-of-sight (LOS) velocity with a spatial resolution of 6 mas (~14 au). The radius of the faintest outer wind regions is ~26 mas (~60 au). At several negative LOS velocities, the primary stellar wind is less extended to the northwest than in other directions. This asymmetry is most likely caused by the WWC. Therefore, we see both the velocity field of the undisturbed primary wind and the WWC cavity. In continuum spectral channels, the primary star wind is more compact than in line channels. A fit of the observed continuum visibilities with the visibilities of a stellar wind CMFGEN model (CMFGEN is an atmosphere code developed to model the spectra of a variety of objects) provides a full width at half maximum fit diameter of the primary stellar wind of 2.84 ± 0.06 mas (6.54 ± 0.14 au). We comparethe derived intensity distributions with the CMFGEN stellar wind model and hydrodynamic WWC models.

There have recently been detections of radio emission from low-mass stars, some of which are indicative of star-planet interactions. Motivated by these exciting new results, in this paper we present Alfv\'en wave-driven stellar wind models of the two active planet-hosting M dwarfs Prox Cen and AU Mic. Our models incorporate large-scale photospheric magnetic field maps reconstructed using the Zeeman-Doppler Imaging method. We obtain a mass-loss rate of $0.25~\dot{M}_{\odot}$ for the wind of Prox Cen. For the young dwarf AU Mic, we explore two cases: a low and high mass-loss rate. Depending on the properties of the Alfv\'en waves which heat the corona in our wind models, we obtain mass-loss rates of $27$ and $590~\dot{M}_{\odot}$ for AU Mic. We use our stellar wind models to assess the generation of electron cyclotron maser instability emission in both systems, through a mechanism analogous to the sub-Alfv\'enic Jupiter-Io interaction. For Prox Cen we do not find any feasible scenario where the planet can induce radio emission in the star's corona, as the planet orbits too far from the star in the super-Alfv\'enic regime. However, in the case that AU Mic has a stellar wind mass-loss rate of $27~\dot{M}_{\odot}$, we find that both planets b and c in the system can induce radio emission from $\sim10$ MHz to 3 GHz in the corona of the host star for the majority of their orbits, with peak flux densities of $\sim10$ mJy. Detection of such radio emission would allow us to place an upper limit on the mass-loss rate of the star., Comment: 8 pages, 5 figures. Accepted for publication in MNRAS

We present a case study of using a novel method to identify red supergiant (RSG) candidates in NGC 6822, based on their 1.6 $\mu$m "H-bump". We collected 32 bands of photometric data for NGC 6822 ranging from optical to MIR. By using the theoretical spectra from MARCS, we demonstrate that there is a prominent difference around 1.6 $\mu$m ("H-bump") between low-surface-gravity (LSG) and high-surface-gravity (HSG) targets. Taking advantage of this feature, we identify efficient color-color diagrams (CCDs) of rzH and rzK to separate HSG and LSG targets from crossmatching of optical and NIR data. Moreover, synthetic photometry from ATLAS9 also give similar results. Further separating RSG candidates from the rest of the LSG candidates is done by using semi-empirical criteria on NIR CMDs and resulted in 323 RSG candidates. Meanwhile, the simulation of foreground stars from Besan\c{c}on models also indicates that our selection criteria is largely free from the contamination of Galactic giants. In addition to the "H-bump" method, we also use the traditional BVR method as a comparison and/or supplement, by applying a slightly aggressive cut to select as much as possible RSG candidates (358 targets). Furthermore, the Gaia astrometric solution is used to constrain the sample, where 181 and 193 targets were selected from the "H-bump" and BVR method, respectively. The percentages of selected targets in both methods are similar as $\sim$60\%, indicating the comparable accuracy of the two methods. In total, there are 234 RSG candidates after combining targets from both methods with 140 ($\sim$60\%) of them in common. The final RSG candidates are in the expected locations on the MIR CMDs, while the spatial distribution is also coincident with the FUV-selected star formation regions, suggesting the selection is reasonable and reliable., Comment: A&A accepted

Jets are a commonly observed phenomenon in post-asymptotic giant branch (post-AGB) binaries. Due to the orbital motion of the binary, the jet causes variable absorption in the Balmer profiles. In previous work, we have developed spatio-kinematic and radiative transfer models to reproduce the observed Balmer line variability and derive the spatio-kinematic structure of the jet and its mass-loss rate. Here, we apply our jet model to five post-AGB binaries with distinct H{\alpha} line variability and diverse orbital properties. Our models fit the H{\alpha} line variations very well. We estimate jet mass-loss rates between 10-8 Mdot yr-1 and 10-4 Mdot yr-1 , from which we deduce accretion rates onto the companion between 10-7 Mdot yr-1 and 10-3 Mdot yr-1 . These accretion rates are somewhat higher than can be comfortably explained with reasonable sources of accretion, but we argue that the circumbinary disc in these systems is most-likely the source feeding the accretion, although accretion from the post-AGB star cannot be ruled out. The diversity of the variability in the five objects is due to their wide ejection cones combined with a range of viewing angles, rather than inherent differences between the objects. The nature of the observations does not let us easily distinguish which jet launching model (stellar jet, disc wind, or X-wind) should be favoured. In conclusion, we show that our jet model includes the physical parameters to successfully reproduce the H{\alpha} line variations and retrieve the structure and mass-loss rates of the jet for all five objects that are representative of the diverse sample of Galactic post-AGB binaries., Comment: 20 pages, 9 figures, accepted for publication in MNRAS

Accurate mass-loss rates and terminal velocities from massive stars winds are essential to obtain synthetic spectra from radiative transfer calculations and to determine the evolutionary path of massive stars. From a theoretical point of view, analytical expressions for the wind parameters and velocity profile would have many advantages over numerical calculations that solve the complex non-linear set of hydrodynamic equations. In a previous work, we obtained an analytical description for the fast wind regime. Now, we propose an approximate expression for the line-force in terms of new parameters and obtain a velocity profile closed-form solution (in terms of the Lambert $W$ function) for the $\delta$-slow regime. Using this analytical velocity profile, we were able to obtain the mass-loss rates based on the m-CAK theory. Moreover, we established a relation between this new set of line-force parameters with the known stellar and m-CAK line-force parameters. To this purpose, we calculated a grid of numerical hydrodynamical models and performed a multivariate multiple regression. The numerical and our descriptions lead to good agreement between their values., Comment: 8 pages, 1 figure, accepted for publication in MNRAS

[Abridged] Ensemble studies of red-giant stars with exquisite asteroseismic, spectroscopic, and astrometric constraints offer a novel opportunity to recast and address long-standing questions concerning the evolution of stars and of the Galaxy. Here, we infer masses and ages for nearly 5400 giants with available Kepler light curves and APOGEE spectra, and discuss some of the systematics that may affect the accuracy of the inferred stellar properties. First, we look at age-chemical-abundances relations. We find a dearth of young, metal-rich stars, and the existence of a significant population of old (8-9 Gyr), low-[$\alpha$/Fe], super-solar metallicity stars, reminiscent of the age and metallicity of the well-studied open cluster NGC6791. The age-chemo-kinematic properties of these stars indicate that efficient radial migration happens in the thin disk. We find that ages and masses of the nearly 400 $\alpha$-element-rich red-giant-branch (RGB) stars in our sample are compatible with those of an old (~11 Gyr), nearly coeval, chemical-thick disk population. Using a statistical model, we show that 95% of the population was born within ~1.5 Gyr. Moreover, we find a difference in the vertical velocity dispersion between low- and high-[$\alpha$/Fe] populations, confirming their different chemo-dynamical histories. We then exploit the almost coeval $\alpha$-rich population to gain insight into processes that may have altered the mass of a star along its evolution, which are key to improve the mapping of the observed stellar mass to age. We find evidence for a mean integrated RGB mass loss = 0.10 $\pm$ 0.02 Msun and that the occurrence of massive (M $\gtrsim$ 1.1 Msun) $\alpha$-rich stars is of the order of 5% on the RGB, and significantly higher in the RC, supporting the scenario in which most of these stars had undergone interaction with a companion., Comment: accepted for publication in A&A, 26 pages, 24 figures, catalogue available via cds

Whilst ``slingshot'' prominences have been observed on M-dwarfs, most if not all theoretical studies have focused on solar-like stars. We present an investigation into stellar prominences around rapidly rotating young M-dwarfs. We have extrapolated the magnetic field in the corona from Zeeman-Doppler maps and determined the sites of mechanical stability where prominences may form. We analyse the prominence mass that could be supported and the latitude range over which this material is distributed. We find that for these maps, much of this prominence mass may be invisible to observation - typically, 13 pages, 10 figures

During the first three observing runs of the Advanced gravitational-wave detector network, the LIGO/Virgo collaboration detected several black hole binary (BHBH) mergers. As the population of detected BHBH mergers grows, it will become possible to constrain different channels for their formation. Here we consider the chemically homogeneous evolution (CHE) channel in close binaries, by performing population synthesis simulations that combine realistic binary models with detailed cosmological calculations of the chemical and star-formation history of the Universe. This allows us to constrain population properties, as well as cosmological and aLIGO/aVirgo detection rates of BHBH mergers formed through this pathway. We predict a BHBH merger rate at redshift zero of $5.8 \, \textrm{Gpc}^{-3} \textrm{yr}^{-1}$ through the CHE channel, to be compared with aLIGO/aVirgo's measured rate of ${53.2}_{-28.2}^{+55.8} \, \text{Gpc}^{-3} \text{yr}^{-1}$, and find that eventual merger systems have BH masses in the range $17 - 43 \, \textrm{M}_{\odot}$ below the pair-instability supernova (PISN) gap, and $>124 \, \textrm{M}_{\odot}$ above the PISN gap. We investigate effects of momentum kicks during black hole formation, and calculate cosmological and magnitude limited PISN rates. We also study the effects of high-redshift deviations in the star formation rate. We find that momentum kicks tend to increase delay times of BHBH systems, and our magnitude limited PISN rate estimates indicate that current deep surveys should be able to detect such events. Lastly, we find that our cosmological merger rate estimates change by at most $\sim 8\%$ for mild deviations of the star formation rate in the early Universe, and by up to $\sim 40\%$ for extreme deviations., 20 pages, 18 figures, 4 tables. References added, text clarifications added, results and figures unchanged

28 pags., 10 figs., 7 tabs., Context. This is the first publication from the DEATHSTAR project. The overall goal of the project is to reduce the uncertainties of the observational estimates of mass-loss rates from evolved stars on the Asymptotic Giant Branch (AGB). Aim. The aim in this first publication is to constrain the sizes of the 12CO emitting region from the circumstellar envelopes around 42 mostly southern AGB stars, of which 21 are M-Type and 21 are C-Type, using the Atacama Compact Array (ACA) at the Atacama Large Millimeter/submillimeter Array. The symmetry of the outflows is also investigated. Methods. Line emission from 12CO J = 2→1 and 3→2 from all of the sources were mapped using the ACA. In this initial analysis, the emission distribution was fit to a Gaussian distribution in the uv-plane. A detailed radiative transfer analysis will be presented in a future publication. The major and minor axis of the best-fit Gaussian at the line center velocity of the 12CO J = 2→1 emission gives a first indication of the size of the emitting region. Furthermore, the fitting results, such as the Gaussian major and minor axis, center position, and the goodness of fit across both lines, constrain the symmetry of the emission distribution. For a subsample of sources, the measured emission distribution is compared to predictions from previous best-fit radiative transfer modeling results. Results. We find that the CO envelope sizes are, in general, larger for C-Type than for M-Type AGB stars, which is as expected if the CO/H2 ratio is larger in C-Type stars. Furthermore, the measurements show a relation between the measured (Gaussian) 12CO J = 2→1 size and circumstellar density that, while in broad agreement with photodissociation calculations, reveals large scatter and some systematic differences between the different stellar types. For lower mass-loss-rate irregular and semi-regular variables of both M-and C-Type AGB stars, the 12CO J = 2→1 size appears to be independent of the ratio of the mass-loss rate to outflow velocity, which is a measure of circumstellar density. For the higher mass-loss-rate Mira stars, the 12CO J = 2→1 size clearly increases with circumstellar density, with larger sizes for the higher CO-Abundance C-Type stars. The M-Type stars appear to be consistently smaller than predicted from photodissociation theory. The majority of the sources have CO envelope sizes that are consistent with a spherically symmetric, smooth outflow, at least on larger scales. For about a third of the sources, indications of strong asymmetries are detected. This is consistent with what was found in previous interferometric investigations of northern sources. Smaller scale asymmetries are found in a larger fraction of sources. Conclusions. These results for CO envelope radii and shapes can be used to constrain detailed radiative transfer modeling of the same stars so as to determine mass-loss rates that are independent of photodissociation models. For a large fraction of the sources, observations at higher spatial resolution will be necessary to deduce the nature and origin of the complex circumstellar dynamics revealed by our ACA observations., This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.2.00025.S and #2017.1.00595.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www. cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/ consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

Stars with initial masses between $\sim0.8$ and 8~$M_\odot$ present copious mass loss during the asymptotic giant branch (AGB) at the end of their lives. The accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the extended AGB stellar atmospheres. The details of this process are not yet well understood, however. Using the extreme-adaptive-optics imager and polarimeter SPHERE/ZIMPOL, we observed light polarised by grains around W\,Hya, SW\,Vir, and R\,Crt, which have mass-loss rates between 10$^{-7}$ and 10$^{-6}~M_\odot~{\rm yr^{-1}}$. We find the distribution of dust to be asymmetric around the three targets. A biconical morphology is seen for R Crt, with a position angle that is very similar to those inferred from interferometric observations of maser emission and of mid-infrared continuum emission. The cause of the biconical outflow cannot be directly inferred from the ZIMPOL data. The dust grains polarise light more efficiently at 0.65~$\mu$m for R\,Crt and SW\,Vir and at 0.82~$\mu$m for W\,Hya. This indicates that at the time of the observations, the grains around SW\,Vir and R\,Crt had sizes $< 0.1~\mu$m, while those around W\,Hya were larger, with sizes $\gtrsim 0.1~\mu$m. The asymmetric distribution of dust around R\,Crt makes the interpretation more uncertain for this star, however. We find that polarised light is produced already from within the visible photosphere of W~Hya, which we reproduce using models with an inner dust shell that is optically thick to scattering. The radial profile of the polarised light observed around W\,Hya reveal a steep dust density profile. We find the wind-acceleration region of W\,Hya to extend to at least $\sim 7~R_\star$, in agreement with theoretical predictions of acceleration up to $\sim 12~R_\star$., Comment: 16 pages, 9 figures, 3 tables

High-mass X-ray binaries (HMXBs) are exceptional astrophysical laboratories that offer a rare glimpse into the physical processes that govern accretion on compact objects, massive-star winds, and stellar evolution. In a subset of the HMXBs, the compact objects accrete matter solely from winds of massive donor stars. These so-called wind-fed HMXBs are divided in persistent HMXBs and supergiant fast X-ray transients (SFXTs) according to their X-ray properties. While it has been suggested that this dichotomy depends on the characteristics of stellar winds, they have been poorly studied. With this investigation, we aim to remedy this situation by systematically analyzing donor stars of wind-fed HMXBs that are observable in the UV, concentrating on those with neutron star (NS) companions. We obtained Swift X-ray data, HST UV spectra, and additional optical spectra for all our targets. Our multi-wavelength approach allows us to provide stellar and wind parameters for six donor stars (four wind-fed systems and two OBe X-ray binaries). The wind properties are in line with the predictions of the line-driven wind theory. Three of the donor stars are in an advanced evolutionary stage, while for some of the stars, the abundance pattern indicates that processed material might have been accreted. When passing by the NS in its tight orbit, the donor star wind has not yet reached its terminal velocity but it is still significantly slower; its speed is comparable with the orbital velocity of the NS companion. There are no systematic differences between the two types of wind-fed HMXBs (persistent versus transients) with respect to the donor stars. For the SFXTs in our sample, the orbital eccentricity is decisive for their transient X-ray nature. Based on the orbital parameters and the further evolution of the donor stars, the investigated HMXBs will presumably form Thorne-\.Zytkow objects in the future., Comment: 44 pages, 13 figures, 15 tables; accepted for publication in A&A

We investigated continuum and molecular line emission of four species (CO, HCN, H$^{13}$CN, and HCO$^+$) at 0.8 mm in the inner arcsec around Eta Car, using ALMA archival observations at a resolution better than 0.2 arcsec. We report the discovery of an asymmetric extended structure northwest of the star, independent from the continuum point source. The structure is only traced by continuum and HCO$^+$, and not detected in the other lines. Kinematics of this structure reveal that the HCO$^+$ gas likely arises from ejecta expelled in the 1890s eruption. The ejecta is propagating outward within the cavity produced by the current wind-wind interaction of Eta Car A and its companion. Chemical analysis of the ejecta reveals an apparent lack of CO and nitrogen-bearing species. We explore possible explanations for this peculiar chemistry, that differentiates this structure from the ejecta of the Great Eruption, rich in HCN and H$^{13}$CN. We also report an absorption component near the continuum point source, only traced by HCN and H$^{13}$CN in their vibrational-ground and vibrationally-excited states. This absorbing gas is attributed to a hot bullet of N-enriched material expelled at a projected velocity of 40 km s$^{-1}$., Accepted for publication in MNRAS. 12 pages, 8 figures

Context. The condensation of inorganic dust grains in the winds of evolved stars is poorly understood. As of today, it is not yet known which molecular clusters form the first dust grains in oxygen-rich (C/O < 1) asymptotic giant branch (AGB) winds. Aluminium oxides and iron-free silicates are often put forward as promising candidates for the first dust seeds. Aims. We aim to constrain the dust formation histories in the winds of oxygen-rich AGB stars. Methods. We obtained Atacama Large Millimeter/sub-millimeter array (ALMA) observations with a spatial resolution of 120 150 mas tracing the dust formation region of the low mass-loss rate AGB star, R Dor, and the high mass-loss rate AGB star, IK Tau. We detected emission line profiles of AlO, AlOH, and AlCl in the ALMA data and used these line profiles to derive a lower limit of atomic aluminium incorporated in molecules. This constrains the aluminium budget that can condense into grains. Results. Radiative transfer models constrain the fractional abundances of AlO, AlOH, and AlCl in IK Tau and R Dor. We show that the gas-phase aluminium chemistry is completely different in both stars with a remarkable difference in the AlO and AlOH abundance stratification. The amount of aluminium locked up in these three molecules is small, 1:110 w.r.t. H2, for both stars, i.e. only 2%of the total aluminium budget. An important result is that AlO and AlOH, which are the direct precursors of alumina (Al2O3) grains, are detected well beyond the onset of the dust condensation, which proves that the aluminium oxide condensation cycle is not fully effcient. The ALMA observations allow us to quantitatively assess the current generation of theoretical dynamical-chemical models for AGB winds. We discuss how the current proposed scenario of aluminium dust condensation for low mass-loss rate AGB stars within a few stellar radii from the star, in particular for R Dor and W Hya, poses a challenge if one wishes to explain both the dust spectral features in the spectral energy distribution (SED) in interferometric data and in the polarized light signal. In particular, the estimated grain temperature of Al2O3 is too high for the grains to retain their amorphous structure. We advocate that large gas-phase (Al2O3)n clusters (n > 34) can be the potential agents of the broad 11 m feature in the SED and in the interferometric data and we propose potential formation mechanisms for these large clusters.

We present a clean, magnitude-limited (IRAC1 or WISE1 $\leq$ 15.0 mag) multiwavelength source catalog for the SMC with 45,466 targets in total, with the purpose of building an anchor for future studies, especially for the massive star populations at low-metallicity. The catalog contains data in 50 different bands including 21 optical and 29 infrared bands, ranging from the ultraviolet to the far-infrared. Additionally, radial velocities and spectral classifications were collected from the literature, as well as infrared and optical variability statistics were retrieved from different projects. The catalog was essentially built upon a $1''$ crossmatching and a $3''$ deblending between the SEIP source list and Gaia DR2 photometric data. Further constraints on the proper motions and parallaxes from Gaia DR2 allowed us to remove the foreground contamination. We estimated that about 99.5\% of the targets in our catalog were most likely genuine members of the SMC. By using the evolutionary tracks and synthetic photometry from MIST and the theoretical $\rm J-K_S$ color cuts, we identified 1,405 RSG, 217 YSG and 1,369 BSG candidates in the SMC in five different CMDs, where attention should also be paid to the incompleteness of our sample. We ranked the candidates based on the intersection of different CMDs. A comparison between the models and observational data shows that the lower limit of initial mass for the RSGs population may be as low as 7 or even 6 $M_{\odot}$ and the RSG is well separated from the AGB population even at faint magnitude, making RSGs a unique population connecting the evolved massive and intermediate stars, since stars with initial mass around 6 to 8 $M_{\odot}$ are thought to go through a second dredge-up to become AGBs. We encourage the interested reader to further exploit the potential of our catalog., 19 pages, 17 figures, A&A accepted, catalogs are included in the ancillary files

SN2010jl was a luminous Type IIn supernova (SN), detected in radio, optical, X-ray and hard X-rays. Here we report on its six year R- and g-band light curves obtained using the Palomar Transient Factory. The light curve was generated using a pipeline based on the proper image subtraction method and we discuss the algorithm performances. As noted before, the R-band light curve, up to about 300 days after maximum light is well described by a power-law decline with a power-law index of about -0.5. Between day 300 and day 2300 after maximum light, it is consistent with a power-law decline, with a power-law index of about -3.4. The longevity of the light curve suggests that the massive circum-stellar material around the progenitor was ejected on time scales of at least tens of years prior to the progenitor explosion., 6 pages, PASP in press

The Small Magellanic Cloud multiple system HD 5980 contains a luminous blue variable (LBV) that underwent a major eruption in 1994, and whose current spectrum is that of a hydrogen-rich Wolf–Rayet (WR) star. Since the eruption, the wind mass-loss rate has been declining while wind speeds have been steadily increasing. Observations obtained in 2014 when Star A (the LBV) eclipses Star B indicate that the fitted mass-loss rate and luminosity have reached the lowest values ever determined for such spectra: Ṁ = 4.5 × 10⁻⁵ M⊙ yr⁻¹, L = 1.7 × 10⁶ L⊙. In addition, the radius of the LBV’s continuum-emitting region is similar to that derived from the eclipse light curves of the late 1970s. Hence, it appears to have attained a similar ‘low’ state to that of the late 1970s. While a good fit to the emission spectrum is obtained using a cmfgen model, there are discrepancies in the UV. In particular, the extent of the observed absorption profiles is ∼1000 km s⁻¹ greater than predicted by the emission-line intensities. Further, HST UV observations obtained in 2016, when Star A is eclipsed by Star B, show unusual P Cygni profiles that are not easily explained. Surprisingly the 2016 emission-line spectrum is similar to that at the opposite eclipse obtained in 2014. The complex UV profiles are likely to arise as a consequence of the dynamics of the wind–wind collision and radiative braking, both of which will cause significant departures from spherical symmetry, and have a strong orbital phase dependence. However, other scenarios, such as intrinsically aspherical winds, cannot be ruled out., Facultad de Ciencias Astronómicas y Geofísicas

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

We follow-up on a previous finding that AGB Mira variables containing the 3DUP indicator technetium (Tc) in their atmosphere form a different sequence of K-[22] colour as a function of pulsation period than Miras without Tc. A near- to mid-infrared colour such as K-[22] is a good probe for the dust mass-loss rate of the stars. Contrary to what might be expected, Tc-poor Miras show redder K-[22] colours (i.e. higher dust mass-loss rates) than Tc-rich Miras at a given period. Here, the previous sample is extended and the analysis is expanded towards other colours and dust spectra. The most important aim is to investigate if the same two sequences can be revealed in the gas mass-loss rate. We analysed new optical spectra and expanded the sample by including more stars from the literature. Near- and mid-IR photometry and ISO dust spectra of our stars were investigated. Literature data of gas mass-loss rates of Miras and semi-regular variables were collected and analysed. Our results show that Tc-poor Miras are redder than Tc-rich Miras in a broad range of the mid-IR, suggesting that the previous finding based on the K-[22] colour is not due to a specific dust feature in the 22 micron band. We establish a linear relation between K-[22] and the gas mass-loss rate. We also find that the 13 micron feature disappears above K-[22]~2.17 mag, corresponding to $\dot{M}_{\rm g}\sim2.6\times10^{-7}M_{\sun}yr^{-1}$. No similar sequences of Tc-poor and Tc-rich Miras in the gas mass-loss rate vs. period diagram are found, most probably owing to limitations in the available data. Different hypotheses to explain the observation of two sequences in the P vs. K-[22] diagram are discussed and tested, but so far none of them convincingly explains the observations. Nevertheless, we might have found an hitherto unknown but potentially important process influencing mass loss on the TP-AGB., 16 pages, 15 figures, 2 online tables, accepted for publication in A&A

Despite their cosmological utility, the progenitors of Type Ia supernovae (SNe Ia) are still unknown, with many efforts focused on whether accretion from a nondegenerate companion can grow a carbon-oxygen white dwarf to near the Chandrasekhar mass. The association of SNe Ia resembling SN 1991T ("91T-like") with circumstellar interaction may be evidence for this "single-degenerate" channel. However, the observed circumstellar medium (CSM) in these interacting systems is unlike a stellar wind -- of particular interest, it is sometimes detached from the stellar surface, residing at $\sim 10^{16}~{\rm cm}$. A Hubble Space Telescope (HST) program to discover detached CSM around 91T-like SNe Ia successfully discovered interaction nearly two years after explosion in SN 2015cp (Graham et al., 2018). In this work, we present radio and X-ray follow-up observations of SN 2015cp and analyze them in the framework of Harris, Nugent, & Kasen (2016) to limit the properties of a constant-density CSM shell in this system. Assuming the HST detection was shortly after the shock crossed the CSM, we constrain the total CSM mass in this system to be $< 0.5~{\rm M_\odot}$. This limit is comparable to the CSM mass of supernova PTF11kx, but does not rule out lower masses predicted for recurrent novae. From lessons learned modeling PTF11kx and SN 2015cp, we suggest a strategy for future observations of these events to increase the sample of known interacting SNe Ia., 13 pages, 5 figures, published in The Astrophysical Journal

Comprehensive spectral analyses of the Galactic Wolf-Rayet stars of the nitrogen sequence (i.e.\ the WN subclass) have been performed in a previous paper. However, the distances of these objects were poorly known. Distances have a direct impact on the "absolute" parameters, such as luminosities and mass-loss rates. The recent Gaia Data Release (DR2) of trigonometric parallaxes includes nearly all WN stars of our Galactic sample. In the present paper, we apply the new distances to the previously analyzed Galactic WN stars and rescale the results accordingly. On this basis, we present a revised catalog of 55 Galactic WN stars with their stellar and wind parameters. The correlations between mass-loss rate and luminosity show a large scatter, for the hydrogen-free WN stars as well as for those with detectable hydrogen. The slopes of the $\log L - \log \dot{M}$ correlations are shallower than found previously. The empirical Hertzsprung-Russell diagram (HRD) still shows the previously established dichotomy between the hydrogen-free early WN subtypes that are located on the hot side of the zero-age main sequence (ZAMS), and the late WN subtypes, which show hydrogen and reside mostly at cooler temperatures than the ZAMS (with few exceptions). However, with the new distances, the distribution of stellar luminosities became more continuous than obtained previously. The hydrogen-showing stars of late WN subtype are still found to be typically more luminous than the hydrogen-free early subtypes, but there is a range of luminosities where both subclasses overlap. The empirical HRD of the Galactic single WN stars is compared with recent evolutionary tracks. Neither these single-star evolutionary models nor binary scenarios can provide a fully satisfactory explanation for the parameters of these objects and their location in the HRD., Comment: 11 pages; accepted for publication in Astronomy & Astrophysics

13 pags., 9 figs., 2 tabs., Context. Asymptotic giant branch (AGB) stars go through a process of strong mass loss that involves pulsations of the atmosphere, which extends to a region in which the conditions are adequate for dust grains to form. Radiation pressure acts on these grains which, coupled to the gas, drive a massive outflow. The details of this process are not clear, including which molecules are involved in the condensation of dust grains. Aims. We seek to study the role of the SiO molecule in the process of dust formation and mass loss in M-type AGB stars. Methods. Using the IRAM NOEMA interferometer we observed the SiO and SiO J = 3-2, v = 0 emission from the inner circumstellar envelope of the evolved stars IK Tau and IRC+10011. We computed azimuthally averaged emission profiles to compare the observations to models using a molecular excitation and ray-tracing code for SiO thermal emission. Results. We observe circular symmetry in the emission distribution. We also find that the source diameter varies only marginally with radial velocity, which is not the expected behaviour for envelopes expanding at an almost constant velocity. The adopted density, velocity, and abundance laws, together with the mass-loss rate, which best fit the observations, give us information concerning the chemical behaviour of the SiO molecule and its role in the dust formation process. Conclusions. The results indicate that there is a strong coupling between the depletion of gas-phase SiO and gas acceleration in the inner envelope. This could be explained by the condensation of SiO into dust grains., The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013) / ERC Grant Agreement n. 610256 NANOCOSMOS and by the Spanish MINECO, Grant FIS2012-32096. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

We link the onset of pulsation-enhanced, dust-driven winds from asymptotic giant branch (AGB) stars in the Magellanic Clouds to the star's transition between period--luminosity sequences (from B to C'). This transition occurs at ~60 days for solar-mass stars, which represent the bulk of the AGB population: this is the same period at which copious dust production starts in solar-neighbourhood AGB stars. It is contemporaneous with the onset of long-secondary period (LSP) variability on sequence D. The combined amplitude of the first-overtone (B+C') and fundamental (C) modes and (perhaps) long-secondary period (D; LSP) variability appears to drive a sudden increase in mass-loss rate to a stable plateau, previously identified to be a few x 10^-7 solar masses per year. We cite this as evidence that pulsations are necessary to initiate mass loss from AGB stars and that these pulsations are significant in controlling stars' mass-loss rates. We also show evidence that LSPs may evolve from long to short periods as the star evolves, counter to the other period-luminosity sequences., Comment: 6 pages, accepted MNRAS

We report the first detection of circumstellar CO in a globular cluster. Observations with ALMA have detected the CO J=3-2 and SiO v=1 J=8-7 transitions at 345 and 344 GHz, respectively, around V3 in 47 Tucanae (NGC 104; [Fe/H] = -0.72 dex), a star on the asymptotic giant branch. The CO line is detected at 7 sigma at a rest velocity v_LSR = -40.6 km/s and expansion velocity of 3.2 +/- ~0.4 km/s. The brighter, asymmetric SiO line may indicate a circumstellar maser. The stellar wind is slow compared to similar Galactic stars, but the dust opacity remains similar to Galactic comparisons. We suggest that the mass-loss rate is set by the levitation of material into the circumstellar environment by pulsations, but that the terminal wind-expansion velocity is determined by radiation pressure on the dust: a pulsation-enhanced dust-driven wind. We suggest the metal-poor nature of the star decreases the grain size, slowing the wind and increasing its density and opacity. Metallic alloys at high altitudes above the photosphere could also provide an opacity increase. The CO line is weaker than expected from Galactic AGB stars, but its strength confirms a model that includes CO dissociation by the strong interstellar radiation field present inside globular clusters., Comment: 5 pages, accepted MNRAS Letters

It is well known that ordinary supernovae (SNe) are powered by 56Ni cascade decay. Broad-lined type Ic SNe (SNe Ic-BL) are a subclass of SNe that are not all exclusively powered by 56Ni decay. It was suggested that some SNe Ic-BL are powered by magnetar spin-down. iPTF16asu is a peculiar broad-lined type Ic supernova discovered by the intermediate Palomar Transient Factory. With a rest-frame rise time of only 4 d, iPTF16asu challenges the existing popular models, for example, the radioactive heating (56Ni-only) and the magnetar +56Ni models. Here we show that this rapid rise could be attributed to interaction between the SN ejecta and a pre-existing circumstellar medium ejected by the progenitor during its final stages of evolution, while the late-time light curve can be better explained by energy input from a rapidly spinning magnetar. This model is a natural extension to the previous magnetar model. The mass-loss rate of the progenitor and ejecta mass are consistent with a progenitor that experienced a common envelope evolution in a binary. An alternative model for the early rapid rise of the light curve is the cooling of a shock propagating into an extended envelope of the progenitor. It is difficult at this stage to tell which model (interaction+magnetar + 56Ni or cooling+magnetar + 56Ni) is better for iPTF16asu. However, it is worth noting that the inferred envelope mass in the cooling+magnetar + 56Ni is very high.© 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society., This work is supported by the National Program on Key Research and Development Project of China (Grant No. 2016YFA0400801 and 2017YFA0402600), and the National Natural Science Foundation of China (Grant Nos. 11573014,11533033, 11673006, and 11833003). XFW is supported by the National Natural Science Foundation of China (NSFC grants 11325313 and 11633002), and the National Program on Key Research and Development Project (grant no. 2016YFA0400803). SQW and LDL are also supported by China Scholarship Program to conduct research at U.C. Berkeley and UNLV, respectively.

We present the highest resolution single-dish submillimetre observations of the detached shell source U Antliae to date. The observations were obtained at $450~\micron$ and $850~\micron$ with SCUBA-2 instrument on the James Clerk Maxwell Telescope as part of the Nearby Evolved Stars Survey. The emission at $850~\micron$ peaks at $40\arcsec$ with hints of a second peak seen at $\sim 20\arcsec$. The emission can be traced out to a radius of $56\arcsec$ at a $3\sigma$ level. The outer peak observed at $850~\micron$ aligns well with the peak observed at Herschel/PACS wavelengths. With the help of spectral energy distribution fitting and radiative transfer calculations of multiple-shell models for the circumstellar envelope, we explore the various shell structures and the variation of grain sizes along the in the circumstellar envelope. We determine a total shell dust mass of $(2.0 \pm 0.3) \times 10^{-5}$ M$_{\odot}$ and established that the thermal pulse which gave rise to the detached shell occurred 3500 $\pm$ 500 years ago., Comment: Accepted for publication in MNRAS

We present detailed observations of ZTF18abukavn (SN2018gep), discovered in high-cadence data from the Zwicky Transient Facility as a rapidly rising ($1.4\pm0.1$ mag/hr) and luminous ($M_{g,\mathrm{peak}}=-20$ mag) transient. It is spectroscopically classified as a broad-lined stripped-envelope supernova (Ic-BL SN). The high peak luminosity ($L_{\mathrm{bol}} \gtrsim 3 \times 10^{44}$ erg $\mathrm{sec}^{-1}$), the short rise time ($t_{\mathrm{rise}}= 3$ days in $g$-band), and the blue colors at peak ($g-r\sim-0.4$) all resemble the high-redshift Ic-BL iPTF16asu, as well as several other unclassified fast transients. The early discovery of SN2018gep (within an hour of shock breakout) enabled an intensive spectroscopic campaign, including the highest-temperature ($T_{\mathrm{eff}}\gtrsim40,000$ K) spectra of a stripped-envelope SN. A retrospective search revealed luminous ($M_g \sim M_r \approx -14\,$mag) emission in the days to weeks before explosion, the first definitive detection of precursor emission for a Ic-BL. We find a limit on the isotropic gamma-ray energy release $E_\mathrm{��,iso}10$ days) light curve requires an additional energy source, which could be the radioactive decay of Ni-56., Accepted to ApJ on 1 Oct 2019. In this version, we made minor corrections and removed extraneous references

5 pags., 3 figs., 1 tab., We present new interferometer molecular observations of R Leo taken at 1.2 mm with the Atacama Large Millimeter Array with an angular resolution up to 0'. 026. These observations permitted us to resolve the innermost envelope of this star, which revealed a complex structure that involves extended continuum emission and molecular emission showing a non-radial gas velocity distribution. This molecular emission displays prominent red-shifted absorptions located immediately in front of the star, which are typical footprints of material infall. This emission also shows lateral gas motions compatible with a torus-like structure., We thank the European Research Council (ERC Grant 610256: NANOCOSMOS) and the Spanish MINECO/MICINN for funding support through grant AYA2016-75066-C-1-P and the ASTROMOL Consolider project CSD2009-00038. MA acknowledges funding support from the Ramón y Cajal program of Spanish MINECO (RyC-2014-16277).