It has been shown recently that the infrared emission of Cepheids, which is constant over the pulsation cycle, might be due to a pulsating shell of ionized gas of about 15\% of the stellar radius, which could be attributed to the chromospheric activity of Cepheids. The aim of this paper is to investigate the dynamical structure of the chromosphere of Cepheids along the pulsation cycle and quantify its size. We present H$\alpha$ and Calcium Near InfraRed triplet (Ca IR) profile variations using high-resolution spectroscopy with the UVES spectrograph of a sample of 24 Cepheids with a good period coverage from $\approx$ 3 to 60 days. After a qualitative analysis of the spectral lines profiles, we quantify the Van Hoof effect (velocity gradient between the H$\alpha$ and Ca IR) as a function of the period of the Cepheids. Then, we use the Schwarzschild mechanism (a line doubling due to a shock wave) to quantify the size of the chromosphere. We find a significant Van Hoof effect for Cepheids with period larger than $P=10$ days, in particular H$\alpha$ lines are delayed with a velocity gradient up to $\Delta v \approx$30 km/s compared to Ca IR. We find that the size of the chromosphere of long-period Cepheids is of at least $\approx$ 50\% of the stellar radius, which is consistent at first order with the size of the shell made of ionized gas previously found from the analysis of infrared excess. Last, for most of the long-period Cepheids in the sample, we report a motionless absorption feature in the H$\alpha$ line that we attribute to a circumstellar envelope that surrounds the chromosphere. Analyzing the Ca~IR lines of Cepheids is of importance to potentially unbias the period-luminosity relation from their infrared excess, particularly in the context of forthcoming observations from the Radial Velocity Spectrometer (RVS) on board \textit{Gaia}, that could be sensitive to their chromosphere.

Context. Asymptotic giant branch stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius and compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 micron, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths and recovered the radius using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations.

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

Context. Red supergiant stars possess surface features and extended molecular atmospheres. Photospheric convection may be a crucial factor of the levitation of the outer atmospheric layers. However, the mechanism responsible is still poorly understood. Aims. We image the stellar surface of V602 Carinae (V602 Car) to constrain the morphology and contrast of the surface features and of the extended atmospheric layers. Methods. We observed V602 Car with the Very Large Telescope Interferometer (VLTI) PIONIER instrument (1.53-1.78 $\mathrm{\mu}$m) between May and July 2016, and April and July 2019 with different telescope configurations. We compared the image reconstructions with 81 temporal snapshots of 3D radiative-hydrodynamics (RHD) CO$^5$BOLD simulations in terms of contrast and morphology, using the Structural Similarity Index. Results. The interferometric data are compatible with an overall spherical disk of angular diameter 4.4$\pm$0.2 mas, and an extended molecular layer. In 2016, the reconstructed image reveals a bright arc-like feature toward the northern rim of the photospheric surface. In 2019, an arc-like feature is seen at a different orientation and a new peak of emission is detected on the opposite side. The contrasts of the reconstructed surface images are 11%$\pm$2% and 9%$\pm$2% for 2016 and 2019, respectively. The morphology and contrast of the two images are consistent with 3D RHD simulations, within our achieved spatial resolution and dynamic range. The extended molecular layer contributes 10--13% of the total flux with an angular diameter of 6--8 mas. It is present but not clearly visible in the reconstructed images because it is close to the limits of the achieved dynamic range. The presence of the molecular layer is not reproduced by the 3D RHD simulations. Conclusions...

Over the past 20 years, KASCADE and its extension KASCADE-Grande were dedicated to measure high-energy cosmic rays with primary energies of 100 TeV to 1 EeV. The data accumulation was fully completed and all experimental components were dismantled, though the analysis of the high-quality data is still continued. E.g., we investigated the validity of the hadronic interaction model of the new SIBYLL version 2.3c. We also published a new result of a search for large-scale anisotropies performed with the KASCADE-Grande data. Investigation of the attenuation length of the muon in the atmosphere is also updated with the predictions of the SIBYLL 2.3 interaction model. We investigated, in addition, the muon content of high-energy air showers and compared them to all post-LHC interaction models. In this contribution, the new and updated results from KASCADE-Grande will be presented. An update of the web-based data center KCDC offering the original scientific data from KASCADE-Grande to the public will be briefly discussed as well.

Red supergiants are cool massive stars and are the largest and the most luminous stars in the universe. They are characterized by irregular or semi-regular photometric variations, the physics of which is not clearly understood. The paper aims at deriving the velocity field in the red supergiant star $\mu$ Cep and relating it to the photometric variability with the help of the tomographic method. The tomographic method allows to recover the line-of-sight velocity distribution over the stellar disk and within different optical-depth slices. The method is applied to a series of high-resolution spectra of $\mu$ Cep, and these results are compared to those obtained from 3D radiative-hydrodynamics CO5BOLD simulations of red supergiants. Fluctuations in the velocity field are compared with photometric and spectroscopic variations, the latter being derived from the TiO band strength and serving (at least partly) a proxy of the variations in effective temperature. The tomographic method reveals a phase shift between the velocity and spectroscopic/photometric variations. This phase shift results in a hysteresis loop in the temperature - velocity plane, with a timescale of a few hundred days, similar to the photometric one. The similarity between the hysteresis loop timescale measured in $\mu$ Cep and the timescale of acoustic waves disturbing the convective pattern suggests that such waves play an important role in triggering the hysteresis loops.

