Deciphering the distribution of metals throughout galaxies is fundamental in our understanding of galaxy evolution. Nearby, low-metallicity, star-forming dwarf galaxies in particular can offer detailed insight into the metal-dependent processes that may have occurred within galaxies in the early Universe. Here we present VLT/MUSE observations of one such system, JKB 18, a blue diffuse dwarf galaxy with a metallicity of only 12+log(O/H)=7.6$\pm$0.2 (~0.08 Zsol). Using high spatial-resolution integral field spectroscopy of the entire system, we calculate chemical abundances for individual HII regions using the direct method and derive oxygen abundance maps using strong-line metallicity diagnostics. With large-scale dispersions in O/H, N/H and N/O of ~0.5-0.6 dex and regions harbouring chemical abundances outside this 1$\sigma$ distribution, we deem JKB 18 to be chemically inhomogeneous. We explore this finding in the context of other chemically inhomogeneous dwarf galaxies and conclude that neither the accretion of metal-poor gas, short mixing timescales, or self-enrichment from Wolf-Rayet stars are accountable. Using a galaxy-scale, multi-phase, hydrodynamical simulation of a low-mass dwarf galaxy, we find that chemical inhomogeneities of this level may be attributable to the removal of gas via supernovae and the specific timing of the observations with respect to star-formation activity. This study not only draws attention to the fact that dwarf galaxies can be chemically inhomogeneous, but also that the methods used in the assessment of this characteristic can be subject to bias.

We present the first spectroscopic measurements of the ATLAS and Aliqa Uma streams from the Southern Stellar Stream Spectroscopic Survey ($S^5$), in combination with the photometric data from the Dark Energy Survey and astrometric data from $Gaia$. From the coherence of spectroscopic members in radial velocity and proper motion, we find out that these two systems are physically one stream with discontinuity in morphology and density on the sky (the "kink" feature). We refer to this entire stream as ATLAS-Aliqa Uma stream, or AAU stream. We perform a comprehensive exploration of the effect of baryonic substructures and find that only an encounter with the Sagittarius dwarf $\sim 0.5$ Gyr ago can create a feature similar to the observed "kink". In addition, we also identify two gaps in the ATLAS component associated with the broadening in the stream width (the "broadening" feature). These gaps have likely been created by small mass perturbers, such as dark matter halos, as AAU stream is the most distant cold stream known with severe variations in both the stream surface density and the stream track on the sky. With the stream track, stream distance and kinematic information, we determine the orbit of the AAU stream and find that it has been affected by the Large Magellanic Cloud, resulting in a misalignment between the proper motion and stream track. Together with the Orphan-Chenab Stream, AAU is the second stream pair that has been found to be a single stream separated into two segments by external perturbation.

In the era of large spectroscopic surveys, massive databases of high-quality spectra provide tools to outline a new picture of our Galaxy. In this framework, an important piece of information is provided by our ability to infer stellar ages. We aim to provide empirical relations between stellar ages and abundance ratios for a sample of solar-like stars. We investigate the dependence on metallicity, and we apply our relations to Gaia-ESO samples of open clusters and field stars. We analyse high-resolution and high-S/N HARPS spectra of a sample of solar-like stars to obtain precise determinations of their atmospheric parameters and abundances through differential spectral analysis and age through isochrone fitting. We investigate the relations between ages and abundance ratios. For the abundance ratios with a steeper dependence on age, we perform multivariate linear regressions, including the dependence on metallicity. We apply our relations to a sample of open clusters located in 4<R$_{GC}$<16 kpc. Using them, we are able to recover the literature ages only for clusters located at R$_{GC}$>7 kpc. In these clusters, the content of s-elements is lower than expected from chemical evolution models, and consequently the [s/$\alpha$] are lower than in clusters of the same age located in the solar neighbourhood. With our chemical evolution model and a set of empirical yields, we suggest that a strong dependence on the star formation history and metallicity-dependent yields of s-elements can substantially modify the slope of the [s/$\alpha$]-[Fe/H]-age relation in different regions of the Galaxy. Our results point towards a non-universal relation [s/$\alpha$]-[Fe/H]-age, indicating the existence of relations at different R$_{GC}$ or for different star formation history. A better understanding of the s-process at high metallicity is necessary to fully understand the origin of these variations.

We report the result of searching for globular clusters (GCs) around 56 Milky Way satellite dwarf galaxies within the distance of 450 kpc from the Galactic Center except for the Large and Small Magellanic Clouds and the Sagittarius dwarf. For each dwarf, we analyze the stellar distribution of sources in Gaia DR2, selected by magnitude, proper motion, and source morphology. Using the kernel density estimation of stellar number counts, we identify 61 possible GC candidates. Nine out of the 61 are known objects and only Fornax GC 1 -- 6 among them are associated with the dwarf galaxy whereas there is no strong evidence that the other 52 are GCs based on the analysis of existing imaging data. Using simulated GCs, we calculate the GC detection limit $M_{\rm V}^{\rm lim}$ that spans the range from $M_{\rm V}^{\rm lim} = -7$ for more distant dwarfs to $M_{\rm V}^{\rm lim} \sim 0$ for nearby systems. Assuming a Gaussian GC luminosity function, we compute that the completeness of the GC search is above 90 percent for most dwarf galaxies. We construct the 90 percent credible intervals/upper limits on the GC specific frequency $s$ of the MW dwarf galaxies: $12 < s < 47$ for the Fornax, $s < 20$ for the dwarfs with $-12 < M_{\rm V} < -10$, $s < 30$ for the dwarfs with $-10 < M_{\rm V} < -7$, and $s < 90$ for the dwarfs with $M_{\rm V} > -7$.

The recent discoveries of gravitational wave events and in one case also its electromagnetic (EM) counterpart allow us to study the Universe in a novel way. The increased sensitivity of the LIGO and Virgo detectors has opened the possibility for regular detections of EM transient events from mergers of stellar remnants. Gravitational wave sources are expected to have sky localisation up to a few hundred square degrees, thus Gaia as an all-sky multi-epoch photometric survey has the potential to be a good tool to search for the EM counterparts. In this paper we study the possibility of detecting EM counterparts to gravitational wave sources using the Gaia Science Alerts system. We develop an extension to current used algorithms to find transients and test its capabilities in discovering candidate transients on a sample of events from the observation periods O1 and O2 of LIGO and Virgo. For the gravitational wave events from the current run O3 we expect that about 16 (25) per cent should fall in sky regions observed by Gaia 7 (10) days after gravitational wave. The new algorithm will provide about 10 candidates per day from the whole sky.

Most binaries are undetected. Astrometric reductions of a system using the assumption that the object moves like a single point mass can be biased by unresolved binary stars. The discrepancy between the centre of mass of the system (which moves like a point mass) and the centre of light (which is what we observe) introduces additional motion. We explore the extent to which binary systems affect single object models fit to astrometric data. This tells us how observations are diluted by binaries and which systems cause the largest discrepancies - but also allows us to make inferences about the binarity of populations based on observed astrometric error. By examining a sample of mock observations, we show that binaries with periods close to one year can mimic parallax and thus bias distance measurements, whilst long period binaries can introduce significant apparent proper motion. Whilst these changes can soak up some of the error introduced by the binary, the total deviation from the best fitting model can be translated into a lower limit on the on-sky separation of the pair. Throughout we link these predictions to data from the Gaia satellite, whilst leaving the conclusions generalizable to other surveys.

For stars with unresolved companions, motions of the centre of light and that of mass decouple, causing a single-source astrometric model to perform poorly. We show that such stars can be easily detected with the reduced chi2 statistic, or RUWE, provided as part of Gaia DR2. We convert RUWE into the amplitude of the image centroid wobble, which, if scaled by the source distance, is proportional to the physical separation between companions (for periods up to several years). We test this idea on a sample of known spectroscopic binaries and demonstrate that the amplitude of the centroid perturbation scales with the binary period and the mass ratio as expected. We apply this technique to the Gaia DR2 data and show how the binary fraction evolves across the Hertzsprung--Russell diagram. The observed incidence of unresolved companions is high for massive young stars and drops steadily with stellar mass, reaching its lowest levels for white dwarfs. We highlight the elevated binary fraction for the nearby Blue Stragglers and Blue Horizontal Branch stars. We also illustrate how unresolved hierarchical triples inflate the relative velocity signal in wide binaries. Finally, we point out a hint of evidence for the existence of additional companions to the hosts of extrasolar hot jupiters.

We present the serendipitous discovery of a low luminosity nova occurring in a symbiotic binary star system in the Milky Way. We lay out the extensive archival data alongside new follow-up observations related to the stellar object V$^*$ CN Cha in the constellation of Chamaeleon. The object had long period ($\sim\! 250\,$day), high amplitude ($\sim\! 3\,$mag) optical variability in its recent past, preceding an increase in optical brightness by $\sim\! 8\,$magnitudes and a persistence at this luminosity for about 3 years, followed by a period of $\sim\! 1.4\,{\rm mag}\,{\rm yr}^{-1}$ dimming. The object's current optical luminosity seems to be dominated by H$\alpha$ emission, which also exhibits blue-shifted absorption (a P-Cygni-like profile). After consideration of a number of theories to explain these myriad observations, we determine that V$^*$ CN Cha is most likely a symbiotic (an evolved star-white dwarf binary) system which has undergone a long-duration, low luminosity, nova. Interpreted in this way, the outburst in V$^*$ CN Cha is among the lowest luminosity novae ever observed.

Using the data from the Sloan Digital Sky Survey and the Gaia satellite, we assemble a pure sample of $\sim$3000 Blue Horizontal Branch (BHB) stars with 7-D information, including positions, velocities and metallicities. We demonstrate that, as traced with BHBs, the Milky Way's stellar halo is largely unmixed and can not be well represented with a conventional Gaussian velocity distribution. A single-component model fails because the inner portions of the halo are swamped with metal-rich tidal debris from an ancient, head-on collision, known as the "Gaia Sausage". Motivated by the data, we build a flexible mixture model which allows us to track the evolution of the halo make-up across a wide range of radii. It is built from two components, one representing the radially anisotropic Sausage stars with their lobed velocity distribution, the other representing a more metal-poor and more isotropic component built up from minor mergers. We show that inside 25 kpc the "Sausage" contributes at least 50 % of the Galactic halo. The fraction of "Sausage" stars diminishes sharply beyond 30 kpc, which is the long-established break radius of the classical stellar halo.

We report the discovery of a young (${\rm age} \sim 130~{\rm Myr}$), low-mass ($M \sim 1200~{\rm M}_\odot$), metal-poor ($[{\rm Fe}/{\rm H}] \sim -1.1$) stellar association at a heliocentric distance $D \approx 29~{\rm kpc}$, placing it far into the Milky Way halo. At its present Galactocentric position $(R, z) \sim (23, 15)~{\rm kpc}$, the association is (on the sky) near the leading arm of the gas stream emanating from the Magellanic cloud system, but is located $\approx 60^\circ$ from the Large Magellanic Cloud (LMC) center on the other side of the Milky Way disk. If we assume that the cluster is co-located with HI gas in the stream, we directly measure the distance to the leading arm of the Magellanic stream. The measured distance is inconsistent with stream predictions from models of the LMC/SMC interaction and infall into the Milky Way that do not account for ram pressure and gas interaction with Milky Way disk. The estimated age of the cluster is consistent with the time of last passage of the leading arm gas through the Galactic midplane, and we therefore speculate that this star-formation event was triggered by its the last disk midplane passage. Most details of this idea remain a puzzle: the Magellanic stream has low column density, the Milky Way disk at this large radius has low gas density, and the relative velocity of the leading arm gas and Milky Way gas is large. However it formed, the discovery of a young stellar cluster in the Milky Way halo presents an interesting opportunity for study. This cluster was discovered with Gaia astrometry and photometry alone, but folow-up DECam photometry was crucial for measuring its properties.

Open clusters (OCs) are recognised as excellent tracers of Galactic thin-disc properties. At variance with intermediate-age and old OCs, for which a significant number of studies is now available, clusters younger than 150 Myr have been mostly overlooked in terms of their chemical composition, with few exceptions. On the other hand, previous investigations seem to indicate an anomalous behaviour of young clusters, which includes slightly sub-solar iron (Fe) abundances and extreme, unexpectedly high barium (Ba) enhancements. In a series of papers, we plan to expand our understanding of this topic and investigate whether these chemical peculiarities are instead related to abundance analysis techniques. We present a new determination of the atmospheric parameters for 23 dwarf stars observed by the Gaia-ESO survey in five young OCs (younger than 150 Myr) and one star-forming region (NGC 2264). We exploit a new method based on titanium (Ti) lines to derive the spectroscopic surface gravity, and most importantly, the microturbulence parameter. A combination of Ti I and Fe I lines is used to obtain effective temperatures. We also infer the abundances of Fe II, Ti II, Na I, Mg I, Al I, Si I, Ca I, Cr I and Ni I. Our findings are in fair agreement with Gaia-ESO iDR5 results for effective temperatures and surface gravities, but suggest that for very young stars, the microturbulence parameter is over-estimated when Fe lines are employed. This affects the derived chemical composition and causes the metal content of very young clusters to be under-estimated. Our clusters display a metallicity [Fe/H] between +0.04 and +0.12; they are not more metal poor than the Sun. Although based on a relatively small sample size, our explorative study suggests that we may not need to call for ad hoc explanations to reconcile the chemical composition of young OCs with Galactic chemical evolution models.

We present the first detailed elemental abundances in the ultra-faint Magellanic satellite galaxies Carina II (Car II) and Carina III (Car III). With high-resolution Magellan/MIKE spectroscopy, we determined abundances of nine stars in Car II including the first abundances of an RR Lyrae star in an ultra-faint dwarf galaxy; and two stars in Car III. The chemical abundances demonstrate that both systems are clearly galaxies and not globular clusters. The stars in these galaxies mostly display abundance trends matching those of other similarly faint dwarf galaxies: enhanced but declining [alpha/Fe] ratios, iron-peak elements matching the stellar halo, and unusually low neutron-capture element abundances. One star displays a low outlying [Sc/Fe] = -1.0. We detect a large Ba scatter in Car II, likely due to inhomogeneous enrichment by low-mass AGB star winds. The most striking abundance trend is for [Mg/Ca] in Car II, which decreases from +0.4 to -0.4 and indicates clear variation in the initial progenitor masses of enriching core-collapse supernovae. So far, the only ultra-faint dwarf galaxies displaying a similar [Mg/Ca] trend are likely satellites of the Large Magellanic Cloud. We find two stars with [Fe/H] < -3.5, whose abundances likely trace the first generation of metal-free Population III stars and are well-fit by Population III core-collapse supernova yields. An appendix describes our new abundance uncertainty analysis that propagates line-by-line stellar parameter uncertainties.

We introduce the Southern Stellar Stream Spectroscopy Survey (${S}^5$), an on-going program to map the kinematics and chemistry of stellar streams in the Southern Hemisphere. The initial focus of ${S}^5$ has been spectroscopic observations of recently identified streams within the footprint of the Dark Energy Survey (DES), with the eventual goal of surveying streams across the entire southern sky. Stellar streams are composed of material that has been tidally striped from dwarf galaxies and globular clusters and hence are excellent dynamical probes of the gravitational potential of the Milky Way, as well as providing a detailed snapshot of its accretion history. Observing with the 3.9-m Anglo-Australian Telescope's 2-degree-Field fibre positioner and AAOmega spectrograph, and combining the precise photometry of DES DR1 with the superb proper motions from $Gaia$ DR2, allows us to conduct an efficient spectroscopic survey to map these stellar streams. So far ${S}^5$ has mapped 9 DES streams and 3 streams outside of DES; the former are the first spectroscopic observations of these recently discovered streams. In addition to the stream survey, we use spare fibres to undertake a Milky Way halo survey and a low-redshift galaxy survey. This paper presents an overview of the ${S}^5$ program, describing the scientific motivation for the survey, target selection, observation strategy, data reduction and survey validation. Finally, we describe early science results on stellar streams and Milky Way halo stars drawn from the survey. Updates on ${S}^5$, including future public data release, can be found at \url{http://s5collab.github.io}.

We present the serendipitous discovery of the fastest Main Sequence hyper-velocity star (HVS) by the Southern Stellar Stream Spectroscopic Survey (S5). The star S5-HVS1 is a $\sim 2.35$ M$_\odot$ A-type star located at a distance of $\sim 9$ kpc from the Sun and has a heliocentric radial velocity of $1017\pm 2.7$ km/s without any signature of velocity variability. The current 3-D velocity of the star in the Galactic frame is $1755\pm50$ km/s. When integrated backwards in time, the orbit of the star points unambiguously to the Galactic Centre, implying that S5-HVS1 was kicked away from Sgr A* with a velocity of $\sim 1800$ km/s and travelled for $4.8$ Myr to its current location. This is so far the only HVS confidently associated with the Galactic Centre. S5-HVS1 is also the first hyper-velocity star to provide constraints on the geometry and kinematics of the Galaxy, such as the Solar motion $V_{y,\odot}= 246.1\pm 5.3$ km/s or position $R_0=8.12\pm 0.23$ kpc. The ejection trajectory and transit time of S5-HVS1 coincide with the orbital plane and age of the annular disk of young stars at the Galactic centre, and thus may be linked to its formation. With the S5-HVS1 ejection velocity being almost twice the velocity of other hyper-velocity stars previously associated with the Galactic Centre, we question whether they have been generated by the same mechanism or whether the ejection velocity distribution has been constant over time.

We present two catalogues of active galactic nucleus (AGN) candidates selected from the latest data of two all-sky surveys -- Data Release 2 (DR2) of the \emph{Gaia} mission and the unWISE catalogue of the \emph{Wide-field Infrared Survey Explorer} (\emph{WISE}). We train a random forest classifier to predict the probability of each source in the \emph{Gaia}-unWISE joint sample being an AGN, $P_{\rm RF}$, based on \emph{Gaia} astrometric and photometric measurements and unWISE photometry. The two catalogues, which we designate C75 and R85, are constructed by applying different $P_{\rm RF}$ threshold cuts to achieve an overall completeness of 75\% ($\approx$90\% at \emph{Gaia} $G\leq20$ mag) and reliability of 85\% respectively. The C75 (R85) catalogue contains 2,734,464 (2,182,193) AGN candidates across the effective 36,000 deg$^2$ sky, of which $\approx$0.91 (0.52) million are new discoveries. Photometric redshifts of the AGN candidates are derived by a random forest regressor using \emph{Gaia} and \emph{WISE} magnitudes and colours. The estimated overall photometric redshift accuracy is 0.11. Cross-matching the AGN candidates with a sample of known bright cluster galaxies, we identify a high-probability strongly-lensed AGN candidate system, SDSS\,J1326$+$4806, with a large image separation of 21\farcs06. All the AGN candidates in our catalogues will have $\sim$5-year long light curves from \emph{Gaia} by the end of the mission, and thus will be a great resource for AGN variability studies. Our AGN catalogues will also be helpful in AGN target selections for future spectroscopic surveys, especially ones in the southern hemisphere. The C75 catalogue can be downloaded at https://www.ast.cam.ac.uk/~ypshu/AGN_Catalogues.html.

Using Gaia DR2, we trace the Anticenter Stream (ACS) in various stellar populations across the sky and find that it is kinematically and spatially decoupled from the Monoceros Ring. Using stars from {\sc lamost} and {\sc segue}, we show that the ACS is systematically more metal-poor than Monoceros by $0.1$ dex with indications of a narrower metallicity spread. Furthermore, the ACS is predominantly populated of old stars ($\sim 10\,\rm{Gyr}$), whereas Monoceros has a pronounced tail of younger stars ($6-10\, \rm{Gyr}$) as revealed by their cumulative age distributions. Put togehter, all of this evidence support predictions from simulations of the interaction of the Sagittarius dwarf with the Milky Way, which argue that the Anticenter Stream (ACS) is the remains of a tidal tail of the Galaxy excited during Sgr's first pericentric passage after it crossed the virial radius, whereas Monoceros consists of the composite stellar populations excited during the more extended phases of the interaction. We suggest that the ACS can be used to constrain the Galactic potential, particularly its flattening, setting strong limits on the existence of a dark disc. Importantly, the ACS can be viewed as a stand-alone fossil of the chemical enrichment history of the Galactic disc.

We cross-match high-precision astrometric data from Gaia DR2 with accurate multi-band photometry from the Dark Energy Survey (DES) DR1 to confidently measure proper motions for nine stellar streams in the DES footprint: Aliqa Uma, ATLAS, Chenab, Elqui, Indus, Jhelum, Phoenix, Tucana III, and Turranburra. We determine low-confidence proper motion measurements for four additional stellar streams: Ravi, Wambelong, Willka Yaku, and Turbio. We find evidence for a misalignment between stream tracks and the systemic proper motion of streams that may suggest a systematic gravitational influence from the Large Magellanic Cloud. These proper motions, when combined with radial velocity measurements, will allow for detailed orbit modeling which can be used to constrain properties of the LMC and its affect on nearby streams, as well as global properties of the Milky Way's gravitational potential.

In the era of large high-resolution spectroscopic surveys, high-quality spectra can contribute to our understanding of the Galactic chemical evolution, providing chemical abundances belonging to the different nucleosynthesis channels, and also providing constraints to stellar age. Some abundance ratios have been proven to be excellent indicators of stellar ages. We aim at providing an empirical relationship between stellar ages and [C/N] using, as calibrators, open star clusters observed by both the Gaia-ESO and APOGEE surveys. We use stellar parameters and abundances from the Gaia-ESO and APOGEE of the Galactic field and open cluster stars. Ages of star clusters are retrieved from the literature sources and validated using a common set of isochrones. We use the same isochrones to determine, for each age and metallicity, the surface gravity at which the first dredge-up and red giant branch bump occur. We study the effect of extra-mixing processes in our sample of giant stars, and we derive the mean [C/N] in evolved stars, including only stars without evidence of extra-mixing. Combining the Gaia-ESO and APOGEE samples of open clusters, we derive a linear relationship between [C/N] and logarithmic cluster ages. We apply our relationship to selected giant field stars in both Gaia-ESO and APOGEE. We find an age separation between thin and thick disc stars and age trends within their populations, with an increasing age towards lower metallicity populations. With such empirical relationship, we are able to provide an age estimate for giant stars in which C and N abundances are measured. Isochrone fitting is less sensitive for giant than dwarf stars at the turn off. The present method can be thus considered as an additional tool to give an independent estimate of the age of giant stars, with uncertainties in their ages comparable to those obtained using isochrone fitting for dwarf stars.

