$\eta$~Car is one of the most massive stars in the Galaxy. It underwent a massive eruption in the 19th century, which produced the impressive bipolar Homunculus nebula now surrounding it. The central star is an eccentric binary with a period of 5.54\,years. Although the companion has not been detected directly, it causes time-variable ionization and colliding-wind X-ray emission. By characterizing the complex structure and kinematics of the ejecta close to the star, we aim to constrain past and present mass loss of $\eta$~Car. $\eta$~Car is observed with the extreme adaptive optics instrument SPHERE at the Very Large Telescope, using its polarimetric mode in the optical with the ZIMPOL camera. A spatial resolution of 20\,mas was achieved, i.e. very close to the presumed 13 mas apastron separation of the companion star. We detect new structures within the inner arcsecond to the star (2\,300\,au at a 2.3\,kpc distance). We can relate these structures to the eruption near 1890 by tracking their proper motions derived from our new images and historical images over a 30\,years time span. Besides, we find a fan-shaped structure in the inner 200~au to the star in the H$\alpha$ line, that could potentially be associated with the wind collision zone of the two stars.

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.

Cool, evolved stars are the main source of chemical enrichment of the interstellar medium (ISM), and understanding their mass loss and structure offers a unique opportunity to study the cycle of matter in the Universe. Pulsation, convection, and other dynamic processes in cool evolved stars create an atmosphere where molecules and dust can form, including those necessary to the formation of life (e.g.~Carbon-bearing molecules). Understanding the structure and composition of these stars is thus vital to several aspects of stellar astrophysics, ranging from ISM studies to modeling young galaxies and exoplanet research. Recent modeling efforts and increasingly precise observations now reveal that our understanding of cool stars photospheric, chromospheric, and atmospheric structures is limited by inadequate knowledge of the dynamic and chemical processes at work. Here we outline promising scientific opportunities for the next decade. We identify and discuss the following main opportunities: (1) identify and model the physical processes that must be included in current 1D and 3D atmosphere models of cool, evolved stars; (2) refine our understanding of photospheric, chromospheric, and outer atmospheric regions of cool evolved stars, their properties and parameters, through high-resolution spectroscopic observations, and interferometric observations at high angular resolution; (3) include the neglected role of chromospheric activity in the mass loss process of red giant branch and red super giant stars, and understand the role played by their magnetic fields; (4) identify the important shaping mechanisms for planetary nebulae and their relation with the parent asymptotic giant branch stars.

We assess the ability to image the photosphere of red supergiants and, in particular Betelgeuse, through the modelling of the observed linear polarization in atomic spectral lines. We also aim to analyse the resulting images over time, to measure the size and dynamics of the convective structures in these stars. Rayleigh scattering polarizes the continuum and spectral lines depolarize it.This depolarization is seen as a linear polarization signal parallel to the radial direction on the stellar disk. Integrated over the disk, it would result in a null signal, except if brightness asymmetries/inhomogeneities are present. This is the basic concept behind our imaging technique. The several tests and comparisons performed prove that our technique reliably retrieves the salient brightness structures in the photosphere of Betelgeuse, and should be relevant to other red supergiants. We demonstrate that these structures are convective cells, with a characteristic size of more than 60% of the stellar radius. We also derive the characteristic upflow and downflow speeds, 22 and 10 km/s respectively. We find weak magnetic fields concentrated in the downflow lanes in between granules. We follow those convective structures in time. Changes happen on timescales of one week, but individual structures can be tracked over four years of observations.