Context. Spectroscopy of exoplanet atmospheres at high resolving powers is rapidly gaining popularity to measure the presence of atomic and molecular species. While this technique is robust against contaminant absorption in the Earth's atmosphere, the non stationary stellar spectrum creates a non-negligible source of noise that can alter or even prevent detection. Aims. We aim at using three-dimensional stellar simulations to directly remove the signature of the star from observations, and prior to cross correlation with templates for the planet's atmosphere Methods. We compute synthetic spectra from 3D simulations of stellar convection and we couple them with an analytical model reproducing the correct geometry of a transiting exoplanet. We apply the method to HD 189733, and analyse transmission and emission spectroscopy of its hosted exoplanet. In addition, we analyse emission spectroscopy of the non transiting exoplanet 51 Pegasi b. Results. We find a significant improvement in the planet detectability when removing the stellar spectrum with our method. We show that the method is superior to a simple parametrisation of the stellar line profile or to the use of one-dimensional stellar models. This is due to the intrinsic treatment of convection, which allows us to correctly reproduce asymmetric spectral lines, and model center-to-limb variation and Rossiter-McLaughlin effect. In the case of 51 Pegasi b, we succeed at confirming a previous tentative detection of the planet's K-band spectrum due to the improved suppression of stellar residuals. Future high-resolution observations will benefit from the synergy with stellar spectroscopy, and can be used to test the correct modelling of physical processes in stellar atmospheres. We highlight key improvements in modelling techniques and knowledge of opacity sources to extend this work to shorter wavelengths and later-type stars.

The detection of gravitational waves and neutrinos from astrophysical sources with gamma-ray counterparts officially started the era of Multi-Messenger Astronomy. Their transient and extreme nature implies that monitoring the VHE sky is fundamental to investigate the non-electromagnetic signals. However, the limited effective area of space-borne instruments prevents observations above a few hundred GeV, while the small field of view and low duty cycle of IACTs make them unsuited for extensive monitoring activities and prompt response to transients. Extensive Air Shower arrays (EAS) can provide a large field of view, a wide effective area and a very high duty cycle. Their main difficulty is the distinction between gamma-ray and cosmic-ray initiated air showers, especially below the TeV range. Here we present some case studies stressing the importance that a new EAS array in the Southern Hemisphere will be able to survey the sky from below 100 GeV up to several TeV. In the energy domain between 100 and 400 GeV we expect the strongest electromagnetic signatures of the acceleration of ultra-relativistic particles in sources like SNRs, blazar jets and gamma-ray bursts, as recently proved by IACT observations. This spectral window is also crucial to understand the Universe opacity to high energy radiation, thus providing constraints on the cosmological parameters. We will discuss the implications of VHE radiation on the mechanisms at work and we will focus on the advantages resulting from the ability to monitor the energy window lying between the domain of space-borne detectors and ground-based facilities.

Core-collapse supernovae exploding in dense winds are favorable sites for cosmic-ray (CR) acceleration to very high energies. We present our CR-radiation-hydrodynamics simulations of the explosion of a red supergiant. We study the evolution of the shock wave during the first day following core collapse, and estimate the time at which CR acceleration can start. We then calculate the maximum CR energy at the forward shock as a function of time, and show that it may already exceed 100 TeV only a few hours after shock breakout from the surface of the star.

Theory surrounding the origin of the dust-laden winds from evolved stars remains mired in controversy. Characterizing the formation loci and the dust distribution within approximately the first stellar radius above the surface is crucial for understanding the physics that underlie the mass-loss phenomenon. By exploiting interferometric polarimetry, we derive the fundamental parameters that govern the dust structure at the wind base of a red supergiant. We present near-infrared aperture-masking observations of Betelgeuse in polarimetric mode obtained with the NACO/SAMPol instrument. We used both parametric models and radiative transfer simulations to predict polarimetric differential visibility data and compared them to SPHERE/ZIMPOL measurements. Using a thin dust shell model, we report the discovery of a dust halo that is located at only 0.5 R$_{\star}$ above the photosphere (i.e. an inner radius of the dust halo of 1.5 R$_{\star}$). By fitting the data under the assumption of Mie scattering, we estimate the grain size and density for various dust species. By extrapolating to the visible wavelengths using radiative transfer simulations, we compare our model with SPHERE/ZIMPOL data and find that models based on dust mixtures that are dominated by forsterite are most favored. Such a close dusty atmosphere has profound implications for the dust formation mechanisms around red supergiants.

We describe plans for the development of the Southern Wide-field Gamma-ray Observatory (SWGO), a next-generation instrument with sensitivity to the very-high-energy (VHE) band to be constructed in the Southern Hemisphere. SWGO will provide wide-field coverage of a large portion of the southern sky, effectively complementing current and future instruments in the global multi-messenger effort to understand extreme astrophysical phenomena throughout the universe. A detailed description of science topics addressed by SWGO is available in the science case white paper [1]. The development of SWGO will draw on extensive experience within the community in designing, constructing, and successfully operating wide-field instruments using observations of extensive air showers. The detector will consist of a compact inner array of particle detection units surrounded by a sparser outer array. A key advantage of the design of SWGO is that it can be constructed using current, already proven technology. We estimate a construction cost of 54M USD and a cost of 7.5M USD for 5 years of operation, with an anticipated US contribution of 20M USD ensuring that the US will be a driving force for the SWGO effort. The recently formed SWGO collaboration will conduct site selection and detector optimization studies prior to construction, with full operations foreseen to begin in 2026. Throughout this document, references to science white papers submitted to the Astro2020 Decadal Survey with particular relevance to the key science goals of SWGO, which include unveiling Galactic particle accelerators [2-10], exploring the dynamic universe [11-21], and probing physics beyond the Standard Model [22-25], are highlighted in red boldface.