We report the discovery of a Milky Way satellite in the constellation of Antlia. The Antlia 2 dwarf galaxy is located behind the Galactic disc at a latitude of $b\sim 11^{\circ}$ and spans $1.26$ degrees, which corresponds to $\sim2.9$ kpc at its distance of 130 kpc. While similar in spatial extent to the Large Magellanic Cloud, Antlia 2 is orders of magnitude fainter at $M_V=-9$ mag, making it by far the lowest surface brightness system known (at $\sim31.9$ mag/arcsec$^2$), $\sim100$ times more diffuse than the so-called ultra diffuse galaxies. The satellite was identified using a combination of astrometry, photometry and variability data from \textit{Gaia} Data Release 2, and its nature confirmed with deep archival DECam imaging, which revealed a conspicuous BHB signal. We have also obtained follow-up spectroscopy using AAOmega on the AAT, identifying 159 member stars, and we used them to measure the dwarf's systemic velocity, $290.9\pm0.5$km/s, its velocity dispersion, $5.7\pm1.1$ km/s, and mean metallicity, [Fe/H]$=-1.4$. From these properties we conclude that Antlia 2 inhabits one of the least dense Dark Matter (DM) halos probed to date. Dynamical modelling and tidal-disruption simulations suggest that a combination of a cored DM profile and strong tidal stripping may explain the observed properties of this satellite. The origin of this core may be consistent with aggressive feedback, or may even require alternatives to cold dark matter (such as ultra-light bosons).

In a companion paper by Koposov et al., RR Lyrae from \textit{Gaia} Data Release 2 are used to demonstrate that stars in the Orphan stream have velocity vectors significantly misaligned with the stream track, suggesting that it has received a large gravitational perturbation from a satellite of the Milky Way. We argue that such a mismatch cannot arise due to any realistic static Milky Way potential and then explore the perturbative effects of the Large Magellanic Cloud (LMC). We find that the LMC can produce precisely the observed motion-track mismatch and we therefore use the Orphan stream to measure the mass of the Cloud. We simultaneously fit the Milky Way and LMC potentials and infer that a total LMC mass of $1.38^{+0.27}_{-0.24} \times10^{11}\,\rm{M_\odot}$ is required to bend the Orphan Stream, showing for the first time that the LMC has a large and measurable effect on structures orbiting the Milky Way. This has far-reaching consequences for any technique which assumes that tracers are orbiting a static Milky Way. Furthermore, we measure the Milky Way mass within 50 kpc to be $3.80^{+0.14}_{-0.11}\times10^{11} M_\odot$. Finally, we use these results to predict that, due to the reflex motion of the Milky Way in response to the LMC, the outskirts of the Milky Way's stellar halo should exhibit a bulk, upwards motion.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Using RR Lyrae stars in the Gaia Data Release 2 and Pan-STARRS1 we study the properties of the Pisces Overdensity, a diffuse sub-structure in the outer halo of the Milky Way. We show that along the line of sight, Pisces appears as a broad and long plume of stars stretching from 40 to 110 kpc with a steep distance gradient. On the sky Pisces's elongated shape is aligned with the Magellanic Stream. Using follow-up VLT FORS2 spectroscopy, we have measured the velocity distribution of the Pisces candidate member stars and have shown it to be as broad as that of the Galactic halo but offset to negative velocities. Using a suite of numerical simulations, we demonstrate that the structure has many properties in common with the predicted behaviour of the Magellanic wake, i.e. the Galactic halo overdensity induced by the in-fall of the Magellanic Clouds.

Variable stars with well-calibrated period-luminosity relationships provide accurate distance measurements to nearby galaxies and are therefore a vital tool for cosmology and astrophysics. While these measurements typically rely on samples of Cepheid and RR-Lyrae stars, abundant populations of luminous variable stars with longer periods of $10 - 1000$ days remain largely unused. We apply machine learning to derive a mapping between lightcurve features of these variable stars and their magnitude to extend the traditional period-luminosity (PL) relation commonly used for Cepheid samples. Using photometric data for long period variable stars in the Large Magellanic cloud (LMC), we demonstrate that our predictions produce residual errors comparable to those obtained on the corresponding Cepheid population. We show that our model generalizes well to other samples by performing a blind test on photometric data from the Small Magellanic Cloud (SMC). Our predictions on the SMC again show small residual errors and biases, comparable to results that employ PL relations fitted on Cepheid samples. The residual biases are complementary between the long period variable and Cepheid fits, which provides exciting prospects to better control sources of systematic error in cosmological distance measurements. We finally show that the proposed methodology can be used to optimize samples of variable stars as standard candles independent of any prior variable star classification.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The combination of precise radial velocities from multi-object spectroscopy and highly accurate proper motions from Gaia DR2 opens up the possibility for detailed 3D kinematic studies of young star forming regions and clusters. Here, we perform such an analysis by combining Gaia-ESO Survey spectroscopy with Gaia astrometry for ~900 members of the Lagoon Nebula cluster, NGC 6530. We measure the 3D velocity dispersion of the region to be $5.35^{+0.39}_{-0.34}$~km~s$^{-1}$, which is large enough to suggest the region is gravitationally unbound. The velocity ellipsoid is anisotropic, implying that the region is not sufficiently dynamically evolved to achieve isotropy, though the central part of NGC 6530 does exhibit velocity isotropy that suggests sufficient mixing has occurred in this denser part. We find strong evidence that the stellar population is expanding, though this is preferentially occurring in the declination direction and there is very little evidence for expansion in the right ascension direction. This argues against a simple radial expansion pattern, as predicted by models of residual gas expulsion. We discuss these findings in the context of cluster formation, evolution and disruption theories.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Since first noticed by Shapley in 1939, a faint object coincident with the Fornax dwarf spheroidal has long been discussed as a possible sixth globular cluster system. However, debate has continued over whether this overdensity is a statistical artifact or a blended galaxy group. In this Letter we demonstrate, using deep DECam imaging data, that this object is well resolved into stars and is a bona fide star cluster. The stellar overdensity of this cluster is statistically significant at the level of ~ 6 - 6.7 sigma in several different photometric catalogs including Gaia. Therefore, it is highly unlikely to be caused by random fluctuation. We show that Fornax 6 is a star cluster with a peculiarly low surface brightness and irregular shape, which may indicate a strong tidal influence from its host galaxy. The Hess diagram of Fornax 6 is largely consistent with that of Fornax field stars, but it appears to be slightly bluer. However, it is still likely more metal-rich than most of the globular clusters in the system. Faint clusters like Fornax 6 that orbit and potentially get disrupted in the centers of dwarf galaxies can prove crucial for constraining the dark matter distribution in Milky Way satellites.

We use astrometry, broad-band photometry and variability information from the Data Release 2 of ESA's Gaia mission (GDR2) to identify members of the Orphan Stream (OS) across the whole sky. The stream is traced above and below the celestial equator and in both Galactic hemispheres, thus increasing its visible length to ~ 210 degrees equivalent to ~150 kpc in physical extent. Taking advantage of the large number of RR Lyrae stars in the OS, we extract accurate distances and proper motions across the entire stretch of the tidal debris studied. As delineated by the GDR2 RR Lyrae, the stream exhibits two prominent twists in its shape on the sky which are accompanied by changes in the tangential motion. We complement the RR Lyrae maps with those created using GDR2 Red Giants and the DECam Legacy Survey Main Sequence Turn-Off stars. The behavior of the OS track on the sky is consistent across all three tracers employed. We detect a strong non-zero motion in the across-stream direction for a substantial portion of the stream. Such a misalignment between the debris track and the streaming velocity cannot be reproduced in a static gravitational potential and signals an interaction with a massive perturber.

In very young clusters, stellar age distribution is the empirical proof of the duration of star formation (SF) and of the physical mechanisms involved in the process. We derived accurate stellar ages for the cluster NGC6530, associated with the Lagoon Nebula to infer its SF history. We use the Gaia-ESO survey observations and Gaia DR2 data, to derive cluster membership and fundamental stellar parameters. We identified 652 confirmed and 9 probable members. The reddening inferred for members and non-members allows us to distinguish MS stars and giants, in agreement with the distances inferred from Gaia DR2 data. The foreground and background stars show a spatial pattern that traces the 3D structure of the nebular dust component. We derive stellar ages for 382 confirmed cluster members and we find that the gravity-sensitive gamma index distribution for M stars is correlated with stellar age. For all members with Teff<5500 K, the mean logarithmic age is 5.84 (units of years) with a dispersion of 0.36 dex. The age distribution of stars with accretion and/or disk (CTTSe) is similar to that of stars without accretion and without disk (WTTSp). We interpret this dispersion as evidence of a real age spread since the total uncertainties on age determinations, derived from Monte Carlo simulations, are significantly smaller than the observed spread. This conclusion is supported by the evidence of a decreasing of the gravity-sensitive gamma index as a function of stellar ages. The presence of the age spread is also supported by the spatial distribution and the kinematics of old and young members. In particular, members with accretion and/or disk, formed in the last 1 Myr, show evidence of subclustering around the cluster center, in the Hourglass Nebula and in the M8-E region, suggesting a possible triggering of star formation events by the O-type star ionization fronts.

Gaia DR2 5932173855446728064 was recently proposed to be unbound from the Milky Way based on the $-614.3\pm2.5\;\mathrm{km}\;\mathrm{s}^{-1}$ median radial velocity given in Gaia DR2. We obtained eight epochs of spectroscopic follow-up and find a very different median radial velocity of $-56.5 \pm 5.3\;\mathrm{km}\;\mathrm{s}^{-1}$. If this difference were to be explained by binarity, then the unseen companion would be an intermediate-mass black hole; we therefore argue that the Gaia DR2 radial velocity must be in error. We find it likely that the spectra obtained by Gaia were dominated by the light from a star $4.3\;\mathrm{arcsec}$ away, and that, due to the slitless, time delay integration nature of Gaia spectroscopy, this angular offset corresponded to a spurious $620\;\mathrm{km}\;\mathrm{s}^{-1}$ shift in the calcium triplet of the second star. We argue that such unanticipated alignments between stars may account for 105 of the 202 stars with radial velocities faster than $500\;\mathrm{km}\;\mathrm{s}^{-1}$ in Gaia DR2 and propose a quality cut to exclude stars that are susceptible. We propose further cuts to remove stars where the colour photometry is suspect and stars where the radial velocity measurement is based on fewer than four transits, and thus produce an unprecedentedly clean selection of Gaia RVS stars for use in studies of Galactic dynamics.

We present the results of a systematic Milky Way satellite search performed across an array of publicly available wide-area photometric surveys. Our aim is to complement previous searches by widening the parameter space covered. Specifically, we focus on objects smaller than $1'$ and include old, young, metal poor and metal rich stellar population masks. As a result we find 9 new likely genuine stellar systems in data from GAIA, DES, and Pan-STARRS, which were picked from the candidate list because of conspicuous counterparts in the cut-out images. The presented systems are all very compact ($r_h<1'$) and faint ($M_V\gtrsim-3$), and are associated either with the Galactic disk, or the Magellanic Clouds. While most of the stellar systems look like Open Clusters, their exact classification is, as of today, unclear. With these discoveries, we extend the parameter space occupied by star clusters to sizes and luminosities previously unexplored and demonstrate that rather than two distinct classes of Globular and Open clusters, there appears to be a continuity of objects, unmarked by a clear decision boundary.

We carried out a spectroscopic follow-up program of the four new stellar stream candidates detected by Belokurov & Koposov (2016) in the outskirts of the Large Magellanic Cloud (LMC) using FORS2 (VLT). The medium-resolution spectra were used to measure the line-of-sight velocities, estimate stellar metallicities and to classify stars into Blue Horizontal Branch (BHB) and Blue Straggler (BS) stars. Using the 4-D phase-space information, we attribute approximately one half of our sample to the Magellanic Clouds, while the rest is part of the Galactic foreground. Only two of the four stream candidates are confirmed kinematically. While it is impossible to estimate the exact levels of MW contamination, the phase-space distribution of the entire sample of our Magellanic stars matches the expected velocity gradient for the LMC halo and extends as far as 33 deg (angular separation) or 29 kpc from the LMC center. Our detections reinforce the idea that the halo of the LMC seems to be larger than previously expected, and its debris can be spread in the sky out to very large separations from the LMC center. Finally, we provide some kinematic evidence that many of the stars analysed here have likely come from the Small Magellanic Cloud.

Using deep wide-field photometry three-year data (Y3) from the Dark Energy Survey (DES), we present a panoramic study of the Fornax dwarf spheroidal galaxy. The data presented here -- a small subset of the full survey -- uniformly covers a region of 25 square degrees centered on the galaxy to a depth of g ~ 23.5. We use this data to study the structural properties of Fornax, overall stellar population, and its member stars in different evolutionary phases. We also search for possible signs of tidal disturbance. Fornax is found to be significantly more spatially extended than what early studies suggested. No statistically significant distortions or signs of tidal disturbances were found down to a surface brightness limit of ~32.1 mag/arcsec^2. However, there are hints of shell-like features located ~30' - 40' from the center of Fornax that may be stellar debris from past merger events. We also find that intermediate age and young main-sequence populations show different orientation at the galaxy center and have many substructures. The deep DES Y3 data allows us to characterize the age of those substructures with great accuracy, both those previously known and those newly discovered in this work, on the basis of their color-magnitude diagram morphology. We find that the youngest overdensities are all found on the Eastern side of Fornax, where the Fornax field population itself is slightly younger than in the West. The high quality DES Y3 data reveals that Fornax has many rich structures, and provides insights into its complex formation history.

The high spatial resolution and precise astrometry and photometry of the Gaia mission should make it particularly apt at discovering and resolving transients occurring in, or near, the centres of galaxies. Indeed, some nuclear transients are reported by the Gaia Science Alerts (GSA) team, but not a single confirmed Tidal Disruption Event has been published. In order to explore the sensitivity of GSA, we performed an independent and systematic search for nuclear transients using Gaia observations. Our transient search is driven from an input galaxy catalogue (derived from the Sloan Digital Sky Survey Release 12). We present a candidate detection metric which is independent from the existing GSA methodology, to see if Gaia Alerts are biased against the discovery of nuclear transients, and in particular which steps may have an impact. Our technique does require significant manual vetting of candidates, making implementation in the GSA system impractical for daily operations, although it could be run weekly, which for month-to-year long transients would make a scientifically valuable addition. Our search yielded ~480 nuclear transients, 5 of which were alerted and published by GSA. The list of (in some cases ongoing) transients includes candidates for events related to enhanced accretion onto a super-massive black hole and TDEs. An implementation of the detection methodology and criteria used in this paper as an extension of GSA could open up the possibility for Gaia to fulfil the role as a main tool to find transient nuclear activity as predicted in the literature.

We report the discovery of a thin stellar stream - which we name the Jet stream - crossing the constellations of Hydra and Pyxis. The discovery was made in data from the SLAMS survey, which comprises deep $g$ and $r$ imaging for a $650$ square degree region above the Galactic disc performed by the CTIO Blanco + DECam. SLAMS photometric catalogues will be made publicly available. The stream is approximately 0.18 degrees wide and 10 degrees long, though it is truncated by the survey footprint. Its colour-magnitude diagram is consistent with an old, metal-poor stellar population at a heliocentric distance of approximately 29 kpc. We corroborate this measurement by identifying a spatially coincident overdensity of likely blue horizontal branch stars at the same distance. There is no obvious candidate for a surviving stream progenitor.

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.

We measure the orbital properties of halo stars using 7-dimensional information provided by Gaia and the Sloan Digital Sky Survey. A metal-rich population of stars, present in both local main sequence stars and more distant blue horizontal branch stars, have very radial orbits (eccentricity ~ 0.9) and apocenters that coincide with the stellar halo "break radius" at galactocentric distance r ~ 20 kpc. Previous work has shown that the stellar halo density falls off much more rapidly beyond this break radius. We argue that the correspondence between the apocenters of high metallicity, high eccentricity stars and the broken density profile is caused by the build-up of stars at the apocenter of a common dwarf progenitor. Although the radially biased stars are likely present down to metallicities of [Fe/H] ~ -2 the increasing dominance at higher metallicities suggests a massive dwarf progenitor, which is at least as massive as the Fornax and Sagittarius dwarf galaxies, and is likely the dominant progenitor of the inner stellar halo.

We use a sample of ~350 RR Lyrae stars with radial velocities and Gaia DR2 proper motions to study orbital properties of the Hercules-Aquila Cloud (HAC) and Virgo Over-density (VOD). We demonstrate that both structures are dominated by stars on highly eccentric orbits, with peri-centres around ~1 kpc and apo-centres between 15 and 25 kpc from the Galactic centre. Given that the stars in the HAC and the VOD occupy very similar regions in the space spanned by integrals of motion, we conclude that these diffuse debris clouds are part of the same accretion event. More precisely, these inner halo sub-structures likely represent two complementary not-fully-mixed portions of an ancient massive merger, also known as the "sausage" event.

The Gaia Sausage is an elongated structure in velocity space discovered by Belokurov et al. (2018) using the kinematics of metal-rich halo stars. It was created by a massive dwarf galaxy ($\sim 5 \times 10^{10} M_\odot$) on a strongly radial orbit that merged with the Milky Way at a redshift $z\lesssim 3$. We search forthe associated Sausage Globular Clusters by analysing the structure of 91 Milky Way globular clusters (GCs) in action space using the Gaia Data Release 2 catalogue, complemented with Hubble Space Telescope proper motions. There is a characteristic energy $E_{\rm crit}$ which separates the in situ objects, such as the bulge/disc clusters, from the accreted objects, such as the young halo clusters. There are 15 old halo GCs that have $E > E_{\rm crit}$. Eight of the high energy, old halo GCs are strongly clumped in azimuthal and vertical action, yet strung out like beads on a chain at extreme radial action. They are very radially anisotropic ($\beta \sim 0.95$) and move on orbits that are all highly eccentric ($e \gtrsim 0.80$). They also form a track in the age-metallicity plane distinct from the bulk of the Milky Way GCs and compatible with a dwarf spheroidal origin. These properties are consistent with GCs associated with the merger event that gave rise to the Gaia Sausage.

We use the SDSS-Gaia catalogue to search for substructure in the stellar halo. The sample comprises 62\,133 halo stars with full phase space coordinates and extends out to heliocentric distances of $\sim 10$ kpc. As actions are conserved under slow changes of the potential, they permit identification of groups of stars with a common accretion history. We devise a method to identify halo substructures based on their clustering in action space, using metallicity as a secondary check. This is validated against smooth models and numerical constructed stellar halos from the Aquarius simulations. We identify 21 substructures in the SDSS-Gaia catalogue, including 7 high significance, high energy and retrograde ones. We investigate whether the retrograde substructures may be material stripped off the atypical globular cluster $\omega$~Centauri. Using a simple model of the accretion of the progenitor of the $\omega$~Centauri, we tentatively argue for the possible association of up to 5 of our new substructures (labelled Rg1, Rg3, Rg4, Rg6 and Rg7) with this event. This sets a minimum mass of $5 \times 10^8 M_\odot$ for the progenitor, so as to bring $\omega$~Centauri to its current location in action -- energy space. Our proposal can be tested by high resolution spectroscopy of the candidates to look for the unusual abundance patterns possessed by $\omega$~Centauri stars.

We present a high resolution synthetic spectral library, INTRIGOSS, designed for studying FGK stars. The library is based on atmosphere models computed with specified individual element abundances via ATLAS12 code. Normalized SPectra (NSP) and surface Flux SPectra (FSP), in the 4830-5400 A, wavelength range, were computed with the SPECTRUM code. INTRIGOSS uses the solar composition by Grevesse et al. 2007 and four [alpha/Fe] abundance ratios and consists of 15,232 spectra. The synthetic spectra are computed with astrophysical gf-values derived by comparing synthetic predictions with a very high SNR solar spectrum and the UVES-U580 spectra of five cool giants. The validity of the NSPs is assessed by using the UVES-U580 spectra of 2212 stars observed in the framework of the Gaia-ESO Survey and characterized by homogeneous and accurate atmospheric parameter values and by detailed chemical compositions. The greater accuracy of NSPs with respect to spectra from the AMBRE, GES_Grid, PHOENIX, C14, and B17 synthetic spectral libraries is demonstrated by evaluating the consistency of the predictions of the different libraries for the UVES-U580 sample stars. The validity of the FSPs is checked by comparing their prediction with both observed spectral energy distribution and spectral indices. The comparison of FSPs with SEDs derived from ELODIE, INDO--U.S., and MILES libraries indicates that the former reproduce the observed flux distributions within a few percent and without any systematic trend. The good agreement between observational and synthetic Lick/SDSS indices shows that the predicted blanketing of FSPs well reproduces the observed one, thus confirming the reliability of INTRIGOSS FSPs.

Several works have found an increase of the abundances of the s-process neutron-capture elements in the youngest Galactic stellar populations, giving important constraints to stellar and Galactic evolution. We aim to trace the abundance patterns and the time-evolution of five s-process elements in the first peak, Y and Zr, and in the second peak, Ba, La and Ce using the Gaia-ESO idr5 results. From the UVES spectra of cluster member stars, we determined the average composition of clusters with ages >0.1 Gyr. We derived statistical ages and distances of field stars, and we separated them in thin and thick disc populations. We studied the time evolution and dependence on metallicity of abundance ratios using open clusters and field stars. Using our large and homogeneous sample of open clusters, thin and thick disc stars, spanning an age range larger than 10 Gyr, we confirm an increase towards young ages of s-process abundances in the Solar neighbourhood. These trends are well defined for open clusters and stars located nearby the solar position and they may be explained by a late enrichment due to significant contribution to the production of these elements from long-living low-mass stars. At the same time, we found a strong dependence of the s-process abundance ratios with the Galactocentric distance and with the metallicity of the clusters and field stars. Our results, derived from the largest and homogeneous sample of s-process abundances in the literature, confirm the growth with decreasing stellar ages of the s-process abundances in both field and open cluster stars. At the same time, taking advantage of the abundances of open clusters located in a wide Galactocentric range, they open a new view on the dependence of the s-process evolution on the metallicity and star formation history, pointing to different behaviours at various Galactocentric distances.

