The Nancy Grace Roman Space Telescope (Roman) will perform a Galactic Exoplanet Survey (RGES) to discover bound exoplanets with semi-major axes greater than 1 au using gravitational microlensing. Roman will even be sensitive to planetary mass objects that are not gravitationally bound to any host star. Such free-floating planetary mass objects (FFPs) will be detected as isolated microlensing events with timescales shorter than a few days. A measurement of the abundance and mass function of FFPs is a powerful diagnostic of the formation and evolution of planetary systems, as well as the physics of the formation of isolated objects via direct collapse. We show that Roman will be sensitive to FFP lenses that have masses from that of Mars ($0.1 M_\oplus$) to gas giants ($M\gtrsim100M_\oplus$) as isolated lensing events with timescales from a few hours to several tens of days, respectively. We investigate the impact of the detection criteria on the survey, especially in the presence of finite-source effects for low-mass lenses. The number of detections will depend on the abundance of such FFPs as a function of mass, which is at present poorly constrained. Assuming that FFPs follow the fiducial mass function of cold, bound planets adapted from Cassan et al. (2012), we estimate that Roman will detect $\sim250$ FFPs with masses down to that of Mars (including $\sim 60$ with masses $\le M_\oplus$). We also predict that Roman will improve the upper limits on FFP populations by at least an order of magnitude compared to currently-existing constraints.

Galactic microlensing datasets now comprise in excess of $10^4$ events, and with the advent of next generation microlensing surveys that may be undertaken with facilities such as the Vera Rubin Observatory (formerly LSST) and Nancy Grace Roman Space Telescope (formerly WFIRST), this number will increase significantly. So too will the fraction of events with measurable higher order information such as finite source effects and lens-source relative proper motion. Analysing such data requires a more sophisticated Galactic microlens modeling approach. We present a new second-generation Manchester-Besan\c{c}on Microlensing Simulator (MaB$\mu$lS-2), which uses a version of the Besan\c{c}on population synthesis Galactic model that provides good agreement with stellar kinematics observed by HST towards the bulge. MaB$\mu$lS-2 provides high-fidelity signal-to-noise limited maps of the microlensing optical depth, rate and average timescale towards a 400 sq. degree region of the Galactic bulge in several optical to near-infrared pass-bands. The maps take full account of the unresolved stellar background as well as limb-darkened source profiles. Comparing MaB$\mu$lS-2 to the efficiency-corrected OGLE-IV 8,000 event sample shows a much improved agreement over the previous version of MaB$\mu$lS, and succeeds in matching even small-scale structural features in the OGLE-IV event rate map. However, there remains evidence for a small under-prediction in the event rate per source and over-prediction in timescale. MaB$\mu$lS-2 is available online (\<www.mabuls.net>) to provide on-the-fly maps for user supplied cuts in survey magnitude, event timescale and relative proper motion.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We report the discovery of a unique collection of metal-poor giant-stars, that exhibit anomalously high levels of $^{28}$Si, clearly above typical Galactic levels. Our sample spans a narrow range of metallicities, peaking at $-1.07\pm 0.06$, and exhibit abundance ratios of [Si,Al/Fe] that are as extreme as those observed in Galactic globular clusters (GCs), and Mg is slightly less overabundant. In almost all the sources we used, the elemental abundances were re-determined from high-resolution spectra, which were re-analyzed assuming LTE. Thus, we compiled the main element families, namely the light elements (C, N), $\alpha-$elements (O, Mg, Si), iron-peak element (Fe), $\textit{s}-$process elements (Ce, Nd), and the light odd-Z element (Al). We also provide dynamical evidence that most of these stars lie on tight (inner)halo-like and retrograde orbits passing through the bulge. Such kinds of objects have been found in present-day halo GCs, providing the clearest chemical signature of past accretion events in the (inner) stellar halo of the Galaxy, formed possibly as the result of dissolved halo GCs. Their chemical composition is, in general, similar to that of typical GCs population, although several differences exist.

We report the identification of possible extended star debris candidates beyond the cluster tidal radius of NGC 6362 based on the second {\it Gaia} data release ({\it Gaia} DR2). We found 259 objects possibly associated with the cluster lying in the vicinity of the giant branch and 1--2 magnitudes fainter/brighter than the main-sequence turn-off in the cluster color-magnitude diagram and which cover an area on the sky of $\sim$4.1 deg$^{2}$ centered on the cluster. We traced back the orbit of NGC 6362 in a realistic Milky-Way potential, using the \texttt{GravPot16} package, for 3 Gyrs. The orbit shows that the cluster shares similar orbital properties as the inner disk, having peri-/apo-galactic distances, and maximum vertical excursion from the Galactic plane inside the corotation radius (CR), moving inwards from CR radius to visit the inner regions of the Milky Way. The dynamical history of the cluster reveals that it has crossed the Galactic disk several times in its lifetime and has recently undergone a gravitational shock, $\sim 15.9$ Myr ago, suggesting that less than 0.1\% of its mass has been lost during the current disk-shocking event. Based on the cluster's orbit and position in the Galaxy, we conclude that the possible extended star debris candidates are a combined effect of the shocks from the Galactic disk and evaporation from the cluster. Lastly, the evolution of the vertical component of the angular momentum shows that the cluster is strongly affected dynamically by the Galactic bar potential.

We report the serendipitous discovery of a nitrogen-rich, mildly metal-poor ([Fe/H]=-1.08) giant star in a single-lined spectroscopic binary system found in the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) survey, Data Release 14 (DR14). Previous work has assumed that the two percent of halo giants with unusual elemental abundances have been evaporated from globular clusters, but other origins for their abundance signatures, including binary mass transfer, must also be explored. We present the results of an abundance re-analysis of the APOGEE-2 high-resolution near-infrared spectrum of 2M12451043+1217401 with the Brussels Automatic Stellar Parameter (BACCHUS) automated spectral analysis code, and re-derive manually the main element families, namely the light elements (C, N), elements (O, Mg, Si), iron-peak element (Fe), \textit{s}-process element (Ce), and the light odd-Z element (Al). Our analysis confirm the N-rich nature of 2M12451043+1217401, which has a [N/Fe] ratio of $+0.69$, and shows that the abundances of C and Al are slightly discrepant from that of a typical mildly metal-poor RGB star, but exhibit Mg, Si, O and \textit{s}-process abundances (Ce) of typical field stars. We also detect a particularly large variability in its radial velocity over the period of the APOGEE-2 observations, and the most likely orbit fit to the radial velocity data has a period of $730.89\pm106.86$ days, a velocity semi-amplitude of $9.92 \pm 0.14$ km s$^{-1}$, and an eccentricity of $\sim 0.1276 \pm0.1174$, which support the hypothesis of a binary companion, and that has probably been polluted by a now-extinct AGB star.

We use Gaia DR2 magnitudes, colours and parallaxes for stars with G<12 to explore a 15-dimensional space that includes simultaneously the initial mass function (IMF) and a non-parametric star formation history (SFH) for the Galactic disc. This inference is performed by combining the Besancon Galaxy Model fast approximate simulations (BGM FASt) and an approximate Bayesian computation algorithm. We find in Gaia DR2 data an imprint of a star formation burst 2-3 Gyr ago, in the Galactic thin disc domain, and a present star formation rate (SFR) of about 1 Msun. Our results show a decreasing trend of the SFR from 9-10 Gyr to 6-7 Gyr ago. This is consistent with the cosmological star formation quenching observed at redshifts z<1.8. This decreasing trend is followed by a SFR enhancement starting at about 5 Gyr ago and continuing until about 1 Gyr ago which is detected with high statistical significance by discarding the null hypothesis of an exponential SFH with a p-value=0.002. We estimate, from our best fit model, that about the 50% of the mass used to generate stars, along the thin disc life, was expended in the period from 5 to 1 Gyr ago. The timescale and the amount of stellar mass generated during the SFR enhancement event lead us to hypothesise that its origin, currently under investigation, is not intrinsic to the disc. Thus, an external perturbation is needed for its explanation. Additionally we find, for the thin disc, a slope of the IMF of $\alpha_3 \approx 2$ for masses M>1.53 Msun and $\alpha_2 \approx 1.3$ for the mass range between 0,5 and 1,53 Msun. This is the first time that we consider a non-parametric SFH for the thin disc in the Besancon Galaxy Model. This new step, together with the capabilities of the Gaia DR2 parallaxes to break degeneracies between different stellar populations, allow us to better constrain the SFH and the IMF.

Connecting star and planet properties in a single model is not straightforward. Stellar population synthesis models are key to explore combined statistical constraints from stars and planets observations. The Besan\c{c}on stellar population synthesis model (Robin et al. 2003, Lagarde et al. 2017) includes now the stellar evolutionary tracks computed with the stellar evolution code STAREVOL (Lagarde et al. 2012, Amard et al. 2016). It provides the global (M, R, Teff, etc) and chemical properties of stars for 54 chemical species. It enables to study the different galactic populations of the Milky Way (the halo, the bulge, the thin and thick disc) and a specific observational survey. Here, we couple the Besan\c{c}on model with a simple stoichiometric model (Santos et al. 2017) in order to determine the expected composition of the planet building blocks (PBB). We investigate the trends and correlations of the expected chemical abundances of PBB in the different stellar populations of the Milky Way (Cabral et al. 2018).

The GES survey using FLAMES at the VLT has obtained high-resolution UVES spectra for a large number of giant stars, allowing a determination of the abundances of the key chemical elements C and N at their surface. The surface abundances of these chemical species are well-known to change in stars during their evolution on the red giant branch after the first dredge-up episod, as a result of extra-mixing phenomena. We investigate the effects of thermohaline mixing on C and N abundances using the first comparison between the GES [C/N] determinations with simulations of the observed fields using a model of stellar population synthesis. We explore the effects of thermohaline mixing on the chemical properties of giants through stellar evolutionary models computed with the stellar evolution code STAREVOL. We include these stellar evolution models in the Besan\c{c}on Galaxy model to simulate the [C/N] distributions determined from the UVES spectra of the GES and compare them with the observations. Theoretical predictions including the effect of thermohaline mixing are in good agreement with the observations. However, the field stars in the GES with C and N-abundance measurements have a metallicity close to solar, where the efficiency of thermohaline mixing is not very large. The C and N abundances derived by the GES in open and globular clusters clearly show the impact of thermohaline mixing at low-metallicity, allowing to explain the [C/N] ratio observed in lower-mass and older giant stars. Using independent observations of carbon isotopic ratio in clump field stars and open clusters, we also confirm that thermohaline mixing should be taken into account to explain the behavior of 12C/13C ratio as a function of stellar age. Overall the current model including thermohaline mixing is able to reproduce very well the C- and N-abundances over the whole metallicity range investigated by the GES data.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

Future space missions will improve considerably our understanding of the formation and history of planetary systems. Currently, observations show that the presence of planetary companions is closely linked to the metallicity and the chemical abundances of the host stars. We aim to build an integrated tool to predict the planet building blocks composition as a function of the stellar populations, for the interpretation of the ongoing and future large surveys. We synthesize stellar populations with the Besan\c{c}on Galaxy model (BGM) which includes stellar evolutionary tracks computed with the stellar evolution code STAREVOL. We integrate to the BGM a simple stoichiometric model to determine the expected composition of the planet building blocks. We determine the expected PBB composition around FGK stars, for the four galactic populations within the Milky Way. Our solar neighborhood simulations are in good agreement with the recent results obtained with the HARPS survey for f_iron, f_w and the heavy mass fraction f_Z. We present evidence of the clear dependence of f_iron and f_w with the initial alpha abundances [\alpha/Fe] of the host star. We find that the different initial [\alpha/Fe] distributions in the different galactic populations lead to a bimodal distribution of PBB composition and to an iron/water valley separating PBB with high and low iron/water mass fractions. We linked host star abundances and expected PBB composition in an integrated model of the Galaxy. Derived trends are an important step for statistical analyses of expected planet properties. In particular, internal structure models may use these results to derive statistical trends of rocky planets properties, constrain habitability and prepare interpretation of on-going and future large scale surveys of exoplanet search.

Context: The Initial Mass Function (IMF) plays a crucial role on galaxy evolution and its implications on star formation theory make it a milestone for the next decade. It is in the intermediate and high mass ranges where the uncertainties of the IMF are larger. This is a major subject of debate and analysis both for Galactic and extragalactic science. Aims: Our goal is to constrain the IMF of the Galactic thin disc population using both Galactic Classical Cepheids and Tycho-2 data. Methods: For the first time the Besan\c{c}on Galaxy Model (BGM) has been used to characterise the Galactic population of the Classical Cepheids. We have modified the age configuration in the youngest populations of the BGM thin disc model to avoid artificial discontinuities in the age distribution of the simulated Cepheids. Three statistical methods, optimized for different mass ranges, have been developed and applied to search for the best IMF that fits the observations. This strategy allows us to quantify variations in the Star Formation History (SFH), the stellar density at Sun position and the thin disc radial scale length. A rigorous treatment of unresolved multiple stellar systems has been undertaken adopting a spatial resolution according to the catalogues used. Results: For intermediate masses, our study favours a composite field-star IMF slope of $\alpha=3.2$ for the local thin disc, excluding flatter values such as the Salpeter IMF ($\alpha=2.35$). Moreover, a constant Star Formation History is definitively excluded, the three statistical methods considered here show that it is inconsistent with the observational data. Conclusions: Using field stars and Galactic Classical Cepheids, we have found, above $1M_\odot$, an IMF steeper than the canonical stellar IMF of associations and young clusters. This result is consistent with the predictions of the Integrated Galactic IMF.

We develop a new theoretical framework to generate Besan\c{c}on Galaxy Model fast approximate simulations (BGM FASt) to address fundamental questions of the Galactic structure and evolution performing multi-parameter inference. As a first application of our strategy we simultaneously infer the IMF, the star formation history and the stellar mass density in the Solar Neighbourhood. The BGM FASt strategy is based on a reweighing scheme, that uses a specific pre-sampled simulation, and on the assumption that the distribution function of the generated stars in the Galaxy can be described by an analytical expression. To validate BGM FASt we execute a set of tests. Finally, we use BGM FASt with an approximate Bayesian computation algorithm to obtain the posterior PDF of the inferred parameters, by comparing synthetic versus Tycho-2 colour-magnitude diagrams. Results: The validation shows a very good agreement between BGM FASt and the standard BGM, with BGM FASt being $\approx 10^4$ times faster. By analysing Tycho-2 data we obtain a thin disc star formation history decreasing in time and a present rate of $1.2 \pm 0.2 M_\odot/yr$. The resulting total stellar mass density in the Solar Neighbourhood is $0.051_{-0.005}^{+0.002} M_\odot/pc^3$ and the local dark matter density is $0.012 \pm 0.001 M_\odot/pc^3$. For the composite IMF we obtain a slope of $\alpha_2={2.1}_{-0.3}^{+0.1}$ in the mass range between $0.5 M_\odot$ and $1.53M_\odot$. The results of the slope at the high mass range are trustable up to $4M_\odot$ and highly depend on the choice of the extinction map (obtaining $\alpha_3={2.9}_{-0.2}^{+0.2}$ and $\alpha_3={3.7}_{-0.2}^{+0.2}$ respectively, for two different extinction maps). Systematic uncertainties are not included. Conclusions: The good performance of BGM FASt demonstrates that it is a very valuable tool to perform multi-parameter inference using Gaia data releases.

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