Cool evolved stars are the main source of chemical enrichment of the interstellar medium. Understanding their mass loss offers a unique opportunity to study the cycle of matter. We discuss interferometric studies and their comparison to latest state-of-the-art dynamic model atmospheres. They show broad agreement for asymptotic giant branch stars. For red supergiants, however, current models cannot explain observed extensions by far, pointing to missing physical processes in their models, and uncertainties in our general understanding of mass loss. We present ongoing imaging and time-series observations that may provide the strongest constraint and may help to identify missing dynamic processes. VLTI studies will remain the highest spatial resolution observations at ESO into the ELT era, complemented by ALMA observations. We discuss crucial improvements in both instrumental and operational areas for surface imaging of cool evolved stars in the era of the ELT.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

The signature of wind patterns caused by the interplay of rotation and energy redistribution in hot Jupiters is detectable at high spectral resolution, yet no direct comparison has been attempted between predictions from general circulation models (GCMs) and observed high-resolution spectra. We present the first of such comparisons on near-infrared transmission spectra of the hot Jupiter HD 189733b. Exploring twelve rotation rates and two chemical regimes, we have created model spectra from 3-D GCMs and cross-correlated them with the observed spectra. Comparing our models against those of HD 189733b, we obtain three key results: (1) we confirm CO and H$_2$O in the planet's atmosphere at a detection significance of 8.2$\sigma$; (2) we recover the signature of $\sim$km/s day-to-night winds at $\sim$mbar pressures; and (3) we constrain the rotation period of the planet to between 1.2 and 4.69 days (synchronous rotation (2.2 days) remains consistent with existing observations). Our results do not suffer from the shortcomings of 1-D models as cross correlation templates - mainly that these models tend to over-constrain the slower rotation rates and show evidence for anomalous blue shifts. Our 3-D models instead match the observed line-of-sight velocity of this planet by self-consistently including the effects of high-altitude day-to-night winds. Overall, we find a high degree of consistency between HD 189733b observations and our GCM-based spectra, implying that the physics and chemistry are adequately described in current 3-D forward models for the purpose of interpreting observations at high spectral resolution.

The paper is a script of a lecture given at the ISAPP-Baikal summer school in 2018. The lecture gives an overview of the Tunka Advanced Instrument for cosmic rays and Gamma Astronomy (TAIGA) facility including historical introduction, description of existing and future setups, and outreach and open data activities.

The polarization state of starlight reflected by a planetary atmosphere uniquely reveals coverage, particle size, and composition of aerosols as well as changing cloud patterns. It is not possible to obtain a comparable level of detailed from flux-only observations. Furthermore, polarization observations can probe the atmosphere of planets independently of the orbital geometry (i.e., transiting and non-transiting planets). We show that a high-resolution spectropolarimeter with a broad wavelength coverage, particularly if attached to a large space telescope, would enable simultaneous study of the polarimetric exoplanet properties of the continuum and to look for and characterize the polarimetric signal due to scattering from single molecules.

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.

In this manuscript I summarise the work I have carried out in the last 10 years. My research projects are articulated around various fields in astrophysics having, as unifying theme, the 3D modelling and the observation of stellar atmospheres. In particular, my field of research covers all late type stars, from the Main Sequence stars to the advanced stages (red giant and supergiant stars). My work is based on the quantitative analysis of images, synthetic spectra, bolometric luminosity, planetary transits of different types of stars and on the modeling of their atmospheres. They include the use and development of 3D simunlations that are necessary for the interpretation of the data obtained with the new instruments equipping the space borne and ground telescopes. Interferometry and spectroscopy hold a privileged place in my research. My research themes concern interdisciplinary aspects by addressing questions related to atmospheric physics in relation to 3D hydrodynamic models. The predictions thus obtained are used to drive a vast campaign of interferometric, astrometric, photometric, spectroscopic and imaging observations. Finally, 3D simulations are also directly used to drive the development of new generation instrument at the VLTI, as well as very large future observatories and are directly linked to the ESA Gaia and PLATO missions.

In this paper we present the most promising science cases for a new generation visible instrument on the VLTI and the conceptual idea for the instrumental configuration. We also present a statistical study of the potential targets that may be accessible for the different classes of objects and for the required spectral resolutions.

A tomographic method, aiming at probing velocity fields at depth in stellar atmospheres, is applied to the red supergiant star {\mu} Cep and to snapshots of 3D radiative-hydrodynamics simulation in order to constrain atmospheric motions and relate them to photometric variability.

GALAH is a large-scale magnitude-limited southern stellar spectroscopic survey. Its second data release (GALAH DR2) provides values of stellar parameters and abundances of 23 elements for 342,682 stars (Buder et al.). Here we add a description of the public release of radial velocities with a typical accuracy of 0.1 km/s for 336,215 of these stars, achievable due to the large wavelength coverage, high resolving power and good signal to noise ratio of the observed spectra, but also because convective motions in stellar atmosphere and gravitational redshift from the star to the observer are taken into account. In the process we derive medians of observed spectra which are nearly noiseless, as they are obtained from between 100 and 1116 observed spectra belonging to the same bin with a width of 50 K in temperature, 0.2 dex in gravity, and 0.1 dex in metallicity. Publicly released 1181 median spectra have a resolving power of 28,000 and trace the well-populated stellar types with metallicities between -0.6 and +0.3. Note that radial velocities from GALAH are an excellent match to the accuracy of velocity components along the sky plane derived by Gaia for the same stars. The level of accuracy achieved here is adequate for studies of dynamics within stellar clusters, associations and streams in the Galaxy. So it may be relevant for studies of the distribution of dark matter.

