Context. While Gaia enables to probe in great detail the extended local neighbourhood, the thin disk structure at larger distances remains sparsely explored. Aims. We aim here to build a non-parametric 3D model of the thin disc structures handling both the extinction and the stellar density simultaneously. Methods. We developed a Bayesian deconvolution method in two dimensions: extinction and distance. It uses a reference catalogue which completeness information defines the selection function. It is designed so that any complementary information from other catalogues can be added. It has also been designed to be robust to outliers, frequent in crowded fields, and differential extinction. The prior information is designed to be minimal: only a reference H-R diagram. We derived for this an empirical H-R diagram of the thin disk using Gaia DR2 data and synthetic isochrone-based H-R diagrams can also be used. Results. We validated the method on simulations and real fields using 2MASS and UKIDSS data complemented by Gaia DR2 photometry and parallaxes. We detail the results of two test fields: a 2MASS field centred around the NGC 4815 open cluster which shows an over-density of both extinction and stellar density at the cluster distance, and a UKIDSS field at $l=10^{\circ}$ where we recover the position of the Galactic bar.

Aims. We aim to map the 3D distribution of the interstellar extinction of the Milky Way disk up to distances larger than those probed with the Gaia parallax alone. Methods. We apply the FEDReD (Field Extinction-Distance Relation Deconvolver) algorithm to the 2MASS near-infrared photometry together with the Gaia DR2 astrometry and photometry. This algorithm uses a Bayesian deconvolution approach, based on an empirical HR-diagram representative of the local thin disk, in order to map the extinction as a function of distance of various fields of view. Results. We analysed more than 5.6 million stars to obtain an extinction map of the entire Galactic disk within $|b|<0.24^{\circ}$. This map provides information up to $5~\mathrm{kpc}$ in the direction of the Galactic centre and at more than $7~\mathrm{kpc}$ in the direction of the anticentre. This map reveals the complete shape of structures known locally, such as the Vela complex or the split of the local arm. Furthermore our extinction map shows many large "clean bubbles" especially one in the Sagittarius -- Carina complex, and four others which define a structure that we nickname the butterfly.

The nature of Milky Way dwarf spheroidals (MW dSphs) has been questioned, in particular whether they are dominated by dark matter (DM). Here we investigate an alternative scenario, for which tidal shocks are exerted by the MW to DM-free dSphs after a first infall of their gas-rich progenitors, and for which theoretical calculations have been verified by pure N-body simulations. Whether or not the dSphs are on their first infall cannot be resolved on the sole basis of their star formation history. In fact, gas removal may cause complex gravitational instabilities and near-pericenter passages can give rise to tidal disruptive processes. Advanced precision with the Gaia satellite in determining both their past orbital motions and the MW velocity curve is, however, providing crucial results. First, tidal shocks explain why DM-free dSphs are found preferentially near their pericenter, where they are in a destructive process, while their chance to be long-lived satellites is associated with a very low probability P~ 2 10^-7, which is at odds with the current DM-dominated dSph scenario. Second, most dSph binding energies are consistent with a first infall. Third, the MW tidal shocks that predict the observed dSph velocity dispersions are themselves predicted in amplitude by the most accurate MW velocity curve. Fourth, tidal shocks accurately predict the forces or accelerations exerted at half-light radius of dSphs, including the MW and the Magellanic System gravitational attractions. The above is suggestive of dSphs that are DM-free and tidally shocked near their pericenters, which may provoke a significant quake in our understanding of near-field cosmology.

The Upper-Scorpius association (5-11 Myr) contains a unique population of low-mass (M<30 MJup) brown-dwarfs. The detailed relative characterization of their physical properties (mass, radius, temperature, composition, ongoing accretion) offers the opportunity to potentially explore their origin and their mechanisms of formation. We aim at characterizing the chemical and physical properties of three young, late-M brown-dwarfs claimed to be companions of the Upper-Scorpius stars USco161031.9-16191305, HIP77900, and HIP78530 using medium resolution spectroscopy at UV (R~3300), optical (R~5400), and near-infrared (R~4300) wavelengths. The spectra of six free-floating analogues from the same association are analyzed for comparison and to explore the potential physical differences between these substellar objects. We also aim at analyzing emission lines at UV and optical wavelengths to investigate the presence of ongoing accretion processes. The X-Shooter spectrograph at VLT was used to obtain the spectra of our targets over the 0.3-2.5 $\mu$m range simultaneously. Performing a forward modelling of the observed spectra with the ForMoSA code (Nested Sampling Bayesian inference method in using BT-SETTL15 models), we infer the Teff, log(g), and radius of our objects. Our results are compatible with physical parameters found for free floating analogues of the association and with evolutionary-model predictions. However the final accuracy on the Teff estimates is strongly limited by non-reproducibilities of the BT-SETTL15 models for late-M brown-dwarfs. We identified emission lines in the spectrum of several objects attributed to chromospheric activity except for USco1608-2315 for which they are indicative of active accretion. We confirm the x4 over-luminosity of USco161031.9-16191305B. It could be explained by the object activity and if the companion is an unresolved multiple system.

We have found that the high velocity dispersions of dwarf spheroidal galaxies (dSphs) can be well explained by Milky Way (MW) tidal shocks, which reproduce precisely the gravitational acceleration previously attributed to dark matter (DM). Here we summarize the main results of Hammer et al. (2019) who studied the main scaling relations of dSphs and show how dark-matter free galaxies in departure from equilibrium reproduce them well, while they appear to be challenging for the DM model. These results are consistent with our most recent knowledge about dSph past histories, including their orbits, their past star formation history and their progenitors, which are likely tiny dwarf irregular galaxies.

High contrast imaging enables the determination of orbital parameters for substellar companions (planets, brown dwarfs) from the observed relative astrometry and the estimation of model and age-dependent masses from their observed magnitudes or spectra. Combining astrometric positions with radial velocity gives direct constraints on the orbit and on the dynamical masses of companions. A brown dwarf was discovered with the VLT/SPHERE instrument in 2017, which orbits at $\sim$ 11 au around HD 206893. Its mass was estimated between 12 and 50 $M_{Jup}$ from evolutionary models and its photometry. However, given the significant uncertainty on the age of the system and the peculiar spectrophotometric properties of the companion, this mass is not well constrained. We aim at constraining the orbit and dynamical mass of HD 206893 B. We combined radial velocity data obtained with HARPS spectra and astrometric data obtained with the high contrast imaging VLT/SPHERE and VLT/NaCo instruments, with a time baseline less than three years. We then combined those data with astrometry data obtained by Hipparcos and Gaia with a time baseline of 24 years. We used a MCMC approach to estimate the orbital parameters and dynamical mass of the brown dwarf from those data. We infer a period between 21 and 33{\deg} and an inclination in the range 20-41{\deg} from pole-on from HD 206893 B relative astrometry. The RV data show a significant RV drift over 1.6 yrs. We show that HD 206893 B cannot be the source of this observed RV drift as it would lead to a dynamical mass inconsistent with its photometry and spectra and with Hipparcos and Gaia data. An additional inner (semimajor axis in the range 1.4-2.6 au) and massive ($\sim$ 15 $M_{Jup}$) companion is needed to explain the RV drift, which is compatible with the available astrometric data of the star, as well as with the VLT/SPHERE and VLT/NaCo nondetection.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

Gaia data are revolutionizing our knowledge of the evolutionary history of the Milky Way. 3D maps of the interstellar dust provide complementary information and are a tool for a wide range of uses. We aimed at building 3D maps of the dust in the Local arm and surrounding regions. To do so, Gaia DR2 photometric data were combined with 2MASS measurements to derive extinction towards stars that possess accurate photometry and relative uncertainties on DR2 parallaxes smaller than 20\%. We applied to the extinctions a new hierarchical inversion algorithm adapted to large datasets and to a inhomogeneous target distribution. Each step associates regularized Bayesian inversions along radial directions and a subsequent inversion in 3D of their results. Each inverted distribution serves as a prior for the subsequent step and the spatial resolution is progressively increased. We present the resulting 3D distribution of the dust in a 6 x 6 x 0.8 kpc3 volume around the Sun. Its main features are found to be elongated along different directions that vary from below to above the mid-plane: the outer part of Carina-Sagittarius, mainly located above the mid-plane, the Local arm/Cygnus Rift around and above the mid-plane and the fragmented Perseus arm are oriented close to the direction of circular motion. The long spur (nicknamed the split) that extends between the Local Arm and Carina-Sagittarius, the compact near side of Carina-Sagittarius and the Cygnus Rift below the Plane are oriented along l=40 to 55 deg. Dust density images in vertical planes reveal in some regions a wavy pattern and show that the solar neighborhood within 500 pc remains atypical by its extent above and below the Plane. We show several comparisons with the locations of molecular clouds, HII regions, O stars and masers. The link between the dust concentration and these tracers is markedly different from one region to the other.