Betelgeuse is an M-type supergiant that presents a circularly polarized (Stokes V) signal in its line profiles, interpreted in terms of a surface magnetic field. The weak circular polarization signal has been monitored over 7.5 years in order to follow its evolution on different timescales, and eventually to determine its physical origin. Linear polarization measurements have also been obtained regularly in the last few years. We used both the ESPaDOnS and Narval spectropolarimeters to obtain high signal-to-noise ratio (S/N) spectra, which were processed by means of the least-squares deconvolution (LSD) method. In order to ensure the reality of the very weak circular polarization, special care has been taken to limit instrumental effects. In addition, several tests were performed on the Stokes V signal to establish its stellar and Zeeman origin. We confirm the magnetic nature of the circular polarization, pointing to a surface magnetic field of the order of 1G. The Stokes V profiles present variations over different timescales, the most prominent one being close to the long secondary period (LSP; around 2000d for Betelgeuse) often invoked in red evolved stars. This long period is also dominant for all the other Stokes parameters. The circular polarization is tentatively modeled by means of magnetic field concentrations mimicking spots, showing in particular that the velocity associated with each "spot" also follows the long timescale, and that this signal is nearly always slightly redshifted. From the coupled variations of both linear and circular polarization signatures in amplitude, velocity and timescale, we favour giant convection cells as the main engine at the origin of polarization signatures and variations in all the Stokes parameters. This strengthens support for the hypothesis that large convective cells are at the origin of the LSP.

Context. Red supergiant stars are one of the latest stages in the evolution of massive stars. Their photospheric convection may play an important role in the launching mechanism of their mass loss; however, its characteristics and dynamics are still poorly constrained. Aims. By observing red supergiant stars with near infrared interferometry at different epochs, we expect to reveal the evolution of bright convective features on their stellar surface. Methods. We observed the M2Iab-Ib red supergiant star CE Tau with the VLTI/PIONIER instrument in the H band at two different epochs separated by one month. Results. We derive the angular diameter of the star and basic stellar parameters, and reconstruct two reliable images of its H-band photosphere. The contrast of the convective pattern of the reconstructed images is $5 \pm 1 \%$ and $6 \pm 1 \%$ for our two epochs of observation. Conclusions. The stellar photosphere shows few changes between the two epochs. The contrast of the convective pattern is below the average contrast variations obtained on 30 randomly chosen snapshots of the best matching 3D radiative hydrodynamics simulation: $23 \pm 1~\%$ for the original simulation images and $16 \pm 1~\%$ for the maps degraded to the reconstruction resolution. We offer two hypotheses to explain this observation. CE Tau may be experiencing a quiet convective activity episode or it could be a consequence of its warmer effective temperature (hence its smaller radius) compared to the simulation.

The objective of the SPHERE Data Center is to optimize the scientific return of SPHERE at the VLT, by providing optimized reduction procedures, services to users and publicly available reduced data. This paper describes our motivation, the implementation of the service (partners, infrastructure and developments), services, description of the on-line data, and future developments. The SPHERE Data Center is operational and has already provided reduced data with a good reactivity to many observers. The first public reduced data have been made available in 2017. The SPHERE Data Center is gathering a strong expertise on SPHERE data and is in a very good position to propose new reduced data in the future, as well as improved reduction procedures.

We observed Betelgeuse using ALMA's extended configuration in band 7 (f~340 GHz, {\lambda}~0.88 mm), resulting in a very high angular resolution of 18 mas. Using a solid body rotation model of the 28SiO(v=2,J=8-7) line emission, we show that the supergiant is rotating with a projected equatorial velocity of v_eq sin i = 5.47 +/- 0.25 km/s at the equivalent continuum angular radius R_star = 29.50 +/- 0.14 mas. This corresponds to an angular rotation velocity of {\omega} sin i = (5.6 +/- 1.3) x 10^(-9) rad/s. The position angle of its north pole is PA = 48.0 +/- 3.5{\deg}. The rotation period of Betelgeuse is estimated to P/sin i = 36 +/- 8 years. The combination of our velocity measurement with previous observations in the ultraviolet shows that the chromosphere is co-rotating with the star up to a radius of ~10 au (45 mas or 1.5x the ALMA continuum radius). The coincidence of the position angle of the polar axis of Betelgeuse with that of the major ALMA continuum hot spot, a molecular plume, and a partial dust shell (from previous observations) suggests that focused mass loss is currently taking place in the polar region of the star. We propose that this hot spot corresponds to the location of a particularly strong "rogue" convection cell, which emits a focused molecular plume that subsequently condenses into dust at a few stellar radii. Rogue convection cells therefore appear to be an important factor shaping the anisotropic mass loss of red supergiants.