Context. The contribution of quiet-Sun regions to the solar irradiance variability is currently unclear. Some solar-cycle variations of the quiet-Sun physical structure, such as the temperature gradient, might affect the irradiance. The synoptic measurement of this quantity along the activity cycle would improve our understanding of long-term irradiance variations. Aims. We intend to test a method previously introduced for measuring the photospheric temperature gradient from high-resolution spectroscopic observation and to study its feasibility with ground-based instruments with and without adaptative optics. Methods. We used synthetic profiles of the FeI 630.15 nm obtained from realistic three-dimensional hydrodynamical simulations of the photospheric granulation and line radiative transfer computations under local thermodynamical equilibrium conditions. Synthetic granulation images at different levels in the line are obtained by convolution with the instrumental point spread function (PSF) under various conditions of atmospheric turbulence, with and without correction by an adaptative optics (AO) system. The PSF are obtained with the PAOLA software, and the AO performances are inspired by the system that will be operating on the Daniel K. Inouye Solar Telescope. Results. We consider two different conditions of atmospheric turbulence, with Fried parameters of 7 cm and 5 cm, respectively. We show that the degraded images lead to both a bias and a loss of precision in the temperature-gradient measurement, and that the correction with the AO system allows us to drastically improve the measurement quality. Conclusions. Long-term synoptic observations of the temperature gradient in the solar photosphere can be undertaken by implement- ing this method on ground-based solar telescopes that are equipped with an AO correction system.

The dynamical structure of the atmosphere of Cepheids has been well studied in the optical. Several authors have found very interesting spectral features in the J band, but little data have been secured beyond 1.6 um. However, such observations can probe different radial velocities and line asymmetry regimes, and are able to provide crucial insights into stellar physics. Our goal was to investigate the infrared line-forming region in the K-band domain, and its impact on the projection factor and the k-term of Cepheids. We secured CRIRES observations for the long-period Cepheid l Car, with a focus on the unblended spectral line NaI2208.969 nm. We measured the corresponding radial velocities (by using the first moment method) and the line asymmetries (by using the bi-Gaussian method). These quantities are compared to the HARPS visible spectra we previously obtained on l Car. The optical and infrared radial velocity curves show the same amplitude (only about 3% of difference), with a slight radial velocity shift of about 0.5 +/- 0.3 km s^-1 between the two curves. Around the minimum radius (phase ~ 0.9) the visible radial velocity curve is found in advance compared to the infrared one (phase lag), which is consistent with an infrared line forming higher in the atmosphere (compared to the visible line) and with a compression wave moving from the bottom to the top of the atmosphere during maximum outward velocity. The asymmetry of the K-band line is also found to be significantly different from that of the optical line.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Context. Asymptotic Giant Branch stars are characterised by complex stellar surface dynamics that affect the measurements and amplify the uncertainties on stellar parameters. The resulting motion of the stellar photo-center could have adverse effects on the parallax determination with Gaia. Aims. We explore the impact of the convection-related surface structure in AGBs on the photocentric variability. We quantify these effects to characterise the observed parallax errors and estimate fundamental stellar parameters and dynamical properties. Methods. We use 3D radiative-hydrodynamics simulations of convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the Gaia G band. From those maps, we calculate the intensity-weighted mean of all emitting points tiling the visible stellar surface and evaluate its motion as a function of time. We extract the parallax error from Gaia DR2 for a sample of semiregular variables in the solar neighbourhood and compare it to the synthetic predictions of photo-center displacements. Results. AGB stars show a complex surface morphology characterised by the presence of few large scale long-lived convective cells accompanied by short-lived and small scale structures. As a consequence, the position of the photo-center displays temporal excursions between 0.077 to 0.198 AU (5 to 11% of the corresponding stellar radius), depending on the simulation considered. We show that the convection-related variability accounts for a substantial part to the Gaia DR2 parallax error of our sample of semiregular variables. Finally, we put in evidence for a correlation between the mean photo-center displacement and the stellar fundamental parameters: surface gravity, radius, and pulsation. We denote that parallax variations could be exploited quantitatively using appropriate RHD simulations corresponding to the observed star.