The assembly process of our Galaxy can be retrieved using the motions and chemistry of individual stars. Chemo-dynamical studies of the nearby halo have long hinted at the presence of multiple components such as streams, clumps, duality and correlations between the stars' chemical abundances and orbital parameters. More recently, the analysis of two large stellar surveys have revealed the presence of a well-populated chemical elemental abundance sequence, of two distinct sequences in the colour-magnitude diagram, and of a prominent slightly retrograde kinematic structure all in the nearby halo, which may trace an important accretion event experienced by the Galaxy. Here report an analysis of the kinematics, chemistry, age and spatial distribution of stars in a relatively large volume around the Sun that are mainly linked to two major Galactic components, the thick disk and the stellar halo. We demonstrate that the inner halo is dominated by debris from an object which at infall was slightly more massive than the Small Magellanic Cloud, and which we refer to as Gaia-Enceladus. The stars originating in Gaia-Enceladus cover nearly the full sky, their motions reveal the presence of streams and slightly retrograde and elongated trajectories. Hundreds of RR Lyrae stars and thirteen globular clusters following a consistent age-metallicity relation can be associated to Gaia-Enceladus on the basis of their orbits. With an estimated 4:1 mass-ratio, the merger with Gaia-Enceladus must have led to the dynamical heating of the precursor of the Galactic thick disk and therefore contributed to the formation of this component approximately 10 Gyr ago. These findings are in line with simulations of galaxy formation, which predict that the inner stellar halo should be dominated by debris from just a few massive progenitors.

The evolution of the Milky Way disk, which contains most of the stars in the Galaxy, is affected by several phenomena. For example, the bar and the spiral arms of the Milky Way induce radial migration of stars and can trap or scatter stars close to orbital resonances. External perturbations from satellite galaxies can also have a role, causing dynamical heating of the Galaxy, ring-like structures in the disk and correlations between different components of the stellar velocity. These perturbations can also cause 'phase wrapping' signatures in the disk, such as arched velocity structures in the motions of stars in the Galactic plane. Some manifestations of these dynamical processes have already been detected, including kinematic substructure in samples of nearby stars, density asymmetries and velocities across the Galactic disk that differ from the axisymmetric and equilibrium expectations, especially in the vertical direction, and signatures of incomplete phase mixing in the disk. Here we report an analysis of the motions of six million stars in the Milky Way disk. We show that the phase-space distribution contains different substructures with various morphologies, such as snail shells and ridges, when spatial and velocity coordinates are combined. We infer that the disk must have been perturbed between 300 million and 900 million years ago, consistent with estimates of the previous pericentric passage of the Sagittarius dwarf galaxy. Our findings show that the Galactic disk is dynamically young and that modelling it as time-independent and axisymmetric is incorrect.

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.

Using Gaia DR2 astrometry, we map the kinematic signature of the Galactic stellar warp out to a distance of 7 kpc from the Sun. Combining Gaia DR2 and 2MASS photometry, we identify, via a probabilistic approach, 599 494 upper main sequence stars and 12 616 068 giants without the need for individual extinction estimates. The spatial distribution of the upper main sequence stars clearly shows segments of the nearest spiral arms. The large-scale kinematics of both the upper main sequence and giant populations show a clear signature of the warp of the Milky Way, apparent as a gradient of 5-6 km/s in the vertical velocities from 8 to 14 kpc in Galactic radius. The presence of the signal in both samples, which have different typical ages, suggests that the warp is a gravitationally induced phenomenon.

Dwarf galaxies populating the Galactic halo are assumed to host the largest fractions of dark matter, as calculated from their velocity dispersions. Their major axes are preferentially aligned with the Vast Polar Structure (VPOS) that is perpendicular to the Galactic disk, and we find their velocity gradients aligned as well. It suggests that tidal forces exerted by the Milky Way are distorting dwarf galaxies. Here we demonstrate on the basis of the impulse approximation that the Galactic gravitational acceleration induces the dwarf line-of-sight velocity dispersion, which is also evidenced by strong dependences between both quantities. Since this result is valid for any dwarf mass value, it implies that dark matter estimate in Milky Way dwarfs cannot be deduced from the product of their radius to the square of their line-of-sight velocity dispersion. This questions the high dark-matter fractions reported for these evanescent systems, and the universally adopted total-to-stellar mass relationship in the dwarf regime. It suggests that many dwarfs are at their first passage and are dissolving into the Galactic halo. It opens a promising way to estimate the Milky Way total mass profile at large distances.

The second Gaia data release (GDR2) provides precise five-parameter astrometric data (positions, proper motions and parallaxes) for an unprecedented amount of sources (more than $1.3$ billion, mostly stars). The use of this wealth of astrometric data comes with a specific challenge: how does one properly infer from these data the astrophysical parameters of interest? The main - but not only - focus of this paper is the issue of the estimation of distances from parallaxes, possibly combined with other information. We start with a critical review of the methods traditionally used to obtain distances from parallaxes and their shortcomings. Then we provide guidelines on how to use parallaxes more efficiently to estimate distances by using Bayesian methods. In particular also we show that negative parallaxes, or parallaxes with relatively larger uncertainties still contain valuable information. Finally, we provide examples that show more generally how to use astrometric data for parameter estimation, including the combination of proper motions and parallaxes and the handling of covariances in the uncertainties. The paper contains examples based on simulated Gaia data to illustrate the problems and the solutions proposed. Furthermore, the developments and methods proposed in the paper are linked to a set of tutorials included in the Gaia archive documentation that provide practical examples and a good starting point for the application of the recommendations to actual problems. In all cases the source code for the analysis methods is provided. Our main recommendation is to always treat the derivation of (astro-) physical parameters from astrometric data, in particular when parallaxes are involved, as an inference problem which should preferably be handled with a full Bayesian approach.

We present the second Gaia data release, Gaia DR2, consisting of astrometry, photometry, radial velocities, and information on astrophysical parameters and variability, for sources brighter than magnitude 21. In addition epoch astrometry and photometry are provided for a modest sample of minor planets in the solar system. A summary of the contents of Gaia DR2 is presented, accompanied by a discussion on the differences with respect to Gaia DR1 and an overview of the main limitations which are still present in the survey. Recommendations are made on the responsible use of Gaia DR2 results. Gaia DR2 contains celestial positions and the apparent brightness in G for approximately 1.7 billion sources. For 1.3 billion of those sources, parallaxes and proper motions are in addition available. The sample of sources for which variability information is provided is expanded to 0.5 million stars. This data release contains four new elements: broad-band colour information in the form of the apparent brightness in the $G_\mathrm{BP}$ (330--680 nm) and $G_\mathrm{RP}$ (630--1050 nm) bands is available for 1.4 billion sources; median radial velocities for some 7 million sources are presented; for between 77 and 161 million sources estimates are provided of the stellar effective temperature, extinction, reddening, and radius and luminosity; and for a pre-selected list of 14000 minor planets in the solar system epoch astrometry and photometry are presented. Finally, Gaia DR2 also represents a new materialisation of the celestial reference frame in the optical, the Gaia-CRF2, which is the first optical reference frame based solely on extragalactic sources. There are notable changes in the photometric system and the catalogue source list with respect to Gaia DR1, and we stress the need to consider the two data releases as independent.