We present results of the first dynamical stream fits to the recently discovered Tucana III stream. These fits assume a fixed Milky Way potential and give proper motion predictions, which can be tested with the upcoming Gaia Data Release 2. These fits reveal that Tucana III is on an eccentric orbit around the Milky Way and, more interestingly, that Tucana III passed within 15 kpc of the Large Magellanic Cloud (LMC) approximately 75 Myr ago. Given this close passage, we fit the Tucana III stream in the combined presence of the Milky Way and the LMC. We find that the predicted proper motions depend on the assumed mass of the LMC and that the LMC can induce a substantial proper motion perpendicular to the stream track. A detection of this misalignment will directly probe the extent of the LMC's influence on our Galaxy, and has implications for nearly all methods which attempt to constraint the Milky Way potential. Such a measurement will be possible with the upcoming Gaia DR2, allowing for a measurement of the LMC's mass.

We use a new panoramic imaging survey, conducted with the Dark Energy Camera, to map the stellar fringes of the Large and Small Magellanic Clouds to extremely low surface brightness V $\gtrsim$ 32 mag arcsec$^{-2}$. Our results starkly illustrate the closely interacting nature of the LMC-SMC pair. We show that the outer LMC disk is strongly distorted, exhibiting an irregular shape, evidence for warping, and significant truncation on the side facing the SMC. Large diffuse stellar substructures are present both to the north and south of the LMC, and in the inter-Cloud region. At least one of these features appears co-spatial with the bridge of RR Lyrae stars that connects the Clouds. The SMC is highly disturbed -- we confirm the presence of tidal tails, as well as a large line-of-sight depth on the side closest to the LMC. Young, intermediate-age, and ancient stellar populations in the SMC exhibit strikingly different spatial distributions. In particular, those with ages $\sim$ 1.5-4 Gyr exhibit a spheroidal distribution with a centroid offset from that of the oldest stars by several degrees towards the LMC. We speculate that the gravitational influence of the LMC may already have been perturbing the gaseous component of the SMC several Gyr ago. With careful modeling, the variety of substructures and tidal distortions evident in the Magellanic periphery should tightly constrain the interaction history of the Clouds.

Investigating the chemical homogeneity of stars born from the same molecular cloud at virtually the same time is very important for our understanding of the chemical enrichment of the interstellar medium and with it the chemical evolution of the Galaxy. One major cause of inhomogeneities in the abundances of open clusters is stellar evolution of the cluster members. In this work, we investigate variations in the surface chemical composition of member stars of the old open cluster M67 as a possible consequence of atomic diffusion effects taking place during the main-sequence phase. The abundances used are obtained from high-resolution UVES/FLAMES spectra within the framework of the Gaia-ESO Survey. We find that the surface abundances of stars on the main sequence decrease with increasing mass reaching a minimum at the turn-off. After deepening of the convective envelope in sub-giant branch stars, the initial surface abundances are restored. We found the measured abundances to be consistent with the predictions of stellar evolutionary models for a cluster with the age and metallicity of M67. Our findings indicate that atomic diffusion poses a non-negligible constraint on the achievable precision of chemical tagging methods.

We report the discovery of a nearby dwarf galaxy in the constellation of Hydrus, between the Large and the Small Magellanic Clouds. Hydrus 1 is a mildy elliptical ultra-faint system with luminosity $M_V\sim$ -4.7 and size $\sim$ 50 pc, located 28 kpc from the Sun and 24 kpc from the LMC. From spectroscopy of $\sim$ 30 member stars, we measure a velocity dispersion of 2.7 km/s and find tentative evidence for a radial velocity gradient consistent with 3 km/s rotation. Hydrus 1's velocity dispersion indicates that the system is dark matter dominated, but its dynamical mass-to-light ratio M/L $\sim$ 66 is significantly smaller than typical for ultra-faint dwarfs at similar luminosity. The kinematics and spatial position of Hydrus~1 make it a very plausible member of the family of satellites brought into the Milky Way by the Magellanic Clouds. While Hydrus 1's proximity and well-measured kinematics make it a promising target for dark matter annihilation searches, we find no evidence for significant gamma-ray emission from Hydrus 1. The new dwarf is a metal-poor galaxy with a mean metallicity [Fe/H]=-2.5 and [Fe/H] spread of 0.4 dex, similar to other systems of similar luminosity. Alpha-abundances of Hyi 1 members indicate that star-formation was extended, lasting between 0.1 and 1 Gyr, with self-enrichment dominated by SN Ia. The dwarf also hosts a highly carbon-enhanced extremely metal-poor star with [Fe/H] $\sim$ -3.2 and [C/Fe] $\sim$ +3.0.

We analyse the structure of the local stellar halo of the Milky Way using $\sim$ 60000 stars with full phase space coordinates extracted from the SDSS--{\it Gaia} catalogue. We display stars in action space as a function of metallicity in a realistic axisymmetric potential for the Milky Way Galaxy. The metal-rich population is more distended towards high radial action $J_R$ as compared to azimuthal or vertical action, $J_\phi$ or $J_z$. It has a mild prograde rotation $(\langle v_\phi \rangle \approx 25$ km s$^{-1}$), is radially anisotropic and highly flattened with axis ratio $q \approx 0.6 - 0.7$. The metal-poor population is more evenly distributed in all three actions. It has larger prograde rotation $(\langle v_\phi \rangle \approx 50$ km s$^{-1}$), a mild radial anisotropy and a roundish morphology ($q\approx 0.9$). We identify two further components of the halo in action space. There is a high energy, retrograde component that is only present in the metal-rich stars. This is suggestive of an origin in a retrograde encounter, possibly the one that created the stripped dwarf galaxy nucleus, $\omega$Centauri. Also visible as a distinct entity in action space is a resonant component, which is flattened and prograde. It extends over a range of metallicities down to [Fe/H] $\approx -3$. It has a net outward radial velocity $\langle v_R \rangle \approx 12$ km s$^{-1}$ within the Solar circle at $|z| <3.5$ kpc. The existence of resonant stars at such extremely low metallicities has not been seen before.

We designed a follow-up program to find the spectroscopic properties of the Hercules-Aquila Cloud (HAC) and test scenarios for its formation. We measured the radial velocities (RVs) of 45 RR Lyrae in the southern portion of the HAC using the facilities at the MDM observatory, producing the first large sample of velocities in the HAC. We found a double-peaked distribution in RVs, skewed slightly to negative velocities. We compared both the morphology of HAC projected onto the plane of the sky and the distribution of velocities in this structure outlined by RR Lyrae and other tracer populations at different distances to N-body simulations. We found that the behaviour is characteristic of an old, well-mixed accretion event with small apo-galactic radius. We cannot yet rule out other formation mechanisms for the HAC. However, if our interpretation is correct, HAC represents just a small portion of a much larger debris structure spread throughout the inner Galaxy whose distinct kinematic structure should be apparent in RV studies along many lines of sight.

Context. The origin and dynamical evolution of star clusters is an important topic in stellar astrophysics. Several models have been proposed to understand the formation of bound and unbound clusters and their evolution, and these can be tested by examining the kinematical and dynamical properties of clusters over a wide range of ages and masses. Aims. We use the Gaia-ESO Survey products to study four open clusters (IC 2602, IC 2391, IC 4665, and NGC 2547) that lie in the age range between 20 and 50 Myr. Methods. We employ the gravity index $\gamma$ and the equivalent width of the lithium line at 6708 $\AA$, together with effective temperature $\rm{T_{eff}}$, and the metallicity of the stars in order to discard observed contaminant stars. Then, we derive the cluster radial velocity dispersions $\sigma_c$, the total cluster mass $\rm{M}_{tot}$, and the half mass radius $r_{hm}$. Using the $Gaia$-DR1 TGAS catalogue, we independently derive the intrinsic velocity dispersion of the clusters from the astrometric parameters of cluster members. Results. The intrinsic radial velocity dispersions derived by the spectroscopic data are larger than those derived from the TGAS data, possibly due to the different masses of the considered stars. Using $\rm{M}_{tot}$ and $r_{hm}$ we derive the virial velocity dispersion $\sigma_{vir}$ and we find that three out of four clusters are supervirial. This result is in agreement with the hypothesis that these clusters are dispersing, as predicted by the "residual gas expulsion" scenario. However, recent simulations show that the virial ratio of young star clusters may be overestimated if it is determined using the global velocity dispersion, since the clusters are not fully relaxed.

Using a large sample of Main Sequence stars with 7-D measurements supplied by Gaia and SDSS, we study the kinematic properties of the local (within ~10 kpc from the Sun) stellar halo. We demonstrate that the halo's velocity ellipsoid evolves strongly with metallicity. At the low [Fe/H] end, the orbital anisotropy (the amount of motion in the radial direction compared to the tangential one) is mildly radial with 0.2<beta<0.4. However, for stars with [Fe/H]>-1.7 we measure extreme values of beta~0.9. Across the metallicity range considered, i.e. -3<[Fe/H]<-1, the stellar halo's spin is minimal, at the level of 20<v_theta (km/s) <30. Using a suite of cosmological zoom-in simulations of halo formation, we deduce that the observed acute anisotropy is inconsistent with the continuous accretion of dwarf satellites. Instead, we argue, the stellar debris in the inner halo were deposited in a major accretion event by a satellite with Mvir>10^10 Msun around the epoch of the Galactic disc formation, i.e. between 8 and 11 Gyr ago. The radical halo anisotropy is the result of the dramatic radialisation of the massive progenitor's orbit, amplified by the action of the growing disc.

Using deep photometric data from CFHT/Megacam, we study the morphology and density of the GD-1 stream, one of the longest and coldest stellar streams in the Milky Way. Our deep data recovers the lower main sequence of the stream with unprecedented quality, clearly separating it from Milky Way foreground and background stars. An analysis of the distance to different parts of the stream shows that GD-1 lies at a heliocentric distance between 8 and 10 kpc, with only a shallow gradient across 45 deg on the sky. Matched filter maps of the stream density show clear density variations, such as deviations from a single orbital track and tentative evidence for stream fanning. We also detect a clear under-density in the middle of the stream track at $\varphi_{1}$=-45 deg surrounded by overdense stream segments on either side. This location is a promising candidate for the elusive missing progenitor of the GD-1 stream. We conclude that the GD-1 stream has clearly been disturbed by interactions with the Milky Way disk or other sub-halos.

We report the discovery of two ultra-faint satellites in the vicinity of the Large Magellanic Cloud (LMC) in data from the Magellanic Satellites Survey (MagLiteS). Situated 18$^{\circ}$ ($\sim 20$ kpc) from the LMC and separated from each other by only $18^\prime$, Carina~II and III form an intriguing pair. By simultaneously modeling the spatial and the color-magnitude stellar distributions, we find that both Carina~II and Carina~III are likely dwarf galaxies, although this is less clear for Carina~III. There are in fact several obvious differences between the two satellites. While both are well described by an old and metal poor population, Carina~II is located at $\sim 36$ kpc from the Sun, with $M_V\sim-4.5$ and $r_h\sim 90$ pc, and it is further confirmed by the discovery of 3 RR Lyrae at the right distance. In contrast, Carina~III is much more elongated, measured to be fainter ($M_V\sim-2.4$), significantly more compact ($r_h\sim30$ pc), and closer to the Sun, at $\sim 28$ kpc, placing it only 8 kpc away from Car~II. Together with several other systems detected by the Dark Energy Camera, Carina~II and III form a strongly anisotropic cloud of satellites in the vicinity of the Magellanic Clouds.

We use the SDSS-Gaia Catalogue to identify six new pieces of halo substructure. SDSS-Gaia is an astrometric catalogue that exploits SDSS data release 9 to provide first epoch photometry for objects in the Gaia source catalogue. We use a version of the catalogue containing $245\,316$ stars with all phase space coordinates within a heliocentric distance of $\sim 10$ kpc. We devise a method to assess the significance of halo substructures based on their clustering in velocity space. The two most substantial structures are multiple wraps of a stream which has undergone considerable phase mixing (S1, with 94 members) and a kinematically cold stream (S2, with 61 members). The member stars of S1 have a median position of ($X,Y,Z$) = ($8.12, -0.22, 2.75$) kpc and a median metallicity of [Fe/H] $= -1.78$. The stars of S2 have median coordinates ($X,Y,Z$) = ($8.66, 0.30, 0.77$) kpc and a median metallicity of [Fe/H] $= -1.91$. They lie in velocity space close to some of the stars in the stream reported by Helmi et al. (1999). By modelling, we estimate that both structures had progenitors with virial masses $\approx 10^{10} M_\odot$ and infall times $\gtrsim 9$ Gyr ago. Using abundance matching, these correspond to stellar masses between $10^6$ and $10^7 M_\odot$. These are somewhat larger than the masses inferred through the mass-metallicity relation by factors of 5 to 15. Additionally, we identify two further substructures (S3 and S4 with 55 and 40 members) and two clusters or moving groups (C1 and C2 with 24 and 12) members. In all 6 cases, clustering in kinematics is found to correspond to clustering in both configuration space and metallicity, adding credence to the reliability of our detections.

We exploit the exquisite, deep Hyper Suprime-Cam (HSC) imaging data to probe the Galactic halo out to 200 kpc. Using the ~100 square degree, multi-band photometry of the first HSC Wide survey data release, we identify blue horizontal branch (BHB) stars beyond 50 kpc in the halo. The presence of the Sagittarius (Sgr) stream in the HSC fields produces a notable excess of stars at the apocentre of the leading arm (~50-60 kpc). For fields excluding Sgr, the BHB counts are consistent with a continuation of a -4 power-law from the inner halo. However, we find that the majority of the non-Sgr BHB stars beyond 50 kpc reside in one 27 square degree HSC field called "VVDS". Curiously, this field is located close to the Magellanic plane, and we hypothesize that the excess of stars between 50 and 200 kpc could be associated with distant Magellanic debris. Indeed, without the VVDS, there are very few BHBs in the remaining portions of the Galaxy probed by the HSC. Accordingly, this scarcity of tracers is consistent with a significant decline in stellar density beyond 50 kpc, with a power-law of -4 or steeper.

In this paper we describe the routine photometric calibration of data taken with the VIRCAM instrument on the ESO VISTA telescope. The broadband ZYJHKs data are directly calibrated from 2MASS point sources visible in every VISTA image. We present the empirical transformations between the 2MASS and VISTA, and WFCAM and VISTA, photometric systems for regions of low reddening. We investigate the long-term performance of VISTA+VIRCAM. An investigation of the dependence of the photometric calibration on interstellar reddening leads to these conclusions: (1) For all broadband filters, a linear colour-dependent correction compensates the gross effects of reddening where $E(B-V)<5.0$. (2) For $Z$ and $Y$, there is a significantly larger scatter above E(B-V)=5.0, and insufficient measurements to adequately constrain the relation beyond this value. (3) The $JHK\!s$ filters can be corrected to a few percent up to E(B-V)=10.0. We analyse spatial systematics over month-long timescales, both inter- and intra-detector and show that these are present only at very low levels in VISTA. We monitor and remove residual detector-to-detector offsets. We compare the calibration of the main pipeline products: pawprints and tiles. We show how variable seeing and transparency affect the final calibration accuracy of VISTA tiles, and discuss a technique, {\it grouting}, for mitigating these effects. Comparison between repeated reference fields is used to demonstrate that the VISTA photometry is precise to better than $\simeq2\%$ for the $Y$$J$$H$$Ks$ bands and $3\%$ for the $Z$ bands. Finally we present empirically determined offsets to transform VISTA magnitudes into a true Vega system.

We describe the methodologies that, taking advantage of Gaia-DR1 and the Gaia-ESO Survey data, enable the comparison of observed open star cluster sequences with stellar evolutionary models. The final, long-term goal is the exploitation of open clusters as age calibrators. We perform a homogeneous analysis of eight open clusters using the Gaia-DR1 TGAS catalogue for bright members, and information from the Gaia-ESO Survey for fainter stars. Cluster membership probabilities for the Gaia-ESO Survey targets are derived based on several spectroscopic tracers. The Gaia-ESO Survey also provides the cluster chemical composition. We obtain cluster parallaxes using two methods. The first one relies on the astrometric selection of a sample of bona fide members, while the other one fits the parallax distribution of a larger sample of TGAS sources. Ages and reddening values are recovered through a Bayesian analysis using the 2MASS magnitudes and three sets of standard models. Lithium depletion boundary (LDB) ages are also determined using literature observations and the same models employed for the Bayesian analysis. For all but one cluster, parallaxes derived by us agree with those presented in Gaia Collaboration et al. (2017), while a discrepancy is found for NGC 2516; we provide evidence supporting our own determination. Inferred cluster ages are robust against models and are generally consistent with literature values. The systematic parallax errors inherent in the Gaia DR1 data presently limit the precision of our results. Nevertheless, we have been able to place these eight clusters onto the same age scale for the first time, with good agreement between isochronal and LDB ages where there is overlap. Our approach appears promising and demonstrates the potential of combining Gaia and ground-based spectroscopic datasets.

We show that tagging RR Lyrae stars according to their location in the period-amplitude diagram can be used to shed light on the genesis of the Galactic stellar halo. The mixture of RR Lyrae of ab type, separated into classes along the lines suggested by Oosterhoff, displays a strong and coherent evolution with Galactocentric radius. The change in the RR Lyrae composition appears to coincide with the break in the halo's radial density profile at ~25 kpc. Using simple models of the stellar halo, we establish that at least three different types of accretion events are necessary to explain the observed RRab behavior. Given that there exists a correlation between the RRab class fraction and the total stellar content of a dwarf satellite, we hypothesize that the field halo RRab composition is controlled by the mass of the progenitor contributing the bulk of the stellar debris at the given radius. This idea is tested against a suite of cosmological zoom-in simulations of Milky Way-like stellar halo formation. Finally, we study some of the most prominent stellar streams in the Milky Way halo and demonstrate that their RRab class fractions follow the trends established previously.

We attempt to determine the relative fraction of stars that have undergone significant radial migration by studying the orbital properties of metal-rich ([Fe/H]$>0.1$) stars within 2 kpc of the Sun using a sample of more than 3,000 stars selected from iDR4 of the Gaia-ESO Survey. We investigate the kinematic properties, such as velocity dispersion and orbital parameters, of stellar populations near the sun as a function of [Mg/Fe] and [Fe/H], which could show evidence of a major merger in the past history of the Milky Way. This was done using the stellar parameters from the Gaia-ESO Survey along with proper motions from PPMXL to determine distances, kinematics, and orbital properties for these stars to analyze the chemodynamic properties of stellar populations near the Sun. Analyzing the kinematics of the most metal-rich stars ([Fe/H]$>0.1$), we find that more than half have small eccentricities ($e<0.2$) or are on nearly circular orbits. Slightly more than 20\% of the metal-rich stars have perigalacticons $R_p>7$ kpc. We find that the highest [Mg/Fe], metal-poor populations have lower vertical and radial velocity dispersions compared to lower [Mg/Fe] populations of similar metallicity by $\sim10$ km s$^{-1}$. The median eccentricity increases linearly with [Mg/Fe] across all metallicities, while the perigalacticon decreases with increasing [Mg/Fe] for all metallicities. Finally, the most [Mg/Fe]-rich stars are found to have significant asymmetric drift and rotate more than 40 km s$^{-1}$ slower than stars with lower [Mg/Fe] ratios. While our results cannot constrain how far stars have migrated, we propose that migration processes are likely to have played an important role in the evolution of the Milky Way, with metal-rich stars migrating from the inner disk toward to solar neighborhood and past mergers potentially driving enhanced migration of older stellar populations in the disk.

We exploit the first \gaia data release to study the properties of the Galactic stellar halo as traced by RR Lyrae. We demonstrate that it is possible to select a pure sample of RR Lyrae using only photometric information available in the Gaia+2MASS catalogue. The final sample contains about 21600 RR Lyrae covering an unprecedented fraction ($\sim60\%$) of the volume of the Galactic inner halo ($\text{R}<28$ kpc). We study the morphology of the stellar halo by analysing the RR Lyrae distribution with parametric and non-parametric techniques. Taking advantage of the uniform all-sky coverage, we test halo models more sophisticated than usually considered in the literature, such as those with varying flattening, tilt and/or offset of the halo with respect to the Galactic disc. A consistent picture emerges: the inner halo is well reproduced by a smooth distribution of stars settled on triaxial ellipsoids. The minor axis is perpendicular to the Milky Ways disc, while the major axis is misaligned by $\sim20^{\circ}$ from the Galactic Y axis. The elongation along the major axis is mild ($\text{p}=1.27$), and the vertical flattening is shown to evolve from a squashed state with $\text{q}\approx0.57$ in the centre to a more spherical $\text{q}\approx0.75$ at the outer edge of our dataset. The density slope is well approximated by a single power-law with exponent $\alpha=-2.96$. Within the range probed, we see no significant evidence for a change of the radial density slope, out of the plane tilt or an offset of the halo with respect to the Galaxy's centre.