The `KASCADE Cosmic ray Data Centre' is a web portal (\url{https://kcdc.ikp.kit.edu}), where the data of the astroparticle physics experiment KASCADE-Grande are made available for the interested public. The KASCADE experiment was a large-area detector for the measurement of high-energy cosmic rays via the detection of extensive air showers. The multi-detector installations KASCADE and its extension KASCADE-Grande stopped the active data acquisition in 2013 of all its components end of 2012 after more than 20 years of data taking. In several updates since our first release in 2013 with KCDC we provide the public measured and reconstructed parameters of more than 433 million air showers. In addition, KCDC provides meta data information and documentation to enable a user outside the community of experts to perform their own data analysis. Simulation data from three different high energy interaction models have been made available as well as a compilation of measured and published spectra from various experiments. In addition, detailed educational examples shall encourage high-school students and early stage researchers to learn about astroparticle physics, cosmic radiation as well as the handling of Big Data and about the sustainable and public provision of scientific data.

Dynamic model atmospheres of Mira stars predict variabilities in the photospheric radius and in atmospheric molecular layers which are not yet strongly constrained by observations. Here we measure the variability of the oxygen-rich Mira star R Peg. We used K band spectro-interferometry with a spectral resolution of about 4000 obtained at four epochs between post-maximum and minimum visual phases employing the newly available GRAVITY beam combiner at the VLTI. Our observations show a continuum radius that is anti-correlated with the visual lightcurve. Uniform disc (UD) angular diameters at a near-continuum wavelength of 2.25 mu are steadily increasing with values of 8.7 mas, 9.4 mas, 9.8 mas, and 9.9 mas at visual phases of 0.15, 0.36, 0,45, 0.53, respectively. UD diameters at a bandpass around 2.05 mu, dominated by water vapour, follow the near-continuum variability at larger UD diameters. UD diameters at the CO 2-0 bandhead, instead, are correlated with the visual lightcurve and anti-correlated with the near-continuum UD diameters. The observed anti-correlation between continuum radius and visual lightcurve is consistent with an earlier study of the oxygen-rich Mira S Lac, and with recent 1D CODEX dynamic model atmosphere predictions. The amplitude of the variation is comparable to the earlier observations of S Lac, and smaller than predicted by CODEX models. The wavelength-dependent visibility variations at our epochs can be reproduced by a set of CODEX models at model phases between 0.3 and 0.6. The anti-correlation of water vapour and CO contributions at our epochs suggests that these molecules undergo different processes along the stellar cycle. The newly available GRAVITY instrument is suited to conducting longer time series observations, which are needed to provide strong constraints on the model-predicted intra- and inter-cycle variability.

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 evolution of the muon content of very high energy air showers (EAS) in the atmosphere is investigated with data of the KASCADE-Grande observatory. For this purpose, the muon attenuation length in the atmosphere is obtained to $\Lambda_\mu = 1256 \, \pm 85 \, ^{+229}_{-232}(\mbox{syst})\, \mbox{g/cm}^2$ from the experimental data for shower energies between $10^{16.3}$ and $10^{17.0} \, \mbox{eV}$. Comparison of this quantity with predictions of the high-energy hadronic interaction models QGSJET-II-02, SIBYLL 2.1, QGSJET-II-04 and EPOS-LHC reveals that the attenuation of the muon content of measured EAS in the atmosphere is lower than predicted. Deviations are, however, less significant with the post-LHC models. The presence of such deviations seems to be related to a difference between the simulated and the measured zenith angle evolutions of the lateral muon density distributions of EAS, which also causes a discrepancy between the measured absorption lengths of the density of shower muons and the predicted ones at large distances from the EAS core. The studied deficiencies show that all four considered hadronic interaction models fail to describe consistently the zenith angle evolution of the muon content of EAS in the aforesaid energy regime.

Large stellar surveys of the Milky Way require validation with reference to a set of "benchmark" stars whose fundamental properties are well-determined. For metal-poor benchmark stars, disagreement between spectroscopic and interferometric effective temperatures has called the reliability of the temperature scale into question. We present new interferometric measurements of three metal-poor benchmark stars, HD 140283, HD 122563, and HD 103095, from which we determine their effective temperatures. The angular sizes of all the stars were determined from observations with the PAVO beam combiner at visible wavelengths at the CHARA array, with additional observations of HD 103095 made with the VEGA instrument, also at the CHARA array. Together with photometrically derived bolometric fluxes, the angular diameters give a direct measurement of the effective temperature. For HD 140283 we find {\theta}_LD = 0.324+/-0.005 mas, Teff = 5787+/-48 K; for HD 122563, {\theta}_LD = 0.926+/-0.011 mas, Teff = 4636+/-37 K; and for HD 103095 {\theta}_LD = 0.595+/-0.007 mas, Teff = 5140+/-49 K. Our temperatures for HD 140283 and HD 103095 are hotter than the previous interferometric measurements by 253 K and 322 K, respectively. We find good agreement between our temperatures and recent spectroscopic and photometric estimates. We conclude some previous interferometric measurements have been affected by systematic uncertainties larger than their quoted errors.

Context. The surface structures and dynamics of cool stars are characterized by the presence of convective motions and turbulent flows which shape the emergent spectrum. Aims. We used realistic three-dimensional radiative hydrodynamical simulations from the Stagger-grid to calculate synthetic spectra with the radiative transfer code Optim3D for stars with different stellar parameters to predict photometric colors and convective velocity shifts. Methods. We calculated spectra from 1000 to 200 000 A with a constant resolving power of 20 000 and from 8470 and 8710 A (Gaia Radial Velocity Spectrometer - RVS - spectral range), with a constant resolving power of 300 000. Results. We used synthetic spectra to compute theoretical colors in the Johnson-Cousins UBV(RI), SDSS, 2MASS, Gaia, SkyMapper, Stromgren systems, and HST-WFC3. We showed that 1D versus 3D differences are limited to a small percent except for the narrow filters that span the optical and UV region of the spectrum. In addition, we derived the effect of the convective velocity fields on selected Fe I lines. We found the overall convective shift for 3D simulations with respect to the reference 1D hydrostatic models, revealing line shifts of between -0.235 and +0.361 km/s. We showed a net correlation of the convective shifts with the effective temperature: lower effective temperatures denote redshifts and higher effective temperatures denote blueshifts. We conclude that the extraction of accurate radial velocities from RVS spectra need an appropriate wavelength correction from convection shifts. Conclusions. The use of realistic 3D hydrodynamical stellar atmosphere simulations has a small but significant impact on the predicted photometry compared with classical 1D hydrostatic models for late-type stars. We make all the spectra publicly available for the community through the POLLUX database.

Cool giant and supergiant star atmospheres are characterized by complex velocity fields originating from convection and pulsation processes which are not fully understood yet. The velocity fields impact the formation of spectral lines, which thus contain information on the dynamics of stellar atmospheres. The tomographic method allows to recover the distribution of the component of the velocity field projected on the line of sight at different optical depths in the stellar atmosphere. The computation of the contribution function to the line depression aims at correctly identifying the depth of formation of spectral lines in order to construct numerical masks probing spectral lines forming at different optical depths. The tomographic method is applied to 1D model atmospheres and to a realistic 3D radiative hydrodynamics simulation performed with CO5BOLD in order to compare their spectral line formation depths and velocity fields. In 1D model atmospheres, each spectral line forms in a restricted range of optical depths. On the other hand, in 3D simulations, the line formation depths are spread in the atmosphere mainly because of temperature and density inhomogeneities. Comparison of CCF profiles obtained from 3D synthetic spectra with velocities from the 3D simulation shows that the tomographic method correctly recovers the distribution of the velocity component projected on the line of sight in the atmosphere.

Cosmic rays provide an unique approach to study hadronic interactions at high energies in the kinematic forward direction. The KASCADE air shower experiment was the first to conduct quantitative tests of hadronic interactions with air shower data. A brief overview is given on results from KASCADE and its extension KASCADE-Grande with respect to investigations of hadronic interactions and the properties of cosmic rays.

KASCADE and KASCADE-Grande were multi-detector installations to measure individual air showers of cosmic rays at ultra-high energy. Based on data sets measured by KASCADE and KASCADE-Grande, 90% C.L. upper limits to the flux of gamma-rays in the primary cosmic ray flux are determined in an energy range of ${10}^{14} - {10}^{18}$ eV. The analysis is performed by selecting air showers with a low muon content as expected for gamma-ray-induced showers compared to air showers induced by energetic nuclei. The best upper limit of the fraction of gamma-rays to the total cosmic ray flux is obtained at $3.7 \times {10}^{15}$ eV with $1.1 \times {10}^{-5}$. Translated to an absolute gamma-ray flux this sets constraints on some fundamental astrophysical models, such as the distance of sources for at least one of the IceCube neutrino excess models.

AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory designed to extend its energy range of detection and to directly measure the muon content of the cosmic ray primary particle showers. The array will be formed by an infill of surface water-Cherenkov detectors associated with buried scintillation counters employed for muon counting. Each counter is composed of three scintillation modules, with a 10 m$^2$ detection area per module. In this paper, a new generation of detectors, replacing the current multi-pixel photomultiplier tube (PMT) with silicon photo sensors (aka. SiPMs), is proposed. The selection of the new device and its front-end electronics is explained. A method to calibrate the counting system that ensures the performance of the detector is detailed. This method has the advantage of being able to be carried out in a remote place such as the one where the detectors are deployed. High efficiency results, i.e. 98 % efficiency for the highest tested overvoltage, combined with a low probability of accidental counting ($\sim$2 %), show a promising performance for this new system.

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

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.

We obtained 16 VLT/X-shooter observations of GX 339-4 in quiescence in the period May - September 2016 and detected absorption lines from the donor star in its NIR spectrum. This allows us to measure the radial velocity curve and projected rotational velocity of the donor for the first time. We confirm the 1.76 day orbital period and we find that $K_2$ = $219 \pm 3$ km s$^{-1}$, $\gamma = 26 \pm 2$ km s$^{-1}$ and $v \sin i = 64 \pm 8$ km s$^{-1}$. From these values we compute a mass function $f(M) =1.91 \pm 0.08~M_{\odot}$, a factor $\sim 3$ lower than previously reported, and a mass ratio $q = 0.18 \pm 0.05$. We confirm the donor is a K-type star and estimate that it contributes $\sim 45-50\%$ of the light in the $J$- and H-band. We constrain the binary inclination to $37^\circ < i < 78^\circ$ and the black hole mass to $2.3~M_{\odot} < M_\mathrm{BH} < 9.5~M_{\odot}$. GX 339-4 may therefore be the first black hole to fall in the 'mass-gap' of $2-5~M_{\odot}$.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

High spatial resolution is the key for the understanding various astrophysical phenomena. But even with the future E-ELT, single dish instruments are limited to a spatial resolution of about 4 mas in the visible. For the closest objects within our Galaxy most of the stellar photosphere remains smaller than 1 mas. With the success of long baseline interferometry these limitations were soom overcome. Today low and high resolution interferometric instruments on the VLTI and CHARA offer an immense range of astrophysical studies. Combining more telescopes and moving to visible wavelengths broadens the science cases even more. With the idea of developing strong science cases for a future visible interferometer, we organized a science group around the following topics: pre-main sequence and main sequence stars, fundamental parameters, asteroseismology and classical pulsating stars, evolved stars, massive stars, active galactic nuclei (AGNs) and imaging techniques. A meeting was organized on the 15th and 16th of January, 2015 in Nice with the support of the Action Specific in Haute Resolution Angulaire (ASHRA), the Programme National en Physique Stellaire (PNPS), the Lagrange Laboratory and the Observatoire de la Cote d'Azur, in order to present these cases and to discuss them further for future visible interferometers. This White Paper presents the outcome of the exchanges. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41.

The Extreme Energy Events Project is a synchronous sparse array of 52 tracking detectors for studying High Energy Cosmic Rays (HECR) and Cosmic Rays-related phenomena. The observatory is also meant to address Long Distance Correlation (LDC) phenomena: the network is deployed over a broad area covering 10 degrees in latitude and 11 in longitude. An overview of a set of preliminary results is given, extending from the study of local muon flux dependance on solar activity to the investigation of the upward-going component of muon flux traversing the EEE stations; from the search for anisotropies at the sub-TeV scale to the hints for observations of km-scale Extensive Air Shower (EAS).

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.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Experimental complex (EC) NEVOD includes a number of unique experimental facilities for studies of main components of cosmic rays on the Earth's surface. The complex is used for the basic research of CR flux characteristics and their interactions in the energy range 10^15 - 10^19 eV, and for applied investigations directed to the development of methods of the muon diagnostics of the atmosphere and the Earth's magnetosphere and near-terrestrial space. To extend the experimental capabilities and raising the status of the installation to the Mega Science level, nowadays new large-scale detectors: array for the EAS registration - NEVOD-EAS, detector of atmospheric neutrons - URAN, and large-area coordinate-tracking detector - TREK, are being deployed around EC NEVOD. The description of new detectors and a common trigger system to ensure the joint operation together with other detectors of EC NEVOD are presented.

A new cluster type shower array NEVOD-EAS is designed for estimating the size, axis position and arrival direction of extensive air showers registered by the Cherenkov water detector NEVOD and coordinate detector DECOR (Moscow, Russia). In 2015-2016, the central part of the array was deployed and started its operation. It includes 4 independent clusters of scintillation detector stations located around the NEVOD-DECOR experimental complex on the area of about $10^{4} m^{2}$. This article presents the results of studying amplitude and timing characteristics of the array clusters which are critical for EAS parameters reconstruction, as well as the examples of registered events.

We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the electromagnetic and muonic components. It allows an accurate determination of the spread of primary masses in the cosmic-ray flux. Up till now, constraints on the spread of primary masses have been dominated by systematic uncertainties. The present correlation measurement is not affected by systematics in the measurement of the depth of shower maximum or the signal in the water Cherenkov stations. The analysis relies on general characteristics of air showers and is thus robust also with respect to uncertainties in hadronic event generators. The observed correlation in the energy range around the `ankle' at $\lg(E/{\rm eV})=18.5-19.0$ differs significantly from expectations for pure primary cosmic-ray compositions. A light composition made up of proton and helium only is equally inconsistent with observations. The data are explained well by a mixed composition including nuclei with mass $A > 4$. Scenarios such as the proton dip model, with almost pure compositions, are thus disfavoured as the sole explanation of the ultrahigh-energy cosmic-ray flux at Earth.

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