Context. Convection is a candidate to explain the trigger of red supergiant star mass loss. Owing to the small size of the convective cells on the photosphere, few of their characteristics are known. Aims. Using near infrared interferometry, we intend to resolve the photosphere of red supergiant stars and to bring new constraints on their modeling. Methods. We observed the nearby red supergiant Antares using the 4-telescope instrument VLTI/PIONIER. We collected data on the three available configurations of the 1.8m telescopes in the H band. Results. We obtained unprecedented angular resolution on the disk of a star (6% of the star angular diameter) that limits the mean size of convective cells, brings new constraints on numerical simulations. Using an analytical model with a distribution of bright spots we determine their effect on the visibility signal. Conclusion. We determine that the interferometric signal on Antares is compatible with convective cells with various sizes from 45% to 5% of the angular diameter. We also conclude that convective cells can strongly affect the angular diameter and limb-darkening measurements. In particular, the apparent angular diameter becomes dependent on the sampled position angles.

The mechanisms responsible for heating the extended atmospheres of early-M spectral-type supergiants are poorly understood. So too is the subsequent role these mechanisms play in driving the large mass-loss rates of these stars. Here we present ALMA long (i.e., $\sim$16 km) baseline 338 GHz (0.89 mm) continuum observations of the free-free emission in the extended atmosphere of the M2 spectral-type supergiant Betelgeuse. The spatial resolution of 14 mas exquisitely resolves the atmosphere, revealing it to have a mean temperature of 2760 K at $\sim$1.3 R$_{\star}$, which is below both the photospheric effective temperature ($T_{\textrm{eff}} = 3690$ K) and the temperatures at $\sim$$2 $R$_{\star}$. This is unambiguous proof for the existence of an inversion of the mean temperature in the atmosphere of a red supergiant. The emission is clearly not spherically symmetric with two notable deviations from a uniform disk detected in both the images and visibilities. The most prominent asymmetry is located in the north-east quadrant of the disk and is spatially resolved showing it to be highly elongated with an axis-ratio of 2.4 and occupying $\sim$$5\%$ of the disk projected area. Its temperature is approximately 1000 K above the measured mean temperature at 1.3 R$_{\star}$. The other main asymmetry is located on the disk limb almost due east of the disk center and occupies $\sim$$3\%$ of the disk projected area. Both emission asymmetries are clear evidence for localized heating taking place in the atmosphere of Betelgeuse. We suggest that the detected localized heating is related to magnetic activity generated by large-scale photospheric convection.

Context. Red giant branch (RGB) stars are very bright objects in galaxies and are often used as standard candles. Interferometry is the ideal tool to characterize the dynamics and morphology of their atmospheres. Aims. We aim at precisely characterising the surface dynamics of a sample of RGB stars. Methods. We obtained interferometric observations for three RGB stars with the MIRC instrument mounted at the CHARA interfer- ometer. We looked for asymmetries on the stellar surfaces using limb-darkening models. Results. We measured the apparent diameters of HD197989 (Epsilon Cyg) = 4.61+-0.02 mas, HD189276 (HR7633) = 2.95+-0.01 mas, and HD161096 (Beta Oph) = 4.43+-0.01 mas. We detected departures from the centrosymmetric case for all three stars with the tendency of a greater effect for lower logg of the sample. We explored the causes of this signal and conclude that a possible explanation to the interferometric signal is the convection-related and/or the magnetic-related surface activity. However, it is necessary to monitor these stars with new observations, possibly coupled with spectroscopy, in order to firmly establish the cause.