We report the discovery, spectroscopic confirmation, and mass modelling of the gravitationally lensed quasar system PS J0630-1201. The lens was discovered by matching a photometric quasar catalogue compiled from Pan-STARRS and WISE photometry to the Gaia DR1 catalogue, exploiting the high spatial resolution of the latter (FWHM $\sim $0.1") to identify the three brightest components of the lens. Follow-up spectroscopic observations with the WHT confirm the multiple objects are quasars at redshift $z_{q}=3.34$. Further follow-up with Keck AO high-resolution imaging reveals that the system is composed of two lensing galaxies and the quasar is lensed into a $\sim$2.8" separation four-image cusp configuration with a fifth image clearly visible, and a 1.0" arc due to the lensed quasar host galaxy. The system is well-modelled with two singular isothermal ellipsoids, reproducing the position of the fifth image. We discuss future prospects for measuring time delays between the images and constraining any offset between mass and light using the faintly detected Einstein arcs associated with the quasar host galaxy.

Gaia's exceptional resolution (FWHM $\sim$ 0.1$^{\prime\prime}$) allows identification and cataloguing of the multiple images of gravitationally lensed quasars. We investigate a sample of 49 known lensed quasars in the SDSS footprint, with image separations less than 2$^{\prime\prime}$, and find that 8 are detected with multiple components in the first Gaia data release. In the case of the 41 single Gaia detections, we generally are able to distinguish these lensed quasars from single quasars when comparing Gaia flux and position measurements to those of Pan-STARRS and SDSS. This is because the multiple images of these lensed quasars are typically blended in ground-based imaging and therefore the total flux and a flux-weighted centroid are measured, which can differ significantly from the fluxes and centroids of the individual components detected by Gaia. We compare the fluxes through an empirical fit of Pan-STARRS griz photometry to the wide optical Gaia bandpass values using a sample of isolated quasars. The positional offsets are calculated from a recalibrated astrometric SDSS catalogue. Applying flux and centroid difference criteria to spectroscopically confirmed quasars, we discover 4 new sub-arcsecond-separation lensed quasar candidates which have two distinct components of similar colour in archival CFHT or HSC data. Our method based on single Gaia detections can be used to identify the $\sim$ 1400 lensed quasars with image separation above 0.5$^{\prime\prime}$, expected to have only one image bright enough to be detected by Gaia.

We study the kinematic properties of the Eastern Banded Structure (EBS) and Hydra I overdensity using exquisite proper motions derived from the Sloan Digital Sky Survey (SDSS) and Gaia source catalog. Main sequence turn-off stars in the vicinity of the EBS are identified from SDSS photometry; we use the proper motions and, where applicable, spectroscopic measurements of these stars to probe the kinematics of this apparent stream. We find that the EBS and Hydra I share common kinematic and chemical properties with the nearby Monoceros Ring. In particular, the proper motions of the EBS, like Monoceros, are indicative of prograde rotation (V_phi ~ 180-220 km/s), which is similar to the Galactic thick disc. The kinematic structure of stars in the vicinity of the EBS suggest that it is not a distinct stellar stream, but rather marks the "edge" of the Monoceros Ring. The EBS and Hydra I are the latest substructures to be linked with Monoceros, leaving the Galactic anti-centre a mess of interlinked overdensities which likely share a unified, Galactic disc origin.

The typical methodology for comparing simulated galaxies with observational surveys is usually to apply a spatial selection to the simulation to mimic the region of interest covered by a comparable observational survey sample. In this work we compare this approach with a more sophisticated post-processing in which the observational uncertainties and selection effects (photometric, surface gravity and effective temperature) are taken into account. We compare a `solar neighbourhood analogue' region in a model Milky Way-like galaxy simulated with RAMSES-CH with fourth release Gaia-ESO survey data. We find that a simple spatial cut alone is insufficient and that observational uncertainties must be accounted for in the comparison. This is particularly true when the scale of uncertainty is large compared to the dynamic range of the data, e.g. in our comparison, the [Mg/Fe] distribution is affected much more than the more accurately determined [Fe/H] distribution. Despite clear differences in the underlying distributions of elemental abundances between simulation and observation, incorporating scatter to our simulation results to mimic observational uncertainty produces reasonable agreement. The quite complete nature of the Gaia-ESO survey means that the selection function has minimal impact on the distribution of observed age and metal abundances but this would become increasingly more important for surveys with narrower selection functions.

We use a new contiguous imaging survey conducted using the Dark Energy Camera to investigate the distribution and properties of young stellar populations in the Magellanic inter-Cloud region. These young stars are strongly spatially clustered, forming a narrow chain of low-mass associations that trace the densest HI gas in the Magellanic Bridge and extend, in projection, from the SMC to the outer disk of the LMC. The associations in our survey footprint have ages $\lesssim 30$ Myr, masses in the range $\sim 100-1200\,{\rm M}_\odot$, and very diffuse structures with half-light radii of up to $\sim 100$ pc. The two most populous are strongly elliptical, and aligned to $\approx 10^{{\rm o}}$ with the axis joining the centres of the LMC and SMC. These observations strongly suggest that the young inter-Cloud populations formed in situ, likely due to the compression of gas stripped during the most recent close LMC-SMC encounter. The associations lie at distances intermediate between the two Clouds, and we find no evidence for a substantial distance gradient across the imaged area. Finally, we identify a vast shell of young stars surrounding a central association, that is spatially coincident with a low column density bubble in the HI distribution. The properties of this structure are consistent with a scenario where stellar winds and supernova explosions from massive stars in the central cluster swept up the ambient gas into a shell, triggering a new burst of star formation. This is a prime location for studying stellar feedback in a relatively isolated environment.

We study the structure of the inner Milky Way using the latest data release of the Vista Variables in Via Lactea (VVV) survey. The VVV is a deep near-infrared, multi-colour photometric survey with a coverage of 300 square degrees towards the Bulge/Bar. We use Red Clump (RC) stars to produce a high-resolution dust map of the VVV's field of view. From de-reddened colour-magnitude diagrams we select Red Giant Branch stars to investigate their 3D density distribution within the central 4 kpc. We demonstrate that our best-fit parametric model of the Bulge density provides a good description of the VVV data, with a median percentage residual of 5$\%$ over the fitted region. The strongest of the otherwise low-level residuals are overdensities associated with a low-latitude structure as well as the so-called X-shape previously identified using the split RC. These additional components contribute only $\sim5\%$ and $\sim7\%$ respectively to the Bulge mass budget. The best-fit Bulge is `boxy' with an axis ratio of [1:0.44:0.31] and is rotated with respect to the Sun-Galactic Centre line by at least $20^{\circ}$. We provide an estimate of the total, full sky, mass of the Bulge of $M_\mathrm{Bulge}^{\mathrm{Chabrier}} = 2.36 \times 10^{10} M_{\odot}$ for a Chabrier initial mass function. We show there exists a strong degeneracy between the viewing angle and the dispersion of the Red Clump absolute magnitude distribution. The value of the latter is strongly dependent on the assumptions made about the intrinsic luminosity function of the Bulge.

The Gaia-ESO Survey (GES) is a large spectroscopic survey that provides a unique opportunity to study the distribution of spectroscopic multiple systems among different populations of the Galaxy. We aim at detecting binarity/multiplicity for stars targeted by the GES from the analysis of the cross-correlation functions (CCFs) of the GES spectra with spectral templates. We develop a method based on the computation of the CCF successive derivatives to detect multiple peaks and determine their radial velocities, even when the peaks are strongly blended. The parameters of the detection of extrema (DOE) code have been optimized for each GES GIRAFFE and UVES setup to maximize detection. This code therefore allows to automatically detect multiple line spectroscopic binaries (SBn, n>1). We apply this method on the fourth GES internal data release and detect 354 SBn candidates (342 SB2, 11 SB3 and even one SB4), including only 9 SB2 known in the literature. This implies that about 98% of these SBn candidates are new (because of their faint visual magnitude that can reach V=19). Visual inspection of the SBn candidate spectra reveals that the most probable candidates have indeed a composite spectrum. Among SB2 candidates, an orbital solution could be computed for two previously unknown binaries: 06404608+0949173 (known as V642 Mon) in NGC 2264 and 19013257-0027338 in Berkeley 81. A detailed analysis of the unique SB4 (four peaks in the CCF) reveals that HD 74438 in the open cluster IC 2391 is a physically bound stellar quadruple system. The SB candidates belonging to stellar clusters are reviewed in detail to discard false detections. We warn against the use of atmospheric parameters for these system components rather than by SB-specific pipelines. Our implementation of an automatic detection of spectroscopic binaries within the GES has allowed an efficient discovery of many new multiple systems.

We present the first detailed study of the behaviour of the stellar proper motion across the entire Galactic Anti-centre area visible in the Sloan Digital Sky Survey data. We use recalibrated SDSS astrometry in combination with positions from {\it Gaia} DR1 to provide tangential motion measurements with a systematic uncertainty $<$5 kms$^{-1}$ for the Main Sequence stars at the distance of the Monoceros Ring. We demonstrate that Monoceros members rotate around the Galaxy with azimuthal speeds of ~230 kms$^{-1}$, only slightly lower than that of the Sun. Additionally, both vertical and azimuthal components of their motion are shown to vary considerably but gradually as a function of Galactic longitude and latitude. The stellar over-density in the Anti-centre region can be split into two components, the narrow, stream-like ACS and the smooth Ring. According to our analysis, these two structures show very similar but clearly distinct kinematic trends, which can be summarised as follows: the amplitude of the velocity variation in $v_{\phi}$ and $v_z$ in the ACS is higher compared to the Ring, whose velocity gradients appear to be flatter. Currently, no model available can explain the entirety of the data in this area of the sky. However, the new accurate kinematic map introduced here should provide strong constraints on the genesis of the Monoceros Ring and the associated sub-structure.

Only in the Milky Way is it possible to conduct an experiment which uses stellar streams to detect low-mass dark matter subhaloes. In smooth and static host potentials, tidal tails of disrupting satellites appear highly symmetric. However, perturbations from dark subhaloes, as well as from GMCs and the Milky Way bar, can induce density fluctuations that destroy this symmetry. Motivated by the recent release of unprecedentedly deep and wide imaging data around the Pal~5 stellar stream, we develop a new probabilistic, adaptive and non-parametric technique which allows us to bring the cluster's tidal tails into clear focus. Strikingly, we uncover a stream whose density exhibits visible changes on a variety of angular scales. We detect significant bumps and dips, both narrow and broad: two peaks on either side of the progenitor, each only a fraction of a degree across, and two gaps, $\sim2^{\circ}$ and $\sim9^{\circ}$ wide, the latter accompanied by a gargantuan lump of debris. This largest density feature results in a pronounced inter-tail asymmetry which cannot be made consistent with an unperturbed stream according to a suite of simulations we have produced. We conjecture that the sharp peaks around Pal 5 are epicyclic overdensities, while the two dips are consistent with impacts by subhaloes. Assuming an age of 3.4 Gyr for Pal 5, these two gaps would correspond to the characteristic size of gaps created by subhaloes in the mass range of $10^6-10^7 M_\odot$ and $10^7-10^8 M_\odot$ respectively. In addition to dark substructure, we find that the bar of the Milky Way can plausibly produce the asymmetric density seen in Pal 5 and that GMCs could cause the smaller gap.

We combine Gaia data release 1 astrometry with Sloan Digital Sky Survey (SDSS) images taken some ~10-15 years earlier, to measure proper motions of stars in the halo of our Galaxy. The SDSS-Gaia proper motions have typical statistical errors of 2 mas/yr down to r ~ 20 mag, and are robust to variations with magnitude and colour. Armed with this exquisite set of halo proper motions, we identify RR Lyrae, blue horizontal branch (BHB), and K giant stars in the halo, and measure their net rotation with respect to the Galactic disc. We find evidence for a gently rotating prograde signal (< $V_\phi$ > ~ 5-25 km/s) in the halo stars, which shows little variation with Galactocentric radius out to 50 kpc. The average rotation signal for the three populations is < $V_\phi$ > = 14 +/- 2 +/- 10 (syst.) km/s. There is also tentative evidence for a kinematic correlation with metallicity, whereby the metal richer BHB and K giant stars have slightly stronger prograde rotation than the metal poorer stars. Using the Auriga simulation suite we find that the old (T >10 Gyr) stars in the simulated halos exhibit mild prograde rotation, with little dependence on radius or metallicity, in general agreement with the observations. The weak halo rotation suggests that the Milky Way has a minor in situ halo component, and has undergone a relatively quiet accretion history.

We present the results of the very first search for faint Milky Way satellites in the Gaia data. Using stellar positions only, we are able to re-discover objects detected in much deeper data as recently as the last couple of years. While we do not identify new prominent ultra-faint dwarf galaxies, we report the discovery of two new star clusters, Gaia 1 and Gaia 2. Gaia 1 is particularly curious, as it is a massive (2.2$\times$10$^4$ M$_\odot$), large ($\sim$9 pc) and nearby (4.6 kpc) cluster, situated 10' away from the brightest star on the sky, Sirius! Even though this satellite is detected at significance in excess of 10, it was missed by previous sky surveys. We conclude that Gaia possesses powerful and unique capabilities for satellite detection thanks to its unrivalled angular resolution and highly efficient object classification.

We report the discovery of tidal tails around the Galactic globular cluster NGC 7492, based on the Data Release 1 of the Pan-STARRS 1 survey. The tails were detected with a version of the matched filter technique applied to the $(g-r,r)$ and $(g-i,i)$ color-magnitude diagrams. Tidal tails emerging from the cluster extend at least $\sim$3.5 degrees in the North-East to South-East direction, equivalent to $\sim1.5$ kpc in projected length.

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.

We identify a halo substructure in the Tycho Gaia Astrometric Solution (TGAS) dataset, cross-matched with the RAVE-on data release. After quality cuts, the stars with large radial action ($J_R > 800$ kms$^{-1}$ kpc) are extracted. A subset of these stars is clustered in longitude and velocity and can be selected with further cuts. The 14 stars are centered on $(X,Y,Z) \approx (9.0,-1.0,-0.6)$ kpc and form a coherently moving structure in the halo with median $(v_R,v_\phi,v_z) = (167.33,0.86,-94.85)$ kms$^{-1}$. They are all metal-poor giants with median [Fe/H] $=-0.83$. To guard against the effects of distance errors, we compute spectrophotometric distances for the 8 out of 14 stars where this is possible. We find that 6 of the stars are still comoving. These 6 stars also have a much tighter [Fe/H] distribution $\sim -0.7$ with one exception ([Fe/H] = -2.12). We conclude that the existence of the comoving cluster is stable against changes in distance estimation and conjecture that this is the dissolving remnant of a long-ago accreted globular cluster.

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.

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.

The spatial distribution of elemental abundances in the disc of our Galaxy gives insights both on its assembly process and subsequent evolution, and on the stellar nucleogenesis of the different elements. Gradients can be traced using several types of objects as, for instance, (young and old) stars, open clusters, HII regions, planetary nebulae. We aim at tracing the radial distributions of abundances of elements produced through different nucleosynthetic channels -the alpha-elements O, Mg, Si, Ca and Ti, and the iron-peak elements Fe, Cr, Ni and Sc - by using the Gaia-ESO idr4 results of open clusters and young field stars. From the UVES spectra of member stars, we determine the average composition of clusters with ages >0.1 Gyr. We derive statistical ages and distances of field stars. We trace the abundance gradients using the cluster and field populations and we compare them with a chemo-dynamical Galactic evolutionary model. Results. The adopted chemo-dynamical model, with the new generation of metallicity-dependent stellar yields for massive stars, is able to reproduce the observed spatial distributions of abundance ratios, in particular the abundance ratios of [O/Fe] and [Mg/Fe] in the inner disc (5 kpc<RGC <7 kpc), with their differences, that were usually poorly explained by chemical evolution models. Often, oxygen and magnesium are considered as equivalent in tracing alpha-element abundances and in deducing, e.g., the formation time-scales of different Galactic stellar populations. In addition, often [alpha/Fe] is computed combining several alpha-elements. Our results indicate, as expected, a complex and diverse nucleosynthesis of the various alpha-elements, in particular in the high metallicity regimes, pointing towards a different origin of these elements and highlighting the risk of considering them as a single class with common features.

We present the first extensive spectroscopic study of the global population in star clusters Trumpler~16, Trumpler~14 and Collinder~232 in the Carina Nebula, using data from the Gaia-ESO Survey, down to solar-mass stars. In addition to the standard homogeneous Survey data reduction, a special processing was applied here because of the bright nebulosity surrounding Carina stars. We find about four hundred good candidate members ranging from OB types down to slightly sub-solar masses. About one-hundred heavily-reddened early-type Carina members found here were previously unrecognized or poorly classified, including two candidate O stars and several candidate Herbig Ae/Be stars. Their large brightness makes them useful tracers of the obscured Carina population. The spectroscopically-derived temperatures for nearly 300 low-mass members allows the inference of individual extinction values, and the study of the relative placement of stars along the line of sight. We find a complex spatial structure, with definite clustering of low-mass members around the most massive stars, and spatially-variable extinction. By combining the new data with existing X-ray data we obtain a more complete picture of the three-dimensional spatial structure of the Carina clusters, and of their connection to bright and dark nebulosity, and UV sources. The identification of tens of background giants enables us also to determine the total optical depth of the Carina nebula along many sightlines. We are also able to put constraints on the star-formation history of the region, with Trumpler~14 stars found to be systematically younger than stars in other sub-clusters. We find a large percentage of fast-rotating stars among Carina solar-mass members, which provide new constraints on the rotational evolution of pre-main-sequence stars in this mass range.

The radial metallicity distribution in the Galactic thin disc represents a crucial constraint for modelling disc formation and evolution. Open clusters allow us to derive both the radial metallicity distribution and its evolution over time. In this paper we perform the first investigation of the present-day radial metallicity distribution based on [Fe/H] determinations in late type members of pre-main-sequence clusters. Because of their youth, these clusters are therefore essential for tracing the current inter-stellar medium metallicity. We used the products of the Gaia-ESO Survey analysis of 12 young regions (age<100 Myr), covering Galactocentric distances from 6.67 to 8.70 kpc. For the first time, we derived the metal content of star forming regions farther than 500 pc from the Sun. Median metallicities were determined through samples of reliable cluster members. For ten clusters the membership analysis is discussed in the present paper, while for other two clusters (Chamaeleon I and Gamma Velorum) we adopted the members identified in our previous works. All the pre-main-sequence clusters considered in this paper have close-to-solar or slightly sub-solar metallicities. The radial metallicity distribution traced by these clusters is almost flat, with the innermost star forming regions having [Fe/H] values that are 0.10-0.15 dex lower than the majority of the older clusters located at similar Galactocentric radii. This homogeneous study of the present-day radial metallicity distribution in the Galactic thin disc favours models that predict a flattening of the radial gradient over time. On the other hand, the decrease of the average [Fe/H] at young ages is not easily explained by the models. Our results reveal a complex interplay of several processes (e.g. star formation activity, initial mass function, supernova yields, gas flows) that controlled the recent evolution of the Milky Way.

Milky Way open clusters are very diverse in terms of age, chemical composition, and kinematic properties. Intermediate-age and old open clusters are less common, and it is even harder to find them inside the solar Galactocentric radius, due to the high mortality rate and strong extinction inside this region. NGC 6802 is one of the inner disk open clusters (IOCs) observed by the $Gaia$-ESO survey (GES). This cluster is an important target for calibrating the abundances derived in the survey due to the kinematic and chemical homogeneity of the members in open clusters. Using the measurements from $Gaia$-ESO internal data release 4 (iDR4), we identify 95 main-sequence dwarfs as cluster members from the GIRAFFE target list, and eight giants as cluster members from the UVES target list. The dwarf cluster members have a median radial velocity of $13.6\pm1.9$ km s$^{-1}$, while the giant cluster members have a median radial velocity of $12.0\pm0.9$ km s$^{-1}$ and a median [Fe/H] of $0.10\pm0.02$ dex. The color-magnitude diagram of these cluster members suggests an age of $0.9\pm0.1$ Gyr, with $(m-M)_0=11.4$ and $E(B-V)=0.86$. We perform the first detailed chemical abundance analysis of NGC 6802, including 27 elemental species. To gain a more general picture about IOCs, the measurements of NGC 6802 are compared with those of other IOCs previously studied by GES, that is, NGC 4815, Trumpler 20, NGC 6705, and Berkeley 81. NGC 6802 shows similar C, N, Na, and Al abundances as other IOCs. These elements are compared with nucleosynthetic models as a function of cluster turn-off mass. The $\alpha$, iron-peak, and neutron-capture elements are also explored in a self-consistent way.

We exploit the first data release from the Gaia mission to identify candidate Mira variables in the outskirts of the Magellanic Clouds. The repeated observations of sources during the initial phase of the Gaia mission is used to identify stars that show signs of variability. This variability information, combined with infrared photometry from 2MASS and WISE, allows us to select a clean sample of giants in the periphery of the LMC. We find evidence for Miras surrounding the LMC out to ~20 deg in all directions, apart from the North-West quadrant. Our sample does not generally follow the gas distribution of the Magellanic system; Miras are notably absent in the gaseous bridge between the LMC and SMC, but they are likely related to the stellar RR Lyrae bridge reported by Belokurov et al. (2016). The stellar stream discovered by Mackey et al. (2016) to the North of the LMC is almost perfectly delineated by our Mira variables, and likely extends further East toward the Galactic plane. The presence of an intermediate-age population in this stream advocates an LMC disc origin. We also find a significant excess of Miras to the East of the LMC; these more diffusely distributed stars are likely stripped SMC stars due to interactions with the LMC. Miras are also identified in regions of the sky away from the Clouds; we locate stars likely associated with known massive substructures, and also find potential associations with stripped SMC debris above the Galactic plane.

We present the discovery and spectroscopic confirmation with the ESO NTT and Gemini South telescopes of eight new 6.0 < z < 6.5 quasars with z$_{AB}$ < 21.0. These quasars were photometrically selected without any star-galaxy morphological criteria from 1533 deg$^{2}$ using SED model fitting to photometric data from the Dark Energy Survey (g, r, i, z, Y), the VISTA Hemisphere Survey (J, H, K) and the Wide-Field Infrared Survey Explorer (W1, W2). The photometric data was fitted with a grid of quasar model SEDs with redshift dependent Lyman-{\alpha} forest absorption and a range of intrinsic reddening as well as a series of low mass cool star models. Candidates were ranked using on a SED-model based $\chi^{2}$-statistic, which is extendable to other future imaging surveys (e.g. LSST, Euclid). Our spectral confirmation success rate is 100% without the need for follow-up photometric observations as used in other studies of this type. Combined with automatic removal of the main types of non-astrophysical contaminants the method allows large data sets to be processed without human intervention and without being over run by spurious false candidates. We also present a robust parametric redshift estimating technique that gives comparable accuracy to MgII and CO based redshift estimators. We find two z $\sim$ 6.2 quasars with HII near zone sizes < 3 proper Mpc which could indicate that these quasars may be young with ages < 10$^6$ - 10$^7$ years or lie in over dense regions of the IGM. The z = 6.5 quasar VDESJ0224-4711 has J$_{AB}$ = 19.75 is the second most luminous quasar known with z > 6.5.

The young (~2 Myr) cluster Chamaeleon I is one of the closest laboratories to study the early stages of star cluster dynamics in a low-density environment. We studied its structural and kinematical properties combining parameters from the high-resolution spectroscopic survey Gaia-ESO with data from the literature. Our main result is the evidence of a large discrepancy between the velocity dispersion (sigma = 1.14 \pm 0.35 km s^{-1}) of the stellar population and the dispersion of the pre-stellar cores (~0.3 km s^{-1}) derived from submillimeter observations. The origin of this discrepancy, which has been observed in other young star clusters is not clear. It may be due to either the effect of the magnetic field on the protostars and the filaments, or to the dynamical evolution of stars driven by two-body interactions. Furthermore, the analysis of the kinematic properties of the stellar population put in evidence a significant velocity shift (~1 km s^{-1}) between the two sub-clusters located around the North and South main clouds. This result further supports a scenario, where clusters form from the evolution of multiple substructures rather than from a monolithic collapse. Using three independent spectroscopic indicators (the gravity indicator $\gamma$, the equivalent width of the Li line, and the H_alpha 10\% width), we performed a new membership selection. We found six new cluster members located in the outer region of the cluster. Starting from the positions and masses of the cluster members, we derived the level of substructure Q, the surface density \Sigma and the level of mass segregation $\Lambda_{MSR}$ of the cluster. The comparison between these structural properties and the results of N-body simulations suggests that the cluster formed in a low density environment, in virial equilibrium or supervirial, and highly substructured.

We take advantage of the deep and wide coverage of the VST ATLAS survey to study the line-of-sight structure of the Sagittarius stellar stream in the Southern hemisphere, only ~40{\deg} away from the progenitor. We use photometrically selected Sub-Giant Branch (SGB) stars to reveal a complex debris morphology of the trailing arm and detect at least two clear peaks in the SGB distance modulus distribution. The separation between the two line-of-sight components is at least 5 kpc at the edge of the VST ATLAS footprint, but appears to change along the stream, which allows us to conclude that these detections correspond to two physically independent stellar structures, rather than a mix of co-distant stellar populations within a single stream. Our discovery of a fork in the Sgr trailing arm is verified using Blue Horizontal Branch stars and our distance measurements are calibrated using RR Lyrae stars from the Catalina Real-Time Transient Survey. Comparing with numerical simulations of the Sgr dwarf disruption, the more distant of the two components in the fork matches perfectly with the track of the trailing debris. However, no obvious counterpart exists in the simulation for the closer line-of-sight component.

We present results from MMT/Hectochelle spectroscopy of red giant candidate stars along the line of sight to the recently-discovered Galactic satellite Crater 2. Modelling the joint distribution of stellar positions, velocities and metallicities as a mixture of Crater 2 and Galactic foreground populations, we identify 62 members of Crater 2, for which we resolve line-of-sight velocity dispersion 2.7 +/- 0.3 km/s about mean velocity of 87.5 +/- 0.4 km/s. We also resolve a metallicity dispersion 0.22 about a mean of [Fe/H]=-1.98 +/- 0.1 that is 0.28 +/- 0.14 poorer than is estimated from photometry. Despite Crater 2's relatively large size (projected halflight radius R(h)=1 kpc) and intermediate luminosity (M_V =-8), its velocity dispersion is the coldest that has been resolved for any dwarf galaxy. These properties make Crater 2 the most extreme low-density outlier in dynamical as well as structural scaling relations among the Milky Way's dwarf spheroidals. Even so, under assumptions of dynamical equilibrium and negligible contamination by unresolved binary stars, the observed velocity distribution implies a gravitationally dominant dark matter halo, with dynamical mass 4.4 x 10^6 Msun and mass-to-light ratio M/L=53 enclosed within a radius of 1 kpc, where the equivalent circular velocity is 4.3 km/s.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

Context. This paper presents an overview of the photometric data that are part of the first Gaia data release. Aims. The principles of the processing and the main characteristics of the Gaia photometric data are presented. Methods. The calibration strategy is outlined briefly and the main properties of the resulting photometry are presented. Results. Relations with other broadband photometric systems are provided. The overall precision for the Gaia photometry is shown to be at the milli-magnitude level and has a clear potential to improve further in future releases.

The search for chemically unevolved galaxies remains prevalent in the nearby Universe, mostly because these systems provide excellent proxies for exploring in detail the physics of high-z systems. The most promising candidates are extremely metal-poor galaxies (XMPs), i.e., galaxies with <1/10 solar metallicity. However, due to the bright emission line based search criteria traditionally used to find XMPs, we may not be sampling the full XMP population. In 2014 we reoriented this search using only morphological properties and uncovered a population of ~150 `blue diffuse dwarf (BDD) galaxies', and published a sub-sample of 12 BDD spectra. Here we present optical spectroscopic observations of a larger sample of 51 BDDs, along with their SDSS photometric properties. With our improved statistics, we use direct-method abundances to confirm that BDDs are chemically unevolved (7.43<12+log(O/H)<8.01), with ~20% of our sample classified as being XMP galaxies, and find they are actively forming stars at rates of 1-33x10^-2 M_sol/yr in HII regions randomly embedded in a blue, low-surface brightness continuum. Stellar masses are calculated from population synthesis models and estimated to be in the range log(M_star/M_sol) ~5-9. Unlike other low-metallicity star-forming galaxies, BDDs are in agreement with the mass-metallicity relation at low masses, suggesting they are not accreting large amounts of pristine gas relative to their stellar mass. BDD galaxies appear to be a population of actively star-forming dwarf irregular (dIrr) galaxies who fall within the class of low-surface brightness dIrr galaxies. Their ongoing star-formation and irregular morphology make them excellent analogues for galaxies in the early Universe.

We present the discovery and preliminary characterization of a gravitationally lensed quasar with a source redshift $z_{s}=2.74$ and image separation of $2.9"$ lensed by a foreground $z_{l}=0.40$ elliptical galaxy. Since the images of gravitationally lensed quasars are the superposition of multiple point sources and a foreground lensing galaxy, we have developed a morphology independent multi-wavelength approach to the photometric selection of lensed quasar candidates based on Gaussian Mixture Models (GMM) supervised machine learning. Using this technique and $gi$ multicolour photometric observations from the Dark Energy Survey (DES), near IR $JK$ photometry from the VISTA Hemisphere Survey (VHS) and WISE mid IR photometry, we have identified a candidate system with two catalogue components with $i_{AB}=18.61$ and $i_{AB}=20.44$ comprised of an elliptical galaxy and two blue point sources. Spectroscopic follow-up with NTT and the use of an archival AAT spectrum show that the point sources can be identified as a lensed quasar with an emission line redshift of $z=2.739\pm0.003$ and a foreground early type galaxy with $z=0.400\pm0.002$. We model the system as a single isothermal ellipsoid and find the Einstein radius $\theta_E \sim 1.47"$, enclosed mass $M_{enc} \sim 4 \times 10^{11}$M$_{\odot}$ and a time delay of $\sim$52 days. The relatively wide separation, month scale time delay duration and high redshift make this an ideal system for constraining the expansion rate beyond a redshift of 1.

We present the discovery of stellar tidal tails around the Large and the Small Magellanic Clouds in the Gaia DR1 data. In between the Clouds, their tidal arms are stretched towards each other to form an almost continuous stellar bridge. Our analysis relies on the exquisite quality of the Gaia's photometric catalogue to build detailed star-count maps of the Clouds. We demonstrate that the Gaia DR1 data can be used to detect variable stars across the whole sky, and in particular, RR Lyrae stars in and around the LMC and the SMC. Additionally, we use a combination of Gaia and Gale to follow the distribution of Young Main Sequence stars in the Magellanic System. Viewed by Gaia, the Clouds show unmistakable signs of interaction. Around the LMC, clumps of RR Lyrae are observable as far as ~20 degrees, in agreement with the most recent map of Mira-like stars reported in Deason et al (2016). The SMC's outer stellar density contours show a characteristic S-shape, symptomatic of the on-set of tidal stripping. Beyond several degrees from the center of the dwarf, the Gaia RR Lyrae stars trace the Cloud's trailing arm, extending towards the LMC. This stellar tidal tail mapped with RR Lyrae is not aligned with the gaseous Magellanic Bridge, and is shifted by some ~5 degrees from the Young Main Sequence bridge. We use the offset between the bridges to put constraints on the density of the hot gaseous corona of the Milky Way.

We show that non-magnetic models for the evolution of pre-main-sequence (PMS) stars *cannot* simultaneously describe the colour-magnitude diagram (CMD) and the pattern of lithium depletion seen in the cluster of young, low-mass stars surrounding $\gamma^2$ Velorum. The age of 7.5+/-1 Myr inferred from the CMD is much younger than that implied by the strong Li depletion seen in the cluster M-dwarfs and the Li depletion occurs at much redder colours than predicted. The epoch at which a star of a given mass depletes its Li and the surface temperature of that star are both dependent on its radius. We demonstrate that if the low-mass stars have radii ~10 per cent larger at a given mass and age, then both the CMD and Li depletion pattern of the Gamma Vel cluster are explained at a common age of 18-21 Myr. This radius inflation could be produced by some combination of magnetic suppression of convection and extensive cool starspots. Models that incorporate radius inflation suggest that PMS stars similar to those in the Gamma Vel cluster, in the range 0.2<M/Msun<0.7, are at least a factor of two older and ~7 per cent cooler than previously thought and that their masses are much larger (by >30 per cent) than inferred from conventional, non-magnetic models in the Hertzsprung-Russell diagram. Systematic changes of this size may be of great importance in understanding the evolution of young stars, disc lifetimes and the formation of planetary systems.

The orbits, atmospheric parameters, chemical abundances, and ages of individual stars in the Milky Way provide the most comprehensive illustration of galaxy formation available. The Tycho-Gaia Astrometric Solution (TGAS) will deliver astrometric parameters for the largest ever sample of Milky Way stars, though its full potential cannot be realized without the addition of complementary spectroscopy. Among existing spectroscopic surveys, the RAdial Velocity Experiment (RAVE) has the largest overlap with TGAS ($\gtrsim$200,000 stars). We present a data-driven re-analysis of 520,781 RAVE spectra using The Cannon. For red giants, we build our model using high-fidelity APOGEE stellar parameters and abundances for stars that overlap with RAVE. For main-sequence and sub-giant stars, our model uses stellar parameters from the K2/EPIC. We derive and validate effective temperature $T_{\rm eff}$, surface gravity $\log{g}$, and chemical abundances of up to seven elements (O, Mg, Al, Si, Ca, Fe, Ni). We report a total of 1,685,851 elemental abundances with a typical precision of 0.07 dex, a substantial improvement over previous RAVE data releases. The synthesis of RAVE-on and TGAS is the most powerful data set for chemo-dynamic analyses of the Milky Way ever produced.

We announce the discovery of the Aquarius~2 dwarf galaxy, a new distant satellite of the Milky Way, detected on the fringes of the VST ATLAS and the SDSS surveys. The object was originally identified as an overdensity of Red Giant Branch stars, but chosen for subsequent follow-up based on the presence of a strong Blue Horizontal Branch, which was also used to measure its distance of $\sim 110$ kpc. Using deeper imaging from the IMACS camera on the 6.5m Baade and spectroscopy with DEIMOS on Keck, we measured the satellite's half-light radius $5.1\pm 0.8$ arcmin, or $\sim 160$ pc at this distance, and its stellar velocity dispersion of $5.4^{+3.4}_{-0.9}$ km s$^{-1}$. With $\mu=30.2$ mag arcsec$^{-2}$ and $M_V=-4.36$, the new satellite lies close to two important detection limits: one in surface brightness; and one in luminosity at a given distance, thereby making Aquarius~2 one of the hardest dwarfs to find.

We announce the discovery of the Crater 2 dwarf galaxy, identified in imaging data of the VST ATLAS survey. Given its half-light radius of ~1100 pc, Crater 2 is the fourth largest dwarf in the Milky Way, surpassed only by the LMC, SMC and the Sgr dwarf. With a total luminosity of $M_V\approx-8$, this satellite galaxy is also one of the lowest surface brightness dwarfs. Falling under the nominal detection boundary of 30 mag arcsec$^{-2}$, it compares in nebulosity to the recently discovered Tuc 2 and Tuc IV and UMa II. Crater 2 is located ~120 kpc from the Sun and appears to be aligned in 3-D with the enigmatic globular cluster Crater, the pair of ultra-faint dwarfs Leo IV and Leo V and the classical dwarf Leo II. We argue that such arrangement is probably not accidental and, in fact, can be viewed as the evidence for the accretion of the Crater-Leo group.

The discovery of lithium-rich giants contradicts expectations from canonical stellar evolution. Here we report on the serendipitous discovery of 20 Li-rich giants observed during the Gaia-ESO Survey, which includes the first nine Li-rich giant stars known towards the CoRoT fields. Most of our Li-rich giants have near-solar metallicities, and stellar parameters consistent with being before the luminosity bump. This is difficult to reconcile with deep mixing models proposed to explain lithium enrichment, because these models can only operate at later evolutionary stages: at or past the luminosity bump. In an effort to shed light on the Li-rich phenomenon, we highlight recent evidence of the tidal destruction of close-in hot Jupiters at the sub-giant phase. We note that when coupled with models of planet accretion, the observed destruction of hot Jupiters actually predicts the existence of Li-rich giant stars, and suggests Li-rich stars should be found early on the giant branch and occur more frequently with increasing metallicity. A comprehensive review of all known Li-rich giant stars reveals that this scenario is consistent with the data. However more evolved or metal-poor stars are less likely to host close-in giant planets, implying that their Li-rich origin requires an alternative explanation, likely related to mixing scenarios rather than external phenomena.

Red giant stars are perhaps the most important type of stars for Galactic and extra-galactic archaeology: they are luminous, occur in all stellar populations, and their surface temperatures allow precise abundance determinations for many different chemical elements. Yet, the full star formation and enrichment history of a galaxy can be traced directly only if two key observables can be determined for large stellar samples - age and chemical composition. While spectroscopy is a powerful method to analyse the detailed abundances of stars, stellar ages are the "missing link in the chain", since they are not a direct observable. However, spectroscopy should be able to estimate stellar masses, which for red giants directly infer ages provided their chemical composition is known. Here we establish a new empirical relation between the shape of the hydrogen line in the observed spectra of red giants and stellar mass determined from asteroseismology. The relation allows to determine stellar masses and ages with the accuracy of 10-15%. The method can be used with confidence for stars in the following range of stellar parameters: 4000 < Teff < 5000 K, 0.5 < log g < 3.5, -2.0 < [Fe/H] < 0.3, and luminosities log L/LSun < 2.5. Our analysis provides observational evidence that the Halpha spectral characteristics of red giant stars are tightly correlated with their mass and therefore their age. We also show that the method samples well all stellar populations with ages above 1 Gyr. Targeting bright giants, the method allows to obtain simultaneous age and chemical abundance information far deeper than would be possible with asteroseismology, extending the possible survey volume to remote regions of the Milky Way and even to neighbouring galaxies like Andromeda or the Magellanic Clouds already with present instrumentation, like VLT and Keck facilities.

We analyse the kinematics of $\sim 2000$ giant stars in the direction of the Galactic bulge, extracted from the Gaia-ESO survey in the region $-10^\circ \lesssim \ell \lesssim 10^\circ$ and $-11^\circ \lesssim b \lesssim -3^\circ$. We find distinct kinematic trends in the metal rich ($\mathrm{[M/H]}>0$) and metal poor ($\mathrm{[M/H]}<0$) stars in the data. The velocity dispersion of the metal-rich stars drops steeply with latitude, compared to a flat profile in the metal-poor stars, as has been seen previously. We argue that the metal-rich stars in this region are mostly on orbits that support the boxy-peanut shape of the bulge, which naturally explains the drop in their velocity dispersion profile with latitude. The metal rich stars also exhibit peaky features in their line-of-sight velocity histograms, particularly along the minor axis of the bulge. We propose that these features are due to stars on resonant orbits supporting the boxy-peanut bulge. This conjecture is strengthened through the comparison of the minor axis data with the velocity histograms of resonant orbits generated in simulations of buckled bars. The 'banana' or 2:1:2 orbits provide strongly bimodal histograms with narrow velocity peaks that resemble the Gaia-ESO metal-rich data.

The nature of the metallicity gradient inside the solar circle (R_GC < 8 kpc) is poorly understood, but studies of Cepheids and a small sample of open clusters suggest that it steepens in the inner disk. We investigate the metallicity gradient of the inner disk using a sample of inner disk open clusters that is three times larger than has previously been studied in the literature to better characterize the gradient in this part of the disk. We used the Gaia-ESO Survey (GES) [Fe/H] values and stellar parameters for stars in 12 open clusters in the inner disk from GES-UVES data. Cluster mean [Fe/H] values were determined based on a membership analysis for each cluster. Where necessary, distances and ages to clusters were determined via comparison to theoretical isochrones. The GES open clusters exhibit a radial metallicity gradient of -0.10+-0.02 dex/kpc, consistent with the gradient measured by other literature studies of field red giant stars and open clusters in the range R_GC ~ 6-12 kpc. We also measure a trend of increasing [Fe/H] with increasing cluster age, as has also been found in the literature. We find no evidence for a steepening of the inner disk metallicity gradient inside the solar circle as earlier studies indicated. The age-metallicity relation shown by the clusters is consistent with that predicted by chemical evolution models that include the effects of radial migration, but a more detailed comparison between cluster observations and models would be premature.

We report the discovery of two Li-rich giants, with A(Li) ~ 1.50, in an analysis of a sample of 40 giants of the open cluster Trumpler 20 (with turnoff mass ~ 1.8 Msun). The cluster was observed in the context of the Gaia-ESO Survey. The atmospheric parameters and Li abundances were derived using high-resolution UVES spectra. The Li abundances were corrected for nonlocal thermodynamical equilibrium (non-LTE) effects. Only upper limits of the Li abundance could be determined for the majority of the sample. Two giants with detected Li turned out to be Li rich: star MG 340 has A(Li) non-LTE = 1.54 \pm 0.21 dex and star MG 591 has A(Li) non-LTE = 1.60 \pm 0.21 dex. Star MG 340 is on average ~ 0.30 dex more rich in Li than stars of similar temperature, while for star MG 591 this difference is on average ~ 0.80 dex. Carbon and nitrogen abundances indicate that all stars in the sample have completed the first dredge-up. The Li abundances in this unique sample of 40 giants in one open cluster clearly show that extra mixing is the norm in this mass range. Giants with Li abundances in agreement with the predictions of standard models are the exception. To explain the two Li-rich giants, we suggest that all events of extra mixing have been inhibited. This includes rotation-induced mixing during the main sequence and the extra mixing at the red giant branch luminosity bump. Such inhibition has been suggested in the literature to occur because of fossil magnetic fields in red giants that are descendants of main-sequence Ap-type stars.

The Gaia-ESO Survey was designed to target all major Galactic components (i.e., bulge, thin and thick discs, halo and clusters), with the goal of constraining the chemical and dynamical evolution of the Milky Way. This paper presents the methodology and considerations that drive the selection of the targeted, allocated and successfully observed Milky Way field stars. The detailed understanding of the survey construction, specifically the influence of target selection criteria on observed Milky Way field stars is required in order to analyse and interpret the survey data correctly. We present the target selection process for the Milky Way field stars observed with VLT/FLAMES and provide the weights that characterise the survey target selection. The weights can be used to account for the selection effects in the Gaia-ESO Survey data for scientific studies. We provide a couple of simple examples to highlight the necessity of including such information in studies of the stellar populations in the Milky Way.

We investigate here the lithium content of young low-mass stars in the 5~Myr-old star forming region NGC~2264 and its relationship with rotation. We combine lithium equivalent width measurements, EW(Li), from the Gaia-ESO Survey with the determination of rotational periods from the CSI~2264 survey. We consider only bona fide non accreting cluster members in order to minimize uncertainties on EW(Li). We report the existence of a relationship between lithium content and rotation in NGC~2264 at an age of 5~Myr. The Li-rotation connection is seen over a restricted temperature range, Teff=3800--4400K, where fast rotators are Li-rich compared to slow ones. This correlation is similar to, albeit of lower amplitude than, the Li-rotation connection previously reported for K dwarfs in the 125 Myr-old Pleiades cluster. We investigate whether the non-standard pre-main sequence models developed so far to explain the Pleiades results, which are based on episodic accretion, pre-main sequence core-envelope decoupling, and/or radius inflation due to enhanced magnetic activity, can account for an early development of the Li-rotation connection. While radius inflation appears to be the most promising possibility, each of these models has issues. We therefore also discuss external causes that might operate during the first few Myr of pre-main sequence evolution, such as planet engulfment and/or steady disk accretion, as possible candidates for the common origin for Li-excess and fast rotation in young low-mass pre-main sequence stars. The emergence of a connection between lithium content and rotation rate at such an early age as 5~Myr suggests a complex link between accretion processes, early angular momentum evolution, and possibly planet formation, which likely impacts early stellar evolution and still is to be fully deciphered.

We present observations from the Gaia-ESO Survey in the lines of H$\alpha$, [N II], [S II] and He I of nebular emission in the central part of the Carina Nebula. We investigate the properties of the two already known kinematic components (approaching and receding, respectively), which account for the bulk of emission. Moreover, we investigate the features of the much less known low-intensity high-velocity (absolute RV $>$50 km/s) gas emission. We show that gas giving rise to H$\alpha$ and He I emission is dynamically well correlated, but not identical, to gas seen through forbidden-line emission. Gas temperatures are derived from line-width ratios, and densities from [S II] doublet ratios. The spatial variation of N ionization is also studied, and found to differ between the approaching and receding components. The main result is that the bulk of the emission lines in the central part of Carina arises from several distinct shell-like expanding regions, the most evident found around $\eta$ Car, the Trumpler 14 core, and the star WR25. Such "shells" are non-spherical, and show distortions probably caused by collisions with other shells or colder, higher-density gas. Part of them is also obscured by foreground dust lanes, while only very little dust is found in their interior. Preferential directions, parallel to the dark dust lanes, are found in the shell geometries and physical properties, probably related to strong density gradients in the studied region. We also find evidence that the ionizing flux emerging from $\eta$ Car and the surrounding Homunculus nebula varies with polar angle. The high-velocity components in the wings of H$\alpha$ are found to arise from expanding dust reflecting the $\eta$ Car spectrum.

We present an update of a spectroscopic follow-up survey at low-resolution of a large number of RR Lyrae halo overdensity candidates found in the southern sky. The substructure candidates were identified in the RR Lyrae catalog of Torrealba et al. (2015) using Catalina Real-time Transient Survey (CRTS) data. Radial velocities and mean metallicities have been estimated for target stars in almost half of the original overdensities to assess their potential membership to coherent halo features.

We report the discovery of a substantial stellar overdensity in the periphery of the Large Magellanic Cloud (LMC), found using public imaging from the first year of the Dark Energy Survey. The structure appears to emanate from the edge of the outer LMC disk at a radius $\approx 13.5$ degrees due north of its centre, and stretches more than $10$ kpc towards the east. It is roughly $1.5$ kpc wide and has an integrated $V$-band luminosity of at least $M_V = -7.4$. The stellar populations in the feature are indistinguishable from those in the outer LMC disk. We attempt to quantify the geometry of the outer disk using simple planar models, and find that only a disk with mild intrinsic ellipticity can simultaneously explain the observed stellar density on the sky and the azimuthal line-of-sight distance profile. We also see possible non-planar behaviour in the outer disk that may reflect a warp and/or flare, as well as deviations that resemble a ring-like structure between $\sim9-12$ degrees from the LMC centre. Based on all these observations, we conclude that our remote, stream-like feature is likely comprised of material that has been stripped from the outskirts of the LMC disk, although we cannot rule out that it represents a transient overdensity in the disk itself. We conduct a simple $N$-body simulation to show that either type of structure could plausibly arise due to the tidal force of the Milky Way; however we also recognize that a recent close interaction between the LMC and the SMC may be the source of the stripping or perturbation. Finally, we observe evidence for extremely diffuse LMC populations extending to radii of $\sim 18.5$ kpc in the disk plane ($\approx 20$ degrees on the sky), corroborating previous spectroscopic detections at comparable distances.

Stellar evolution models predict that internal mixing should cause some sodium overabundance at the surface of red giants more massive than ~ 1.5--2.0 Msun. The surface aluminium abundance should not be affected. Nevertheless, observational results disagree about the presence and/or the degree of the Na and Al overabundances. In addition, Galactic chemical evolution models adopting different stellar yields lead to quite different predictions for the behavior of [Na/Fe] and [Al/Fe] versus [Fe/H]. Overall, the observed trends of these abundances with metallicity are not well reproduced. We readdress both issues, using new Na and Al abundances determined within the Gaia-ESO Survey, using two samples: i) more than 600 dwarfs of the solar neighborhood and of open clusters and ii) low- and intermediate-mass clump giants in six open clusters. Abundances of Na in giants with mass below ~2.0 Msun, and of Al in giants below ~3.0 Msun, seem to be unaffected by internal mixing processes. For more massive giants, the Na overabundance increases with stellar mass. This trend agrees well with predictions of stellar evolutionary models. Chemical evolution models that are able to fit well the observed [Na/Fe] vs. [Fe/H] trend in solar neighborhood dwarfs can not simultaneously explain the run of [Al/Fe] with [Fe/H], and viceversa. The comparison with stellar ages is hampered by severe uncertainties. Indeed, reliable age estimates are available for only a half of the stars of the sample. We conclude that Al is underproduced by the models, except for stellar ages younger than about 7 Gyr. In addition, some significant source of late Na production seems to be missing in the models. Either current Na and Al yields are affected by large uncertainties, and/or some important Galactic source(s) of these elements has not been taken into account up to now. [abridged]

We study the chemodynamical properties of the Carina dwarf spheroidal by combining an intermediate spectroscopic resolution dataset of more than 900 red giant and red clump stars, with high-precision photometry to derive the atmospheric parameters, metallicities and age estimates for our targets. Within the red giant branch population, we find evidence for the presence of three distinct stellar sub-populations with different metallicities, spatial distributions, kinematics and ages. As in the Fornax and Sculptor dwarf spheroidals, the subpopulation with the lowest average metallicity is more extended and kinematically hotter than all other populations. However, we identify an inversion in the parallel ordering of metallicity, kinematics and characteristic length scale in the two most metal rich subpopulations, which therefore do not contribute to a global negative chemical gradient. Contrary to common trends in the chemical properties with radius, the metal richest population is more extended and mildly kinematically hotter than the main component of intermediate metallicity. More investigations are required to ascertain the nature of this inversion, but we comment on the mechanisms that might have caused it.

We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of $147$ stellar targets along the line of sight to the newly-discovered `ultrafaint' stellar systems Tucana 2 (Tuc 2) and Grus 1 (Gru 1). Based on simultaneous estimates of line-of-sight velocity and stellar-atmospheric parameters, we identify 8 and 7 stars as probable members of Tuc 2 and and Gru 1, respectively. Our sample for Tuc 2 is sufficient to resolve an internal velocity dispersion of $8.6_{-2.7}^{+4.4}$ km s$^{-1}$ about a mean of $-129.1_{-3.5}^{+3.5}$ km s$^{-1}$ (solar rest frame), and to estimate a mean metallicity of [Fe/H]= $-2.23_{-0.12}^{+0.18}$. These results place Tuc 2 on chemodynamical scaling relations followed by dwarf galaxies, suggesting a dominant dark matter component with dynamical mass $2.7_{-1.3}^{+3.1}\times 10^6$ $\mathrm{M}_{\odot}$ enclosed within the central $\sim 160$ pc, and dynamical mass-to-light ratio $1900_{-900}^{+2200}$ $\mathrm{M}_{\odot}/L_{V,\odot}$. For Gru 1 we estimate a mean velocity of $-140.5_{-1.6}^{+2.4}$ km s$^{-1}$ and a mean metallicity of [Fe/H]=$-1.42_{-0.42}^{+0.55}$, but our sample does not resolve Gru 1's velocity dispersion. The radial coordinates of Tuc 2 and Gru 1 in Galactic phase space suggest that their orbits are among the most energetic within distance $\leq 300$ kpc. Moreover, their proximity to each other in this space arises naturally if both objects are trailing the Large Magellanic Cloud.

(Abridged) Recent spectroscopic surveys have begun to explore the Galactic disk system outside the solar neighborhood on the basis of large data samples. In this way, they provide valuable information for testing spatial and temporal variations of disk structure kinematics and chemical evolution. We used a Gaussian mixture model algorithm, as a rigurous mathematical approach, to separate in the [Mg/Fe] vs. [Fe/H] plane a clean disk star subsample from the Gaia-ESO survey internal data release 2. We find that the sample is separated into five groups associated with major Galactic components; the metal-rich end of the halo, the thick disk, and three subgroups for the thin disk sequence. This is confirmed with a sample of red clump stars from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. The two metal-intermediate and metal-rich groups of the thin disk decomposition ([Fe/H]>-0.25 dex) highlight a change in the slope at solar metallicity. This holds true at different radial regions. The distribution of Galactocentric radial distances of the metal-poor part of the thin disk ([Fe/H]<-0.25 dex) is shifted to larger distances than those of the more metal-rich parts. Moreover, the metal-poor part of the thin disk presents indications of a scale height intermediate between those of the thick and the rest of the thin disk, and it displays higher azimuthal velocities than the latter. These stars might have formed and evolved in parallel and/or dissociated from the inside-out formation taking place in the internal thin disk. Their enhancement levels might be due to their origin from gas pre-enriched by outflows from the thick disk or the inner halo. The smooth trends of their properties (their spatial distribution with respect to the plane, in particular) with [Fe/H] and [Mg/Fe] suggested by the data indicates a quiet dynamical evolution, with no relevant merger events.

Using Blue Horizontal Branch stars identified in the Dark Energy Survey Year 1 data, we report the detection of an extended and lumpy stellar debris distribution around the Magellanic Clouds. At the heliocentric distance of the Clouds, overdensities of BHBs are seen to reach at least to ~30 degrees, and perhaps as far as ~50 degrees from the LMC. In 3D, the stellar halo is traceable to between 25 and 50 kpc from the LMC. We catalogue the most significant of the stellar sub-structures revealed, and announce the discovery of a number of narrow streams and diffuse debris clouds. Two narrow streams appear approximately aligned with the Magellanic Clouds' proper motion. Moreover, one of these overlaps with the gaseous Magellanic Stream on the sky. Curiously, two diffuse BHB agglomerations seem coincident with several of the recently discovered DES satellites. Given the enormous size and the conspicuous lumpiness of the LMC's stellar halo, we speculate that the dwarf could easily have been more massive than previously had been assumed.

We analyze the resolved stellar populations of the faint stellar system, Crater, based on deep optical imaging taken with the Hubble Space Telescope. The HST/ACS-based color-magnitude diagram (CMD) of Crater extends $\sim$4 magnitudes below the oldest main sequence turnoff, providing excellent leverage on Crater's physical properties. Structurally, Crater has a half-light radius of $\sim$20 pc and shows no evidence for tidal distortions. Crater is well-described by a simple stellar population with an age of $\sim$7.5 Gyr, [M/H]$\sim-1.65$, a M$_{\star}\sim10^4$ M$_{\odot}$, M$_{\rm V}\sim -5.3$, located at a distance of d$\sim$ 145 kpc, with modest uncertainties in these properties due to differences in the underlying stellar evolution models. The sparse sampling of stars above the turnoff and sub-giant branch are likely to be 1.0-1.4 M$_{\odot}$ binary star systems (blue stragglers) and their evolved descendants, as opposed to intermediate age main sequence stars. Confusion of these populations highlights a substantial challenge in accurately characterizing sparsely populated stellar systems. Our analysis shows that Crater is not a dwarf galaxy, but instead is an unusually young cluster given its location in the Milky Way's very outer stellar halo. Crater is similar to SMC cluster Lindsay 38, and its position and velocity are in good agreement with observations and models of the Magellanic stream debris, suggesting it may have accreted from the Magellanic Clouds. However, its age and metallicity are also in agreement with the age-metallicity relationships of lower mass dwarf galaxies such as Leo I or Carina. Despite uncertainty over its progenitor system, Crater appears to have been incorporated into the Galaxy more recently than $z\sim1$ (8 Gyr ago), providing an important new constraint on the accretion history of the Milky Way. [abridged]

We present a study of the detectability of transient events associated with galaxies for the Gaia European Space Agency astrometric mission. We simulated the on-board detections, and on-ground processing for a mock galaxy catalogue to establish the properties required for the discovery of transient events by Gaia, specifically tidal disruption events (TDEs) and supernovae (SNe). Transients may either be discovered by the on-board detection of a new source or by the brightening of a previously known source. We show that Gaia transients can be identified as new detections on-board for offsets from the host galaxy nucleus of 0.1--0.5,arcsec, depending on magnitude and scanning angle. The Gaia detection system shows no significant loss of SNe at close radial distances to the nucleus. We used the detection efficiencies to predict the number of transients events discovered by Gaia. For a limiting magnitude of 19, we expect around 1300 SNe per year: 65% SN Ia, 28% SN II and 7% SN Ibc, and ~20 TDEs per year.

Understanding the history and the evolution of the Milky Way disc is one of the main goals of modern astrophysics. We study the velocity dispersion behaviour of Galactic disc stars as a function of the [Mg/Fe] ratio, which can be used as a proxy of relative age. This key relation is essential to constrain the formation mechanisms of the disc stellar populations as well as the cooling processes. We used the recommended parameters and chemical abundances of 7800 FGK Milky Way field stars from the second internal data release of the Gaia-ESO Survey. These stars were observed with the GIRAFFE spectrograph, and cover a large spatial volume (6<R<10kpc and |Z|<2kpc). Based on the [Mg/Fe] and [Fe/H] ratios, we separated the thin- from the thick-disc sequence. From analysing the Galactocentric velocity of the stars for the thin disc, we find a weak positive correlation between Vphi and [Fe/H], due to a slowly rotating Fe-poor tail. For the thick disc, a strong correlation with [Fe/H] and [Mg/Fe] is established. We have detected an inversion of the radial velocity dispersion with [Mg/Fe] for thick-disc stars with [Fe/H]<-0.1dex and [Mg/Fe]>+0.2dex. First, the velocity dispersion increases with [Mg/Fe] at all [Fe/H] ratios for the thin-disc stars, and then it decreases for the thick-disc at the highest [Mg/Fe] abundances. Similar trends are observed within the errors for the azimuthal velocity dispersion, while a continuous increase with [Mg/Fe] is observed for the vertical velocity dispersion. The velocity dispersion decrease agrees with previous measurements of the RAVE survey, although it is observed here for a greater metallicity interval and a larger spatial volume. We confirm the existence of [Mg/Fe]-rich thick-disc stars with cool kinematics in the generally turbulent context of the primitive Galactic disc. This is discussed in the framework of the different disc formation scenarios.

It is well known that newly formed planetary systems undergo processes of orbital reconfiguration and planetary migration. As a result, planets or protoplanetary objects may accrete onto the central star, being fused and mixed into its external layers. If the accreted mass is sufficiently high and the star has a sufficiently thin convective envelope, such events may result in a modification of the chemical composition of the stellar photosphere in an observable way, enhancing it with elements that were abundant in the accreted mass. The recent Gaia-ESO Survey observations of the 10-20 Myr old Gamma Velorum cluster have enabled identifying a star that is significantly enriched in iron with respect to other cluster members. In this Letter we further investigate the abundance pattern of this star, showing that its abundance anomaly is not limited to iron, but is also present in the refractory elements, whose overabundances are correlated with the condensation temperature. This finding strongly supports the hypothesis of a recent accretion of rocky material.

SCUSS is a u-band photometric survey covering about 4000 square degree of the South Galactic Cap, reaching depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this advantage we combine SCUSS u band with SDSS DR9 gri photometric bands to identify BHB stars and explore halo substructures. We confirm the existence of the Pisces overdensity, which is a structure in the outer halo (at around 80 kpc) that was discovered using RR Lyrae stars. For the first time we are able to determine its spatial extent, finding that it appears to be part of a stream with a clear distance gradient. The stream, which is ~5 degrees wide and stretches along ~25 degrees, consists of 20-30 BHBs with a total significance of around 6sigma over the background. Assuming we have detected the entire stream and that the progenitor has fully disrupted, then the number of BHBs suggests the original system was similar to smaller classical or a larger ultra-faint dwarf galaxy. On the other hand, if the progenitor still exists, it can be hunted for by reconstructing its orbit from the distance gradient of the stream. This new picture of the Pisces overdensity sheds new light on the origin of this intriguing system.

We investigate, using the Gaia-ESO Survey internal Data-Release 2, the properties of the double sequence of the Milky Way discs (defined chemically as the high-alpha and low-alpha populations), and discuss their compatibility with discs defined by other means such as metallicity, kinematics or positions. This investigation uses two different approaches: in velocity space for stars located in the extended Solar neighbourhood, and in chemical space for stars at different ranges of Galactocentric radii and heights from the plane. The separation we find in velocity space allows us to investigate, in a novel manner, the extent in metallicity of each of the two sequences, identifying them with the two discs, without making any assumption about the shape of their metallicity distribution functions. Then, using the separation in chemical space, we characterise the spatial variation of the slopes of the [alpha/Fe] - [Fe/H] sequences for the thick and thin discs and the way in which the relative proportions of the two discs change across the Galaxy. We find that the thick disc (high-alpha sequence), extends up to [Fe/H]~ +0.2 and the thin disc (low-alpha sequence), at least down to [Fe/H]~ -0.8. Radial and vertical gradients in alpha-abundances are found for the thin disc, with mild spatial variations in its [alpha/Fe] - [Fe/H] paths, whereas for the thick disc we do not detect any such spatial variations. The small variations in the spatial [alpha/Fe] - [Fe/H] paths of the thin disc do not allow us to distinguish between formation models of this structure. On the other hand, the lack of radial gradients and [alpha/Fe] - [Fe/H] variations for the thick disc indicate that the mechanism responsible for the mixing of the metals in the young Galaxy (e.g. radial stellar migration or turbulent gaseous disc) was more efficient before the (present) thin disc started forming.

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events.

We present a search for gamma-ray emission from the direction of the newly discovered dwarf galaxy Reticulum II. Using Fermi-LAT data, we detect a signal that exceeds expected backgrounds between ~2-10 GeV and is consistent with annihilation of dark matter for particle masses less than a few x 10^2 GeV. Modeling the background as a Poisson process based on Fermi-LAT diffuse models, and taking into account trials factors, we detect emission with p-value less than 9.8 x 10^-5 (>3.7 sigma). An alternative, model-independent treatment of background reduces the significance, raising the p-value to 9.7 x 10^-3 (2.3 sigma). Even in this case, however, Reticulum II has the most significant gamma-ray signal of any known dwarf galaxy. If Reticulum II has a dark matter halo that is similar to those inferred for other nearby dwarfs, the signal is consistent with the s-wave relic abundance cross section for annihilation.

We report five Local Volume dwarf galaxies (two of which are presented here for the first time) uncovered during a comprehensive archival search for optical counterparts to ultra-compact high velocity clouds (UCHVCs). The UCHVC population of HI clouds are thought to be candidate gas-rich, low mass halos at the edge of the Local Group and beyond, but no comprehensive search for stellar counterparts to these systems has been presented. Careful visual inspection of all publicly available optical and ultraviolet imaging at the position of the UCHVCs revealed six blue, diffuse counterparts with a morphology consistent with a faint dwarf galaxy beyond the Local Group. Optical spectroscopy of all six candidate dwarf counterparts show that five have an H$\alpha$-derived velocity consistent with the coincident HI cloud, confirming their association, the sixth diffuse counterpart is likely a background object. The size and luminosity of the UCHVC dwarfs is consistent with other known Local Volume dwarf irregular galaxies. The gas fraction ($M_{HI}/M_{star}$) of the five dwarfs are generally consistent with that of dwarf irregular galaxies in the Local Volume, although ALFALFA-Dw1 (associated with ALFALFA UCHVC HVC274.68+74.70$-$123) has a very high $M_{HI}/M_{star}$$\sim$40. Despite the heterogenous nature of our search, we demonstrate that the current dwarf companions to UCHVCs are at the edge of detectability due to their low surface brightness, and that deeper searches are likely to find more stellar systems. If more sensitive searches do not reveal further stellar counterparts to UCHVCs, then the dearth of such systems around the Local Group may be in conflict with $\Lambda$CDM simulations.

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.

We present the first detailed quantitative study of the stellar populations of the Sagittarius (Sgr) streams within the Stripe 82 region, using photometric and spectroscopic observations from the Sloan Digital Sky Survey (SDSS). The star formation history (SFH) is determined separately for the bright and faint Sgr streams, to establish whether both components consist of a similar stellar population mix or have a distinct origin. Best fit SFH solutions are characterised by a well-defined, tight sequence in age-metallicity space, indicating that star formation occurred within a well-mixed, homogeneously enriched medium. Star formation rates dropped sharply at an age of ~5-7 Gyr, possibly related to the accretion of Sgr by the MW. Finally, the Sgr sequence displays a change of slope in age-metallicity space at an age between 11-13 Gyr consistent with the Sgr alpha-element knee, indicating that supernovae type Ia started contributing to the abundance pattern ~1-3 Gyr after the start of star formation. Results for both streams are consistent with being drawn from the parent Sgr population mix, but at different epochs. The SFH of the bright stream starts from old, metal-poor populations and extends to a metallicity of [Fe/H]~-0.7, with peaks at ~7 and 11 Gyr. The faint SFH samples the older, more metal-poor part of the Sgr sequence, with a peak at ancient ages and stars mostly with [Fe/H]<-1.3 and age>9 Gyr. Therefore, we argue in favour of a scenario where the faint stream consists of material stripped i) earlier, and ii) from the outskirts of the Sgr dwarf.