The radio technique is a promising method for detection of cosmic-ray air showers of energies around $100\,$PeV and higher with an array of radio antennas. Since the amplitude of the radio signal can be measured absolutely and increases with the shower energy, radio measurements can be used to determine the air-shower energy on an absolute scale. We show that calibrated measurements of radio detectors operated in coincidence with host experiments measuring air showers based on other techniques can be used for comparing the energy scales of these host experiments. Using two approaches, first via direct amplitude measurements, and second via comparison of measurements with air shower simulations, we compare the energy scales of the air-shower experiments Tunka-133 and KASCADE-Grande, using their radio extensions, Tunka-Rex and LOPES, respectively. Due to the consistent amplitude calibration for Tunka-Rex and LOPES achieved by using the same reference source, this comparison reaches an accuracy of approximately $10\,\%$ - limited by some shortcomings of LOPES, which was a prototype experiment for the digital radio technique for air showers. In particular we show that the energy scales of cosmic-ray measurements by the independently calibrated experiments KASCADE-Grande and Tunka-133 are consistent with each other on this level.

We add four warmer late-type supergiants to our previous spectro-interferometric studies of red giants and supergiants. V766 Cen (=HR 5171 A) is found to be a high-luminosity log(L/L_sun)=5.8+-0.4 source of Teff 4290+-760 K and radius 1490+-540 Rsun located close to both the Hayashi and Eddington limits; this source is consistent with a 40 Msun evolutionary track without rotation and current mass 27-36 Msun. It exhibits NaI in emission arising from a shell of radius 1.5 Rphot and a photocenter displacement of about 0.1 Rphot. V766 Cen shows strong extended molecular (CO) layers and a dusty circumstellar background component. This suggest an optically thick pseudo-photosphere at about 1.5 Rphot at the onset of the wind. V766 Cen is a red supergiant located close to the Hayashi limit instead of a yellow hypergiant already evolving back toward warmer Teff as previously discussed. The stars sigma Oph, BM Sco, and HD 206859 are found to have lower luminosities of about log(L/Lsun)=3.4-3.5 and Teff of 3900-5300 K, corresponding to 5-9 Msun tracks. They do not show extended molecular layers as observed for higher luminosity red supergiants of our sample. BM Sco shows an unusually strong contribution by an over-resolved circumstellar dust component. These stars are more likely high-mass red giants instead of red supergiants. This leaves us with an unsampled locus in the HR diagram corresponding to luminosities log(L/Lsun)~3.8-4.8 or masses 10-13 Msun, possibly corresponding to the mass region where stars explode as type II-P supernovae during the RSG stage. Our previously found relation of increasing strength of extended molecular layers with increasing luminosities is now confirmed to extend to double our previous luminosities and up to the Eddington limit. This might further point to steadily increasing radiative winds with increasing luminosity. [Abridged]