The Gaia Data Release 2 contains the 1st release of radial velocities complementing the kinematic data of a sample of about 7 million relatively bright, late-type stars. Aims: This paper provides a detailed description of the Gaia spectroscopic data processing pipeline, and of the approach adopted to derive the radial velocities presented in DR2. Methods: The pipeline must perform four main tasks: (i) clean and reduce the spectra observed with the Radial Velocity Spectrometer (RVS); (ii) calibrate the RVS instrument, including wavelength, straylight, line-spread function, bias non-uniformity, and photometric zeropoint; (iii) extract the radial velocities; and (iv) verify the accuracy and precision of the results. The radial velocity of a star is obtained through a fit of the RVS spectrum relative to an appropriate synthetic template spectrum. An additional task of the spectroscopic pipeline was to provide 1st-order estimates of the stellar atmospheric parameters required to select such template spectra. We describe the pipeline features and present the detailed calibration algorithms and software solutions we used to produce the radial velocities published in DR2. Results: The spectroscopic processing pipeline produced median radial velocities for Gaia stars with narrow-band near-IR magnitude Grvs < 12 (i.e. brighter than V~13). Stars identified as double-lined spectroscopic binaries were removed from the pipeline, while variable stars, single-lined, and non-detected double-lined spectroscopic binaries were treated as single stars. The scatter in radial velocity among different observations of a same star, also published in DR2, provides information about radial velocity variability. For the hottest (Teff > 7000 K) and coolest (Teff < 3500 K) stars, the accuracy and precision of the stellar parameter estimates are not sufficient to allow selection of appropriate templates. [Abridged]

To illustrate the potential of GDR2, we provide a first look at the kinematics of the Milky Way disc, within a radius of several kiloparsecs around the Sun. We benefit for the first time from a sample of 6.4 million F-G-K stars with full 6D phase-space coordinates, precise parallaxes, and precise Galactic cylindrical velocities . From this sample, we extracted a sub-sample of 3.2 million giant stars to map the velocity field of the Galactic disc from $\sim$5~kpc to $\sim$13~kpc from the Galactic centre and up to 2~kpc above and below the plane. We also study the distribution of 0.3 million solar neighbourhood stars ($r < 200$~pc), with median velocity uncertainties of 0.4~km/s, in velocity space and use the full sample to examine how the over-densities evolve in more distant regions. GDR2 allows us to draw 3D maps of the Galactocentric median velocities and velocity dispersions with unprecedented accuracy, precision, and spatial resolution. The maps show the complexity and richness of the velocity field of the galactic disc. We observe streaming motions in all the components of the velocities as well as patterns in the velocity dispersions. For example, we confirm the previously reported negative and positive galactocentric radial velocity gradients in the inner and outer disc, respectively. Here, we see them as part of a non-axisymmetric kinematic oscillation, and we map its azimuthal and vertical behaviour. We also witness a new global arrangement of stars in the velocity plane of the solar neighbourhood and in distant regions in which stars are organised in thin substructures with the shape of circular arches that are oriented approximately along the horizontal direction in the $U-V$ plane. Moreover, in distant regions, we see variations in the velocity substructures more clearly than ever before, in particular, variations in the velocity of the Hercules stream. (abridged)

The second Gaia data release (DR2), contains very precise astrometric and photometric properties for more than one billion sources, astrophysical parameters for dozens of millions, radial velocities for millions, variability information for half a million of stellar sources and orbits for thousands of solar system objects. Before the Catalogue publication, these data have undergone dedicated validation processes. The goal of this paper is to describe the validation results in terms of completeness, accuracy and precision of the various Gaia DR2 data. The validation processes include a systematic analysis of the Catalogue content to detect anomalies, either individual errors or statistical properties, using statistical analysis, and comparisons to external data or to models. Although the astrometric, photometric and spectroscopic data are of unprecedented quality and quantity, it is shown that the data cannot be used without a dedicated attention to the limitations described here, in the Catalogue documentation and in accompanying papers. A particular emphasis is put on the caveats for the statistical use of the data in scientific exploitation.

For Gaia DR2 (GDR2), 280 million spectra, collected by the RVS instrument on-board Gaia, were processed and median radial velocities were derived for 9.8 million sources brighter than Grvs = 12 mag. This paper describes the validation and properties of the median radial velocities published in GDR2. Quality tests and filters are applied to select, from the 9.8 million radial velocities, those with the quality to be published in GDR2. The accuracy of the selected sample is assessed with respect to ground-based catalogues. Its precision is estimated using both ground-based catalogues and the distribution of the Gaia radial velocity uncertainties. GDR2 contains median radial velocities for 7 224 631 stars, with Teff in the range [3550, 6900] K, which passed succesfully the quality tests. The published median radial velocities provide a full sky-coverage and have a completness with respect to the astrometric data of 77.2\% (for $G \leq 12.5$ mag). The median radial velocity residuals with respect to the ground-based surveys vary from one catalogue to another, but do not exceed a few 100s m/s. In addition, the Gaia radial velocities show a positive trend as a function of magnitude, which starts around Grvs $\sim 9$ mag and reaches about $+500$ m/s at Grvs $= 11.75$ mag. The overall precision, estimated from the median of the Gaia radial velocity uncertainties, is 1.05 km/s. The radial velocity precision is function of many parameters, in particular the magnitude and effective temperature. For bright stars, Grvs in [4, 8] mag, the precision is in the range 200-350 m/s, which is about 3 to 5 times more precise than the pre-launch specification of 1 km/s. At the faint end, Grvs = 11.75 mag, the precisions for Teff = 5000 K and 6500 K are respectively 1.4 km/s and 3.7 km/s.