We present the stellar density profile of the outer halo of the Galaxy traced over a range of Galactocentric radii from $15< R_{GC} < 220$ kpc by blue horizontal branch (BHB) stars. These stars are identified photometrically using deep $u-$band imaging from the new Canada-France-Imaging-Survey (CFIS) that reaches 24.5 mag. This is combined with $griz$ bands from Pan-STARRS 1 and covers a total of $\sim4000$ deg$^2$ of the northern sky. We present a new method to select BHB stars that has low contamination from blue stragglers and high completeness. We use this sample to measure and parameterize the three dimensional density profile of the outer stellar halo. We fit the profile using (i) a simple power-law with a constant flattening (ii) a flattening that varies as a function of Galactocentric radius (iii) a broken power law profile. We find that outer stellar halo traced by the BHB is well modelled by a broken power law with a constant flattening of $q=0.86 \pm 0.02$, with an inner slope of $\gamma=4.24 \pm 0.08$. This is much steeper than the preferred outer profile that has a slope of $\beta=3.21\pm 0.07$ after a break radius of $r_b=41.4^{+2.5}_{-2.4}$ kpc. The outer profile of the stellar halo trace by BHB stars is shallower than that recently measured using RR Lyrae, a surprising result given the broad similarity of the ages of these stellar populations.

Context. Dynamically self-consistent galactic models are necessary for analysing and interpreting star counts, stellar density distributions, and stellar kinematics in order to understand the formation and the evolution of our Galaxy. Aims. We modify and improve the dynamical self-consistency of the Besan\c{c}on Galaxy model in the case of a stationary and axisymmetric gravitational potential. Methods. Each stellar orbit is modelled by determining a Staeckel approximate integral of motion. Generalised Shu distribution functions with three integrals of motion are used to model the stellar distribution functions. Results. This new version of the Besan\c{c}on model is compared with the previous axisymmetric BGM2014 version and we find that the two versions have similar densities for each stellar component. The dynamically self-consistency is improved and can be tested by recovering the forces and the potential through the Jeans equations applied to each stellar distribution function. Forces are recovered with an accuracy better than one per cent over most of the volume of the Galaxy.

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

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.

The Wide Field InfraRed Survey Telescope (WFIRST) is the next NASA astrophysics flagship mission, to follow the James Webb Space Telescope (JWST). The WFIRST mission was chosen as the top-priority large space mission of the 2010 astronomy and astrophysics decadal survey in order to achieve three primary goals: to study dark energy via a wide-field imaging survey, to study exoplanets via a microlensing survey, and to enable a guest observer program. Here we assess the ability of the several WFIRST designs to achieve the goal of the microlensing survey to discover a large sample of cold, low-mass exoplanets with semimajor axes beyond roughly one AU, which are largely impossible to detect with any other technique. We present the results of a suite of simulations that span the full range of the proposed WFIRST architectures, from the original design envisioned by the decadal survey, to the current design, which utilizes a 2.4-m telescope donated to NASA. By studying such a broad range of architectures, we are able to determine the impact of design trades on the expected yields of detected exoplanets. In estimating the yields we take particular care to ensure that our assumed Galactic model predicts microlensing event rates that match observations, consider the impact that inaccuracies in the Galactic model might have on the yields, and ensure that numerical errors in lightcurve computations do not bias the yields for the smallest mass exoplanets. For the nominal baseline WFIRST design and a fiducial planet mass function, we predict that a total of ${\sim}1400$ bound exoplanets with mass greater than ${\sim}0.1~M_{\oplus}$ should be detected, including ${\sim}200$ with mass ${\lesssim}3~M_{\oplus}$. WFIRST should have sensitivity to planets with mass down to ${\sim}0.02~M_{\oplus}$, or roughly the mass of Ganymede.

IC 166 is an intermediate-age open cluster ($\sim 1$ Gyr) which lies in the transition zone of the metallicity gradient in the outer disc. Its location, combined with our very limited knowledge of its salient features, make it an interesting object of study. We present the first high-resolution spectroscopic and precise kinematical analysis of IC 166, which lies in the outer disc with $R_{GC} \sim 12.7$ kpc. High resolution \textit{H}-band spectra were analyzed using observations from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We made use of the Brussels Automatic Stellar Parameter (BACCHUS) code to provide chemical abundances based on a line-by-line approach for up to eight chemical elements (Mg, Si, Ca, Ti, Al, K, Mn and Fe). The $\alpha-$element (Mg, Si, Ca and whenever available Ti) abundances, and their trends with Fe abundances have been analysed for a total of 13 high-likelihood cluster members. No significant abundance scatter was found in any of the chemical species studied. Combining the positional, heliocentric distance, and kinematic information we derive, for the first time, the probable orbit of IC 166 within a Galactic model including a rotating boxy bar, and found that it is likely that IC 166 formed in the Galactic disc, supporting its nature as an unremarkable Galactic open cluster with an orbit bound to the Galactic plane.

The goal of this paper is to demonstrate the outstanding quality of the second data release of the Gaia mission and its power for constraining many different aspects of the dynamics of the satellites of the Milky Way. We focus here on determining the proper motions of 75 Galactic globular clusters, nine dwarf spheroidal galaxies, one ultra-faint system, and the Large and Small Magellanic Clouds. Using data extracted from the Gaia archive, we derived the proper motions and parallaxes for these systems, as well as their uncertainties. We demonstrate that the errors, statistical and systematic, are relatively well understood. We integrated the orbits of these objects in three different Galactic potentials, and characterised their properties. We present the derived proper motions, space velocities, and characteristic orbital parameters in various tables to facilitate their use by the astronomical community. Our limited and straightforward analyses have allowed us for example to (i) determine absolute and very precise proper motions for globular clusters; (ii) detect clear rotation signatures in the proper motions of at least five globular clusters; (iii) show that the satellites of the Milky Way are all on high-inclination orbits, but that they do not share a single plane of motion; (iv) derive a lower limit for the mass of the Milky Way of 9.8^{+6.7}_{-2.7} x 10^{11} Msun based on the assumption that the Leo I dwarf spheroidal is bound; (v) derive a rotation curve for the Large Magellanic Cloud based solely on proper motions that is competitive with line-of-sight velocity curves, now using many orders of magnitude more sources; and (vi) unveil the dynamical effect of the bar on the motions of stars in the Large Magellanic Cloud. All these results highlight the incredible power of the Gaia astrometric mission, and in particular of its second data release.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

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

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

BANYAN $\Sigma$ is a new Bayesian algorithm to identify members of young stellar associations within 150 pc of the Sun. It includes 27 young associations with ages in the range ~1-800 Myr, modelled with multivariate Gaussians in 6-dimensional XYZUVW space. It is the first such multi-association classification tool to include the nearest sub-groups of the Sco-Cen OB star-forming region, the IC 2602, IC 2391, Pleiades and Platais 8 clusters, and the $\rho$ Ophiuchi, Corona Australis, and Taurus star-formation regions. A model of field stars is built from a mixture of multivariate Gaussians based on the Besan\c{c}on Galactic model. The algorithm can derive membership probabilities for objects with only sky coordinates and proper motion, but can also include parallax and radial velocity measurements, as well as spectrophotometric distance constraints from sequences in color-magnitude or spectral type-magnitude diagrams. BANYAN $\Sigma$ benefits from an analytical solution to the Bayesian marginalization integrals that makes it more accurate and significantly faster than its predecessor BANYAN II. A contamination versus hit rate analysis is presented and demonstrates that BANYAN $\Sigma$ achieves a better classification performance than other moving group classification tools, especially in terms of cross-contamination between young associations. An updated list of bona fide members in the 27 young associations, augmented by the Gaia-DR1 release, are presented. This new tool will make it possible to analyze large data sets such as the upcoming Gaia-DR2 to identify new young stars. IDL and Python versions of BANYAN $\Sigma$ are made available with this publication. (shortened)

Metal-poor globular clusters (GCs) exhibit intriguing Al-Mg anti-correlations and possible Si-Al correlations, which are important clues to decipher the multiple-population phenomenon. NGC 5053 is one of the most metal-poor GCs in the nearby Universe, and has been suggested to be associated with the Sagittarius (Sgr) dwarf galaxy, due to its similarity in location and radial velocity with one of the Sgr arms. In this work, we simulate the orbit of NGC 5053, and argue against a physical connection between Sgr and NGC 5053. On the other hand, the Mg, Al, and Si spectral lines, which are difficult to detect in the optical spectra of NGC 5053 stars, have been detected in the near-infrared APOGEE spectra. We use three different sets of stellar parameters and codes to derive the Mg, Al, and Si abundances. Regardless of which method is adopted, we see a large Al variation, and a substantial Si spread. Along with NGC 5053, metal-poor GCs exhibit different Mg, Al, and Si variations. Moreover, NGC 5053 has the lowest cluster mass among the GCs that have been identified to exhibit an observable Si spread until now.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) provides spectroscopic information of regions of the inner Milky Way inaccessible to optical surveys. We present the first large study of the metallicity distribution of the innermost Galactic regions based on homogeneous measurements from the SDSS Data Release 12 for 7545 red giant stars within 4.5 kpc of the Galactic center, with the goal to shed light on the structure and origin of the Galactic bulge. Stellar metallicities are found, through multiple-Gaussian decompositions, to be distributed in several components indicative of the presence of various stellar populations such as the bar, or the thin and the thick disk. A super-solar ([Fe/H]=+0.32) and a solar ([Fe/H]=+0.00) metallicity components, tentatively associated with the thin disk and the Galactic bar, respectively, seem to be the major contributors near the midplane. The solar-metallicity component extends outwards in the midplane but is not observed in the innermost regions. The central regions (within 3 kpc of the Galactic center) reveal, on the other hand, the presence of a significant metal-poor population ([Fe/H] =-0.46), tentatively associated with the thick disk, and which becomes the dominant component far from the midplane ($|Z| > +0.75 kpc). Varying contributions from these different components produce a transition region at +0.5 kpc < |Z| < +1.0 kpc characterized by a significant vertical metallicity gradient.

By performing an orbital analysis within a Galactic model including a bar, we found that it is plausible that the newly discovered stars that show enhanced Al and N accompanied by Mg underabundances may have formed in the outer halo, or were brought in by satellites field possibly accreted a long time ago. However, another subsample of three N- and Al-rich stars with Mg-deficiency are kinematically consistent with the inner stellar halo. A speculative scenario to explain the origin of the atypical chemical composition of these stars in the inner halo is that they migrated to the inner stellar halo as unbound stars due to the mechanism of bar-induced resonant trapping.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm deg^2}$ of $u$-band photometry taken as part of the Canada-France Imaging Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the Northern sky. By combing the CFHT $u$-band photometry together with SDSS and Pan-STARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably the metallicities of individual stars to $\sim 0.2$ dex, and hence additionally obtain good photometric distance estimates. This survey thus permits the measurement of metallicities and distances of the dominant main-sequence population out to approximately 30 kpc, and provides much higher number of stars at large extraplanar distances than have been available from previous surveys. We develop a non-parametric distance-metallicity decomposition algorithm and apply it to the sky at $30\deg < |b| < 70\deg$ and to the North Galactic Cap. We find that the metallicity-distance distribution is well-represented by three populations whose metallicity distributions do not vary significantly with vertical height above the disk. As traced in main-sequence stars, the stellar halo component shows a vertical density profile that is close to exponential, with a scale height of around 3 kpc. This may indicate that the inner halo was formed partly from disk stars ejected in an ancient minor merger.

The Canada-France Imaging Survey (CFIS) will map the northern high Galactic latitude sky in the $u$-band ("CFIS-u", 10,000$\, {\rm deg^2}$) and in the $r$-band ("CFIS-r", 5,000$\, {\rm deg^2}$), enabling a host of stand-alone science investigations, and providing some of the ground-based data necessary for photometric redshift determination for the Euclid mission. In this first contribution we present the $u$-band component of the survey, describe the observational strategy, and discuss some first highlight results, based on approximately one third of the final area. We show that the Galactic anticenter structure is distributed continuously along the line of sight, out to beyond 20 kpc, and possesses a metallicity distribution that is essentially identical to that of the outer disk sampled by APOGEE. This suggests that it is probably a buckled disk of old metal-rich stars, rather than a stream or a flare. We also discuss the future potential for CFIS-u in discovering star-forming dwarf galaxies around the Local Group, the characterization of the white dwarf and blue straggler population of the Milky Way, as well as its sensitivity to low-surface brightness structures in external galaxies.

We provide new insight on the origin of the cold high-V$_{\rm los}$ peaks ($\sim$200 kms$^{-1}$) in the Milky Way bulge discovered in the APOGEE commissioning data \citep{Nidever2012}. Here we show that such kinematic behaviour present in the field regions towards the Galactic bulge is not likely associated with orbits that build the boxy/peanut (B/P) bulge. To this purpose, a new set of test particle simulations of a kinematically cold stellar disk evolved in a 3D steady-state barred Milky Way galactic potential, has been analysed in detail. Especially bar particles trapped into the bar are identified through the orbital Jacobi energy $E_{J}$, which allows us to identify the building blocks of the B/P feature and investigate their kinematic properties. Finally, we present preliminary results showing that the high-V$_{\rm los}$ features observed towards the Milky Way bulge are a natural consequence of a large-scale \textit{midplane} particle structure, which is unlikely associated with the Galactic bar.

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.

We report the peculiar chemical abundance patterns of eleven atypical Milky Way (MW) field red giant stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). These atypical giants exhibit strong Al and N enhancements accompanied by C and Mg depletions, strikingly similar to those observed in the so-called second-generation (SG) stars of globular clusters (GCs). Remarkably, we find low-Mg abundances ([Mg/Fe]$<$0.0) together with strong Al and N overabundances in the majority (5/7) of the metal-rich ([Fe/H]$\gtrsim - 1.0$) sample stars, which is at odds with actual observations of SG stars in Galactic CGs of similar metallicities. This chemical pattern is unique and unprecedented among MW stars, posing urgent questions about its origin. These atypical stars could be former SG stars of dissolved GCs formed with intrinsically lower abundances of Mg and enriched Al (subsequently self-polluted by massive AGB stars) or the result of exotic binary systems. We speculate that the stars Mg-deficiency as well as the orbital properties suggest that they could have an extragalactic origin. This discovery should guide future dedicated spectroscopic searches of atypical stellar chemical patterns in our Galaxy; a fundamental step forward to understand the Galactic formation and evolution.

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

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We attempt to constrain the kinematics of the thin and thick disks using the Besancon population synthesis model together with RAVE DR4 and Gaia first data release (TGAS). The RAVE fields were simulated by applying a detailed target selection function and the kinematics was computed using velocity ellipsoids depending on age in order to study the secular evolution. We accounted for the asymmetric drift computed from fitting a Stackel potential to orbits. Model parameters such as velocity dispersions, mean motions, and velocity gradients were adjusted using an ABC-MCMC method. We made use of the metallicity to enhance the separation between thin and thick disks. We show that this model is able to reproduce the kinematics of the local disks in great detail. The disk follows the expected secular evolution, in very good agreement with previous studies of the thin disk. The new asymmetric drift formula, fitted to our previously described St\"ackel potential, fairly well reproduces the velocity distribution in a wide solar neighborhood. The U and W components of the solar motion determined with this method agree well with previous studies. However, we find a smaller V component than previously thought, essentially because we include the variation of the asymmetric drift with distance to the plane. The thick disk is represented by a long period of formation (at least 2 Gyr), during which, as we show, the mean velocity increases with time while the scale height and scale length decrease, very consistently with a collapse phase with conservation of angular momentum. This new Galactic dynamical model is able to reproduce the observed velocities in a wide solar neighborhood at the quality level of the TGAS-RAVE sample, allowing us to constrain the thin and thick disk dynamical evolution, as well as determining the solar motion.

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.

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

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