Six billion years from now, while evolving on the asymptotic giant branch (AGB), the Sun will metamorphose from a red giant into a beautiful planetary nebula. This spectacular evolution will impact the Solar System planets, but observational confirmations of the predictions of evolution models are still elusive as no planet orbiting an AGB star has yet been discovered. The nearby AGB red giant L2 Puppis (d = 64 pc) is surrounded by an almost edge-on circumstellar dust disk. We report new observations with ALMA at very high angular resolution (18 x 15 mas) in band 7 (f ~ 350 GHz) that allow us to resolve the velocity profile of the molecular disk. We establish that the gas velocity profile is Keplerian within the central cavity of the dust disk, allowing us to derive the mass of the central star L2 Pup A, mA = 0.659 +/- 0.011 +/- 0.041 Msun (+/- 6.6%). From evolutionary models, we determine that L2 Pup A had a near-solar main sequence mass, and is therefore a close analog of the future Sun in 5 to 6 Gyr. The continuum map reveals the presence of a secondary source (B) at a radius of 2 AU contributing fB/ fA = 1.3 +/- 0.1% of the flux of the AGB star. L2 Pup B is also detected in CO emission lines at a radial velocity of vB = 12.2 +/- 1.0 km/s. The close coincidence of the center of rotation of the gaseous disk with the position of the continuum emission from the AGB star allows us to constrain the mass of the companion to mB = 12 +/- 16 MJup. L2 Pup B is most likely a planet or low mass brown dwarf with an orbital period around 5 years. Its continuum brightness and molecular emission suggest that it may be surrounded by an extended molecular atmosphere or an accretion disk. L2 Pup therefore emerges as a promising vantage point on the distant future of our Solar System.

Betelgeuse is an M supergiant that harbors spots and giant granules at its surface and presents linear polarization of its continuum. We have previously discovered linear polarization signatures associated with individual lines in the spectra of cool and evolved stars. Here, we investigate whether a similar linearly polarized spectrum exists for Betelgeuse. We used the spectropolarimeter Narval, combining multiple polarimetric sequences to obtain high signal-to-noise ratio spectra of individual lines, as well as the least-squares deconvolution (LSD) approach. We have discovered the existence of a linearly polarized spectrum for Betelgeuse, detecting a rather strong signal (at a few times 10$^{-4}$ of the continuum intensity level), both in individual lines and in the LSD profiles. Studying its properties and the signal observed for the resonant \ion{Na}{i}\,D lines, we conclude that we are mainly observing depolarization of the continuum by the absorption lines. The linear polarization of the Betelgeuse continuum is due to the anisotropy of the radiation field induced by brightness spots at the surface and Rayleigh scattering in the atmosphere. We have developed a geometrical model to interpret the observed polarization, from which we infer the presence of two brightness spots and their positions on the surface of Betelgeuse. We show that applying the model to each velocity bin along the Stokes Q and U profiles allows the derivation of a map of the bright spots. We use the Narval linear polarization observations of Betelgeuse obtained over a period of 1.4 years to study the evolution of the spots. Our study of the linearly polarized spectrum of Betelgeuse provides a novel method for studying the evolution of brightness spots at its surface and complements quasi-simultaneous observations obtained with PIONIER at the VLTI.

Context. The mass-loss mechanism of cool massive evolved stars is poorly understood. The proximity of Betelgeuse makes it an appealing target to study its atmosphere, map the shape of its envelope, and follow the structure of its wind from the photosphere out to the interstellar medium. Aims. A link is suspected between the powerful convective motions in Betelgeuse and its mass loss. We aim to constrain the spatial structure and temporal evolution of the convective pattern on the photosphere and to search for evidence of this link. Methods. We report new interferometric observations in the infrared H band using the VLTI/PIONIER instrument. We monitored the photosphere of Betelgeuse between 2012 January and 2014 November to look for evolutions that may trigger the outflow. Results. Our interferometric observations at low spatial frequencies are compatible with the presence of a hot spot on the photosphere that has a characteristic width of one stellar radius. It appears to be superposed on the smaller scale convective pattern. In the higher spatial frequency domain, we observe a significant difference between the observations and the predictions of 3D hydrodynamical simulations. Conclusions. We bring new evidence for the presence of a convective pattern in the photosphere of red supergiants. The inferred hot spot is probably the top of a giant convection cell although an asymmetric extension of the star cannot be excluded by these interferometric observations alone. The properties of the observed surface features show a stronger contrast and inhomogeneity as predicted by 3D radiative hydrodynamical simulations. We propose that the large observed feature is modifying the signature of the convective pattern at the surface of the star in a way that simulations cannot reproduce.