Stellar activity and convection-related surface structures might cause bias in planet detection and characterization that use these transits. Surface convection simulations help to quantify the granulation signal. We used realistic three-dimensional radiative hydrodynamical simulations from the Stagger grid and synthetic images computed with the radiative transfer code Optim3D to model the transits of three prototype planets: a hot Jupiter, a hot Neptune, and a terrestrial planet. We computed intensity maps from RHD simulations of the Sun and a K-dwarf star at different wavelength bands from optical to far-infrared. We modeled the transit using synthetic stellar-disk images and emulated the temporal variation of the granulation intensity. We identified two types of granulation noise that act simultaneously during the planet transit: (i) the intrinsic change in the granulation pattern with timescales smaller than the usual planet transit, and (ii) the fact that the transiting planet occults isolated regions of the photosphere that differ in local surface brightness. Our modeling approach shows that the granulation pattern has a non-negligible effect on the light curve depth during the transit, and, consequentially on the determination of the radius of the planet transiting. The granulation noise appears to be correlated among the different wavelength ranges either in the visible or in the infrared regions. The granulation has to be considered as an intrinsic uncertainty (as a result of stellar variability) on the precise measurements of exoplanet transits of planets. The full characterization of the granulation is essential for determining the degree of uncertainty on the planet parameters. In this context, the use of 3D RHD simulations is important to measure the convection-related fluctuations.