The cornerstone mission Gaia, together with complementary surveys, will revolutionize our understanding of the formation and history of our Galaxy, providing accurate stellar masses, radii, ages, distances, as well as chemical properties for a very large sample of stars across different Galactic stellar populations.Using an improved population synthesis approach and new stellar evolution models we attempt to evaluate the possibility of deriving ages and masses of clump stars from their chemical properties.A new version of the Besancon Galaxy model (BGM) uses new stellar evolutionary tracks computed with STAREVOL.These provide chemical and seismic properties from the PMS to the early-AGB.For the first time, the BGM can explore the effects of an extra-mixing occurring in red-giant stars.In particular we focus on the effects of thermohaline instability on chemical properties as well as on the determination of stellar ages and masses using the surface [C/N] abundance ratio.The impact of extra-mixing on 3He, 12C/13C, N, and [C/N] abundances along the giant branch is quantified.We underline the crucial contribution of asteroseismology to discriminate between evolutionary states of field giants belonging to the Galactic disc.The inclusion of thermohaline instability has a significant impact on 12C/13C,3He as well as on the [C/N] values.We show the efficiency of thermohaline mixing at different metallicities and its influence on the determined stellar mass and age from the observed [C/N] ratio.We then propose simple relations to determine ages and masses from chemical abundances according to these models.We emphasize the usefulness of population synthesis tools to test stellar models and transport processes inside stars.We show that transport processes occurring in red-giant stars should be taken into account in the determination of ages for future Galactic archaeology studies.(abridged)

Context. Galactic structure studies can be used as a path to constrain the scenario of formation and evolution of our Galaxy. The dependence with the age of stellar population parameters would be linked with the history of star formation and dynamical evolution. Aims. We aim to investigate the structures of the outer Galaxy, such as the scale length, disc truncation, warp and flare of the thin disc and study their dependence with age by using 2MASS data and a population synthesis model (the so-called Besan\c{c}on Galaxy Model). Methods. We have used a genetic algorithm to adjust the parameters on the observed colour-magnitude diagrams at longitudes 80 deg <= l <= 280 deg for |b| <= 5.5 deg. We explored parameter degeneracies and uncertainties. Results. We identify a clear dependence of the thin disc scale length, warp and flare shapes with age. The scale length is found to vary between 3.8 kpc for the youngest to about 2 kpc for the oldest. The warp shows a complex structure, clearly asymmetrical with a node angle changing with age from approximately 165 deg for old stars to 195 deg for young stars. The outer disc is also flaring with a scale height that varies by a factor of two between the solar neighbourhood and a Galactocentric distance of 12 kpc. Conclusions. We conclude that the thin disc scale length is in good agreement with the inside-out formation scenario and that the outer disc is not in dynamical equilibrium. The warp deformation with time may provide some clues to its origin.

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

We report on detection, from observations obtained with the APOGEE spectroscopic survey, of a metal-poor ([Fe/H] $= -1.3$ dex) field giant star with an extreme Mg-Al abundance ratio ([Mg/Fe] $= -0.31$ dex; [Al/Fe] $= 1.49$ dex). Such low Mg/Al ratios are seen only among the second-generation population of globular clusters, and are not present among Galactic disk field stars. The light element abundances of this star, 2M16011638-1201525, suggest that it could have been born in a globular cluster. We explore several origin scenarios, in particular studying the orbit of the star to check the probability of it being kinematically related to known globular clusters. We performed simple orbital integrations assuming the estimated distance of 2M16011638-1201525 and the available six-dimensional phase-space coordinates of 63 globular clusters, looking for close encounters in the past with a minimum distance approach within the tidal radius of each cluster. We found a very low probability that 2M16011638-1201525 was ejected from most globular clusters; however, we note that the best progenitor candidate to host this star is globular cluster $\omega$ Centauri (NGC 5139). Our dynamical investigation demonstrates that 2M16011638-1201525 reaches a distance $|Z_{max}| < 3 $ kpc from the Galactic plane and a minimum and maximum approach to the Galactic center of $R_{min}<0.62$ kpc and $R_{max}<7.26$ kpc in an eccentric ($e\sim0.53$) and retrograde orbit. Since the extreme chemical anomaly of 2M16011638-1201525 has also been observed in halo field stars, this object could also be considered a halo contaminant, likely been ejected into the Milky Way disk from the halo. We conclude that, 2M16011638-20152 is also kinematically consistent with the disk but chemically consistent with halo field stars.

Photometric data from the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC) and the Sloan Digital Sky Survey (SDSS) are used to derive the global structure parameters of the smooth components of the Milky Way. The data, which cover nearly 11,000 deg$^2$ sky area and the full range of Galactic latitude, allow us to construct a globally representative Galactic model. The number density distribution of Galactic halo stars is fitted with an oblate spheroid that decays by power law. The best-fit yields an axis ratio and a power law index $\kappa=0.65$ and $p=2.79$, respectively. The $r$-band differential star counts of three dwarf samples are then fitted with a Galactic model. The best-fit model yielded by a Markov Chain Monte Carlo analysis has thin and thick disk scale heights and lengths of $H_{1}=$ 322\,pc and $L_{1}=$2343\,pc, $H_{2}=$794\,pc and $L_{2}=$3638\,pc, a local thick-to-thin disk density ratio of $f_2=$11\,per\,cent, and a local density ratio of the oblate halo to the thin disk of $f_h=$0.16\,per\,cent. The measured star count distribution, which is in good agreement with the above model for most of the sky area, shows a number of statistically significant large scale overdensities, including some of the previously known substructures, such as the Virgo overdensity and the so-called "north near structure", and a new feature between 150\degr $< l < $ 240\degr~and $-1$5\degr $< b < $ $-$5\degr, at an estimated distance between 1.0 and 1.5\,kpc. The Galactic North-South asymmetry in the anticentre is even stronger than previously thought.

Ground-based optical microlensing surveys have provided tantalising, if inconclusive, evidence for a significant population of free-floating planets (FFPs). Both ground and space-based facilities are being used and developed which will be able to probe the distrubution of FFPs with much better sensitivity. It is vital also to develop a high-precision microlensing simulation framework to evaluate the completeness of such surveys. We present the first signal-to-noise limited calculations of the FFP microlensing rate using the Besancon Galactic model. The microlensing distribution towards the Galactic centre is simulated for wide-area ground-based optical surveys such as OGLE or MOA, a wide-area ground-based near-IR survey, and a targeted space-based near-IR survey which could be undertaken with Euclid or WFIRST. We present a calculation framework for the computation of the optical and near-infrared microlensing rate and optical depth for simulated stellar catalogues which are signal-to-noise limited, and take account of extinction, unresolved stellar background light and finite source size effects, which can be significant for FFPs. We find that the global ground-based I-band yield over a central 200 deg^2 region covering the Galactic centre ranges from 20 Earth-mass FFPs year^-1 up to 3,500 year^-1 for Jupiter FFPs in the limit of 100% detection efficiency, and almost an order of magnitude larger for a K-band survey. For ground-based surveys we find that the inclusion of finite source and the unresolved background reveals a mass-dependent variation in the spatial distribution of FFPs. For a space-based H-band covering 2 deg^2, the yield depends on the target field but maximizes close to the Galactic centre with around 76 Earth through to 1,700 Jupiter FFPs year^-1. For near-IR space-based surveys the spatial distribution of FFPs is found to be largely insensitive to the FFP mass scale.

Formation of globular clusters (GCs), the Galactic bulge, or galaxy bulges in general, are important unsolved problems in Galactic astronomy. Homogeneous infrared observations of large samples of stars belonging to GCs and the Galactic bulge field are one of the best ways to study these problems. We report the discovery by APOGEE of a population of field stars in the inner Galaxy with abundances of N, C, and Al that are typically found in GC stars. The newly discovered stars have high [N/Fe], which is correlated with [Al/Fe] and anti-correlated with [C/Fe]. They are homogeneously distributed across, and kinematically indistinguishable from, other field stars in the same volume. Their metallicity distribution is seemingly unimodal, peaking at [Fe/H]~-1, thus being in disagreement with that of the Galactic GC system. Our results can be understood in terms of different scenarios. N-rich stars could be former members of dissolved GCs, in which case the mass in destroyed GCs exceeds that of the surviving GC system by a factor of ~8. In that scenario, the total mass contained in so-called "first-generation" stars cannot be larger than that in "second-generation" stars by more than a factor of ~9 and was certainly smaller. Conversely, our results may imply the absence of a mandatory genetic link between "second generation" stars and GCs. Last, but not least, N-rich stars could be the oldest stars in the Galaxy, the by-products of chemical enrichment by the first stellar generations formed in the heart of the Galaxy.

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

The most accurate 6D phase-space information from the Radial Velocity Experiment (RAVE) was used to integrate the orbits of 105 stars around the galactic globular cluster 47 Tucanae, to look for close encounters between them in the past, with a minimum distance approach less than the cluster tidal radius. The stars are currently over the distance range 3.0 kpc $<$ d $<$ 5.5 kpc. Using the uncertainties in the current position and velocity vector for both, star and cluster, 105 pairs of star-cluster orbits were generated in a Monte Carlo numerical scheme, integrated over 2 Gyr and considering an axisymmetric and non-axisymmetric Milky-Way-like Galactic potential, respectively. In this scheme, we identified 20 potential cluster members that had close encounters with the globular cluster 47 Tucanae, all of which have a relative velocity distribution (V$_{rel}$) less than 200 km s$^{-1}$ at the minimum distance approach. Among these potential members, 9 had close encounters with the cluster with velocities less than the escape velocity of 47 Tucanae, therefore a scenario of tidal stripping seems likely. These stars have been classified with a 93\% confidence level, leading to the identification of extratidal cluster stars. For the other 11 stars, V$_{rel}$ exceeds the escape velocity of the cluster, therefore they were likely ejected or are unassociated interlopers.

Galactic bulge microlensing surveys provide a probe of Galactic structure. We present the first field-by-field comparison between microlensing observations and the Besan\c{c}on population synthesis Galactic model. Using an updated version of the model we provide maps of optical depth, average event duration and event rate for resolved source populations and for difference imaging (DIA) events. We also compare the predicted event timescale distribution to that observed. The simulation follows the selection criteria of the MOA-II survey (Sumi et al. 2013). We modify the Besan\c{c}on model to include M dwarfs and brown dwarfs. Our best fit model requires a brown dwarf mass function slope of $-0.4$. The model provides good agreement with the observed average duration, and respectable consistency with the shape of the timescale distribution (reduced $\chi^2 \simeq 2.2$). The DIA and resolved source limiting yields bracket the observed number of events by MOA-II ($2.17\times$ and $0.83\times$ the number observed, respectively). We perform a 2-dimensional fit to the event spatial distribution to predict the optical depth and event rate across the Galactic bulge. The most serious difficulty for the model is that it provides only $\sim 50\%$ of the measured optical depth and event rate per star at low Galactic latitude around the inner bulge ($|b|<3{^\circ}$). This discrepancy most likely is associated with known under-estimated extinction and star counts in the innermost regions and therefore provides additional support for a missing inner stellar population.

The Besan\c{c}on Galaxy model was used to compare the infrared colour distribution of synthetic stars with those from 2MASS observations taking the selection function of the data into account, in order to study the shape of the stellar halo of the Milky Way, with complemetary spectroscopic data from SDSS-III/APOGEE survey. Furthermore, we compared the generated mock metallicity distribution of the Besan\c{c}on Galaxy model, to the intrinsic metallicity distribution with reliable stellar parameters from the APOGEE Stellar Parameters and Chemical Abundances Pipeline (ASPCAP). The comparison was carried accross a large volume of the inner part of the Galaxy, revealing that a metal-poor population, [M/H]$<-1.2$ dex, could fill an extended component of the inner galactic halo. With this data set, we are able to model a more realistic mass density distribution of the stellar halo component of the Milky Way, assuming a six-parameters double power-law model, and reconstruct the behaviour of the rotation curve in the inner part of the Galaxy.

We present preliminary results of the CIDA Equatorial Variability Survey (CEVS), looking for quasar (hereafter QSO) candidates near the Galactic plane. The CEVS contains photometric data from extended and adjacent regions of the Milky Way disk ($\sim$ 500 sq. deg.). In this work 2.5 square degrees with moderately high temporal sampling in the CEVS were analyzed. The selection of QSO candidates was based on the study of intrinsic optical photometric variability of 14,719 light curves. We studied samples defined by cuts in the variability index (Vindex $>$ 66.5), periodicity index (Q $>$ 2), and the distribution of these sources in the plane (AT , ${\gamma}$), using a slight modification of the first-order of the structure function for the temporal sampling of the survey. Finally, 288 sources were selected as QSO candidates. The results shown in this work are a first attempt to develop a robust method to detect QSO towards the Galactic plane in the era of massive surveys such as VISTA and Gaia.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE), one of the programs in the Sloan Digital Sky Survey III (SDSS-III), has now completed its systematic, homogeneous spectroscopic survey sampling all major populations of the Milky Way. After a three year observing campaign on the Sloan 2.5-m Telescope, APOGEE has collected a half million high resolution (R~22,500), high S/N (>100), infrared (1.51-1.70 microns) spectra for 146,000 stars, with time series information via repeat visits to most of these stars. This paper describes the motivations for the survey and its overall design---hardware, field placement, target selection, operations---and gives an overview of these aspects as well as the data reduction, analysis and products. An index is also given to the complement of technical papers that describe various critical survey components in detail. Finally, we discuss the achieved survey performance and illustrate the variety of potential uses of the data products by way of a number of science demonstrations, which span from time series analysis of stellar spectral variations and radial velocity variations from stellar companions, to spatial maps of kinematics, metallicity and abundance patterns across the Galaxy and as a function of age, to new views of the interstellar medium, the chemistry of star clusters, and the discovery of rare stellar species. As part of SDSS-III Data Release 12, all of the APOGEE data products are now publicly available.