Our team has obtained observations of the photosphere of the two closest red supergiant stars Betelgeuse ($\alpha$ Ori) and Antares ($\alpha$ Sco) using near infrared interferometry. We have been monitoring the photosphere of Betelgeuse with the VLTI/PIONIER instrument for three years. On Antares, we obtained an unprecedented sampling of the visibility function. These data allow us to probe the convective photosphere of massive evolved stars.

Adaptive optics observations in the infrared (VLT/NACO, Kervella et al. 2014) and visible (VLT/SPHERE, Kervella et al. 2015) domains revealed that the nearby AGB star L2 Pup (d=64 pc) is surrounded by a dust disk seen almost edge-on. Thermal emission from a large dust "loop" is detected at 4 microns up to more than 10 AU from the star. We also detect a secondary source at a separation of 32 mas, whose nature is uncertain. L2 Pup is currently a relatively "young" AGB star, so we may witness the formation of a planetary nebula. The mechanism that breaks the spherical symmetry of mass loss is currently uncertain, but we propose that the dust disk and companion are key elements in the shaping of the bipolar structure. L2 Pup emerges as an important system to test this hypothesis.

The bright southern star L2 Pup is a particularly prominent asymptotic giant branch (AGB) star, as its distance of 64 pc makes it the nearest of its type. We report new adaptive optics observations of L2 Pup at visible wavelengths with the SPHERE/ZIMPOL instrument of the VLT that confirm the presence of the circumstellar dust disk at high inclination discovered recently by Kervella et al. (2014b). The signature of the three-dimensional structure of the disk is clearly observed in the map of the degree of linear polarization pL. We identify the inner rim of the disk through its polarimetric signature at a radius of 6 AU from the AGB star. The ZIMPOL intensity images in the V and R bands also reveal a close-in secondary source at a projected separation of 2 AU from the primary. The identification of the spectral type of this companion is uncertain due to the strong reddening from the disk, but its photometry suggests that it is a late K giant, of comparable mass to the AGB star. We present refined physical parameters for the dust disk derived using the RADMC-3D radiative transfer code. We also interpret the pL map using a simple polarization model to infer the three-dimensional structure of the envelope. Interactions between the inner binary system and the disk apparently form spiral structures that propagate along the orthogonal axis to the disk to form streamers. Two dust plumes propagating orthogonally to the disk are also detected. They originate in the inner stellar system, and are possibly related to the interaction of the wind of the two stars with the material in the disk. Based on the morphology of the envelope of L2 Pup, we propose that this star is at an early stage of the formation of a bipolar planetary nebula.

The physical mechanism through which the outgoing material of massive red supergiants is accelerated above the escape velocity is unclear. Thanks to the transparency of its circumstellar envelope, the nearby red supergiant Betelgeuse gives the opportunity to probe the innermost layers of the envelope of a typical red supergiant down to the photosphere, i.e. where the acceleration of the wind is expected to occur. We took advantage of the SPHERE/ZIMPOL adaptive optics imaging polarimeter to resolve the visible photosphere and close envelope of Betelgeuse. We detect an asymmetric gaseous envelope inside a radius of 2 to 3 times the near-infrared photospheric radius of the star (R*), and a significant Halpha emission mostly contained within 3 R*. From the polarimetric signal, we also identify the signature of dust scattering in an asymmetric and incomplete dust shell located at a similar radius. The presence of dust so close to the star may have a significant impact on the wind acceleration through radiative pressure on the grains. The 3 R* radius emerges as a major interface between the hot gaseous and dusty envelopes. The detected asymmetries strengthen previous indications that the mass loss of Betelgeuse is likely tied to the vigorous convective motions in its atmosphere.