We present an overview of the scientific potential of MATISSE, the Multi Aperture mid-Infrared SpectroScopic Experiment for the Very Large Telescope Interferometer. For this purpose we outline selected case studies from various areas, such as star and planet formation, active galactic nuclei, evolved stars, extrasolar planets, and solar system minor bodies and discuss strategies for the planning and analysis of future MATISSE observations. Moreover, the importance of MATISSE observations in the context of complementary high-angular resolution observations at near-infrared and submillimeter/millimeter wavelengths is highlighted.

The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

MATISSE represents a great opportunity to image the environment around massive and evolved stars. This will allow one to put constraints on the circumstellar structure, on the mass ejection of dust and its reorganization , and on the dust-nature and formation processes. MATISSE measurements will often be pivotal for the understanding of large multiwavelength datasets on the same targets collected through many high-angular resolution facilities at ESO like sub-millimeter interferometry (ALMA), near-infrared adaptive optics (NACO, SPHERE), interferometry (PIONIER, GRAVITY), spectroscopy (CRIRES), and mid-infrared imaging (VISIR). Among main sequence and evolved stars, several cases of interest have been identified that we describe in this paper.

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

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.

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec \theta)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec \theta)_\mathrm{max}$.

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

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.

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

We obtained a total of 20 near-infrared K-band spectro-interferometric snapshot observations of the Mira variables o Cet, R Leo, R Aqr, X Hya, W Vel, and R Cnc with a spectral resolution of about 1500. We compared observed flux and visibility spectra with predictions by CODEX 1D dynamic model atmospheres and with azimuthally averaged intensities based on CO5BOLD 3D dynamic model atmospheres including convection. Our visibility data confirm the presence of spatially extended molecular atmospheres located above the continuum radii with large-scale inhomogeneities or clumps that contribute a few percent of the total flux. The detailed structure of the inhomogeneities or clumps show a variability on time scales of 3 months and above. Both modeling attempts provided satisfactory fits to our data. In particular, they are both consistent with the observed decrease in the visibility function at molecular bands of water vapor and CO, indicating a spatially extended molecular atmosphere. Observational variability phases are mostly consistent with those of the best-fit CODEX models, except for near-maximum phases, where data are better described by near-minimum models. Rosseland angular diameters derived from the model fits are broadly consistent between those based on the 1D and the 3D models and with earlier observations. We derived fundamental parameters including absolute radii, effective temperatures, and luminosities for our sources. Our results provide a first observational support for theoretical results that shocks induced by convection and pulsation in the 3D CO5BOLD models of AGB stars are roughly spherically expanding and of similar nature to those of self-excited pulsations in 1D CODEX models. Unlike for red supergiants, the pulsation- and shock-induced dynamics can levitate the molecular atmospheres of Mira variables to extensions that are consistent with observations.

Cosmic ray air showers emit radio pulses at MHz frequencies, which can be measured with radio antenna arrays - like LOPES at the Karlsruhe Institute of Technology in Germany. To improve the understanding of the radio emission, we test theoretical descriptions with measured data. The observables used for these tests are the absolute amplitude of the radio signal, and the shape of the radio lateral distribution. We compare lateral distributions of more than 500 LOPES events with two recent and public Monte Carlo simulation codes, REAS 3.11 and CoREAS (v 1.0). The absolute radio amplitudes predicted by REAS 3.11 are in good agreement with the LOPES measurements. The amplitudes predicted by CoREAS are lower by a factor of two, and marginally compatible with the LOPES measurements within the systematic scale uncertainties. In contrast to any previous versions of REAS, REAS 3.11 and CoREAS now reproduce the shape of the measured lateral distributions correctly. This reflects a remarkable progress compared to the situation a few years ago, and it seems that the main processes for the radio emission of air showers are now understood: The emission is mainly due to the geomagnetic deflection of the electrons and positrons in the shower. Less important but not negligible is the Askaryan effect (net charge variation). Moreover, we confirm that the refractive index of the air plays an important role, since it changes the coherence conditions for the emission: Only the new simulations including the refractive index can reproduce rising lateral distributions which we observe in a few LOPES events. Finally, we show that the lateral distribution is sensitive to the energy and the mass of the primary cosmic ray particles.

LOPES was a digital antenna array detecting the radio emission of cosmic-ray air showers. The calibration of the absolute amplitude scale of the measurements was done using an external, commercial reference source, which emits a frequency comb with defined amplitudes. Recently, we obtained improved reference values by the manufacturer of the reference source, which significantly changed the absolute calibration of LOPES. We reanalyzed previously published LOPES measurements, studying the impact of the changed calibration. The main effect is an overall decrease of the LOPES amplitude scale by a factor of $2.6 \pm 0.2$, affecting all previously published values for measurements of the electric-field strength. This results in a major change in the conclusion of the paper 'Comparing LOPES measurements of air-shower radio emission with REAS 3.11 and CoREAS simulations' published in Astroparticle Physics 50-52 (2013) 76-91: With the revised calibration, LOPES measurements now are compatible with CoREAS simulations, but in tension with REAS 3.11 simulations. Since CoREAS is the latest version of the simulation code incorporating the current state of knowledge on the radio emission of air showers, this new result indicates that the absolute amplitude prediction of current simulations now is in agreement with experimental data.