Using six-dimesional phase-space information from the Fourth Data release of the Radial Velocity Experiment (RAVE) over the range of Galactic longitude 240$^{\circ}< l <$ 360$^{\circ}$ and $V_{LSR} < -239$ kms$^{-1}$, we have computed orbits for 329 RAVE stars that were originally selected as chemically and kinematically related to $\omega$ Centauri. The orbits were integrated in a Milky-Way-like axisymmetric Galactic potential, ignoring the effects of the dynamical evolution of $\omega$ Centauri due to the tidal effects of the Galaxy disk on the cluster along time. We also ignored secular changes in the Milky Way potential over time. In a Monte Carlo scheme, and under the assumption that the stars may have been ejected with velocities greater than the escape velocity ($V_{rel}>V_{esc,0}$) from the cluster, we identified 15 stars as having close encounters with $\omega$ Centauri: (\textit{i}) 8 stars with relative velocities $V_{rel}< 200 $ kms$^{-1}$ may have been ejected $\sim$ 200 Myr ago from $\omega$ Centauri; (\textit{ii}) other group of 7 stars were identified with high relative velocity $V_{rel}> 200 $ kms$^{-1}$ during close encounters, and seems unlikely that they have been ejected from $\omega$ Centauri. We also confirm the link between J131340.4-484714 as potential member of $\omega$ Centauri, and probably ejected $\sim$ 2.0 Myr ago, with a relative velocity $V_{rel}\sim80$ kms$^{-1}$.

We present an estimate of the third integral of motion for axisymmetric three-dimensional potentials. This estimate is based on a Staeckel approximation and is explicitly written as a function of the potential. We tested this scheme for the Besancon Galactic model and two other disc-halo models and find that orbits of disc stars have an accurately conserved third quasi integral. The accuracy ranges from of 0.1% to 1% for heights varying from z = 0~kpc to z= 6 kpc and Galactocentric radii R from 5 to 15kpc. We also tested the usefulness of this quasi integral in analytic distribution functions of disc stellar populations: we show that the distribution function remains approximately stationary and that it allows to recover the potential and forces by applying Jeans equations to its moments.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We analyze a sample of 3,944 low-resolution (R ~ 2000) optical spectra from the Sloan Digital Sky Survey (SDSS), focusing on stars with effective temperatures 5800 < Teff < 6300 K, and distances from the Milky Way plane in excess of 5 kpc, and determine their abundances of Fe, Ca, and Mg. We followed the same methodology as in the previous paper in this series, deriving atmospheric parameters by chi2 minimization, but this time we obtained the abundances of individual elements by fitting their associated spectral lines. Distances were calculated from absolute magnitudes obtained by a statistical comparison of our stellar parameters with stellar-evolution models. The observations reveal a decrease in the abundances of iron, calcium, and magnesium at large distances from the Galactic center. The median abundances for the halo stars analyzed are fairly constant up to a Galactocentric distance r ~ 20 kpc, rapidly decrease between r ~ 20 and r ~ 40 kpc, and flatten out to significantly lower values at larger distances, consistent with previous studies. In addition, we examine the [Ca/Fe] and [Mg/Fe] as a function of Fe/H and Galactocentric distance. Our results show that the most distant parts of the halo show a steeper variation of the [Ca/Fe] and [Mg/Fe] with iron. We found that at the range -1.6 < [Fe/H] < -0.4 [Ca/Fe] decreases with distance, in agreement with earlier results based on local stars. However, the opposite trend is apparent for [Mg/Fe]. Our conclusion that the outer regions of the halo are more metal-poor than the inner regions, based on in situ observations of distant stars, agrees with recent results based on inferences from the kinematics of more local stars, and with predictions of recent galaxy formation simulations for galaxies similar to the Milky Way.

Using a sample of 69,919 red giants from the SDSS-III/APOGEE Data Release 12, we measure the distribution of stars in the [$\alpha$/Fe] vs. [Fe/H] plane and the metallicity distribution functions (MDF) across an unprecedented volume of the Milky Way disk, with radius $3<R<15$ kpc and height $|z|<2$ kpc. Stars in the inner disk ($R<5$ kpc) lie along a single track in [$\alpha$/Fe] vs. [Fe/H], starting with $\alpha$-enhanced, metal-poor stars and ending at [$\alpha$/Fe]$\sim0$ and [Fe/H]$\sim+0.4$. At larger radii we find two distinct sequences in [$\alpha$/Fe] vs. [Fe/H] space, with a roughly solar-$\alpha$ sequence that spans a decade in metallicity and a high-$\alpha$ sequence that merges with the low-$\alpha$ sequence at super-solar [Fe/H]. The location of the high-$\alpha$ sequence is nearly constant across the disk, however there are very few high-$\alpha$ stars at $R>11$ kpc. The peak of the midplane MDF shifts to lower metallicity at larger $R$, reflecting the Galactic metallicity gradient. Most strikingly, the shape of the midplane MDF changes systematically with radius, with a negatively skewed distribution at $3<R<7$ kpc, to a roughly Gaussian distribution at the solar annulus, to a positively skewed shape in the outer Galaxy. For stars with $|z|>1$ kpc or [$\alpha$/Fe]$>0.18$, the MDF shows little dependence on $R$. The positive skewness of the outer disk MDF may be a signature of radial migration; we show that blurring of stellar populations by orbital eccentricities is not enough to explain the reversal of MDF shape but a simple model of radial migration can do so.

The SDSS-III/APOGEE survey operated from 2011-2014 using the APOGEE spectrograph, which collects high-resolution (R~22,500), near-IR (1.51-1.70 microns) spectra with a multiplexing (300 fiber-fed objects) capability. We describe the survey data products that are publicly available, which include catalogs with radial velocity, stellar parameters, and 15 elemental abundances for over 150,000 stars, as well as the more than 500,000 spectra from which these quantities are derived. Calibration relations for the stellar parameters (Teff, log g, [M/H], [alpha/M]) and abundances (C, N, O, Na, Mg, Al, Si, S, K, Ca, Ti, V, Mn, Fe, Ni) are presented and discussed. The internal scatter of the abundances within clusters indicates that abundance precision is generally between 0.05 and 0.09 dex across a broad temperature range; within more limited ranges and at high S/N, it is smaller for some elemental abundances. We assess the accuracy of the abundances using comparison of mean cluster metallicities with literature values, APOGEE observations of the solar spectrum and of Arcturus, comparison of individual star abundances with other measurements, and consideration of the locus of derived parameters and abundances of the entire sample, and find that it is challenging to determine the absolute abundance scale; external accuracy may be good to 0.1-0.2 dex. Uncertainties may be larger at cooler temperatures (Teff<4000K). Access to the public data release and data products is described, and some guidance for using the data products is provided.

We present a survey for RR Lyrae stars in an area of 50 deg$^2$ around the globular cluster $\omega$ Centauri, aimed to detect debris material from the alleged progenitor galaxy of the cluster. We detected 48 RR Lyrae stars of which only 11 have been previously reported. Ten among the eleven previously known stars were found inside the tidal radius of the cluster. The rest were located outside the tidal radius up to distances of $\sim 6$ degrees from the center of the cluster. Several of those stars are located at distances similar to that of $\omega$ Centauri. We investigated the probability that those stars may have been stripped off the cluster by studying their properties (mean periods), calculating the expected halo/thick disk population of RR Lyrae stars in this part of the sky, analyzing the radial velocity of a sub-sample of the RR Lyrae stars, and finally, studying the probable orbits of this sub-sample around the Galaxy. None of these investigations support the scenario that there is significant tidal debris around $\omega$ Centauri, confirming previous studies in the region. It is puzzling that tidal debris have been found elsewhere but not near the cluster itself.

The Exoplanet Euclid Legacy Survey (ExELS) proposes to determine the frequency of cold exoplanets down to Earth mass from host separations of ~1 AU out to the free-floating regime by detecting microlensing events in Galactic Bulge. We show that ExELS can also detect large numbers of hot, transiting exoplanets in the same population. The combined microlensing+transit survey would allow the first self-consistent estimate of the relative frequencies of hot and cold sub-stellar companions, reducing biases in comparing "near-field" radial velocity and transiting exoplanets with "far-field" microlensing exoplanets. The age of the Bulge and its spread in metallicity further allows ExELS to better constrain both the variation of companion frequency with metallicity and statistically explore the strength of star-planet tides. We conservatively estimate that ExELS will detect ~4100 sub-stellar objects, with sensitivity typically reaching down to Neptune-mass planets. Of these, ~600 will be detectable in both Euclid's VIS (optical) channel and NISP H-band imager, with ~90% of detections being hot Jupiters. Likely scenarios predict a range of 2900-7000 for VIS and 400-1600 for H-band. Twice as many can be expected in VIS if the cadence can be increased to match the 20-minute H-band cadence. The separation of planets from brown dwarfs via Doppler boosting or ellipsoidal variability will be possible in a handful of cases. Radial velocity confirmation should be possible in some cases, using 30-metre-class telescopes. We expect secondary eclipses, and reflection and emission from planets to be detectable in up to ~100 systems in both VIS and NISP-H. Transits of ~500 planetary-radius companions will be characterised with two-colour photometry and ~40 with four-colour photometry (VIS,YJH), and the albedo of (and emission from) a large sample of hot Jupiters in the H-band can be explored statistically.

We employ the first two years of data from the near-infrared, high-resolution SDSS-III/APOGEE spectroscopic survey to investigate the distribution of metallicity and alpha-element abundances of stars over a large part of the Milky Way disk. Using a sample of ~10,000 kinematically-unbiased red-clump stars with ~5% distance accuracy as tracers, the [alpha/Fe] vs. [Fe/H] distribution of this sample exhibits a bimodality in [alpha/Fe] at intermediate metallicities, -0.9<[Fe/H]<-0.2, but at higher metallicities ([Fe/H]=+0.2) the two sequences smoothly merge. We investigate the effects of the APOGEE selection function and volume filling fraction and find that these have little qualitative impact on the alpha-element abundance patterns. The described abundance pattern is found throughout the range 5<R<11 kpc and 0<|Z|<2 kpc across the Galaxy. The [alpha/Fe] trend of the high-alpha sequence is surprisingly constant throughout the Galaxy, with little variation from region to region (~10%). Using simple galactic chemical evolution models we derive an average star formation efficiency (SFE) in the high-alpha sequence of ~4.5E-10 1/yr, which is quite close to the nearly-constant value found in molecular-gas-dominated regions of nearby spirals. This result suggests that the early evolution of the Milky Way disk was characterized by stars that shared a similar star formation history and were formed in a well-mixed, turbulent, and molecular-dominated ISM with a gas consumption timescale (1/SFE) of ~2 Gyr. Finally, while the two alpha-element sequences in the inner Galaxy can be explained by a single chemical evolutionary track this cannot hold in the outer Galaxy, requiring instead a mix of two or more populations with distinct enrichment histories.

The Sloan Digital Sky Survey III's Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a high-resolution near-infrared spectroscopic survey covering all of the major components of the Galaxy, including the dust-obscured regions of the inner Milky Way disk and bulge. Here we present a sample of 10,341 likely red-clump stars (RC) from the first two years of APOGEE operations, selected based on their position in color-metallicity-surface-gravity-effective-temperature space using a new method calibrated using stellar-evolution models and high-quality asteroseismology data. The narrowness of the RC locus in color-metallicity-luminosity space allows us to assign distances to the stars with an accuracy of 5 to 10%. The sample extends to typical distances of about 3 kpc from the Sun, with some stars out to 8 kpc, and spans a volume of approximately 100 kpc^3 over 5 kpc <~ R <~ 14 kpc, |Z| <~ 2 kpc, and -15 deg <~ Galactocentric azimuth <~ 30 deg. The APOGEE red-clump (APOGEE-RC) catalog contains photometry from 2MASS, reddening estimates, distances, line-of-sight velocities, stellar parameters and elemental abundances determined from the high-resolution APOGEE spectra, and matches to major proper motion catalogs. We determine the survey selection function for this data set and discuss how the RC selection samples the underlying stellar populations. We use this sample to limit any azimuthal variations in the median metallicity within the ~45 degree-wide azimuthal region covered by the current sample to be <= 0.02 dex, which is more than an order of magnitude smaller than the radial metallicity gradient. This result constrains coherent non-axisymmetric flows within a few kpc from the Sun.

We study the shape of the thick disc using photometric data at high and intermediate latitudes from SDSS and 2MASS surveys. We use the population synthesis approach using an ABC-MCMC method to characterize the thick disc shape, scale height, scale length, local density and flare, and we investigate the extend of the thick disc formation period by simulating several formation episodes. We find that the vertical variation in density is not exponential, but much closer to an hyperbolic secant squared. Assuming a single formation epoch, the thick disc is better fitted with a sech2 scale height of 470 pc and a scale length of 2.3 kpc. However if one simulates two successive formation episodes, mimicking an extended formation period, the older episode has a higher scale height and a larger scale length than the younger episode, indicating a contraction during the collapse phase. The scale height decreases from 800 pc to 340 pc, and the scale length from 3.2 kpc to 2 kpc. The star formation increases from the old episode to the young one. During the fitting process, the halo parameters are also determined. The constraint on the halo shows that a transition between a inner and outer halo, if exists, cannot be at a distance of less than about 30 kpc, which is the limit of our investigation using turnoff halo stars. Finally, we show that extrapolating the thick disc towards the bulge region explains well the stellar populations observed there, that there is no longer need to invoke a classical bulge. To explain these results, the most probable scenario for the thick disc is that it formed while the Galaxy was gravitationally collapsing from well mixed gas-rich giant clumps sustained by high turbulence which awhile prevent a thin disc to form, as proposed by Bournaud et al. (2009). This scenario explains the observations in the thick disc region as well as in the bulge region. (abridged)

We analyze a sample of tens of thousands of spectra of halo turnoff stars, obtained with the optical spectrographs of the Sloan Digital Sky Survey (SDSS), to characterize the stellar halo population "in situ" out to a distance of a few tens of kpc from the Sun. In this paper we describe the derivation of atmospheric parameters. We also derive the overall stellar metallicity distribution based on F-type stars observed as flux calibrators for the Baryonic Oscillations Spectroscopic Survey (BOSS). Our analysis is based on an automated method that determines the set of parameters of a model atmosphere that reproduces each observed spectrum best. We used an optimization algorithm and evaluate model fluxes by means of interpolation in a precomputed grid. In our analysis, we account for the spectrograph's varying resolution as a function of fiber and wavelength. Our results for early SDSS (pre-BOSS upgrade) data compare well with those from the SEGUE Stellar Parameter Pipeline (SSPP), except for stars with logg (cgs units) lower than 2.5. An analysis of stars in the globular cluster M13 reveals a dependence of the inferred metallicity on surface gravity for stars with logg < 2.5, confirming the systematics identified in the comparison with the SSPP. We find that our metallicity estimates are significantly more precise than the SSPP results. We obtain a halo metallicity distribution that is narrower and more asymmetric than in previous studies. The lowest gravity stars in our sample, at tens of kpc from the Sun, indicate a shift of the metallicity distribution to lower abundances, consistent with what is expected from a dual halo system in the Milky Way.

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