Context. Betelgeuse is the closest red supergiant (RSG); therefore, it is well suited for studying the complex processes in its atmosphere that lead to the chemical enrichment of the interstellar medium. Aims. We intend to investigate the shape and composition of the close molecular layer (also known as the MOLsphere) that surrounds the star. This analysis is part of a wider program that aims at understanding the dynamics of the circumstellar envelope of Betelgeuse. Methods. On January and February 2011, Betelgeuse was observed using the Astronomical Multi-BEam combineR (AMBER) instrument of the Very Large Telescope Interferometer (VLTI) in the H and K bands. Using the medium spectral resolution of the instrument (R $\sim$ 1500), we were able to investigate the carbon monoxide band heads and the water-vapor bands. We used two different approaches to analyse our data: a model fit in both the continuum and absorption lines and then a fit with a Radiative HydroDynamics (RHD) simulation. Results. Using the continuum data, we derive a uniform disk diameter of $41.01 \pm 0.41$~mas, a power law type limb-darkened disk diameter of $42.28 \pm 0.43$~mas and a limb-darkening exponent of $0.155 \pm 0.009$. Within the absorption lines, using a single layer model, we obtain parameters of the MOLsphere. Using a RHD simulation, we unveil the convection pattern in the visibilities. Conclusions. We derived a new value of the angular diameter of Betelgeuse in the K band continuum. Our observations in the absorption lines are well reproduced by a molecular layer at 1.2 stellar radii containing both CO and H$_2$O. The visibilities at higher spatial frequencies are matching a convection pattern in a RHD simulation.

As the nearest known AGB star (d=64pc) and one of the brightest (mK-2), L2 Pup is a particularly interesting benchmark object to monitor the final stages of stellar evolution. We report new lucky imaging observations of this star with the VLT/NACO adaptive optics system in twelve narrow band filters covering the 1.0-4.0 microns wavelength range. These diffraction limited images reveal an extended circumstellar dust lane in front of the star, that exhibits a high opacity in the J band and becomes translucent in the H and K bands. In the L band, extended thermal emission from the dust is detected. We reproduce these observations using Monte-Carlo radiative transfer modeling of a dust disk with the RADMC-3D code. We also present new interferometric observations with the VLTI/VINCI and MIDI instruments. We measure in the K band an upper limit to the limb-darkened angular diameter of theta_LD = 17.9 +/- 1.6 mas, converting to a maximum linear radius of R = 123 +/- 14 Rsun. Considering the geometry of the extended K band emission in the NACO images, this upper limit is probably close to the actual angular diameter of the star. The position of L2 Pup in the Herzsprung-Russell diagram indicates that this star has a mass around 2 Msun and is probably experiencing an early stage of the asymptotic giant branch. We do not detect any stellar companion of L2 Pup in our adaptive optics and interferometric observations, and we attribute its apparent astrometric wobble in the Hipparcos data to variable lighting effects on its circumstellar material. We however do not exclude the presence of a binary companion, as the large loop structure extending to more than 10 AU to the North-East of the disk in our L band images may be the result of interaction between the stellar wind of L2 Pup and a hidden secondary object. The geometric configuration that we propose, with a large dust disk seen almost edge-on, appears particularly favorable to test and develop our understanding of the formation of bipolar nebulae.

We present the results of the analysis of our recent interferometric observations of Betelgeuse, using the AMBER instrument of the VLTI. Using the medium spectral resolution mode ($R \sim 1500$) we detected the presence of the water vapour and carbon monoxide (CO) molecules in the H and K bands. We also derived the photospheric angular diameter in the continuum. By analysing the depth of the molecular lines and the interferometric visibilities, we derived the column densities of the molecules, as well as the temperature and the size of the corresponding regions in the atmosphere of Betelgeuse (the MOLsphere) using a single shell model around the photosphere. Our results confirm the findings by Perrin \et al\ (\cite{Perrin2004}) and Ohnaka \et al\ (\cite{Ohnaka2011}) that the H$_2$O and CO molecules are distributed around Betelgeuse in a MOLsphere extending to approximately 1.3 times the star's photospheric radius.