Photometry alone is not sufficient to unambiguously distinguish between ultra-faint star clusters and dwarf galaxies because of their overlap in morphological properties. Here we report on VLT/GIRAFFE spectra of candidate member stars in two recently discovered ultra-faint satellites Reticulum 2 and Horologium 1, obtained as part of the ongoing Gaia-ESO Survey. We identify 18 members in Reticulum 2 and 5 in Horologium 1. We find Reticulum 2 to have a velocity dispersion of ~3.22 km/s, implying a M/L ratio of ~ 500. We have inferred stellar parameters for all candidates and we find Reticulum 2 to have a mean metallicity of [Fe/H] = -2.46+/-0.1, with an intrinsic dispersion of ~ 0.29, and is alpha-enhanced to the level of [alpha/Fe]~0.4. We conclude that Reticulum 2 is a dwarf galaxy. We also report on the serendipitous discovery of four stars in a previously unknown stellar substructure near Reticulum 2 with [Fe/H] ~ -2 and V_hel ~ 220 km/s, far from the systemic velocity of Reticulum 2. For Horologium 1 we infer a velocity dispersion of 4.9^{+2.8}_{-0.9} km/s and a consequent M/L ratio of ~ 600, leading us to conclude that Horologium 1 is also a dwarf galaxy. Horologium 1 is slightly more metal-poor than Reticulum 2 [Fe/H] = -2.76 +/- 0.1 and is similarly alpha-enhanced: [alpha/Fe] ~ 0.3. Despite a large error-bar, we also measure a significant spread of metallicities of 0.17 dex which strengthen the evidence that Horologium 1 is indeed a dwarf galaxy. The line-of-sight velocity of Reticulum 2 is offset by some 100 km/s from the prediction of the orbital velocity of the LMC, thus making its association with the Cloud uncertain. However, at the location of Horologium 1, both the backward integrated LMC's orbit and the LMC's halo are predicted to have radial velocities similar to that of the dwarf. Therefore, it is likely that Horologium 1 is or once was a member of the Magellanic Family.