Context. White dwarfs can be used to study the structure and evolution of the Galaxy by analysing their luminosity function and initial mass function. Among them, the very cool white dwarfs provide the information for the early ages of each population. Because white dwarfs are intrinsically faint only the nearby (about 20 pc) sample is reasonably complete. The Gaia space mission will drastically increase the sample of known white dwarfs through its 5-6 years survey of the whole sky up to magnitude V = 20-25. Aims. We provide a characterisation of Gaia photometry for white dwarfs to better prepare for the analysis of the scientific output of the mission. Transformations between some of the most common photometric systems and Gaia passbands are derived. We also give estimates of the number of white dwarfs of the different galactic populations that will be observed. Methods. Using synthetic spectral energy distributions and the most recent Gaia transmission curves, we computed colours of three different types of white dwarfs (pure hydrogen, pure helium, and mixed composition with H/He= 0.1). With these colours we derived transformations to other common photometric systems (Johnson-Cousins, Sloan Digital Sky Survey, and 2MASS). We also present numbers of white dwarfs predicted to be observed by Gaia. Results. We provide relationships and colour-colour diagrams among different photometric systems to allow the prediction and/or study of the Gaia white dwarf colours. We also include estimates of the number of sources expected in every galactic population and with a maximum parallax error. Gaia will increase the sample of known white dwarfs tenfold to about 200 000. Gaia will be able to observe thousands of very cool white dwarfs for the first time, which will greatly improve our understanding of these stars and early phases of star formation in our Galaxy.

K2 is the mission concept for a repurposed Kepler mission that uses two reaction wheels to maintain the satellite attitude and provide ~85 days of coverage for ten 105-deg^2 fields along the ecliptic in the first 2.5 years of operation. We examine stellar populations based on the improved Besancon model of the Galaxy, comment on the general properties for the entire ecliptic plane, and provide stellar occurrence rates in the first 6 tentative K2 campaigns grouped by spectral type and luminosity class. For each campaign we distinguish between main sequence stars and giants and provide their density profile as a function of galactic latitude. We introduce the crowding metric that serves for optimized target selection across the campaigns. For all main sequence stars we compute the expected planetary occurrence rates for three planet sizes: 2-4 REarth, 4-8 REarth, and 8-32 REarth with orbital periods up to 50 days. In conjunction with Gaia and the upcoming TESS and Plato missions, K2 will become a gold mine for stellar and planetary astrophysics.

GALEX, the first all sky imaging UV satellite, has imaged a large part of the sky providing an excellent opportunity for studying UV star counts. The aim of our study is to investigate in detail the observed UV star counts obtained by GALEX vis-a-vis the model simulated catalogs produced by the Besancon model of stellar population synthesis in various Galactic directions, and to explore the potential for studying the structure of our Galaxy from images in multiple NUV and FUV filters of the forthcoming Ultraviolet Imaging Telescope (UVIT) to be flown onboard ASTROSAT. We have upgraded the Besancon model of stellar population synthesis to include the UV bands of GALEX and UVIT. Depending on the availability of contiguous GALEX, SDSS, WISE and 2MASS overlapping regions, we have chosen a set of 19 GALEX fields which spread over a range of Galactic directions. We cross-matched GALEX sources with the WISE+2MASS and SDSS catalogs and UV stars in the GALEX catalog are identified by choosing a suitable IR colour, J - W1 (W1 is a WISE band at 3.4 microns). The IR colour cut method, which is used for the first time for separation of stars, is discussed in comparison with the GALEX+SDSS star counts method. We present the results of the UV star counts analysis carried out using the data from GALEX. We find that the Besancon model simulations represent the observed star counts of both the GALEX AIS and MIS well within the error bars in various Galactic directions. Based on the model analysis, we separated out white dwarfs of the disc and blue horizontal branch stars of the halo from the observed sample by selecting a suitable FUV - NUV colour. The Besancon model is now ready for further comparisons in the UV domain and will be used for prospective studies for the UVIT instrument to be flown onboard ASTROSAT.

The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allows us to test hypotheses on the star formation history, star evolution, and chemical and dynamical evolution of the Galaxy. The new version of the Besancon Galaxy model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the initial mass function (IMF) and star formation rate (SFR) of the Galactic disc. We present a new strategy for the generation of thin disc stars, which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. The local dynamical self-consistency is maintained in this new scheme. We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. For the first time, a full sky comparison is performed between BGM and data. This strategy allows us to constrain the IMF slope at high masses, which is found to be close to 3.0 and excludes a shallower slope such as Salpeter's one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 Mo/pc^3 implying that there is no compelling evidence for the significant amount of dark matter in the disc. While the model is fitted to Tycho-2 data, which is a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

We investigate the chemo-kinematic properties of the Milky Way disc by exploring the first year of data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE), and compare our results to smaller optical high-resolution samples in the literature, as well as results from lower resolution surveys such as GCS, SEGUE and RAVE. We start by selecting a high-quality sample in terms of chemistry ($\sim$ 20.000 stars) and, after computing distances and orbital parameters for this sample, we employ a number of useful subsets to formulate constraints on Galactic chemical and chemodynamical evolution processes in the Solar neighbourhood and beyond (e.g., metallicity distributions -- MDFs, [$\alpha$/Fe] vs. [Fe/H] diagrams, and abundance gradients). Our red giant sample spans distances as large as 10 kpc from the Sun. We find remarkable agreement between the recently published local (d $<$ 100 pc) high-resolution high-S/N HARPS sample and our local HQ sample (d $<$ 1 kpc). The local MDF peaks slightly below solar metallicity, and exhibits an extended tail towards [Fe/H] $= -$1, whereas a sharper cut-off is seen at larger metallicities. The APOGEE data also confirm the existence of a gap in the [$\alpha$/Fe] vs. [Fe/H] abundance diagram. When expanding our sample to cover three different Galactocentric distance bins, we find the high-[$\alpha$/Fe] stars to be rare towards the outer zones, as previously suggested in the literature. For the gradients in [Fe/H] and [$\alpha$/Fe], measured over a range of 6 $ < $ R $ <$ 11 kpc in Galactocentric distance, we find a good agreement with the gradients traced by the GCS and RAVE dwarf samples. For stars with 1.5 $<$ z $<$ 3 kpc, we find a positive metallicity gradient and a negative gradient in [$\alpha$/Fe].

Using simulations of box/peanut- (B/P-) shaped bulges, we explore the nature of the X-shape of the Milky Way's bulge. An X-shape can be associated with a B/P-shaped bulge driven by a bar. By comparing in detail the simulations and the observations we show that the principal kinematic imprint of the X-shape is a minimum in the difference between the near and far side mean line-of-sight velocity along the minor axis. This minimum occurs at around |b| = 4{\deg}, which is close to the lower limit at which the X-shape can be detected. No coherent signature of an X-shape can be found in Galactocentric azimuthal velocities, vertical velocities, or any of the dispersions. After scaling our simulations, we find that a best fit to the BRAVA data leads to a bar angle of 15{\deg}. We also explore a purely geometric method for determining the distance to the Galactic Centre by tracing the arms of the X-shape. We find that we are able to determine this ill-known distance to an accuracy of about 5% per cent with sufficiently accurate distance measurements for the red clump stars in the arms.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

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

We present the first microlensing candidate for a free-floating exoplanet-exomoon system, MOA-2011-BLG-262, with a primary lens mass of M_host ~ 4 Jupiter masses hosting a sub-Earth mass moon. The data are well fit by this exomoon model, but an alternate star+planet model fits the data almost as well. Nevertheless, these results indicate the potential of microlensing to detect exomoons, albeit ones that are different from the giant planet moons in our solar system. The argument for an exomoon hinges on the system being relatively close to the Sun. The data constrain the product M pi_rel, where M is the lens system mass and pi_rel is the lens-source relative parallax. If the lens system is nearby (large pi_rel), then M is small (a few Jupiter masses) and the companion is a sub-Earth-mass exomoon. The best-fit solution has a large lens-source relative proper motion, mu_rel = 19.6 +- 1.6 mas/yr, which would rule out a distant lens system unless the source star has an unusually high proper motion. However, data from the OGLE collaboration nearly rule out a high source proper motion, so the exoplanet+exomoon model is the favored interpretation for the best fit model. However, the alternate solution has a lower proper motion, which is compatible with a distant (so stellar) host. A Bayesian analysis does not favor the exoplanet+exomoon interpretation, so Occam's razor favors a lens system in the bulge with host and companion masses of M_host = 0.12 (+0.19 -0.06) M_solar and m_comp = 18 (+28 -100 M_earth, at a projected separation of a_perp ~ 0.84 AU. The existence of this degeneracy is an unlucky accident, so current microlensing experiments are in principle sensitive to exomoons. In some circumstances, it will be possible to definitively establish the low mass of such lens systems through the microlensing parallax effect. Future experiments will be sensitive to less extreme exomoons.

Stellar-substellar binary systems are quite rare, and provide interesting benchmarks. They constrain the complex physics of substellar atmospheres, because several physical parameters of the substellar secondary can be fixed from the much better characterized main sequence primary. We report the discovery of CFBDS J111807-064016, a T2 brown dwarf companion to 2MASS J111806.99-064007.8, a low-mass M4.5-M5 star. The brown-dwarf was identified from the Canada France Brown Dwarf Survey. At a distance of 50-120 pc, the 7.7 arcsec angular separation corresponds to projected separations of 390-900 AU. The primary displays no Halpha emission, placing a lower limit on the age of the system of about 6 Gyr. The kinematics is also consistent with membership in the old thin disc. We obtained near-infrared spectra, which together with recent atmosphere models allow us determine the effective temperature and gravity of both components. From these parameters and the age constraint, evolutionary models estimate masses of 0.10 to 0.15 Msol for the M dwarf, and 0.06 to 0.07 Msol for the T dwarf. This system is a particularly valuable benchmark because the brown dwarf is an early T: the cloud-clearing that occurs at the L/T transition is very sensitive to gravity, metallicity, and detailed dust properties, and produces a large scatter in the colours. This T2 dwarf, with its metallicity measured from the primary and its mass and gravity much better constrained than those of younger early-Ts, will anchor our understanding of the colours of L/T transition brown dwarfs. It is also one of the most massive T dwarfs, just below the hydrogen-burning limit, and all this makes it a prime probe a brown dwarf atmosphere and evolution models.

It has recently been discovered that the Galactic bulge is X-shaped, with the two southern arms of the X both crossing the lines of sight at l=0 and |b|>4, hence producing a double red clump (RC) in the bulge CMD. Dynamical models predict the formation of X-shaped bulges, as extreme cases of boxy-peanut bulges. However, since X-shaped bulges were known to be present only in external galaxies, models have never been compared to 3D kinematical data for individual stars. We study the orbital motion of Galactic bulge stars, in the two arms of the X in the southern hemisphere. The goal is to provide observational constraints to bulge formation models that predict the formation of X-shapes through bar dynamical instabilities. Radial velocities have been obtained for a sample of 454 bulge giants, roughly equally distributed between the bright and faint RC, in a field at (l,b)=(0,-6). Proper motions were derived for all RC stars in the same field by combining images from two epochs obtained 11 years apart. The proper motions for the spectroscopic subsample are analyzed taking into account the radial velocities and metallicities measured from near-IR Calcium triplet lines, for a subsample of 352 stars. The radial velocity distribution of stars in the bright RC, tracing the closer overdensity of bulge stars, shows an excess of stars moving towards the Sun. Similarly, an excess of stars receding from the Sun is seen in the far overdensity, traced by faint RC stars. This is explained with the presence of stars on elongated orbits, most likely streaming along the arms of the X-shaped bulge. Proper motions for these stars are consistent with qualitative predictions of dynamical models of peanut-shaped bulges. Surprisingly, stars on elongated orbits have preferentially metal poor metallicities, while the metal rich ones, in both overdensities, are preferentially found in more axysimmetric orbits.

Despite their large number in the Galaxy, M dwarfs remain elusive objects and the modeling of their photospheres has long remained a challenge (molecular opacities, dust cloud formation). Our objectives are to validate the BT-Settl model atmospheres, update the M dwarf T eff -spectral type relation, and find the atmospheric parameters of the stars in our sample. We compare two samples of optical spectra covering the whole M dwarf sequence with the most recent BT-Settl synthetic spectra and use a \c{hi}2 minimization technique to determine Teff . The first sample consists of 97 low-resolution spectra obtained with NTT at La Silla Observatory. The second sample contains 55 mid-resolution spectra obtained at the Siding Spring Observatory (SSO). The spectral typing is realized by comparison with already classified M dwarfs. We show that the BT-Settl synthetic spectra reproduce the slope of the spectral energy distribution and most of its features. Only the CaOH band at 5570{\AA} and AlH and NaH hydrides in the blue part of the spectra are still missing in the models. The Teff-scale obtained with the higher resolved SSO 2.3 m spectra is consistent with that obtained with the NTT spectra. We compare our Teff scale with those of other authors and to published isochrones using the BT-Settl colors. We also present relations between effective temperature, spectral type and colors of the M dwarfs.

Euclid is the second M-class mission of the ESA Cosmic Vision programme, with the principal science goal of studying dark energy. Euclid is also expected to undertake additional Legacy Science programmes. One proposal is the Exoplanet Euclid Legacy Survey (ExELS) which will be the first survey able to measure the abundance of exoplanets down to Earth mass for host separations from ~1AU out to the free-floating (unbound) regime. The cold and free-floating exoplanet regimes represent a crucial discovery space for testing planet formation theories. ExELS will use the gravitational microlensing technique and will detect 1000 microlensing events per month over 1.6 deg^2 of the Galactic bulge. We assess how many of these events will have detectable planetary signatures using a detailed multi-wavelength microlensing simulator (MABuLS) which incorporates the Besancon Galactic model with 3D extinction. MABuLS is the first theoretical simulation of microlensing to treat the effects of point spread function (PSF) blending self-consistently with the underlying Galactic model. We use MABuLS, together with current numerical models for the Euclid PSFs, to explore a number of designs and de-scope options for ExELS, including the exoplanet yield as a function of filter choice and slewing time, and the effect of systematic photometry errors. Using conservative extrapolations of current empirical exoplanet mass functions determined from ground-based microlensing and radial velocity surveys, ExELS can expect to detect a few hundred cold exoplanets around mainly G, K and M-type stellar hosts, including ~45 Earth-mass planets and ~6 Mars-mass planets for an observing programme totalling 10 months. ExELS will be capable of measuring the cold exoplanet mass function down to Earth mass or below, with orbital separations from ~1AU to the free-floating regime. (Abridged)

(Abridged) We analyze chemical and kinematical properties of about 850 FGK solar neighborhood long-lived dwarfs observed with the HARPS high-resolution spectrograph. The stars in the sample have logg > 4 dex, 5000 < Teff < 6500 K, and -1.39 < [Fe/H] < 0.55 dex. We apply a purely chemical analysis approach based on the [alpha/Fe] vs. [Fe/H] plot to separate Galactic stellar populations into the thin disk, thick disk and high-alpha metal-rich (hamr). Our analysis shows a negative gradient of the rotational velocity of the thin disk stars with [Fe/H] (-17 km s^-1 dex^-1), and a steep positive gradient for both the thick disk and hamr stars with the same magnitude of about +42 km s^-1 dex^-1. For the thin disk stars we observed no correlation between orbital eccentricities and metallicity, but observed a steep negative gradient for the thick disk and hamr stars with practically the same magnitude (about -0.18 dex^-1). Our results suggest that radial migration played an important role in the formation and evolution of the thin disk. For the thick disk stars it is not possible to reach a firm conclusion about their origin. Based on the eccentricity distribution of the thick disk stars only their accretion origin can be ruled out, and the heating and migration scenario could explain the positive steep gradient of V_phi with [Fe/H]. Analyzing the hamr stellar population we found that they share properties of both the thin and thick disk population. A comparison of the properties of the hamr stars with that of the subsample of stars from the N-body/SPH simulation using radial migration suggest that they may have originated from the inner Galaxy. Further detailed investigations would help to clarify their exact nature and origin.