Gaia data and stellar surveys open the way to the construction of detailed 3D maps of the Galactic interstellar (IS) dust based on the synthesis of star distances and extinctions. Reliable extinction measurements require very accurate photometric calibrations. We show the first step of an iterative process linking 3D dust maps and photometric calibrations and improving them simultaneously. Our previous 3D map of nearby IS dust was used to select low reddening SDSS/APOGEE-DR14 red giants, and this database served for an empirical effective temperature- and metallicity-dependent photometric calibration in the Gaia G and 2MASS Ks bands. This calibration has been combined with Gaia G-band empirical extinction coefficients recently published, G, J and Ks photometry and APOGEE atmospheric parameters to derive the extinction of a large fraction of the survey targets. Distances were estimated independently using isochrones and the magnitude-independent extinction K(J-Ks). This new dataset has been merged with the one used for the earlier version of dust map. A new Bayesian inversion of distance-extinction pairs has been performed to produce an updated 3D map. We present several properties of the new map. Its comparison with 2D dust emission reveals that all large dust shells seen in emission at mid- and high-latitude are closer than 300pc. The updated distribution constrains the well debated, X-ray bright North Polar Spur to originate beyond 800 pc. We use the Orion region to illustrate additional details and distant clouds. On the large scale the map reveals a complex structure of the Local Arm. 2 to 3 kpc-long chains of clouds appear in planes tilted by 15 deg with respect to the Galactic plane. A series of cavities oriented along a l=60-240deg axis crosses the Arm. (http://stilism.obspm.fr)

The first Gaia data release unlocked the access to the photometric information of 1.1 billion sources in the $G$-band. Yet, given the high level of degeneracy between extinction and spectral energy distribution for large passbands such as the Gaia $G$-band, a correction for the interstellar reddening is needed in order to exploit Gaia data. The purpose of this manuscript is to provide the empirical estimation of the Gaia $G$-band extinction coefficient $k_G$ for both the red giants and main sequence stars, in order to be able to exploit the first data release DR1. We selected two samples of single stars: one for the red giants and one for the main sequence. Both samples are the result of a cross-match between Gaia DR1 and 2MASS catalogues; they consist of high quality photometry in the $G$-, $J$- and $Ks$-bands. These samples were complemented by temperature and metallicity information retrieved from, respectively, APOGEE DR13 and LAMOST DR2 surveys. We implemented a Markov Chain Monte Carlo method where we used $(G-Ks)_0$ vs $T_\mathrm{eff}$ and $(J-Ks)_0$ vs $(G-Ks)_0$ calibration relations to estimate the extinction coefficient $k_G$ and we quantify its corresponding confidence interval via bootstrap resampling method. We tested our method on samples of red giants and main sequence stars, finding consistent solutions. We present here the determination of the Gaia extinction coefficient through a completely empirical method. Furthermore we provide the scientific community a formula for measuring the extinction coefficient as a function of stellar effective temperature, the intrinsic colour $(G-Ks)_0$ and absorption.

Gaia Data Release 1 allows to recalibrate standard candles such as the Red Clump stars. To use those, they first need to be accurately characterised. In particular, colours are needed to derive the interstellar extinction. As no filter is available for the first Gaia data release and to avoid the atmosphere model mismatch, an empirical calibration is unavoidable. The purpose of this work is to provide the first complete and robust photometric empirical calibration of the Gaia Red Clump stars of the solar neighbourhood, through colour-colour, effective temperature-colour and absolute magnitude-colour relations, from the Gaia, Johnson, 2MASS, Hipparcos, Tycho-2, APASS-SLOAN and WISE photometric systems, and the APOGEE DR13 spectroscopic temperatures. We used a 3D extinction map to select low reddening red giants. To calibrate the colour-colour and the effective temperature-colour relations, we developed a MCMC method which accounts for all variable uncertainties and selects the best model for each photometric relation. We estimate the Red Clump absolute magnitude through the mode of a kernel-based distribution function. We provide 20 colour vs G-Ks relations and the first Teff vs G-Ks calibration. We obtained the Red Clump absolute magnitudes for 15 photometric bands with, in particular, M_Ks = -1.606 +/- 0.009 and M_G = (0.495 +/- 0.009) + (1.121 +/- 0.128) (G-Ks-2.1). We present an unreddened Gaia-TGAS HR diagram and use the calibrations to compare its Red Clump and its Red Giant Branch Bump with the Padova isochrones.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

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.

Abridged: We used the fourth internal data release of the Gaia-ESO survey to characterize the bulge chemistry, spatial distribution, kinematics, and to compare it chemically with the thin and thick disks. The sample consist on ~2500 red clump stars in 11 bulge fields ($-10^\circ\leq l\leq+8^\circ$ and $-10^\circ\leq b\leq-4^\circ$), and a set of ~6300 disk stars selected for comparison. The bulge MDF is confirmed to be bimodal across the whole sampled area, with metal-poor stars dominating at high latitudes. The metal-rich stars exhibit bar-like kinematics and display a bimodality in their magnitude distribution, a feature which is tightly associated with the X-shape bulge. They overlap with the metal-rich end of the thin disk sequence in the [Mg/Fe] vs. [Fe/H] plane. Metal-poor bulge stars have a more isotropic hot kinematics and do not participate in the X-shape bulge. With similar Mg-enhancement levels, the position of the metal-poor bulge sequence "knee" is observed at [Fe/H]$_{knee}=-0.37\pm0.09$, being 0.06 dex higher than that of the thick disk. It suggests a higher SFR for the bulge than for the thick disk. Finally, we present a chemical evolution model that suitably fits the whole bulge sequence by assuming a fast ($<1$ Gyr) intense burst of stellar formation at early epochs. We associate metal-rich stars with the B/P bulge formed from the secular evolution of the early thin disk. On the other hand, the metal-poor subpopulation might be the product of an early prompt dissipative collapse dominated by massive stars. Nevertheless, our results do not allow us to firmly rule out the possibility that these stars come from the secular evolution of the early thick disk. This is the first time that an analysis of the bulge MDF and $\alpha$-abundances has been performed in a large area on the basis of a homogeneous, fully spectroscopic analysis of high-resolution, high S/N data.

Due to their volatile nature, when sulfur and zinc are observed in external galaxies, their determined abundances represent the gas-phase abundances in the interstellar medium. This implies that they can be used as tracers of the chemical enrichment of matter in the Universe at high redshift. Comparable observations in stars are more difficult and, until recently, plagued by small number statistics. We wish to exploit the Gaia ESO Survey (GES) data to study the behaviour of sulfur and zinc abundances of a large number of Galactic stars, in a homogeneous way. By using the UVES spectra of the GES sample, we are able to assemble a sample of 1301 Galactic stars, including stars in open and globular clusters in which both sulfur and zinc were measured. We confirm the results from the literature that sulfur behaves as an alpha-element. We find a large scatter in [Zn/Fe] ratios among giant stars around solar metallicity. The lower ratios are observed in giant stars at Galactocentric distances less than 7.5 kpc. No such effect is observed among dwarf stars, since they do not extend to that radius. Given the sample selection, giants and dwarfs are observed at different Galactic locations, and it is plausible, and compatible with simple calculations, that Zn-poor giants trace a younger population more polluted by SN Ia yields. It is necessary to extend observations in order to observe both giants and dwarfs at the same Galactic location. Further theoretical work on the evolution of zinc is also necessary.

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.

We take advantage of the Gaia-ESO Survey iDR4 bulge data to search for abundance anomalies that could shed light on the composite nature of the Milky Way bulge. The alpha-elements (Mg, Si, and whenever available, Ca) abundances, and their trends with Fe abundances have been analysed for a total of 776 bulge stars. In addition, the aluminum abundances and their ratio to Fe and Mg have also been examined. Our analysis reveals the existence of low-alpha element abundance stars with respect to the standard bulge sequence in the [alpha/Fe] vs. [Fe/H] plane. 18 objects present deviations in [alpha/Fe] ranging from 2.1 to 5.3 sigma with respect to the median standard value. Those stars do not show Mg-Al anti-correlation patterns. Incidentally, this sign of the existence of multiple stellar populations is reported firmly for the first time for the bulge globular cluster NGC 6522. The identified low-alpha abundance stars have chemical patterns compatible with those of the thin disc. Their link with massive dwarf galaxies accretion seems unlikely, as larger deviations in alpha abundance and Al would be expected. The vision of a bulge composite nature and a complex formation process is reinforced by our results. The used approach, a multi-method and model-driven analysis of high resolution data seems crucial to reveal this complexity.

Before the publication of the Gaia Catalogue, the contents of the first data release have undergone multiple dedicated validation tests. These tests aim at analysing in-depth the Catalogue content to detect anomalies, individual problems in specific objects or in overall statistical properties, either to filter them before the public release, or to describe the different caveats of the release for an optimal exploitation of the data. Dedicated methods using either Gaia internal data, external catalogues or models have been developed for the validation processes. They are testing normal stars as well as various populations like open or globular clusters, double stars, variable stars, quasars. Properties of coverage, accuracy and precision of the data are provided by the numerous tests presented here and jointly analysed to assess the data release content. This independent validation confirms the quality of the published data, Gaia DR1 being the most precise all-sky astrometric and photometric catalogue to-date. However, several limitations in terms of completeness, astrometric and photometric quality are identified and described. Figures describing the relevant properties of the release are shown and the testing activities carried out validating the user interfaces are also described. A particular emphasis is made on the statistical use of the data in scientific exploitation.

Several recent studies have demonstrated that the Galactic bulge hosts two components with different mean metallicities, and possibly different spatial distribution and kinematics. As a consequence, both the metallicity distribution and the radial velocity of bulge stars vary across different line of sights. We present here the metallicity distribution function of red clump stars in 26 fields spread across a wide area of the bulge, with special emphasis on fields close to Galactic plane, at latitudes b=-2 and b=-1, that were not explored before. This paper includes new metallicities from a sample of about 5000 K giant stars, observed at spectral resolution R=6500, in the Calcium II Triplet region. They are the main dataset of the GIRAFFE Inner Bulge Survey. As part of the same survey we have previously published results for a sample of about 600 K giant stars, at latitude b=-4 , derived from higher resolution spectra (R=22,500). Results. The combined sample allows us to trace and characterize the metal poor and metal rich bulge populations down to the inner bulge. We present a density map for each of the two components. Contrary to the expectations from previous works, we found the metal poor population to be more centrally concentrated than the metal rich one, and with a more axisymmetric spatial distribution. The metal rich population, on the other hand, is arranged in a boxy distribution, consistent with an edge-on bar. By coupling metallicities and radial velocities we show that the metal poor population has a velocity dispersion that varies rather mildly with latitude. On the contrary, the metal rich population has a low velocity dispersion far from the plane (b=-8.5), but it has a steeper gradient with latitude, becoming higher than the metal poor one in the innermost field (b=-1). [abridged]

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

This short document is intended as a companion and introduction to the Detailed Science Case (DSC) for the Maunakea Spectroscopic Explorer. It provides a concise summary of the essential characteristics of MSE from the perspective of the international astronomical community. MSE is a wide field telescope (1.5 square degree field of view) with an aperture of 11.25m. It is dedicated to multi-object spectroscopy at several different spectral resolutions in the range R ~ 2500 - 40000 over a broad wavelength range (0.36 - 1.8{\mu}m). MSE will enable transformational science in areas as diverse as exoplanetary host characterization; stellar monitoring campaigns; tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of the distant Galaxy; connecting galaxies to the large scale structure of the Universe; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars. MSE is the largest ground based optical and near infrared telescope in its class, and it will occupy a unique and critical role in the emerging network of astronomical facilities active in the 2020s. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, eROSITA, SKA, and WFIRST, and is an ideal feeder facility for E-ELT, TMT and GMT.

(Abridged) By using an N-body simulation of a bulge that was formed via a bar instability mechanism, we analyse the imprints of the initial (i.e. before bar formation) location of stars on the bulge kinematics, in particular on the heliocentric radial velocity distribution of bulge stars. Four different latitudes were considered: $b=-4^\circ$, $-6^\circ$, $-8^\circ$, and $-10^\circ$, along the bulge minor axis as well as outside it, at $l=\pm5^\circ$ and $l=\pm10^\circ$. The bulge X-shaped structure comprises stars that formed in the disk at different locations. Stars formed in the outer disk, beyond the end of the bar, which are part of the boxy peanut-bulge structure may show peaks in the velocity distributions at positive and negative heliocentric radial velocities with high absolute values that can be larger than 100 $\rm km$ $\rm s^{-1}$, depending on the observed direction. In some cases the structure of the velocity field is more complex and several peaks are observed. Stars formed in the inner disk, the most numerous, contribute predominantly to the X-shaped structure and present different kinematic characteristics. Our results may enable us to interpret the cold high-velocity peak observed in the APOGEE commissioning data, as well as the excess of high-velocity stars in the near and far arms of the X-shaped structure at $l$=$0^\circ$ and $b$=$-6^\circ$. When compared with real data, the kinematic picture becomes more complex due to the possible presence in the observed samples of classical bulge and/or thick disk stars. Overall, our results point to the existence of complex patterns and structures in the bulge velocity fields, which are generated by the bar. This suggests that caution should be used when interpreting the bulge kinematics: the presence of substructures, peaks and clumps in the velocity fields is not necessarily a sign of past accretion events.

Recent large scale surveys of galactic bulge stars allowed to build a detailed map of the bulge kinematics. The bulge exhibits cylindrical rotation consistent with a disky origin which evolved through bar driven secular evolution. However correlations between metallicity and kinematics complicate this picture. In particular a metal-poor component with distinct kinematic signatures has been detected. Its origin, density profile and link with the other Milky Way stellar populations are currently still poorly constrained.

High resolution (R$\sim$22,500) spectra for 400 red clump giants, in four fields within $\rm -4.8^{\circ} \lesssim b \lesssim -3.4^{\circ}$ and $\rm -10^{\circ} \lesssim l \lesssim +10^{\circ}$, were obtained within the GIRAFFE Inner Bulge Survey (GIBS) project. To this sample we added another $\sim$ 400 stars in Baade's Window, observed with the identical instrumental configuration. We constructed the metallicity distributions for the entire sample, as well as for each field individually, in order to investigate the presence of gradients or field-to-field variations in the shape of the distributions. The metallicity distributions in the five fields are consistent with being drawn from a single parent population, indicating the absence of a gradient along the major axis of the Galactic bar. The global metallicity distribution is well fitted by two Gaussians. The metal poor component is rather broad, with a mean at $\rm <[Fe/H]>=-0.31$ dex and $\sigma=0.31$ dex. The metal-rich one is narrower, with mean $\rm <[Fe/H]>=+0.26$ and $\sigma=0.2$ dex. The [Mg/Fe] ratio follows a tight trend with [Fe/H], with enhancement with respect to solar in the metal-poor regime, similar to the one observed for giant stars in the local thick disc. [Ca/Fe] abundances follow a similar trend, but with a considerably larger scatter than [Mg/Fe]. A decrease in [Mg/Fe] is observed at $\rm [Fe/H]=-0.44$ dex. This \textit{knee} is in agreement with our previous bulge study of K-giants along the minor axis, but is 0.1 dex lower in metallicity than the one reported for the Bulge micro lensed dwarf and sub-giant stars. We found no variation in $\alpha$-element abundance distributions between different fields.

The Gaia-ESO Survey (GES) is a large public spectroscopic survey at the European Southern Observatory Very Large Telescope. A key aim is to provide precise radial velocities (RVs) and projected equatorial velocities (v sin i) for representative samples of Galactic stars, that will complement information obtained by the Gaia astrometry satellite. We present an analysis to empirically quantify the size and distribution of uncertainties in RV and v sin i using spectra from repeated exposures of the same stars. We show that the uncertainties vary as simple scaling functions of signal-to-noise ratio (S/N) and v sin i, that the uncertainties become larger with increasing photospheric temperature, but that the dependence on stellar gravity, metallicity and age is weak. The underlying uncertainty distributions have extended tails that are better represented by Student's t-distributions than by normal distributions. Parametrised results are provided, that enable estimates of the RV precision for almost all GES measurements, and estimates of the v sin i precision for stars in young clusters, as a function of S/N, v sin i and stellar temperature. The precision of individual high S/N GES RV measurements is 0.22-0.26 km/s, dependent on instrumental configuration.

This paper describes the analysis of UVES and GIRAFFE spectra acquired by the Gaia-ESO Public Spectroscopic Survey in the fields of young clusters whose population includes pre-main sequence (PMS) stars. Both methods that have been extensively used in the past and new ones developed in the contest of the Gaia-ESO survey enterprise are available and used. The internal precision of these quantities is estimated by inter-comparing the results obtained by such different methods, while the accuracy is estimated by comparison with independent external data, like effective temperature and surface gravity derived from angular diameter measurements, on a sample of benchmarks stars. Specific strategies are implemented to deal with fast rotation, accretion signatures, chromospheric activity, and veiling. The analysis carried out on spectra acquired in young clusters' fields during the first 18 months of observations, up to June 2013, is presented in preparation of the first release of advanced data products. Stellar parameters obtained with the higher resolution and larger wavelength coverage from UVES are reproduced with comparable accuracy and precision using the smaller wavelength range and lower resolution of the GIRAFFE setup adopted for young stars, which allows us to provide with confidence stellar parameters for the much larger GIRAFFE sample. Precisions are estimated to be $\approx$ 120 K r.m.s. in Teff, $\approx$0.3 dex r.m.s. in logg, and $\approx$0.15 dex r.m.s. in [Fe/H], for both the UVES and GIRAFFE setups.

We study how diffuse interstellar bands (DIBs) measured toward distance-distributed target stars can be used to locate dense interstellar (IS) clouds in the Galaxy and probe a line-of-sight (LOS) kinematical structure, a potential useful tool when gaseous absorption lines are saturated or not available in the spectral range. Cool target stars are numerous enough for this purpose. We have devised automated DIB fitting methods appropriate to cool star spectra and multiple IS components. The data is fitted with a combination of a synthetic stellar spectrum, a synthetic telluric transmission, and empirical DIB profiles. In parallel, stellar distances and extinctions are estimated self-consistently by means of a 2D Bayesian method, from spectroscopically-derived stellar parameters and photometric data. We have analyzed Gaia-ESO Survey (GES) and previously recorded spectra that probe between $\sim$ 2 and 10 kpc long LOS in five different regions of the Milky Way. Depending on the observed spectral intervals, we extracted one or more of the following DIBs: $\lambda\lambda$ 6283.8, 6613.6 and 8620.4. For each field, we compared the DIB strengths with the Bayesian distances and extinctions, and the DIB Doppler velocities with the HI emission spectra. For all fields, the DIB strength and the target extinction are well correlated. In case of targets widely distributed in distance, marked steps in DIBs and extinction radial distance profiles match with each other and broadly correspond to the expected locations of spiral arms. For all fields, the DIB velocity structure agrees with HI emission spectra and all detected DIBs correspond to strong NaI lines. This illustrates how DIBs can be used to locate the Galactic interstellar gas and to study its kinematics at the kpc scale.

The Gaia-ESO Survey is obtaining high-quality spectroscopic data for about 10^5 stars using FLAMES at the VLT. UVES high-resolution spectra are being collected for about 5000 FGK-type stars. These UVES spectra are analyzed in parallel by several state-of-the-art methodologies. Our aim is to present how these analyses were implemented, to discuss their results, and to describe how a final recommended parameter scale is defined. We also discuss the precision (method-to-method dispersion) and accuracy (biases with respect to the reference values) of the final parameters. These results are part of the Gaia-ESO 2nd internal release and will be part of its 1st public release of advanced data products. The final parameter scale is tied to the one defined by the Gaia benchmark stars, a set of stars with fundamental atmospheric parameters. A set of open and globular clusters is used to evaluate the physical soundness of the results. Each methodology is judged against the benchmark stars to define weights in three different regions of the parameter space. The final recommended results are the weighted-medians of those from the individual methods. The recommended results successfully reproduce the benchmark stars atmospheric parameters and the expected Teff-log g relation of the calibrating clusters. Atmospheric parameters and abundances have been determined for 1301 FGK-type stars observed with UVES. The median of the method-to-method dispersion of the atmospheric parameters is 55 K for Teff, 0.13 dex for log g, and 0.07 dex for [Fe/H]. Systematic biases are estimated to be between 50-100 K for Teff, 0.10-0.25 dex for log g, and 0.05-0.10 dex for [Fe/H]. Abundances for 24 elements were derived: C, N, O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ba, Nd, and Eu. The typical method-to-method dispersion of the abundances varies between 0.10 and 0.20 dex.

(Abridged) We analyzed the stellar parameters and radial velocities of ~1200 stars in five bulge fields as determined from the Gaia-ESO survey data (iDR1). We use VISTA Variables in The Via Lactea (VVV) photometry to obtain reddening values by using a semi-empirical T_eff-color calibration. From a Gaussian decomposition of the metallicity distribution functions, we unveil a clear bimodality in all fields, with the relative size of components depending of the specific position on the sky. In agreement with some previous studies, we find a mild gradient along the minor axis (-0.05 dex/deg between b=-6 and b=-10) that arises from the varying proportion of metal-rich and metal-poor components. The number of metal-rich stars fades in favor of the metal-poor stars with increasing b. The K-magnitude distribution of the metal-rich population splits into two peaks for two of the analyzed fields that intersects the near and far branches of the X-shaped bulge structure. In addition, two lateral fields at (l,b)=(7,-9) and (l,b)=(-10,-8) present contrasting characteristics. In the former, the metallicity distribution is dominated by metal-rich stars, while in the latter it presents a mix of a metal-poor population and and a metal-intermediate one, of nearly equal sizes. Finally, we find systematic differences in the velocity dispersion between the metal-rich and the metal-poor components of each field. Our chemo-kinematical analysis is consistent with a varying field-to-field proportion of stars belonging to (i) a metal-rich boxy/peanut X-shaped component, with bar-like kinematics, and (ii) a metal-poor more extended rotating structure with a higher velocity dispersion that dominates far from the Galactic plane. These first GES data allow studying the detailed spatial dependence of the Galactic bulge populations, thanks to the analysis of individual fields with relatively high statistics.

Aims: An effort has been undertaken to simulate the expected Gaia Catalogue, including the effect of observational errors. A statistical analysis of this simulated Gaia data is performed in order to better understand what can be obtained from the Gaia astrometric mission. This catalogue is used in order to investigate the potential yield in astrometric, photometric and spectroscopic information, and the extent and effect of observational errors on the true Gaia Catalogue. This article is a follow-up to Robin et. al. (2012), where the expected Gaia Catalogue content was reviewed but without the simulation of observational errors. Methods: The Gaia Object Generator (GOG) catalogue is analysed using the Gaia Analysis Tool (GAT), producing a number of statistics on the catalogue. Results: A simulated catalogue of one billion objects is presented, with detailed information on the 523 million individual single stars it contains. Detailed information is provided for the expected errors in parallax, position, proper motion, radial velocity, photometry in the four Gaia bands, and physical parameter determination including temperature, metallicity and line of sight extinction.

The Gaia-ESO Survey is a large public spectroscopic survey that aims to derive radial velocities and fundamental parameters of about 10^5 Milky Way stars in the field and in clusters. Observations are carried out with the multi-object optical spectrograph FLAMES, using simultaneously the medium resolution (R~20,000) GIRAFFE spectrograph and the high resolution (R~47,000) UVES spectrograph. In this paper, we describe the methods and the software used for the data reduction, the derivation of the radial velocities, and the quality control of the FLAMES-UVES spectra. Data reduction has been performed using a workflow specifically developed for this project. This workflow runs the ESO public pipeline optimizing the data reduction for the Gaia-ESO Survey, performs automatically sky subtraction, barycentric correction and normalisation, and calculates radial velocities and a first guess of the rotational velocities. The quality control is performed using the output parameters from the ESO pipeline, by a visual inspection of the spectra and by the analysis of the signal-to-noise ratio of the spectra. Using the observations of the first 18 months, specifically targets observed multiple times at different epochs, stars observed with both GIRAFFE and UVES, and observations of radial velocity standards, we estimated the precision and the accuracy of the radial velocities. The statistical error on the radial velocities is sigma~0.4 km s^-1 and is mainly due to uncertainties in the zero point of the wavelength calibration. However, we found a systematic bias with respect to the GIRAFFE spectra (~0.9 km s^-1) and to the radial velocities of the standard stars (~0.5 kms^-1) retrieved from the literature. This bias will be corrected in the future data releases, when a common zero point for all the setups and instruments used for the survey will be established.

The Galactic bulge is a massive, old component of the Milky Way. It is known to host a bar, and it has recently been demonstrated to have a pronounced boxy/peanut structure in its outer region. Several independent studies suggest the presence of more than one stellar populations in the bulge, with different origins and a relative fraction changing across the bulge area. This is the first of a series of papers presenting the results of the Giraffe Inner Bulge Survey, carried out at the ESO-VLT with the multifibre spectrograph FLAMES. Spectra of ~5000 red clump giants in 24 bulge fields have been obtained at resolution R=6500, in the infrared Calcium triplet wavelength region at 8500 {\AA}. They are used to derive radial velocities and metallicities, based on new calibration specifically devised for this project. Radial velocities for another ~1200 bulge red clump giants, obtained from similar archive data, have been added to the sample. Higher resolution spectra have been obtained for 450 additional stars at latitude b=-3.5, with the aim of investigating chemical abundance patterns variations with longitude, across the inner bulge. In total we present here radial velocities for 6392 RC stars. We derive a radial velocity, and velocity dispersion map of the Milky Way bulge, useful to be compared with similar maps of external bulges, and to infer the expected velocities and dispersion at any line of sight. The K-type giants kinematics is consistent with the cylindrical rotation pattern of M-giants from the BRAVA survey. Our sample enables to extend this result to latitude b=-2, closer to the Galactic plane than probed by previous surveys. Finally, we find strong evidence for a velocity dispersion peak at (0,-1) and (0,-2), possibly indicative of a high density peak in the central 250 pc of the bulge

Constraints on the Galactic bulge/bar structure and formation history from stellar kinematics and metallicities mainly come from relatively high-latitude fields (|b|>4) where a complex mix of stellar population is seen. We aim here to constrain the formation history of the Galactic bar by studying the radial velocity and metallicity distributions of stars in-situ (|b|<1). We observed red clump stars in four fields along the bar's major axis (l=10,6,-6 and b=0 plus a field at l=0,b=1) with low-resolution spectroscopy from VLT/FLAMES, observing around the CaII triplet. We developed robust methods for extracting radial velocity and metallicity estimates from these low signal-to-noise spectra. We derived distance probability distributions using Bayesian methods rigorously handling the extinction law. We present radial velocities and metallicity distributions, as well as radial velocity trends with distance. We observe an increase in the radial velocity dispersion near the Galactic plane. We detect the streaming motion of the stars induced by the bar in fields at l=+/-6, the highest velocity components of this bar stream being metal-rich ([Fe/H]~0.2 dex). Our data is consistent with a bar inclined at 26+/-3 from the Sun-Galactic centre line. We observe a significant fraction of metal-poor stars, in particular in the field at l=0,b=1. We confirm the flattening of the metallicity gradient along the minor axis when getting closer to the plane, with a hint that it could actually be inverted. Our stellar kinematics corresponds to the expected behaviour of a bar issued from the secular evolution of the Galactic disc. The mix of several populations, seen further away from the plane, is also seen in the bar in-situ since our metallicity distributions highlight a different spatial distribution between metal-poor and metal-rich stars, the more metal-poor stars being more centrally concentrated.

Interstellar diffuse bands are usually extracted from hot star spectra because they are characterized by smooth continua. It introduces a strong limitation on the number of available targets, and reduces potential studies of the IS matter and the use of absorptions for cloud mapping. We have developed a new automatic fitting method appropriate to interstellar absorptions in spectra of cool stars that possess stellar atmospheric parameters. We applied this method to the extraction of three DIBs in high resolution VLT FLAMES/GIRAFFE spectra of red clump stars from the bulge. By combining all stellar synthetic spectra, HITRAN-LBLRTM atmospheric transmission spectra and diffuse band empirical absorption profiles, we determine the 6196, 6204, and 6284 A DIB strength toward the 219 target stars and discuss the sources of uncertainties. In order to test the sensitivity of the DIB extraction, we intercompare the three results and compare the DIB equivalent widths with the reddening derived from an independent extinction map based on OGLE photometric data. Most stellar spectra could be well reproduced by the composite stellar, atmospheric and interstellar models. Measurement uncertainties on the EWs are smaller for the broad and strong 6284 A DIB, and are of the order of 10-15%. Uncertainties on the two narrow and weaker DIBs are larger, as expected, and found to be highly variable from one target to the other. They strongly depend on the radial velocity of the star . DIB-DIB correlations among the three bands demonstrate that a meaningful signal is extracted. For the 6284 and 6204 A DIBs, the star-to-star variability of the equivalent width (EW) also reflects features of the OGLE extinction map...

Until recently our knowledge of the Galactic Bulge stellar populations was based on the study of a few low extinction windows. Large photometric and spectroscopic surveys are now underway to map large areas of the bulge. They probe several complex structures which are still to be fully characterized as well as their links with the inner disc, the thick disc and the inner halo. I will review our current, rapidly increasing, knowledge of the bulge stellar populations and the new insight expected towards the Gaia era to disentangle the formation history of the Galactic inner regions.

MOONS is a new conceptual design for a Multi-Object Optical and Near-infrared Spectrograph for the Very Large Telescope (VLT), selected by ESO for a Phase A study. The baseline design consists of 1000 fibers deployable over a field of view of 500 square arcmin, the largest patrol field offered by the Nasmyth focus at the VLT. The total wavelength coverage is 0.8um-1.8um and two resolution modes: medium resolution and high resolution. In the medium resolution mode (R=4,000-6,000) the entire wavelength range 0.8um-1.8um is observed simultaneously, while the high resolution mode covers simultaneously three selected spectral regions: one around the CaII triplet (at R=8,000) to measure radial velocities, and two regions at R=20,000 one in the J-band and one in the H-band, for detailed measurements of chemical abundances. The grasp of the 8.2m Very Large Telescope (VLT) combined with the large multiplex and wavelength coverage of MOONS - extending into the near-IR - will provide the observational power necessary to study galaxy formation and evolution over the entire history of the Universe, from our Milky Way, through the redshift desert and up to the epoch of re-ionization at z>8-9. At the same time, the high spectral resolution mode will allow astronomers to study chemical abundances of stars in our Galaxy, in particular in the highly obscured regions of the Bulge, and provide the necessary follow-up of the Gaia mission. Such characteristics and versatility make MOONS the long-awaited workhorse near-IR MOS for the VLT, which will perfectly complement optical spectroscopy performed by FLAMES and VIMOS.

We present a new database of 5780.5 and 6283.8 {\AA} DIB measurements and also study their correlation with the reddening. The database is based on high-resolution, high-quality spectra of early-type nearby stars located in the southern hemisphere at an average distance of 300 pc. Equivalent widths of the two DIBs were determined by means of a realistic continuum fitting and synthetic atmospheric transmissions. For all stars that possess a precise measurement of their color excess, we compare the DIBs and the extinction. We find average linear relationships of the DIBS and the color excess that agree well with those of a previous survey of northern hemisphere stars closer than 550 pc. This similarity shows that there is no significant spatial dependence of the average relationship in the solar neighborhood within $\simeq$ 600 pc. A noticeably different result is our higher degree of correlation of the two DIBs with the extinction. We demonstrate that it is simply due to the lower temperature and intrinsic luminosity of our targets. Using cooler target stars reduces the number of outliers, especially for nearby stars, confirming that the radiation field of UV bright stars has a significant influence on the DIB strength. We have used the cleanest data to compute updated DIB shapes.

Context. This study has been developed in the framework of the computational simulations executed for the preparation of the ESA Gaia astrometric mission. Aims. We focus on describing the objects and characteristics that Gaia will potentially observe without taking into consideration instrumental effects (detection efficiency, observing errors). Methods. The theoretical Universe Model prepared for the Gaia simulation has been statistically analyzed at a given time. Ingredients of the model are described, giving most attention to the stellar content, the double and multiple stars, and variability. Results. In this simulation the errors have not been included yet. Hence we estimate the number of objects and their theoretical photometric, astrometric and spectroscopic characteristics in the case that they are perfectly detected. We show that Gaia will be able to potentially observe 1.1 billion of stars (single or part of multiple star systems) of which about 2% are variable stars, 3% have one or two exoplanets. At the extragalactic level, observations will be potentially composed by several millions of galaxies, half million to 1 million of quasars and about 50,000 supernovas that will occur during the 5 years of mission. The simulated catalogue will be made publicly available by the DPAC on the Gaia portal of the ESA web site http://www.rssd.esa.int/gaia/.

We seek to constrain the formation of the Galactic bulge by means of analysing the detailed chemical composition of a large sample of red clump stars in Baade's window. We measure [Fe/H] in a sample of 219 bulge red clump stars from R=20000 resolution spectra obtained with FLAMES/GIRAFFE at the VLT, using an automatic procedure, differentially to the metal-rich local reference star muLeo. For a subsample of 162 stars, we also derive [Mg/H] from spectral synthesis around the MgI triplet at 6319A. The Fe and Mg metallicity distributions are both asymmetric, with median values of +0.16 and +0.21 respectively. The iron distribution is clearly bimodal, as revealed both by a deconvolution (from observational errors) and a Gaussian decomposition. The decomposition of the observed Fe and Mg metallicity distributions into Gaussian components yields two populations of equal sizes (50% each): a metal-poor component centred around [Fe/H]=-0.30 and [Mg/H]=-0.06 with a large dispersion and a narrow metal-rich component centred around [Fe/H]=+0.32 and [Mg/H]=+0.35. The metal poor component shows high [Mg/Fe] ratios (around 0.3) whereas stars in the metal rich component are found to have near solar ratios. Babusiaux et al. (2010) also find kinematical differences between the two components: the metal poor component shows kinematics compatible with an old spheroid whereas the metal rich component is consistent with a population supporting a bar. In view of their chemical and kinematical properties, we suggest different formation scenarios for the two populations: a rapid formation timescale as an old spheroid for the metal poor component (old bulge) and for the metal rich component, a formation over a longer time scale driven by the evolution of the bar (pseudo-bulge).

We obtained FLAMES-GIRAFFE spectra (R=22,500) at the ESO Very Large Telescope for 650 bulge red giant branch (RGB) stars and performed spectral synthesis to measure Mg, Ca, Ti, and Si abundances. This sample is composed of 474 giant stars observed in 3 fields along the minor axis of the Galactic bulge and at latitudes b=-4, b=-6, b=-12. Another 176 stars belong to a field containing the globular cluster NGC 6553, located at b=-3 and 5 degrees away from the other three fields along the major axis. Our results confirm, with large number statistics, the chemical similarity between the Galactic bulge and thick disk, which are both enhanced in alpha elements when compared to the thin disk. In the same context, we analyze [alpha/Fe] vs. [Fe/H] trends across different bulge regions. The most metal rich stars, showing low [alpha/Fe] ratios at b=-4 disappear at higher Galactic latitudes in agreement with the observed metallicity gradient in the bulge. Metal-poor stars ([Fe/H]<-0.2) show a remarkable homogeneity at different bulge locations. We have obtained further constrains for the formation scenario of the Galactic bulge. A metal-poor component chemically indistinguishable from the thick disk hints for a fast and early formation for both the bulge and the thick disk. Such a component shows no variation, neither in abundances nor kinematics, among different bulge regions. A metal-rich component showing low [alpha/Fe] similar to those of the thin disk disappears at larger latitudes. This allows us to trace a component formed through fast early mergers (classical bulge) and a disk/bar component formed on a more extended timescale.

Context. It has been shown that convection in red supergiant stars gives rise to large granules causing surface inhomogeneities together with shock waves in the photosphere. The resulting motion of the photocenter (on time scales ranging from months to years) could possibly have adverse effects on the parallax determination with Gaia. Aims. We explore the impact of the granulation on the photocentric and photometric variability. We quantify these effects in order to better characterize the error possibly altering the parallax. Methods. We use 3D radiative-hydrodynamics simulations of convection with CO5BOLD and the post-processing radiative transfer code OPTIM3D to compute intensity maps and spectra in the Gaia G band [325-1030 nm]. Results. We provide astrometric and photometric predictions from 3D simulations of RSGs that are used to evaluate the degradation of the astrometric parameters of evolved stars derived by Gaia. We show from RHD simulations that a supergiant like Betelgeuse exhibits a photocentric noise characterised by a standard deviation of the order of 0.1 AU. The number of bright giant and supergiant stars whose Gaia parallaxes will be altered by the photocentric noise ranges from a few tens to several thousandths. The degradation of the astrometric fit due to the presence of this photocentric noise will be noticeable up to about 5 kpc for the brightest supergiants. Moreover, parallaxes of supergiants are affected by a error of the order of a few percents. We show that the photocentric noise, as predicted by the 3D simulation, does account for a substantial part of the supplementary 'cosmic noise' that affects Hipparcos measurements of Betelgeuse and Antares.

We have chosen the name of GYES, one of the mythological giants with one hundred arms, offspring of Gaia and Uranus, for our instrument study of a multifibre spectrograph for the prime focus of the Canada-France-Hawaii Telescope. Such an instrument could provide an excellent ground-based complement for the Gaia mission and a northern complement to the HERMES project on the AAT. The CFHT is well known for providing a stable prime focus environment, with a large field of view, which has hosted several imaging instruments, but has never hosted a multifibre spectrograph. Building upon the experience gained at GEPI with FLAMES-Giraffe and X-Shooter, we are investigating the feasibility of a high multiplex spectrograph (about 500 fibres) over a field of view 1 degree in diameter. We are investigating an instrument with resolution in the range 15000 to 30000, which should provide accurate chemical abundances for stars down to 16th magnitude and radial velocities, accurate to 1 km/s for fainter stars. The study is led by GEPI-Observatoire de Paris with a contribution from Oxford for the study of the positioner. The financing for the study comes from INSU CSAA and Observatoire de Paris. The conceptual study will be delivered to CFHT for review by October 1st 2010.

Two main scenarios for the formation of the Galactic bulge are invoked, the first one through gravitational collapse or hierarchical merging of subclumps, the second through secular evolution of the Galactic disc. We aim to constrain the formation of the Galactic bulge through studies of the correlation between kinematics and metallicities in Baade's Window (l=1, b=-4) and two other fields along the bulge minor axis (l=0, b=-6 and b=-12). We combine the radial velocity and the [Fe/H] measurements obtained with FLAMES/GIRAFFE at the VLT with a spectral resolution of R=20000, plus for the Baade's Window field the OGLE-II proper motions, and compare these with published N-body simulations of the Galactic bulge. We confirm the presence of two distinct populations in Baade's Window found in Hill et al. 2010: the metal-rich population presents bar-like kinematics while the metal-poor population shows kinematics corresponding to an old spheroid or a thick disc one. In this context the metallicity gradient along the bulge minor axis observed by Zoccali et al. (2008), visible also in the kinematics, can be related to a varying mix of these two populations as one moves away from the Galactic plane, alleviating the apparent contradiction between the kinematic evidence of a bar and the existence of a metallicity gradient. We show evidences that the two main scenarios for the bulge formation co-exist within the Milky Way bulge.

Chemical abundances of six extremely metal-poor ([Fe/H]<-2.5) stars in the Sextans dwarf spheroidal galaxy are determined based on high resolution spectroscopy (R=40,000) with the Subaru Telescope High Dispersion Spectrograph. (1) The Fe abundances derived from the high resolution spectra are in good agreement with the metallicity estimated from the Ca triplet lines in low resolution spectra. The lack of stars with [Fe/H]=<-3 in Sextans, found by previous estimates from the Ca triplet, is confirmed by our measurements, although we note that high resolution spectroscopy for a larger sample of stars will be necessary to estimate the true fraction of stars with such low metallicity. (2) While one object shows an overabundance of Mg (similar to Galactic halo stars), the Mg/Fe ratios of the remaining five stars are similar to the solar value. This is the first time that low Mg/Fe ratios at such low metallicities have been found in a dwarf spheroidal galaxy. No evidence for over-abundances of Ca and Ti are found in these five stars, though the measurements for these elements are less certain. Possible mechanisms to produce low Mg/Fe ratios, with respect to that of Galactic halo stars, are discussed. (3) Ba is under-abundant in four objects, while the remaining two stars exhibit large and moderate excesses of this element. The abundance distribution of Ba in this galaxy is similar to that in the Galactic halo, indicating that the enrichment of heavy elements, probably by the r-process, started at metallicities [Fe/H] < -2.5, as found in the Galactic halo.

Hipparcos Intermediate Astrometric Data (IAD) have been used to derive astrometric orbital elements for spectroscopic binaries from the newly released Ninth Catalogue of Spectroscopic Binary Orbits (SB9). Among the 1374 binaries from SB9 which have an HIP entry, 282 have detectable orbital astrometric motion (at the 5% significance level). Among those, only 70 have astrometric orbital elements that are reliably determined (according to specific statistical tests discussed in the paper), and for the first time for 20 systems, representing a 10% increase relative to the 235 DMSA/O systems already present in the Hipparcos Double and Multiple Systems Annex. The detection of the astrometric orbital motion when the Hipparcos IAD are supplemented by the spectroscopic orbital elements is close to 100% for binaries with only one visible component, provided that the period is in the 50 - 1000 d range and the parallax is larger than 5 mas. This result is an interesting testbed to guide the choice of algorithms and statistical tests to be used in the search for astrometric binaries during the forthcoming ESA Gaia mission. Finally, orbital inclinations provided by the present analysis have been used to derive several astrophysical quantities. For instance, 29 among the 70 systems with reliable astrometric orbital elements involve main sequence stars for which the companion mass could be derived. Some interesting conclusions may be drawn from this new set of stellar masses, like the enigmatic nature of the companion to the Hyades F dwarf HIP 20935. This system has a mass ratio of 0.98 but the companion remains elusive.

The least well known structure of the Galaxy is its central region, because of the high extinction, the crowding and the confusion between disk and bulge sources along the line of sight. We show how the structure of the inner Galaxy and the dust distribution can be strongly constrained by using red clump stars as distance indicator. The results of this method applied on the deep near-infrared survey of the inner Galactic bulge made with the Cambridge Infrared Survey Instrument (CIRSI) are presented.

We present here a detailed study of the Sculptor dSph galaxy red giant branch (RGB) and horizontal branch (HB) morphology, combining new near-infrared photometry from CIRSI, with optical data from the ESO WFI. For a Sculptor-like old and generally metal-poor system, the position of RGB stars on the colour-magnitude diagram is mainly metallicity dependent. The advantage of using optical-NIR colours is that the position of the RGB locus is much more sensitive to metallicity than with optical colours alone. In contrast the horizontal branch (HB) morphology is strongly dependent on both metallicity and age. Therefore a detailed study of both the RGB in optical-NIR colours and the HB can help break the age-metallicity degeneracy. Our measured photometric width of the Sculptor giant branch corresponds to a range in metallicity of 0.75 dex. We detect the RGB and AGB bumps in both the NIR and optical luminosity functions, and derive from them a mean metallicity of [M/H] = -1.3 +/- 0.1. From isochrone fitting we derive a mean metallicity of [Fe/H] = -1.42 with a dispersion of 0.2 dex. These photometric estimators are for the first time consistent with individual metallicity measurements derived from spectroscopic observations. No spatial gradient is detected in the RGB morphology within a radius of 13 arcmin, twice the core radius. On the other hand, a significant gradient is observed in the HB morphology index, confirming the `second parameter problem' present in this galaxy. These observations are consistent with an early extended period of star formation continuing in time for a few Gyr. (Abridged)

We present a deep near-infrared wide-angle photometric analysis of the structure of the inner Galactic bar and central disk. The presence of a triaxial structure at the centre of the Galaxy is confirmed, consistent with a bar inclined at 22+/-5.5 deg from the Sun-Galactic centre line, extending to about 2.5 kpc from the Galactic centre and with a rather small axis ratio. A feature at l=-9.8 deg not aligned with this triaxiality suggests the existence of a second structure in the inner Galaxy, a double triaxiality or an inner ring. We argue that this is likely to be the signature of the end of the Galactic bar, at about 2.5-3 kpc, which is circumscribed by an inner pseudo-ring. No thick dust lane preceding the bar is detected, and a hole in the disc's dust distribution inside the bar radius is inferred.