We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of 182 stellar targets along the line of sight to the newly-discovered `ultrafaint' object Reticulum 2 (Ret 2). For 38 of these targets, the spectra are sufficient to provide simultaneous estimates of line-of-sight velocity ($v_{\rm los}$, median random error $\delta_{v_{\rm los}}=1.3$ km s$^{-1}$), effective temperature ($T_{\rm eff}$, $\delta_{T_{\rm eff}}=464$ K), surface gravity ($\log g$, $\delta_{\rm logg}=0.54$ dex) and iron abundance ([Fe/H], $\delta_{\mathrm{[Fe/H]}}=0.45$ dex). We use these results to confirm 18 stars as members of Ret 2. From the member sample we estimate a velocity dispersion of $\sigma_{v_{\rm los}}=3.6_{-0.6}^{+0.9}$ km s$^{-1}$ about a mean of $\langle v_{\rm los}\rangle =64.8_{-1.0}^{+1.1}$ km s$^{-1}$ in the solar rest frame ($\sim -90.9$ km s$^{-1}$ in the Galactic rest frame), and a metallicity dispersion of $\sigma_{\rm [Fe/H]}=0.50_{-0.13}^{+0.17}$ dex about a mean of $\langle \mathrm{[Fe/H]} \rangle =-2.67_{-0.34}^{+0.34}$. These estimates marginalize over possible velocity and metallicity gradients, which are consistent with zero. Our results place Ret 2 on chemodynamical scaling relations followed by the Milky Way's dwarf-galactic satellites. Under assumptions of dynamic equilibrium and negligible contamination from binary stars---both of which must be checked with deeper imaging and repeat spectroscopic observations---the estimated velocity dispersion suggests a dynamical mass of $M(R_{\rm h})\approx 5R_{\rm h}\sigma_{v_{\rm los}}^2/(2G)=2.4_{-0.8}^{+1.3}\times 10^5$ $M_{\odot}$ enclosed within projected halflight radius $R_{\rm h}\sim 32$ pc, with mass-to-light ratio $\approx 2M(R_{\rm h})/L_{\rm V}=462_{-157}^{+264}$ in solar units.

We have used the publicly released Dark Energy Survey data to hunt for new satellites of the Milky Way in the Southern hemisphere. Our search yielded a large number of promising candidates. In this paper, we announce the discovery of 9 new unambiguous ultra-faint objects, whose authenticity can be established with the DES data alone. Based on the morphological properties, three of the new satellites are dwarf galaxies, one of which is located at the very outskirts of the Milky Way, at a distance of 380 kpc. The remaining 6 objects have sizes and luminosities comparable to the Segue~1 satellite and can not be classified straightforwardly without follow-up spectroscopic observations. The satellites we have discovered cluster around the LMC and the SMC. We show that such spatial distribution is unlikely under the assumption of isotropy, and, therefore, conclude that at least some of the new satellites must have been associated with the Magellanic Clouds in the past.

Extremely metal poor (XMP) galaxies are known to be very rare, despite the large numbers of low-mass galaxies predicted by the local galaxy luminosity function. This paper presents a sub-sample of galaxies that were selected via a morphology-based search on SDSS images with the aim of finding these elusive XMP galaxies. By using the recently discovered extremely metal-poor galaxy, Leo P, as a guide, we obtained a collection of faint, blue systems, each with isolated HII regions embedded in a diffuse continuum, that have remained optically undetected until now. Here we show the first results from optical spectroscopic follow-up observations of 12 of ~100 of these blue, diffuse dwarf (BDD) galaxies yielded by our search algorithm. Oxygen abundances were obtained via the direct method for eight galaxies, and found to be in the range 7.45<12+log(O/H)<8.0, with two galaxies being classified as XMPs. All BDDs were found to currently have a young star-forming population (<10 Myr) and relatively high ionisation parameters of their HII regions. Despite their low luminosities (-11< M_B < -18) and low surface brightnesses (~23-25 mag arcsec^-2), the galaxies were found to be actively star-forming, with current star-formation rates between 0.0003 and 0.078 Msol/yr. From our current subsample, BDD galaxies appear to be a population of non-quiescent dwarf irregular (dIrr) galaxies, or the diffuse counterparts to blue compact galaxies (BCDs) and as such may bridge the gap between these two populations. Our search algorithm demonstrates that morphology-based searches are successful in uncovering more diffuse metal-poor star-forming galaxies, which traditional emission-line based searches overlook.

A small fraction of the halo field is made up of stars that share the light element (Z<=13) anomalies characteristic of second generation globular cluster (GC) stars. The ejected stars shed light on the formation of the Galactic halo by tracing the dynamical history of the clusters, which are believed to have once been more massive. Some of these ejected stars are expected to show strong Al enhancement at the expense of shortage of Mg, but until now no such star has been found. We search for outliers in the Mg and Al abundances of the few hundreds of halo field stars observed in the first eighteen months of the Gaia-ESO public spectroscopic survey. One halo star at the base of the red giant branch, here referred to as 22593757-4648029 is found to have [Mg/Fe]=-0.36+-0.04 and [Al/Fe]=0.99+-0.08, which is compatible with the most extreme ratios detected in GCs so far. We compare the orbit of 22593757-4648029 to GCs of similar metallicity and find it unlikely that this star has been tidally stripped with low ejection velocity from any of the clusters. However, both chemical and kinematic arguments render it plausible that the star has been ejected at high velocity from the anomalous GC omega Centauri within the last few billion years. We cannot rule out other progenitor GCs, because some may have disrupted fully, and the abundance and orbital data are inadequate for many of those that are still intact.

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.

The Gaia-ESO survey is a large public spectroscopic survey aimed at investigating the origin and formation history of our Galaxy by collecting spectroscopy of representative samples (about 10^5 Milky Way stars) of all Galactic stellar populations, in the field and in clusters. The survey uses globular clusters as intra- and inter-survey calibrators, deriving stellar atmospheric parameters and abundances of a significant number of stars in clusters, along with radial velocity determinations. We used precise radial velocities of a large number of stars in seven globular clusters (NGC 1851, NGC 2808, NGC 4372, NGC 4833, NGC 5927, NGC 6752, and NGC 7078) to validate pipeline results and to preliminarily investigate the cluster internal kinematics. Radial velocity measurements were extracted from FLAMES/GIRAFFE spectra processed by the survey pipeline as part of the second internal data release of data products to ESO. We complemented our sample with ESO archival data obtained with different instrument configurations. Reliable radial velocity measurements for 1513 bona fide cluster star members were obtained in total. We measured systemic rotation, estimated central velocity dispersions, and present velocity dispersion profiles of all the selected clusters, providing the first velocity dispersion curve and the first estimate of the central velocity dispersion for the cluster NGC~5927. Finally, we explore the possible link between cluster kinematics and other physical parameters. The analysis we present here demonstrates that Gaia-ESO survey data are sufficiently accurate to be used in studies of kinematics of stellar systems and stellar populations in the Milky Way.

We present the results of a deep, wide-area variability survey in the Southern hemisphere, the first of its kind. As part of the Catalina Sky Surveys, the Siding Spring Survey (SSS) has covered $14,800$ square degrees in the declination range of $-75^{\circ}\leq\delta\leq-15^{\circ}$. To mine the enormous SSS dataset efficiently we have developed two algorithms: Automatic Period Selection (APS) and Automatic Fourier Decomposition (AFD), which aim to sharpen the period estimation and produce robust lightcurve models. Armed with the APS and AFD outputs we classify $10,540$ ab-type RR Lyrae (RRab) stars ($\sim$90% of which are new) across the Southern sky. As well as the positional information we supply photometric metallicities, and unreddened distances. For the RRab stars in the halo, a study of the photometric metallicity distribution reveals a nearly Gaussian shape with a mean metallicity of ${\rm [Fe/H]}=-1.4$ dex and a dispersion of $0.3$ dex. A spatial study of the RRab metallicities shows no significant radial gradient in the first $\sim7$ kpc from the Galaxy center. However, further out, a small negative gradient is clearly present. This is complemented by a very obvious correlation of the mean RR Lyrae metallicity with distance above the Galactic plane, $z$. We have also carried out an initial substructure search using the discovered RRab, and present the properties of the candidates with significance greater than $2 \sigma$. Most prominent among these is a southern extension of the Sagittarius dwarf galaxy's stream system, reaching down to declinations $\sim -40\deg$.

We present the combination of optical data from the Science Verification phase of the Dark Energy Survey (DES) with near infrared data from the ESO VISTA Hemisphere Survey (VHS). The deep optical detections from DES are used to extract fluxes and associated errors from the shallower VHS data. Joint 7-band ($grizYJK$) photometric catalogues are produced in a single 3 sq-deg DECam field centred at 02h26m$-$04d36m where the availability of ancillary multi-wavelength photometry and spectroscopy allows us to test the data quality. Dual photometry increases the number of DES galaxies with measured VHS fluxes by a factor of $\sim$4.5 relative to a simple catalogue level matching and results in a $\sim$1.5 mag increase in the 80\% completeness limit of the NIR data. Almost 70\% of DES sources have useful NIR flux measurements in this initial catalogue. Photometric redshifts are estimated for a subset of galaxies with spectroscopic redshifts and initial results, although currently limited by small number statistics, indicate that the VHS data can help reduce the photometric redshift scatter at both $z<0.5$ and $z>1$. We present example DES+VHS colour selection criteria for high redshift Luminous Red Galaxies (LRGs) at $z\sim0.7$ as well as luminous quasars. Using spectroscopic observations in this field we show that the additional VHS fluxes enable a cleaner selection of both populations with $<$10\% contamination from galactic stars in the case of spectroscopically confirmed quasars and $<0.5\%$ contamination from galactic stars in the case of spectroscopically confirmed LRGs. The combined DES+VHS dataset, which will eventually cover almost 5000 sq-deg, will therefore enable a range of new science and be ideally suited for target selection for future wide-field spectroscopic surveys.

We show that a combination of infrared photometry from WISE and 2MASS surveys can yield highly pure samples of M giant stars. We take advantage of the new WISE$\cap$2MASS M giant selection to trace the Sagittarius trailing tail behind the Galactic disk in the direction of the anti-centre. The M giant candidates selected via broad-band photometry are confirmed spectroscopically using AAOmega on the AAT in 3 fields around the extremity of the Sgr trailing tail in the Southern Galactic hemisphere. We demonstrate that at the Sgr longitude $\tilde \Lambda_{\odot} = 204^{\circ}$, the line-of-sight velocity of the trailing tail starts to deviate from the track of the Law & Majewski (2010) model, confirming the prediction of Belokurov et al. (2014). This discovery serves to substantiate the measurement of low differential orbital precession of the Sgr stream which in turn may imply diminished dark matter content within 100 kpc.

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.

We explore the Triangulum-Andromeda (TriAnd) overdensity in the SPLASH (Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo) and SEGUE (the Sloan Extension for Galactic Understanding and Exploration) spectroscopic surveys. Milky Way main sequence turn-off stars in the SPLASH survey reveal that the TriAnd overdensity and the recently discovered PAndAS stream (Martin et al. 2014) share a common heliocentric distance (D ~ 20 kpc), position on the sky, and line-of-sight velocity (V_GSR ~ 50 km/s). Similarly, A-type, giant, and main sequence turn-off stars selected from the SEGUE survey in the vicinity of the Segue 2 satellite show that TriAnd is prevalent in these fields, with a velocity and distance similar to Segue 2. The coincidence of the PAndAS stream and Segue 2 satellite in positional and velocity space to TriAnd suggests that these substructures are all associated, and may be a fossil record of group-infall onto the Milky Way halo. In this scenario, the Segue 2 satellite and PAndAS stream are "satellites of satellites", and the large, metal-rich TriAnd overdensity is the remains of the group central.

In the rest frame of the Local Group (LG), the total momentum of the Milky Way (MW) and Andromeda (M31) should balance to zero. We use this fact to constrain new solutions for the solar motion with respect to the LG centre-of-mass, the total mass of the LG, and the individual masses of M31 and the MW. Using the set of remote LG galaxies at $>350$ kpc from the MW and M31, we find that the solar motion has amplitude $V_{\odot}=299\pm 15 {\rm ~km~s^{-1}}$ in a direction pointing toward galactic longitude $l_{\odot}=98.4^{\circ}\pm 3.6^{\circ}$ and galactic latitude $b_{\odot}=-5.9^{\circ}\pm 3.0^{\circ}$. The velocities of M31 and the MW in this rest frame give a direct measurement of their mass ratio, for which we find $\log_{10} (M_{\rm M31}/M_{\rm MW})=0.36 \pm 0.29$. We combine these measurements with the virial theorem to estimate the total mass within the LG as $M_{\rm LG}=(2.5\pm 0.4)\times 10^{12}~{\rm M}_{\odot}$. Our value for $M_{\rm LG}$ is consistent with the sum of literature values for $M_{\rm MW}$ and $M_{\rm M31}$. This suggests that the mass of the LG is almost entirely located within the two largest galaxies rather than being dispersed on larger scales or in a background medium. The outskirts of the LG are seemingly rather empty. Combining our measurement for $M_{\rm LG}$ and the mass ratio, we estimate the individual masses of the MW and M31 to be $M_{\rm MW}=(0.8\pm 0.5)\times 10^{12}~{\rm M}_{\odot}$ and $M_{\rm M31}=(1.7\pm 0.3)\times 10^{12}~{\rm M}_{\odot}$, respectively. Our analysis favours M31 being more massive than the MW by a factor of $\sim$2.3, and the uncertainties allow only a small probability (9.8%) that the MW is more massive. This is consistent with other properties such as the maximum rotational velocities, total stellar content, and numbers of globular clusters and dwarf satellites, which all suggest that $M_{\rm M31}/M_{\rm MW}>1$.

Knowledge of the abundance distribution of star forming regions and young clusters is critical to investigate a variety of issues, from triggered star formation and chemical enrichment by nearby supernova explosions to the ability to form planetary systems.In spite of this, detailed abundance studies are currently available for relatively few regions. In this context, we present the analysis of the metallicity of the Gamma Velorum cluster, based on the products distributed in the first internal release of the Gaia-ESO Survey. The Gamma Velorum candidate members have been observed with FLAMES, using both UVES and Giraffe, depending on the target brightness and spectral type. In order to derive a solid metallicity determination for the cluster, membership of the observed stars must be first assessed. To this aim, we use several membership criteria including radial velocities, surface gravity estimates, and the detection of the photospheric lithium line. Out of the 80 targets observed with UVES, we identify 14 high-probability members. We find that the metallicity of the cluster is slightly subsolar, with a mean [Fe/H]=-0.057+/-0.018 dex. Although J08095427-4721419 is one of the high-probability members, its metallicity is significantly larger than the cluster average. We speculate about its origin as the result of recent accretion episodes of rocky bodies of ~60 M_Sun hydrogen-depleted material from the circumstellar disk.

We present ~47,000 periodic variables found during the analysis of 5.4 million variable star candidates within a 20,000 square degree region covered by the Catalina Surveys Data Release-1 (CSDR1). Combining these variables with type-ab RR Lyrae from our previous work, we produce an on-line catalog containing periods, amplitudes, and classifications for ~61,000 periodic variables. By cross-matching these variables with those from prior surveys, we find that > 90% of the ~8,000 known periodic variables in the survey region are recovered. For these sources we find excellent agreement between our catalog and prior values of luminosity, period and amplitude, as well as classification. We investigate the rate of confusion between objects classified as contact binaries and type-c RR Lyrae (RRc's) based on periods, colours, amplitudes, metalicities, radial velocities and surface gravities. We find that no more than few percent of these variables in these classes are misidentified. By deriving distances for this clean sample of ~5,500 RRc's, we trace the path of the Sagittarius tidal streams within the Galactic halo. Selecting 146 outer-halo RRc's with SDSS radial velocities, we confirm the presence of a coherent halo structure that is inconsistent with current N-body simulations of the Sagittarius tidal stream. We also find numerous long-period variables that are very likely associated within the Sagittarius tidal streams system. Based on the examination of 31,000 contact binary light curves we find evidence for two subgroups exhibiting irregular lightcurves. One subgroup presents significant variations in mean brightness that are likely due to chromospheric activity. The other subgroup shows stable modulations over more than a thousand days and thereby provides evidence that the O'Connell effect is not due to stellar spots.

We present an algorithm for classifying the nearby transient objects detected by the Gaia satellite. The algorithm will use the low-resolution spectra from the blue and red spectro-photometers on board of the satellite. Taking a Bayesian approach we model the spectra using the newly constructed reference spectral library and literature-driven priors. We find that for magnitudes brighter than 19 in Gaia $G$ magnitude, around 75\% of the transients will be robustly classified. The efficiency of the algorithm for SNe type I is higher than 80\% for magnitudes $G\leq$18, dropping to approximately 60\% at magnitude $G$=19. For SNe type II, the efficiency varies from 75 to 60\% for $G\leq$18, falling to 50\% at $G$=19. The purity of our classifier is around 95\% for SNe type I for all magnitudes. For SNe type II it is over 90\% for objects with $G \leq$19. GS-TEC also estimates the redshifts with errors of $\sigma_z \le$ 0.01 and epochs with uncertainties $\sigma_t \simeq$ 13 and 32 days for type SNe I and SNe II respectively. GS-TEC has been designed to be used on partially calibrated Gaia data. However, the concept could be extended to other kinds of low resolution spectra classification for ongoing surveys.

We report the discovery of a narrow stellar stream crossing the constellations of Sculptor and Fornax in the Southern celestial hemisphere. The portion of the stream detected in the Data Release 1 photometry of the ATLAS survey is at least 12 degrees long, while its width is $\approx$ 0.25 deg. The Color Magnitude Diagram of this halo sub-structure is consistent with a metal-poor [Fe/H] $\lesssim -1.4$ stellar population located at a heliocentric distance of 20 $\pm$ 2 kpc. There are three globular clusters that could tentatively be associated with the stream: NGC 7006, NGC 7078 (M15) and Pyxis, but NGC 7006 and 7078 seem to have proper motions incompatible with the stream orbit.

We use A-type stars selected from Sloan Digital Sky Survey data release 9 photometry to measure the outer slope of the Milky Way stellar halo density profile beyond 50 kpc. A likelihood-based analysis is employed that models the ugr photometry distribution of blue horizontal branch (BHB) and blue straggler (BS) stars. In the magnitude range, 18.5 < g < 20.5, these stellar populations span a heliocentric distance range of: 10 kpc < D_BS < 75 kpc, 40 kpc < D_BHB < 100 kpc. Contributions from contaminants, such as QSOs, and the effect of photometric uncertainties, are also included in our modeling procedure. We find evidence for a very steep outer halo profile, with power-law index alpha ~ 6 beyond Galactocentric radii r=50 kpc, and even steeper slopes favored (alpha ~ 6-10) at larger radii. This result holds true when stars belonging to known overdensities, such as the Sagittarius stream, are included or excluded. We show that, by comparison to numerical simulations, stellar halos with shallower slopes at large distances tend to have more recent accretion activity. Thus, it is likely that the Milky Way has undergone a relatively quiet accretion history over the past several Gyr. Our measurement of the outer stellar halo profile may have important implications for dynamical mass models of the Milky Way, where the tracer density profile is strongly degenerate with total mass-estimates.

We study the relationship between age, metallicity, and alpha-enhancement of FGK stars in the Galactic disk. The results are based upon the analysis of high-resolution UVES spectra from the Gaia-ESO large stellar survey. We explore the limitations of the observed dataset, i.e. the accuracy of stellar parameters and the selection effects that are caused by the photometric target preselection. We find that the colour and magnitude cuts in the survey suppress old metal-rich stars and young metal-poor stars. This suppression may be as high as 97% in some regions of the age-metallicity relationship. The dataset consists of 144 stars with a wide range of ages from 0.5 Gyr to 13.5 Gyr, Galactocentric distances from 6 kpc to 9.5 kpc, and vertical distances from the plane 0 < |Z| < 1.5 kpc. On this basis, we find that i) the observed age-metallicity relation is nearly flat in the range of ages between 0 Gyr and 8 Gyr; ii) at ages older than 9 Gyr, we see a decrease in [Fe/H] and a clear absence of metal-rich stars; this cannot be explained by the survey selection functions; iii) there is a significant scatter of [Fe/H] at any age; and iv) [Mg/Fe] increases with age, but the dispersion of [Mg/Fe] at ages > 9 Gyr is not as small as advocated by some other studies. In agreement with earlier work, we find that radial abundance gradients change as a function of vertical distance from the plane. The [Mg/Fe] gradient steepens and becomes negative. In addition, we show that the inner disk is not only more alpha-rich compared to the outer disk, but also older, as traced independently by the ages and Mg abundances of stars.

We announce the discovery of a new Galactic companion found in data from the ESO VST ATLAS survey, and followed up with deep imaging on the 4m William Herschel Telescope. The satellite is located in the constellation of Crater (the Cup) at a distance of $\sim$ 170 kpc. Its half-light radius is $r_h=30$ pc and its luminosity is $M_V=-5.5$. The bulk of its stellar population is old and metal-poor. We would probably have classified the newly discovered satellite as an extended globular cluster were it not for the presence of a handful of Blue Loop stars and a sparsely populated Red Clump. The existence of the core helium burning population implies that star-formation occurred in Crater perhaps as recently as 400 Myr ago. No globular cluster has ever accomplished the feat of prolonging its star-formation by several Gyrs. Therefore, if our hypothesis that the blue bright stars in Crater are Blue Loop giants is correct, the new satellite should be classified as a dwarf galaxy with unusual properties. Note that only ten degrees to the North of Crater, two ultra-faint galaxies Leo IV and V orbit the Galaxy at approximately the same distance. This hints that all three satellites may once have been closely associated before falling together into the Milky Way halo.

NGC 4815 is a populous ~500 Myr open cluster at a Galactocentric radius of 7 kpc observed in the first six months of the Gaia-ESO Survey. Located in the inner Galactic disk, NGC 4815 is an important tracer of the abundance gradient, where few intermediate age open clusters are found. We use the survey derived radial velocities, stellar atmospheric parameters, metallicity, and elemental abundances for stars targeted as potential members of this cluster to carry out an analysis of cluster properties. The radial velocity distribution of stars in the cluster field is used to define the cluster systemic velocity and derive likely cluster membership for stars observed by the Gaia-ESO Survey. We investigate the distributions of Fe and Fe-peak elements, alpha-elements, and the light elements Na and Al and characterize the cluster's internal chemical homogeneity comparing it to the properties of radial velocity non-member stars. Utilizing these cluster properties, the cluster color-magnitude diagram is analyzed and theoretical isochrones are fit to derive cluster reddening, distance, and age. NGC 4815 is found to have a mean [Fe/H]=+0.03. Elemental abundances of cluster members show typically small internal variation, with internal dispersions of ~0.05 dex. [Ca/Fe] and [Si/Fe] show solar ratios, but [Mg/Fe] is moderately enhanced, while [Ti/Fe] appears slightly deficient. As with many open clusters, [Na/Fe] and [Al/Fe] are enhanced, [Na/Fe] significantly so, although the role of internal mixing and the assumption of LTE in the analysis remain to be investigated. From isochrone fits to color-magnitude diagrams, we find a cluster age of 0.5 to 0.63 Gyr, a reddening of E(B-V) = 0.59 to 0.65, and a true distance modulus (m-M) = 11.95 to 12.20, depending on the choice of theoretical models, leading to a Galactocentric distance of 6.9 kpc.

(Abridged) We have used the atmospheric parameters, [alpha/Fe] abundances and radial velocities, determined from the Gaia-ESO Survey GIRAFFE spectra of FGK-type stars (iDR1), to provide a chemo-kinematical characterisation of the disc stellar populations. We focuss on a subsample of 1016 stars with high quality parameters, covering the volume |Z|<4.5kpc and R in the range 2-13kpc. We have identified a thin to thick disc separation in the [alpha/Fe] vs [M/H] plane, thanks to the presence of a low-density region in the number density distribution. The thick disc stars seem to lie in progressively thinner layers above the Galactic plane, as metallicity increases and [alpha/Fe] decreases. The thin disc population presents a constant value of the mean distance to the plane at all metallicities. Our data confirm the already known correlations between V_phi and [M/H] for the two discs. For the thick disc sequence, a study of the possible contamination by thin disc stars suggests a gradient up to 64km/s/dex. The distributions of V_phi, V_Z, and orbital parameters are analysed for the chemically separated samples. Concerning the gradients with galactocentric radius, we find for the thin disc a flat behaviour of V_phi, a [M/H] gradient of -0.058dex/kpc and a small positive [alpha/Fe] gradient. For the thick disc, flat gradients in [M/H] and [alpha/Fe] are derived. Our chemo-kinematical analysis suggests a picture in which the thick disc seems to have experienced a settling process, during which its rotation increased progressively, and, possibly, the V_phi dispersion decreased. At [M/H]-0.25dex and [alpha/Fe]0.1dex, the mean characteristics of the thick disc in distance to the Galactic plane, V_phi, V_phi dispersion and eccentricity agree with those of the thin disc stars, suggesting a possible connection between these populations at a certain epoch of the disc evolution.

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.

Context: A key science goal of the Gaia-ESO survey (GES) is to use the kinematics of low-mass stars in young clusters to probe their dynamical histories and how they populate the field as they become unbound. The clustering of low-mass stars around the massive W-R binary gamma^2 Velorum was one of the first GES targets. Aims: To empirically determine the radial velocity (RV) precision of GES data, construct a kinematically unbiased sample of cluster members and characterise their dynamical state. Methods: Targets were selected from colour-magnitude diagrams and intermediate resolution spectroscopy used to derive RVs and assess membership from the strength of the Li6708A line. The RV distribution was analysed using a maximum likelihood technique that accounts for unresolved binaries. Results: The GES RV precision is about 0.25km/s and sufficient to resolve velocity structure in the low-mass population around gamma^2 Vel. The structure is well fitted by two kinematic components with roughly equal numbers of stars; the first has an intrinsic dispersion of 0.34+/-0.16km/s, consistent with virial equilibrium. The second has a broader dispersion of 1.60+/-0.37km/s and is offset from the first by ~2km/s. The first population is older by 1-2Myr based on a greater level of Li depletion seen among its M-stars and is probably more centrally concentrated around gamma^2 Vel. Conclusions: We consider several formation scenarios, concluding that the two kinematic components are a bound remnant of the original, denser cluster that formed gamma^2 Vel, and a dispersed population from the wider Vela OB2 association, of which gamma^2 Vel is the most massive member. The apparent youth of gamma^2 Vel compared to the older (>=10Myr) low-mass population surrounding it suggests a scenario where the massive binary formed in a clustered environment after the formation of the bulk of the low-mass stars.[abridged]

We map the large-scale sub-structure in the Galactic stellar halo using accurate 3D positions of ~14,000 RR Lyrae reported by the Catalina Sky Survey. In the heliocentric distance range of 10-25 kpc, in the region of the sky approximately bounded by 30{\deg} < l < 55{\deg} and -45{\deg} < b < -25{\deg}, there appears to be a strong excess of RRab stars. This overdensity, peaking at 18 kpc, is most likely associated with the so-called Hercules-Aquila Cloud, previously detected using Main Sequence tracers at similar distances in the Sloan Digital Sky Survey data. Our analysis of the period-amplitude distribution of RR Lyrae in this region indicates that the HAC is dominated by the Oosterhoff I type population. By comparing the measured RR Lyrae number density to models of a smooth stellar halo, we estimate the significance of the observed excess and provide an updated estimate of the total luminosity of the Cloud's progenitor.

Open clusters are key tools to study the spatial distribution of abundances in the disk and their evolution with time. Using the first release of stellar parameters and abundances of the Gaia-ESO Survey, we analyse the chemical properties of stars in three old/intermediate-age open clusters, namely NGC 6705, NGC 4815, and Trumpler 20, all located in the inner part of the Galactic disk at Galactocentric radius R$_{GC}\sim$7 kpc, aiming at proving their homogeneity and at comparing them with the field population. We study the abundance ratios of elements belonging to two different nucleosynthetic channels: $\alpha$-elements and iron-peak elements. The main results can be summarised as follows: i) cluster members are chemically homogeneous within 3-$\sigma$ in all analysed elements; ii) the three clusters have comparable [El/Fe] patters within $\sim$1-$\sigma$, but they differ in their global metal content [El/H], with NGC 4815 having the lowest metallicity. Their [El/Fe] ratios show differences and analogies with those of the field population, both in the solar neighbourhood and in the bulge/inner disk; iii) comparing the abundance ratios with the results of two chemical evolution models and with field star abundance distributions, we find that the abundance ratios of Mg, Ni, Ca in NGC 6705 might require an inner birthplace, implying a subsequent variation of its R$_{GC}$ during its lifetime, consistent with previous orbit determination. The full dataset of the Gaia-ESO Survey will be a superlative tool to constrain the chemical evolution of our Galaxy by disentangling different formation and evolution scenarios.

Trumpler 20 is an old open cluster (OC) located toward the Galactic centre, at about 3 kpc from the Sun and $\sim$7 kpc from the Galactic centre. Its position makes this cluster particularly interesting in the framework of the chemical properties of the Galactic disc because very few old OCs reside in the inner part of the disc. For this reason it has been selected as a cluster target of the Gaia-ESO Survey, and spectra of many stars in the main sequence and red clump phases are now available. Moreover, although it has been studied by several authors in the past, no consensus on the evolutionary status of Tr 20 has been reached. The heavy contamination of field stars (the line of sight of Tr 20 crosses the Carina spiral arm) complicates a correct interpretation. Another interesting aspect of the cluster is that it shows a broadened main-sequence turn-off and a prominent and extended red-clump, characteristics that are not easily explained by classical evolutionary models. Exploiting both spectroscopic information from the Gaia-ESO Survey (and the ESO archive) and literature photometry, we obtain a detailed and accurate analysis of the properties of the cluster. We make use of the first accurate metallicity measurement ever obtained from several spectra of red clump stars, and of cluster membership determination using radial velocities. According to the evolutionary models adopted, we find that Tr 20 has an age in the range 1.35-1.66 Gyr, an average reddening $E(B-V)$ in the range 0.31-0.35 mag, and a distance modulus $(m-M)_0$ between 12.64 and 12.72 mag. The spectroscopic metallicity is [Fe/H]=+0.17 dex. We discuss the structural properties of the object and constrain possible hypotheses for its broadened upper main sequence by estimating the effect of differential reddening and its extended red clump.

Using a variety of stellar tracers -- blue horizontal branch stars, main-sequence turn-off stars and red giants -- we follow the path of the Sagittarius (Sgr) stream across the sky in Sloan Digital Sky Survey data. Our study presents new Sgr debris detections, accurate distances and line-of-sight velocities that together help to shed new light on the puzzle of the Sgr tails. For both the leading and the trailing tail, we trace the points of their maximal extent, or apo-centric distances, and find that they lie at $R^L$ = 47.8 $\pm$ 0.5 kpc and $R^T$ = 102.5 $\pm$ 2.5 kpc respectively. The angular difference between the apo-centres is 93.2 $\pm$ 3.5 deg, which is smaller than predicted for logarithmic haloes. Such differential orbital precession can be made consistent with models of the Milky Way in which the dark matter density falls more quickly with radius. However, currently, no existing Sgr disruption simulation can explain the entirety of the observational data. Based on its position and radial velocity, we show that the unusually large globular cluster NGC 2419 can be associated with the Sgr trailing stream. We measure the precession of the orbital plane of the Sgr debris in the Milky Way potential and show that, surprisingly, Sgr debris in the primary (brighter) tails evolves differently to the secondary (fainter) tails, both in the North and the South.

We employ a large sample of 20171 optically-selected groups and clusters at 0.15 < z < 0.4 in the SDSS to investigate how the stacked stellar mass fraction varies across a wide range of total mass, $M_{500}$. Our study improves upon previous observational studies in a number of important ways, including having a much larger sample size, an explicit inclusion of the intracluster light (ICL) component, and a thorough examination of the accuracy of our total mass estimates via comparisons to simulations and weak lensing observations. We find that the stellar mass fraction depends only weakly on total mass and that the contribution of ICL to the total stellar mass fraction is significant (typically 20-40 per cent). Both of these findings are in excellent accordance with the predictions of cosmological simulations. Under the assumption of a Chabrier (Salpeter) IMF, the derived star formation efficiency ($f_{star}$/$f_{b}$, where $f_b=\Omega_b/\Omega_m$) is relatively low at 8 per cent (14 per cent) and is consistent with the global star formation efficiency of semi-analytic models that reproduce the galaxy stellar mass function. When our measured stellar mass fractions are combined with the observed relation between hot gas mass fraction and total mass from X-ray observations, our results imply that galaxy groups have significantly lower baryon fractions than massive clusters. Ejection of gas due to energetic AGN feedback (most likely at high redshift) provides a plausible mechanism for explaining the trends we observe.

We investigate the relationship between the halo mass, M_200, and concentration, c, for a sample of 26 group- and cluster-scale strong gravitational lenses. In contrast with previous results, we find that these systems are only ~ 0.1 dex more over-concentrated than similar-mass halos from dark matter simulations; the concentration of a halo with M_200 = 10^14 M_sun is log c = 0.78\pm0.05, while simulations of halos with this mass at similar redshifts (z ~ 0.4) predict log c ~ 0.56 - 0.71. We also find that we are unable to make informative inference on the slope of the M_200-c relation in spite of our large sample size; we note that the steep slopes found in previous studies tend to follow the slope in the covariance between M_200 and c, indicating that these results may be measuring the scatter in the data rather than the intrinsic signal. Furthermore, we conclude that our inability to constrain the M_200-c slope is due to a limited range of halo masses, as determined by explicitly modelling our halo mass distribution, and we suggest that other studies may be producing biased results by using an incorrect distribution for their halo masses.

We present the discovery of a plausible super-luminous supernova (SLSN), found in the archival data of Sloan Digital Sky Survey (SDSS) Stripe 82, called PSN 000123+000504. The supernova peaked at M_g<-21.3 mag in the second half of September 2005, but was missed by the real-time supernova hunt. The observed part of the light curve (17 epochs) showed that the rise to the maximum took over 30 days, while the decline time lasted at least 70 days (observed frame), closely resembling other SLSNe of SN2007bi type. Spectrum of the host galaxy reveals a redshift of z=0.281 and the distance modulus of \mu=40.77 mag. Combining this information with the SDSS photometry, we found the host galaxy to be an LMC-like irregular dwarf galaxy with the absolute magnitude of M_B=-18.2+/-0.2 mag and the oxygen abundance of 12+log[O/H]=8.3+/-0.2. Our SLSN follows the relation for the most energetic/super-luminous SNe exploding in low-metallicity environments, but we found no clear evidence for SLSNe to explode in low-luminosity (dwarf) galaxies only. The available information on the PSN 000123+000504 light curve suggests the magnetar-powered model as a likely scenario of this event. This SLSN is a new addition to a quickly growing family of super-luminous SNe.

Using multi-epoch observations of the Stripe 82 region done by Sloan Digital Sky Survey, we measure precise statistical proper motions of the stars in the Sagittarius stellar stream. The multi-band photometry and SDSS radial velocities allow us to efficiently select Sgr members and thus enhance the proper motion precision to ~0.1 mas/yr. We measure separately the proper motion of a photometrically selected sample of the main sequence turn-off stars, as well as of a spectroscopically selected Sgr giants. The data allow us to determine the proper motion separately for the two Sgr streams in the South found in Koposov et al.(2012). Together with the precise velocities from SDSS, our proper motion provide exquisite constraints of the 3-D motions of the stars in the Sgr streams.

The BVRcIc CCD photometry in the fields of six open clusters toward the Perseus spiral arm is presented. These data along with JHKs magnitudes taken from 2MASS catalog, have been used to determine cluster's ages, distances and color excesses. In addition, the gaps in mass function of Be 97, King 12 and NGC 7788 clusters have been revealed in mass intervals of 1.3-1.5, 1.4-1.6 and 1.5-1.7 solar masses, respectively.

We present the analysis of 12227 type-ab RR Lyrae found among the 200 million public lightcurves in the Catalina Surveys Data Release 1 (CSDR1). These stars span the largest volume of the Milky Way ever surveyed with RR Lyrae, covering ~20,000 square degrees of the sky (0 < RA < 360, -22 < Dec < 65 deg) to heliocentric distances of up to 60kpc. Each of the RR Lyrae are observed between 60 and 419 times over a six-year period. Using period finding and Fourier fitting techniques we determine periods and apparent magnitudes for each source. We find that the periods at generally accurate to sigma = 0.002% by comparison with 2842 previously known RR Lyrae and 100 RR Lyrae observed in overlapping survey fields. We photometrically calibrate the light curves using 445 Landolt standard stars and show that the resulting magnitudes are accurate to ~0.05 mags using SDSS data for ~1000 blue horizontal branch stars and 7788 of the RR Lyrae. By combining Catalina photometry with SDSS spectroscopy, we analyze the radial velocity and metallicity distributions for > 1500 of the RR Lyrae. Using the accurate distances derived for the RR Lyrae, we show the paths of the Sagittarius tidal streams crossing the sky at heliocentric distances from 20 to 60 kpc. By selecting samples of Galactic halo RR Lyrae, we compare their velocity, metallicity, and distance with predictions from a recent detailed N-body model of the Sagittarius system. We find that there are some significant differences between the distances and structures predicted and our observations.

We present the first sample of spectroscopically confirmed heavily reddened broad-line quasars selected using the new near infra-red VISTA Hemisphere Survey and WISE All-Sky Survey. Observations of four candidates with (J-K)>2.5 and K<=16.5 over ~180 sq deg, leads to confirmation that two are highly dust-reddened broad-line Type 1 quasars at z~2. The typical dust extinctions are AV ~ 2-2.5 mags. We measure black-hole masses of ~10^9 M0 and extinction corrected bolometric luminosities of ~10^47 erg/s, making these among the brightest Type 1 quasars currently known. Despite this, these quasars lie well below the detection limits of wide-field optical surveys like the SDSS with iAB > 22. We also present WISE photometry at 3--22um, for our full sample of spectroscopically confirmed reddened quasars including those selected from the UKIDSS Large Area Survey (Banerji et al. 2012a). We demonstrate that the rest-frame infrared SEDs of these reddened quasars are similar to UV-luminous Type 1 quasars with significant hot dust emission and starburst quasar hosts like Mrk231. The average 12um flux density of our reddened quasars is similar to that of the recently discovered HyLIRG WISE1814+3412 (z=2.452) at similar redshifts, with two of our reddened quasars also having comparable 22um flux densities to this extreme HyLIRG. These optically faint, heavily reddened broad-line quasars are therefore among the most mid infrared luminous galaxies at z~2, now being discovered using WISE.

Gaia is the cornerstone mission of the European Space Agency. From late 2013 it will start collecting superb astrometric, photometric and spectroscopic data for around a billion of stars of our Galaxy. While surveying the whole sky down to V=20mag Gaia will be detecting transients and anomalous behaviour of objects, providing near-real-time alerts to the entire astronomical community. Gaia should detected about 6000 supernovae, 1000 microlensing events and many other interesting types of transients. Thanks to its on-board low-dispersion spectrograph the classification of transients will be robust, assuring low false-alert rate. We describe the operation of the Photometric Science Alerts system, outline the scientific possibilities and conclude with an invitation to collaborate in the ground-based follow-up Gaia alerts during the early months of the mission when the outcome of the alerting pipeline needs to be verified.

The aim of the project is to improve our current knowledge of the density of T dwarfs and the shape of the substellar initial mass function by identifying a magnitude-limited sample of T dwarfs in the full southern sky. We present the results of a photometric search aimed at discovering cool brown dwarfs in the Southern sky imaged at infrared wavelengths by the Visible and Infrared Survey Telescope for Astronomy (VISTA) and the Wide Infrared Survey Explorer (WISE) satellite mission. We combined the first data release (DR1) of the VISTA Hemisphere Survey (VHS) and the WISE preliminary data release to extract candidates with red mid-infrared colours and near- to mid-infrared colours characteristics of cool brown dwarfs. The VHS DR1 vs WISE search returned tens of T dwarf candidates, 13 of which are presented here, including two previously published in the literature and five new ones confirmed spectroscopically with spectral types between T4.5 and T8. We estimate that the two T6 dwarfs lie within 16 pc and the T4.5 within 25 pc. The remaining three are 30-50 pc distant. The only T7 dwarf in our sample is the faintest of its spectral class with J=19.28 mag. The other six T dwarf candidates remain without spectroscopic follow-up. We also improve our knowledge on the proper motion accuracy for three bright T dwarfs by combining multi-epoch data from public databases (DENIS, 2MASS, VHS, WISE, Spitzer).

We calculate the energy that baryons must inject in cold dark matter (CDM) haloes in order to remove centrally-divergent DM cusps on scales relevant to observations of dwarf spheroidal galaxies (dSphs). We estimate that the CDM haloes often associated with the Milky Way's dSphs (M_vir/M_\odot \sim 10^{9-10}) require \Delta E/erg \sim 10^{53-55} in order to form cores on scales comparable to the luminous size of these galaxies. While supernova type II (SNeII) explosions can in principle generate this energy, the actual contribution is limited by the low star formation efficiency implied by the abundance of luminous satellites. Considering that CDM's well-known `core/cusp' and `missing satellite' problems place opposing demands on star formation efficiencies, existing observational evidences for large cores in the most luminous dSphs require that CDM models invoke some combination of the following: (i) efficient (of order unity) coupling of SNeII energy into dark matter particles, (ii) star formation histories peaking at unexpectedly high redshifts (z>6), (iii) a top-heavy stellar IMF, and/or (iv) substantial satellite disruption or other stochastic effects to ease the substructure abundance constraints. Our models show that the tension between CDM problems on small scales would increase if cored DM profiles were to be found in fainter dwarves.