Despite its importance for understanding the nature of early stellar generations and for constraining Galactic bulge formation models, at present little is known about the metal-poor stellar content of the central Milky Way. This is a consequence of the great distances involved and intervening dust obscuration, which challenge optical studies. However, the Apache Point Observatory Galactic Evolution Experiment (APOGEE), a wide-area, multifiber, high-resolution spectroscopic survey within Sloan Digital Sky Survey III (SDSS-III), is exploring the chemistry of all Galactic stellar populations at infrared wavelengths, with particular emphasis on the disk and the bulge. An automated spectral analysis of data on 2,403 giant stars in twelve fields in the bulge obtained during APOGEE commissioning yielded five stars with low metallicity([Fe/H]$\le-1.7$), including two that are very metal-poor [Fe/H]$\sim-2.1$ by bulge standards. Luminosity-based distance estimates place the five stars within the outer bulge, where other 1,246 of the analyzed stars may reside. A manual reanalysis of the spectra verifies the low metallicities, and finds these stars to be enhanced in the $\alpha$-elements O, Mg, and Si without significant $\alpha$-pattern differences with other local halo or metal-weak thick-disk stars of similar metallicity, or even with other more metal-rich bulge stars. While neither the kinematics nor chemistry of these stars can yet definitively determine which, if any, are truly bulge members, rather than denizens of other populations co-located with the bulge, the newly-identified stars reveal that the chemistry of metal-poor stars in the central Galaxy resembles that of metal-weak thick-disk stars at similar metallicity.

Studies of the properties of the inner Galactic Bulge depend strongly on the assumptions about the interstellar extinction. Most of the extinction maps available in the literature lack the information about the distance. We combine the observations with the Besancon model of the Galaxy to investigate the variations of extinction along different lines of sight towards the inner Galactic bulge as a function of distance. In addition we study the variations in the extinction law in the Bulge. We construct color-magnitude diagrams with the following sets of colors: H-Ks and J-Ks from the VVV catalogue as well as Ks-[3.6], Ks-[4.5], Ks-[5.8] and Ks-[8.0] from GLIMPSE-II catalogue matched with 2MASS. Using the newly derived temperature-color relation for M giants that match better the observed color-magnitude diagrams we then use the distance-color relations to derive the extinction as a function of distance. The observed colors are shifted to match the intrinsic colors in the Besan\c{c}on model, as a function of distance, iteratively thereby creating an extinction map with three dimensions: two spatial and one distance dimension along each line of sight towards the bulge. Colour excess maps are presented at a resolution of 15' x 15' for 6 different combinations of colors, in distance bins of 1 kpc. The high resolution and depth of the photometry allows to derive extinction maps to 10 kpc distance and up to 35 magnitudes of extinction in Av (3.5 mag in Aks). Integrated maps show the same dust features and consistent values with the other 2D maps. Starting from the color excess between the observations and the model we investigate the extinction law in near-infrared and its variation along different lines of sight.

The presence of two stellar populations in the Milky Way bulge has been reported recently. We aim at studying the abundances and kinematics of stars in the outer bulge, thereby providing additional constraints on models of its formation. Spectra of 401 red giant stars in a field at (l,b)=(0{\deg},-10{\deg}) were obtained with FLAMES at the VLT. Stars of luminosities down to below the two bulge red clumps (RCs) are included. From these spectra we measure general metallicities, abundances of Fe and the alpha-elements, and radial velocities (RV) of the stars. These measurements as well as photometric data are compared to simulations with the Besancon and TRILEGAL models of the Galaxy. We confirm the presence of two populations among our sample stars: i) a metal-rich one at [M/H] ~+0.3, comprising about 30% of the sample, with low RV dispersion and low alpha-abundance, and ii) a metal-poor population at [M/H] ~-0.6 with high RV dispersion and high alpha-abundance. The metal-rich population could be connected to the Galactic bar. We identify this population as the carrier of the double RC feature. We do not find a significant difference in metallicity or RV between the two RCs, a small difference in metallicity being probably due to a selection effect. The RV dispersion agrees well with predictions of the Besancon Galaxy model, but the metallicity of the "thick bulge" model component should be shifted to lower metallicity by 0.2 to 0.3dex to well reproduce the observations. We present evidence that the metallicity distribution function depends on the evolutionary state of the sample stars, suggesting that enhanced mass loss preferentially removes metal-rich stars. We also confirm the decrease of \alpha-element over-abundance with increasing metallicity.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

Commissioning observations with the Apache Point Observatory Galactic Evolution Experiment (APOGEE), part of the Sloan Digital Sky Survey III, have produced radial velocities (RVs) for ~4700 K/M-giant stars in the Milky Way bulge. These high-resolution (R \sim 22,500), high-S/N (>100 per resolution element), near-infrared (1.51-1.70 um; NIR) spectra provide accurate RVs (epsilon_v~0.2 km/s) for the sample of stars in 18 Galactic bulge fields spanning -1<l<20 deg, |b|<20 deg, and dec>-32 deg. This represents the largest NIR high-resolution spectroscopic sample of giant stars ever assembled in this region of the Galaxy. A cold (sigma_v~30 km/s), high-velocity peak (V_GSR \sim +200 km/s) is found to comprise a significant fraction (~10%) of stars in many of these fields. These high RVs have not been detected in previous MW surveys and are not expected for a simple, circularly rotating disk. Preliminary distance estimates rule out an origin from the background Sagittarius tidal stream or a new stream in the MW disk. Comparison to various Galactic models suggests that these high RVs are best explained by stars in orbits of the Galactic bar potential, although some observational features remain unexplained.

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

We re-analyse photometric near-infrared data in order to investigate why it is so hard to get a consensus for the shape and density law of the bulge, as seen from the literature. To solve the problem we use the Besancon Galaxy Model to provide a scheme for parameter fitting of the structural characteristics of the bulge region. The fitting process allows the determination of the global shape of the bulge main structure. We explore various parameters and shape for the bulge/bar structure based on Ferrer's ellipsoids and fit the shape of the inner disc in the same process. The results show that the main structure is a quite standard triaxial boxy bar/bulge with an orientation of about 13 degree with respect to the Sun-centre direction. But the fit is greatly improved when we add a second structure, which is a longer and thicker ellipsoid. We emphasize that our first ellipsoid represent the main boxy bar of the Galaxy, and that the thick bulge could be either a classical bulge slightly flattened by the effect of the bar potential, or a inner thick disc counterpart. We show that the double clump seen at intermediate latitudes can be reproduced by adding a slight flare to the bar. In order to better characterize the populations, we further simulate several fields which have been surveyed in spectroscopy and for which metallicity distribution function (MDF) are available. The model is in good agreement with these MDF along the minor axis if we assume that the main bar has a mean solar metallicity and the second thicker population has a lower metallicity. It then creates naturally a vertical metallicity gradient by the mixing of the two poulations. (abridged)

We present new radial velocity measurements from the Bulge Radial Velocity Assay (BRAVA), a large scale spectroscopic survey of M-type giants in the Galactic bulge/bar region. The sample of ~4500 new radial velocities, mostly in the region -10 deg < l < +10 deg and b ~ -6 deg more than doubles the existent published data set. Our new data extend our rotation curve and velocity dispersion profile to +20 deg, which is ~2.8 kpc from the Galactic Center. The new data confirm the cylindrical rotation observed at -6 deg and -8 deg, and are an excellent fit to the Shen et al. (2010) N-body bar model. We measure the strength of the TiO molecular band as a first step towards a metallicity ranking of the stellar sample, from which we confirm the presence of a vertical abundance gradient. Our survey finds no strong evidence of previously unknown kinematic streams. We also publish our complete catalog of radial velocities, photometry, TiO band strengths, and spectra, which is available at the IRSA archive: http://irsa.ipac.caltech.edu/ as well as at UCLA: http://brava.astro.ucla.edu/.

The European Southern Observatory's (ESO) Visible and Infrared Survey Telescope for Astronomy (VISTA) is a 4-m class survey telescope for wide-field near-infrared imaging. VISTA is currently running a suite of six public surveys, which will shortly deliver their first Europe wide public data releases to ESO. The VISTA Kilo-degree Infrared Galaxy Survey (VIKING) forms a natural intermediate between current wide shallow, and deeper more concentrated surveys, by targeting two patches totalling 1500 sq.deg in the northern and southern hemispheres with measured 5-sigma limiting depths of Z ~ 22.4, Y ~ 21.4, J ~ 20.9, H ~ 19.9 and Ks ~19.3 (Vega). This architecture forms an ideal working parameter space for the discovery of a significant sample of 6.5 <= z <= 7.5 quasars. In the first data release priority has been placed on small areas encompassing a number of fields well sampled at many wavelengths, thereby optimising science gains and synergy whilst ensuring a timely release of the first products. For rare object searches e.g. high-z quasars, this policy is not ideal since photometric selection strategies generally evolve considerably with the acquisition of data. Without a reasonably representative data set sampling many directions on the sky it is not clear how a rare object search can be conducted in a highly complete and efficient manner. In this paper, we alleviate this problem by supplementing initial data with a realistic model of the spatial, luminosity and colour distributions of sources known to heavily contaminate photometric quasar selection spaces, namely dwarf stars of spectral type M, L and T. We use this model along with a subset of available data to investigate contamination of quasar selection space by cool stars and galaxies and lay down a set of benchmark selection constraints that limit contamination to reasonable levels whilst maintaining high completeness...

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

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

We present new results from BRAVA, a large scale radial velocity survey of the Galactic bulge, using M giant stars selected from the Two Micron All Sky Survey catalogue as targets for the Cerro Tololo Inter-American Observatory 4m Hydra multi-object spectrograph. The purpose of this survey is to construct a new generation of self-consistent bar models that conform to these observations. We report the dynamics for fields at the edge of the Galactic bulge at latitudes b=-8 deg. and compare to the dynamics at b=-4 deg. We find that the rotation curve V(r) is the same at b=-8 deg. as at b=-4 deg. That is, the Galactic boxy bulge rotates cylindrically, as do boxy bulges of other galaxies. The summed line of sight velocity distribution at b=-8 deg. is Gaussian, and the binned longitude-velocity plot shows no evidence for either a (disk) population with cold dynamics or for a (classical bulge) population with hot dynamics. The observed kinematics are well modeled by an edge-on N-body bar, in agreement with published structural evidence. Our kinematic observations indicate that the Galactic bulge is a prototypical product of secular evolution in galaxy disks, in contrast with stellar population results that are most easily understood if major mergers were the dominant formation process.

The bulge is a region of the Galaxy of tremendous interest for understanding galaxy formation. However measuring photometry and kinematics in it raises several inherent issues, such as severe crowding and high extinction in the visible. Using the Besancon Galaxy model and a 3D extinction map, we estimate the stellar density as a function of longitude, latitude and apparent magnitude and we deduce the possibility of reaching and measuring bulge stars with Gaia. We also present an ongoing analysis of the bulge using the Canada-France-Hawaii Telescope.

Context. Thanks to recent large scale surveys in the near infrared such as 2MASS, the galactic plane that most suffers from extinction is revealed and its overall structure can be studied. Aims. This work aims at constraining the structure of the Milky Way external disc as seen in 2MASS data, and in particular the warp. Methods. We use the Two Micron All Sky Survey (hereafter 2MASS) along with the Stellar Population Synthesis Model of the Galaxy, developed in Besancon, to constrain the external disc parameters such as its scale length, its cutoff radius, and the slope of the warp. In order to properly interpret the observations, the simulated stars are reddened using a three dimensional extinction map. The shape of the stellar warp is then compared with previous results and with similar structures in gas and dust. Results. We find new constraints on the stellar disc, which is shown to be asymmetrical, similar to observations of HI. The positive longitude side is found to be easily modelled with a S shape warp but with a slope significantly smaller than the slope seen in the HI warp. At negative longitudes, the disc presents peculiarities which are not well reproduced by any simple model. Finally, comparing with the warp seen in the dust, it seems to follow a slope intermediate between the gas and the stars.

We present detailed maps of the microlensing optical depth and event density over an area of 195 sq. deg towards the Galactic bulge. The maps are computed from synthetic stellar catalogues generated from the Besancon Galaxy Model, which comprises four stellar populations and a three-dimensional extinction map calibrated against the Two-Micron All-Sky Survey. The optical depth maps have a resolution of 15 arcminutes, corresponding to the angular resolution of the extinction map. We compute optical depth and event density maps for all resolved sources above I=19, for unresolved (difference image) sources magnified above this limit, and for bright standard candle sources in the bulge. We show that the resulting optical depth contours are dominated by extinction effects, exhibiting fine structure in stark contrast to previous theoretical optical depth maps. Optical depth comparisons between Galactic models and optical microlensing survey measurements cannot safely ignore extinction or assume it to be smooth. We show how the event distribution for hypothetical J and K-band microlensing surveys, using existing ground-based facilities such as VISTA, UKIRT or CFHT, would be much less affected by extinction, especially in the K band. The near infrared provides a substantial sensitivity increase over current I-band surveys and a more faithful tracer of the underlying stellar distribution, something which upcoming variability surveys such as VVV will be able to exploit. Synthetic population models offer a promising way forward to fully exploit large microlensing datasets for Galactic structure studies.

We present the first results of a wide field survey for cool brown dwarfs with the MegaCam camera on the CFHT telescope, the Canada-France Brown Dwarf Survey, hereafter CFBDS. Our objectives are to find ultracool brown dwarfs and to constrain the field-brown dwarf mass function thanks to a larger sample of L and T dwarfs. We identify candidates in CFHT/MegaCam i' and z' images using optimised psf-fitting within Source Extractor, and follow them up with pointed near-infrared imaging on several telescopes. We have so far analysed over 350 square degrees and found 770 brown dwarf candidates brighter than z'{AB}=22.5. We currently have J-band photometry for 220 of these candidates, which confirms 37% as potential L or T dwarfs. Some are among the reddest and farthest brown dwarfs currently known, including an independent identification of the recently published ULAS J003402.77-005206.7 and the discovery of a second brown dwarf later than T8, CFBDS J005910.83-011401.3. Infrared spectra of three T dwarf candidates confirm their nature, and validate the selection process. The completed survey will discover ~100 T dwarfs and ~500 L dwarfs or M dwarfs later than M8, approximately doubling the number of currently known brown dwarfs. The resulting sample will have a very well-defined selection function, and will therefore produce a very clean luminosity function.