We introduce the Minimum Entropy Method, a simple statistical technique for constraining the Milky Way gravitational potential and simultaneously testing different gravity theories directly from 6D phase-space surveys and without adopting dynamical models. We demonstrate that orbital energy distributions that are separable (i.e. independent of position) have an associated entropy that increases under wrong assumptions about the gravitational potential and/or gravity theory. Of known objects, `cold' tidal streams from low-mass progenitors follow orbital distributions that most nearly satisfy the condition of separability. Although the orbits of tidally stripped stars are perturbed by the progenitor's self-gravity, systematic variations of the energy distribution can be quantified in terms of the cross-entropy of individual tails, giving further sensitivity to theoretical biases in the host potential. The feasibility of using the Minimum Entropy Method to test a wide range of gravity theories is illustrated by evolving restricted N-body models in a Newtonian potential and examining the changes in entropy introduced by Dirac, MONDian and f(R) gravity modifications.

By undertaking deep long-slit spectroscopy with the focal reducer SCORPIO of the Russian 6m telescope, we studied stellar population properties and their variation with radius in 15 nearby S0 galaxies sampling a wide range of luminosities and environments. For the large-scale stellar disks of S0s, we have measured SSP-equivalent metallicities ranging from the solar one down to [Z/H]=-0.4 - -0.7, rather high magnesium-to-iron ratios, [Mg/Fe] > +0.2, and mostly old SSP-equivalent ages. Nine of 15 (60%) galaxies have large-scale stellar disks older than 10 Gyr, and among those we find all the galaxies which reside in denser environments. The isolated galaxies may have intermediate-age stellar disks which are 7-9 Gyr old. Only two galaxies of our sample, NGC 4111 and NGC 7332, reveal SSP-equivalent ages of their disks of 2-3 Gyrs. Just these two young disks appear to be thin, while the other, older disks have scale heights typical for thick stellar disks. The stellar populations in the bulges at radii of 0.5r_eff are on the contrary more metal-rich than the solar metallicity, with the ages homogeneously distributed between 2 and 15 Gyr, being almost always younger than the disks. We conclude that S0 galaxies could not form in groups at z=0.4 as is thought now; a new scenario of the general evolution of disk galaxies is proposed instead.