We report the discovery of CFBDS J005910.90-011401.3 (hereafter CFBDS0059), the coolest brown dwarf identified to date. We found CFBDS0059 using i' and z' images from the Canada-France-Hawaii Telescope (CFHT), and present optical and near-infrared photometry, Keck laser guide star adaptive optics imaging, and a complete near-infrared spectrum, from 1.0 to 2.2 $\mu$m. A side to side comparison of the near-infrared spectra of CFBDS0059 and ULAS J003402.77-005206.7 (hereafter ULAS0034), previously the coolest known brown dwarf, indicates that CFBDS0059 is ~50+/-15K cooler. We estimate a temperature of Teff ~ 620K and gravity of log g ~ 4.75. Evolutionary models translate these parameters into an age of 1-5 Gyr and a mass of 15-30 M_Jup. We estimate a photometric distance of ~13pc, which puts CFBDS0059 within easy reach of accurate parallax measurements. Its large proper motion suggests membership in the older population of the thin disk. The spectra of both CFBDS0059 and ULAS J0034 shows probable absorption by a wide ammonia band on the blue side of the $H$-band flux peak. If, as we expect, that feature deepens further for still lower effective temperatures, its appearance will become a natural breakpoint for the transition between the T spectral class and the new Y spectral type. CFBDS0059 and ULAS J0034 would then be the first Y0 dwarfs.

(abridged) By comparing the distribution of dust and gas in the central regions of the Galaxy, we aim to obtain new insights into the properties of the offset dust lanes leading the bar's major axis in the Milky Way. On the one hand, the molecular emission of the dust lanes is extracted from the observed CO l-b-V distribution according to the interpretation of a dynamical model. On the other hand, a three dimensional extinction map of the Galactic central region constructed from near-infrared observations is used as a tracer of the dust itself and clearly reveals dust lanes in its face-on projection. Comparison of the position of both independent detections of the dust lanes is performed in the (l, b) plane. These two completely independent methods are used to provide a coherent picture of the dust lanes in the Milky Way bar. In both the gas and dust distributions, the dust lanes are found to be out of the Galactic plane, appearing at negative latitudes for l > 0 deg and at positive latitudes for l < 0 deg. However, even though there is substantial overlap between the two components, they are offset from one another with the dust appearing to lie closer to the b = 0 deg plane. Two scenarios are proposed to explain the observed offset. The first involves grain destruction by the bar shock and reformation downstream. Due to the decrease in velocity caused by the shock, this occurs at lower z. The second assumes that the gas and dust remain on a common tilted plane, but that the molecular gas decouples from the Milky Way's magnetic field, itself strong enough to resist the shear of the bar's shock. The diffuse gas and dust remain coupled to the field and are carried further downstream. This second scenario has recently been suggested in order to explain observations of the barred galaxy NGC 1097.

The self-binding energy and stability of a galaxy in MOND-based gravity are curiously decreasing functions of its center of mass acceleration towards neighbouring mass concentrations. A tentative indication of this breaking of the Strong Equivalence Principle in field galaxies is the RAVE-observed escape speed in the Milky Way. Another consequence is that satellites of field galaxies will move on nearly Keplerian orbits at large radii (100 - 500 kpc), with a declining speed below the asymptotically constant naive MOND prediction. But consequences of an environment-sensitive gravity are even more severe in clusters, where member galaxies accelerate fast: no more Dark-Halo-like potential is present to support galaxies, meaning that extended axisymmetric disks of gas and stars are likely unstable. These predicted reappearance of asymptotic Keplerian velocity curves and disappearance of "stereotypic galaxies" in clusters are falsifiable with targeted surveys.

With a primary goal of conducting precision weak lensing measurements from space, the COSMOS survey has imaged the largest contiguous area observed by the Hubble Space Telescope (HST) to date using the Advanced Camera for Surveys (ACS). This is the first paper in a series where we describe our strategy for addressing the various technical challenges in the production of weak lensing measurements from the COSMOS data. The COSMOS ACS catalog is constructed from 575 ACS/WFC tiles (1.64 deg^2) and contains a total 1.2x10^6 objects to a limiting magnitude of F814W=26.5. This catalog is made publicly available. The shapes of galaxies have been measured and corrected for the distortion induced by the time varying ACS Point Spread Function and for Charge Transfer Efficiency effects. Next, simulated images are used to derive the shear susceptibility factors that are necessary in order to transform shape measurements into unbiased shear estimators. Finally, for each galaxy, we derive a shape measurement error and utilize this quantity to extract the intrinsic shape noise of the galaxy sample. Interestingly, our results indicate that the intrinsic shape noise varies little with either size, magnitude or redshift. Representing a number density of 66 galaxies per arcmin^2, the final COSMOS weak lensing catalog contains 3.9x10^5 galaxies with accurate shape measurements. The properties of the COSMOS weak lensing catalog described throughout this paper will provide key input numbers for the preparation and design of next-generation wide field space missions.

We report on the stellar content of the COSMOS two degree field, as derived from a rigorous star/galaxy separation approach developed for using stellar sources to define the point spread function variation map used in a study of weak galaxy lensing. The catalog obtained in one filter from the ACS (Advanced Camera for Survey on the Hubble Space Telescope) is cross-identified with ground based multi-wavelength catalogs. The classification is reliable to magnitude $F_{814W}=24$ and the sample is complete even fainter. We construct a color-magnitude diagram and color histograms and compare them with predictions of a standard model of population synthesis. We find features corresponding to the halo subdwarf main sequence turnoff, the thick disk, and the thin disk. This data set provides constraints on the thick disk and spheroid density laws and on the IMF at low mass. We find no evidence of a sharp spheroid edge out to this distance. We identify a blue population of white dwarfs with counts that agree with model predictions. We find a hint for a possible slight stellar overdensity at about 22-34 kpc but the data are not strong enough at present to claim detection of a stream feature in the halo (abridged).

We present spectra of 59 nearby stars candidates, M dwarfs and white dwarfs, previously identified using high proper motion catalogues and the DENIS database. We review the existing spectral classification schemes and spectroscopic parallax calibrations in the near-infrared $J$-band and derive spectral types and distances of the nearby candidates. 42 stars have spectroscopic distances smaller than 25 pc, three of them being white dwarfs. Two targets lie within 10 pc, one M8 star at 10.0 pc (APMPM J0103-3738), and one M4 star at 8.3 pc (LP 225-57). One star, LHS 73, is found to be among the few subdwarfs lying within 20 pc. Furthermore, together with LHS 72, it probably belongs to the closest pair of subdwarfs we know.

The Two Micron All Sky Survey, along with the Stellar Population Synthesis Model of the Galaxy, developed in Besancon, is used to calculate the extinction distribution along different lines of sight. By combining many lines of sight, the large scale distribution of interstellar material can be deduced. The Galaxy model is used to provide the intrinsic colour of stars and their probable distances, so that the near infrared colour excess, and hence the extinction, may be calculated and its distance evaluated. Such a technique is dependent on the model used, however we are able to show that moderate changes in the model parameters result in insignificant changes in the predicted extinction. This technique has now been applied to over 64000 lines of sight, each separated by 15 arcmin, in the inner Galaxy (|l|<=100 deg, |b|<=10 deg). Using our extinction map, we have derived the main characteristics of the large scale structure of the dust distribution: scale height and warp of the ISM disc as well as the angle of the dust in the Galactic Bar. This resulting extinction map will be useful for studies of the inner Galaxy and its stellar populations.

We present a stellar populations analysis of the first release of the CFHTLS (Canada-France-Hawai Telescope Legacy Survey) data. A detailed comparison between the Besan\c{c}on model of the Galaxy and the first data release of the CFHTLS-Deep survey is performed by implementing the MEGACAM photometric system in this model using stellar atmosphere model libraries. The reliability of the theoretical libraries to reproduce the observed colours in the MEGACAM system is investigated. The locations of various stellar species like subdwarfs, white dwarfs, late-type and brown dwarfs, binary systems are identified. The contamination of the stellar sample by quasars and compact galaxies is quantified using spectroscopic data from the VIMOS-VLT Deep Survey (VVDS) as a function of $i'$ magnitude and $r'-i'$ colour. A comparison between simulated counts using the standard IMF at low masses show that the number of very low mass dwarfs may have been underestimated in previous studies. These observations favour a power law IMF following $d(n)/dm \propto m^{-\alpha}$ with $\alpha=2.5$ for $m < $ 0.25 \Msun or $\alpha=3.0$ for $m < $ 0.2 \Msun for single stars. The resulting LF is in agreement with the local LF as measured from the 5 or 25 pc samples. It is in strong disagreement with the Zheng et al (2001) LF measured from deep HST data. We show that this discrepancy can be understood as an indication of a different IMF at low masses at early epochs of the Galaxy compared to the local thin disc IMF.

The bulge is a region of the Galaxy which is of tremendous interest for understanding Galaxy formation. However, measuring photometry and kinematics in it raises several inherent issues, like high extinction in the visible and severe crowding. Here we attempt to estimate the problem of the visibility of the bulge at optical wavelengths, where large CCD mosaics allow to easily cover wide regions from the ground, and where future astrometric missions are planned. Assuming the Besancon Galaxy model and high resolution extinction maps, we estimate the stellar density as a function of longitude, latitude and apparent magnitude and we deduce the possibility of reaching and measuring bulge stars. The method is applied to three Gaia instruments, the BBP and MBP photometers, and the RVS spectrograph. We conclude that, while in the BBP most of the bulge will be accessible, in the MBP there will be a small but significant number of regions where bulge stars will be detected and accurately measured in crowded fields. Assuming that the RVS spectra may be extracted in moderately crowded fields, the bulge will be accessible in most regions apart from the strongly absorbed inner plane regions, because of high extinction, and in low extinction windows like the Baades's window where the crowding is too severe.

We attempt to study the characteristics of the different stellar populations present in the Galactic central region. A Monte Carlo method is used to simultaneously fit 11 thin disc and triaxial outer bulge density parameters on (Ks,J-Ks) star count data in almost 100 windows from the DENIS near infrared large scale survey at -8 deg < l< 12 deg and -4 deg < b < 4 deg. Various bulge density profiles and luminosity functions were tested uing a population synthesis scheme. The best models, selected by a maximum likelihood test, give the following description: the outer bulge is proxy, prolate, and oriented 10.6+/-3 deg with respect to the Sun-center direction. It seems that the main bulge population is not older than 10 Gyr, but this preliminary result needs further work to be confirmed. A significant central hole is found in the middle of the thin disc. We discuss these results in regard to previous findings and the scenario of bulge formation.

Having ruled out the possibility that stellar objects are the main contributor of the dark matter embedding galaxies, microlensing experiments cannot exclude the hypothesis that a significant fraction of the Milky Way dark halo might be made of MACHOs with masses in the range $0.5-0.8 \msun$. Ancient white dwarfs are generally considered the most plausible candidates for such MACHOs. We report the results of a search for such white dwarfs in a proper motion survey covering a 0.16 sqd field at three epochs at high galactic latitude, and 0.938 sqd at two epochs at intermediate galactic latitude (VIRMOS survey), using the CFH telescope. Both surveys are complete to I = 23, with detection efficiency fading to 0 at I = 24.2. Proper motion data are suitable to separate unambiguously halo white dwarfs identified by belonging to a non rotating system. No candidates were found within the colour-magnitude-proper motion volume where such objects can be safely discriminated from any standard population as well as from possible artefacts. In the same volume, we estimate the maximum white dwarf halo fraction compatible with this observation at different significance levels if the halo is at least 14 gigayears old and under different ad hoc initial mass functions. Our data alone rules out a halo fraction greater than 14% at 95% confidence level. Combined with two previous investigations exploring comparable volumes pushes the limit below 4 % (95% confidence level) or below 1.3% (64% confidence), this implies that if baryonic dark matter is present in galaxy halos, it is not, or it is only marginally in the form of faint hydrogen white dwarfs.

We present a systematic search for new neighbourhood stars among the Liverpool-Edinburgh high proper motion survey which we cross-identified with the DENIS survey. Their high proper motions ensure that they are not giant stars. The distances are estimated using DENIS photometry and we found that 100 stars probably lie within 25 parsecs from the Sun. They are mostly M-dwarfs, and 10 are probably white dwarfs. This is the first distance estimate for 84 stars among them. 10 stars, L 170-14A, DENIS J0146291-533931, DENIS J0235219-240038, LP 942-107, DENIS J0320518-635148, LTT 1732, DENIS J0428054-620929, DENIS J2210200-701005, DENIS J2230096-534445, and NLTT 166-3 are estimated to be closer than 15 parsecs. In addition, one star, DENIS J2343155-241047, could also lie within 15 pc if it belongs to the halo.

Since the Hipparcos mission and recent large scale surveys in the optical and the near-infrared, new constraints have been obtained on the structure and evolution history of the Milky Way. The population synthesis approach is a useful tool to interpret such data sets and to test scenarios of evolution of the Galaxy. We present here new constraints on evolution parameters obtained from the Besancon model of population synthesis and analysis of optical and near-infrared star counts. The Galactic potential is computed self-consistently, in agreement with Hipparcos results and the observed rotation curve. Constraints are posed on the outer bulge structure, the warped and flared disc, the thick disc and the spheroid populations. The model is tuned to produce reliable predictions in the visible and the near-infrared in wide photometric bands from U to K. Finally, we describe applications such as photometric and astrometric simulations and a new classification tool based on a Bayesian probability estimator, which could be used in the framework of Virtual Observatories. As examples, samples of simulated star counts at different wavelengths and directions are also given.

The ISOGAL project is an infrared survey of specific regions sampling the Galactic Plane selected to provide information on Galactic structure,stellar populations,stellar mass-loss and the recent star formation history of the inner disk and Bulge of the Galaxy. ISOGAL combines 7 and 15 micron ISOCAM observations - with a resolution of 6'' at worst - with DENIS IJKs data to determine the nature of the sources and theinterstellar extinction. We have observed about 16 square degrees with a sensitivity approaching 10-20mJy, detecting ~10^5 sources,mostly AGB stars,red giants and young stars. The main features of the ISOGAL survey and the observations are summarized in this paper,together with a brief discussion of data processing and quality. The primary ISOGAL products are described briefly (a full description is given in Schuller et al. 2003, astro-ph/0304309): viz. the images and theISOGAL-DENIS five-wavelength point source catalogue. The main scientific results already derived or in progress are summarized. These include astrometrically calibrated 7 and 15um images,determining structures of resolved sources; identification and properties of interstellar dark clouds; quantification of the infrared extinction law and source dereddening; analysis of red giant and (especially) AGB stellar populations in the central Bulge,determining luminosity,presence of circumstellar dust and mass--loss rate,and source classification,supplemented in some cases by ISO/CVF spectroscopy; detection of young stellar objects of diverse types,especially in the inner Bulge with information about the present and recent star formation rate; identification of foreground sources with mid-IR excess. These results are the subject of about 25 refereed papers published or in preparation.

We present new nearby stars extracted from a proper motion catalogue and having a DENIS counterpart. Their distances and spectral type are estimated using the DENIS colours. 107 stars are within 50 pc. 31 stars among them have previously measured distances. In addition, 40 stars may enter within the 50 pc limit depending on which population they belong to. 6 stars among them have already measured distances. 5 objects, LHS5045, LP225-57, LP831-45, LHS1767, and WT792, are probably closer than 15 pc, with LP225-57 at 9.5 pc. Most of these stars are M-type while 4 stars are white dwarfs. 88 M-dwarfs are disc stars, 14 belong to the thick disc and 1 to the spheroid.

We discuss near-infrared star counts at the Galactic pole with a view to guiding the NGST and ground-based NIR cameras. Star counts from deep K-band images from the CFHT are presented, and compared with results from the 2MASS survey and some Galaxy models. With appropriate corrections for detector artifacts and galaxies, the data agree with the models down to K~18, but indicate a larger population of fainter red stars. There is also a significant population of compact galaxies that extend to the observational faint limit of K=20.5. Recent Galaxy models agree well down to K$\sim$19, but diverge at fainter magnitudes.

An automated search for star clusters close to the Galactic plane (|b| < 5 deg) was carried out on the Point Source Catalogue of the DENIS survey. 44% of the Galactic plane have been observed and calibrated. The method allowed to retrieve 22 known star clusters and to identify two new ones, not published yet although previously presented in the 2MASS web site as embedded clusters in HII regions. Extinction in the field and in front of the clusters are estimated using a model of population synthesis. We present the method and give the properties of these clusters.

The interpretation of high proper motion white dwarfs detected by Oppenheimer et al (2001) was the start of a lively controversy. While the discoverers identify a large fraction of their findings as dark halo members, others interpret the same sample as essentially made of disc and/or thick disc stars. We use the comprehensive description of Galactic stellar populations provided by the "Besancon" model to produce a realistic simulation of Oppenheimer et al. data, including all observational selections and calibration biases. The conclusion is unambiguous: Thick disc white dwarfs resulting from ordinary hypotheses on the local density and kinematics are sufficient to explain the observed objects, there is no need for halo white dwarfs. This conclusion is robust to reasonable changes in model ingredients. The main cause of the misinterpretation seems to be that the velocity distribution of a proper motion selected star sample is severely biased in favour of high velocities. This has been neglected in previous analyses. Obviously this does not prove that no such objects like halo white dwarfs can exist, but Oppenheimer et al. observations drive their possible contribution in the dark matter halo down to an extremely low fraction.

Star counts at high and intermediate galactic latitudes, in the visible and the near infrared, are used to determine the density law and the initial mass function of the thick disc population. The combination of shallow fields dominated by stars at the turnoff with deep fields allows the determination of the thick disc mass function in the mass range 0.2-0.8 Msun. Star counts are compared with simulations of a synthesis population model. The fit is based on a maximum likelihood criterion. The best fit model gives a scale height of 800 pc, a scale length of 2500 pc and a local density of 10e-3 stars/pc3 or 7.1e-4 Msun/pc3 for Mv<8. The IMF is found to follow a power law with a slope alpha=-0.5. This is the first determination of the thick disc mass function.

By cross-correlating the results of two recent large-scale surveys, the general properties of a well defined sample of semi-regular variable stars have been determined. ISOGAL mid-infrared photometry and MACHO lightcurves are assembled for approximately 300 stars in the Baade's Windows of low extinction towards the Galactic bulge. These stars are mainly giants of late M spectral type, evolving along the asymptotic giant branch (AGB). They are found to possess a wide and continuous distribution of pulsation periods and to obey an approximate log~period -- bolometric magnitude relation or set of such relations. Approximate mass-loss rates in the range of 1e-8 to 5e-7 M_sun per year are derived from ISOGAL mid-infrared photometry and models of stellar spectra adjusted for the presence of optically-thin circumstellar silicate dust. Mass-loss rates depend on luminosity and pulsation period. Some stars lose mass as rapidly as short-period Miras but do not show Mira-like amplitudes. A period of 70 days or longer is a necessary but not a sufficient condition for mass loss to occur. For AGB stars in the mass-loss ranges that we observe, the functional dependence of mass-loss rate on temperature and luminosity is found to be in agreement with recent theoretical predictions. If we include our mass-loss rates with a sample of extreme mass-losing AGB stars in the Large Magellanic Cloud, we get the general result for AGB stars that mass-loss rate is proportional to luminosity^{2.7}, valid for AGB stars with 10^{-8} to 10^{-4} M_sun per year (Abridged).

We have compiled the near infrared Point Source Catalogue (PSC) towards the Magellanic Clouds (MCs) extracted from the data obtained with the Deep Near Infrared Survey of the Southern Sky - DENIS (Epchtein et al. 1997). The catalogue covers an area of of 19.87* 16 square degrees centered on (RA, DEC)=(5h27m20s, -69o00'00'') for the Large Magellanic Cloud (LMC) and 14.7* 10 square degrees centered on (RA, DEC)=(h02m40s, -73o00'00'' for the Small Magellanic Cloud (SMC) at the epoch J2000. It contains about 1300000 sources towards the LMC and 300000 sources towards the SMC each detected in at least 2 of the 3 photometric bands involved in the survey (I, J, Ks). 70% of the detected sources are true members of the Magellanic Clouds, respectively and consist mainly of red giants, asymptotic giant branch stars and super-giants. The observations have all been made with the same instrument and the data have been calibrated and reduced uniformly. The catalogue provides a homogeneous set of photometric data.

As part of a global analysis of deep star counts to constrain scenarii of galaxy formation and evolution, we investigate possible links between the galactic spheroid and the dark matter halo. A wide set of deep star counts at high and intermediate galactic latitudes is used to determine the large scale density law of the spheroid. Assuming a power density law, the exponent, flattening, local density and IMF slope of this population are estimated. The estimation is checked for robustness against contamination of star counts by the thick disc population. Contamination effects are derived from a model of population synthesis under a broad variety of thick disc parameters. The parameter fit is based on a maximum likelihood criterion. The best fit spheroid density law has a flattening of 0.76, a power index of 2.44. There is a significant degeneracy between these two parameters. The data are also compatible with a slightly less flattened spheroid (c/a = 0.85), in combination with a larger power index (2.75). A flatter spheroid (c/a = 0.6) with a power index of 2 is not excluded either. We also constrain the spheroid IMF slope \alpha to be 1.9 +/- 0.2, leading to a local density of 1.64 10$^{-4}$ stars pc$^{-3}$ and a mass density of 4.15 10$^{-5}$ \Msun pc$^{-3}$. With this slope the expected mass density of brown dwarfs in the halo makes a negligible part of the dark matter halo, as already estimated from microlensing surveys. So, as star count data progresses in depth and extent, the picture of the spheroid star population that comes out points to a shape quite compatible with what we know about the distribution of baryonic dark matter if it is made of stellar remnants, suggesting a common dynamical origin.

A photometric and proper motion survey has been obtained in 2 directions at intermediate latitude: ($l=167.5^\circ$, $b=47.4^\circ$; $\alpha_{2000}=9^h41^m26^s$,$\delta_{2000}=+49^\circ53'27''$) and ($l=278^\circ$, $b=47^\circ$; $\alpha_{2000}=11^h42^m56^s$, $\delta_{2000}=-12^\circ31'42''$). The survey covers 7.13 and 20.84 square degrees, respectively. The limiting magnitude is about 18.5 in V for both directions. We have derived the density laws for stars (M$_{V}$ $\ge$ 3.5) as a function of distance from the galactic plane. The density laws for stars follow a sum of two exponentials with scale heights of 240 pc (thin disk) and 790 pc (thick disk), respectively. The local density of thick disk is found to be 6.1$\pm$3 % relative to the thin disk. The kinematical distribution of stars has been probed to distances up to 3.5 kpc above the galactic plane. New estimates of the parameters of velocity ellipsoid have been derived for the thick disk of the Galaxy. A comparison of our data sets with the Besan\c con model star count predictions has been performed, giving a good agreement in the magnitude range V = 13 to 18.

The ISOGAL project is a survey of the stellar populations, structure, and recent star formation history of the inner disk and bulge of the Galaxy. ISOGAL combines 15 and 7micron ISOCAM observations with DENIS IJKs data to determine the nature of a source and the interstellar extinction. In this paper we report an ISOGAL study of a small field in the inner Galactic Bulge (l=0deg, b=1.0deg, area=0.035 sq. deg) as a prototype of the larger area ISOGAL survey of the inner Galaxy. The five wavelengths of ISOGAL+DENIS, together with the relatively low and constant extinction in front of this specific field, allow reliable determination of the nature of the sources. The primary scientific result of this paper is evidence that the most numerous class of ISOGAL 15micron sources are Red Giants in the Galactic bulge and central disk, with luminosities just above or close to the RGB tip and weak mass-loss rates. They form loose sequences in the magnitude-colour diagrams [15]/Ks-[15] and [15]/[7]-[15]. Their large excesses at 15micron with respect to 2micron and 7micron is due to circumstellar dust produced by mass-loss at low rates. These ISOGAL results are the first systematic evidence and study of dust emission at this early stage (''Intermediate'' AGB), before the onset of the large mass-loss phase. It is thus well established that efficient dust formation is already associated with such low mass-loss rates during this early phase.

We have carried out a sample survey in UBVR photometry and proper motions in different directions in the Galaxy, as part of an investigation of galactic structure and evolution. Three fields in the direction of galactic anticentre, centre and antirotation have been surveyed. Using our new data together with wide-area surveys in other fields available to date, we discuss the radial and vertical structure of the Galaxy. The density laws for stars with M$_{V}\ge$3.5 follow a sum of two exponentials with scale heights of 260$\pm$50 pc (thin disk) and 760$\pm$50 pc, respectively. This second exponential corresponds to a thick disk with a local density of 7.4$^{+2.5}_{-1.5}$ \% relative to the thin disk. The scale lengths for these two populations are respectively 2.3$\pm$0.6 kpc and 3$\pm$1 kpc. The thick disk appears as a kinematically distinct population from the thin disk and shows no vertical gradient. Of the many models that have been proposed for the origin of the thick disk, the evidence at present seems to favour a model in which thick disk formed through the rapid dynamical heating of an early disk by satellite accretion into the disk.

Accurate characterization of thick disc properties from recent kinematic and photometric surveys provides converging evidences that this intermediate population is a sequel of the violent heating of early disc populations by a merging satellite galaxy. The thick disc population is revisited under the light of new data in a number of galactic sample fields. Various thick disc hypotheses are fitted to observational data through a maximum likelihood technique. The resulting characteristics of the thick disc are the following : a scale height of 760 $\pm$ 50 pc, with a local density of 5.6 $\pm$ 1 \% of the thin disc. The scale length is constrained to be 2.8 $\pm$ 0.8 kpc, well in agreement with the disc scale length (2.5 $\pm$ 0.3 kpc). The mean metallicity of the thick disc is found to be -0.7 $\pm$ 0.2 dex, with no {significant} metallicity gradients. These photometric constraints in combination with kinematic data give new constraints on the thick disc formation. We show that thick disc characteristics are {hardly compatible with a top-down formation scenario} but fully compatible with a violent merging event arising at the early thin disc life time as described by Quinn, Hernquist \& Fullagar (1993).

As part of a stellar population sampling program, a series of photometric probes at various field sizes and depths have been obtained in a low extinction window in the galactic anticentre direction. Very deep CCD frames probe the most external parts of the disc, providing strong evidence that the galactic density scale length for the old disc population is rather short (2.5 kpc) and drops abruptly beyond 5.5-6 kpc. Deeper frames in the I band allow to estimate photometric distances and confirm the position of the disc edge. A few stars are found at larger distances. Their number is exactly what we expect if the thick disc does not have any cutoff. We discuss the implications for the formation and evolution of the disc, for the star formation threshold, and for the origin of the thick disc population.

We investigate the kinematics of the intermediate population using the photometric and astrometric sample surveys towards two opposite directions (galactic centre and anticentre) at intermediate latitude. A multivariate discriminant analysis (MDA) is used to distinguish the `Thick Disk' from other populations with the help of the Besan\c con model of population synthesis. The data constrain the asymmetric drift of the intermediate population, which is found to be 40$\pm$10 km sec$^{-1}$ and does not depend on the galactic radius.

We have been doing a photometric and astrometric sample survey programme on Schmidt plates in a few selected directions in the Galaxy with an aim to study galactic stellar populations. In order to calibrate Schmidt plates photometrically, it is necessary to have a number of photometric standards which cover the entire range of magnitudes to be studied on the plates. As Schmidt plates cover a wide area of the sky, the standards should also be spread over the surface of the plate so as to minimise geometrical effects present on the plate. This paper is first in a series of papers aimed at to provide photometric standards in the selected galactic directions to be studied. In this paper we present photometric magnitudes obtained through CCD observation in the galactic direction (l=170deg, b=45deg) centered at alpha(2000) = 9h 41m 20s and delta(2000) = +49deg 54' 20'' and covering an area of 14 sq degs. B and V magnitudes have been obtained for a total of 214 stars out of which for 73 stars U magnitudes also have been obtained. The magnitude range covered is 11<V<21 and 0.18<B-V<1.85.

We have derived a new photographic photometry and proper motions for 20000 stars with completeness to V = 18 in the direction of galactic center at the intermediate latitude (l = 3 deg,b = 47 deg; alpha(1950)=15h 18m, delta(1950) = +02deg 16') for a 15.5 sq deg field. The combination of four glass copies of the Palomar Observatory Sky Survey (i.e. POSS 1402 E&O and POSS 1429 E&O) has been used as a first epoch for proper motion determination. The random error of the proper motions is approximately 0.''3/cen to V=17. The photometric accuracy ranges between 0.07 to 0.10 in the V, B and U bands. We stress the importance of the magnitude and color effects in astrometric surveys of field stars. Using color-magnitude diagrams of a few cluster member stars, a new distance of 6.9+-0.5 kpc is derived for M5 and 20.3+-0.8 kpc for Pal 5 globular clusters. This is in good agreement with other determinations. We have analyzed the components of U+W and V galactic space motions resulting from the accurate proper motions survey. No dependence with z distance is found in the asymmetric drift of the thick disk population. New estimates of the parameters of the velocity ellipsoid have been derived for the thin disk, thick disk and halo populations of the Galaxy.

Models of population synthesis for the Galaxy have been developed in order to understand galactic structure and evolution. They allow to test scenarii of evolution by comparisons between model predictions and observed distributions. Forthcoming near-infrared surveys will enormously increase the amount of information about the distributions of late type stars in the Galaxy. We show here how population synthesis models will help to analyse these surveys and to interprete them in terms of mass function, star formation history and stellar evolution.

As part of a stellar population sampling program, a series of photometric probes at various field sizes and depths have been obtained in a low extinction window in the galactic anticentre direction. Such data set strong constraints on the radial structure of the disc. At the forefront of this "drilling" program, very deep CCD frames probe the most external parts of the disc. Over the whole effective magnitude range (18 to 25), all contributions in the statistics which should be expected from old disc stars beyond 6 kpc vanish, although such stars dominate by far at distances less than 5 kpc. This is the signature of a sharp cut-off in the star density: the edge of the galactic disc between 5.5 and 6 kpc. As a consequence, the galactic radius does not exceed 14 kpc (assuming $R/-/(//sun/)$=8.5). Colours of elliptical galaxies measured in the field rule out the risk of being misled by undetected extinction.

Three colour photometry on CCD frames in the Special Area SpA23 provides a deep probe of the galactic disc in a low absorption window towards the anticenter. Magnitudes to better than 10% at V = 25 and B-V colour down to V = 23 have been obtained. These new data, used in combination with lower magnitude photographic data in a wider field, give a strong evidence that the galactic density scale length is rather short (2.5 kpc) and drops abruptly beyond 6 kpc.