Unveiling the stellar halo with TGAS

The detailed study of the Galactic stellar halo may hold the key to unlocking the assembly history of the Milky Way. Here, we present a machine learning model for selecting metal poor stars from the TGAS catalogue using 5 dimensional phase-space information, coupled with optical and near-IR photometry. We characterise the degree of substructure in our halo sample in the Solar neighbourhood by measuring the velocity correlation function.<br/

Schwarzschild models of the Sculptor dSph galaxy

We have developed a spherically symmetric dynamical model of a dwarf spheroidal galaxy using the Schwarzschild method. This type of modelling yields constraints both on the total mass distribution (e.g. enclosed mass and scale radius) as well as on the orbital structure of the system modelled (e.g. velocity anisotropy). Therefore not only can we derive the dark matter content of these systems, but also explore possible formation scenarios. Here we present preliminary results for the Sculptor dSph. We find that the mass of Sculptor within 1kpc is 8.5\times10^(7\pm0.05) M\odot, its anisotropy profile is tangentially biased and slightly more isotropic near the center. For an NFW profile, the preferred concentration (~15) is compatible with cosmological models. Very cuspy density profiles (steeper than NFW) are strongly disfavoured for Sculptor.Comment: 2 pages, 4 figures, to appear in the proceedings of "Assembling the Puzzle of the Milky Way", Le Grand Bornand (Apr. 17-22, 2011

The properties of the local spiral arms from RAVE data: two-dimensional density wave approach

Using the RAVE survey, we recently brought to light a gradient in the mean galactocentric radial velocity of stars in the extended solar neighbourhood. This gradient likely originates from non-axisymmetric perturbations of the potential, among which a perturbation by spiral arms is a possible explanation. Here, we apply the traditional density wave theory and analytically model the radial component of the two-dimensional velocity field. Provided that the radial velocity gradient is caused by relatively long-lived spiral arms that can affect stars substantially above the plane, this analytic model provides new independent estimates for the parameters of the Milky Way spiral structure. Our analysis favours a two-armed perturbation with the Sun close to the inner ultra-harmonic 4:1 resonance, with a pattern speed \Omega_p=18.6^{+0.3}_{-0.2} km/s/kpc and a small amplitude A=0.55 \pm 0.02% of the background potential (14% of the background density). This model can serve as a basis for numerical simulations in three dimensions, additionally including a possible influence of the galactic bar and/or other non-axisymmetric modes.Comment: 9 pages, 4 figures, accepted for publication in MNRA

A spectroscopic survey of thick disc stars outside the solar neighbourhood

We performed a spectroscopic survey of nearly 700 stars probing the galactic thick disc far from the solar neighbourhood towards the galactic coordinates (l~277, b~47). The derived effective temperatures, surface gravities and overall metallicities were then combined with stellar evolution isochrones, radial velocities and proper motions to derive the distances, kinematics and orbital parameters of the sample stars. The targets belonging to each galactic component (thin disc, thick disc, halo) were selected either on their kinematics or according to their position above the galactic plane, and the vertical gradients were also estimated. We present here atmospheric parameters, distances and kinematics for this sample, and a comparison of our kinematic and metallicity distributions with the Besancon model of the Milky Way. The thick disc far from the solar neighbourhood is found to differ only slightly from the thick disc properties as derived in the solar vicinity. For regions where the thick disc dominates, we measured vertical velocity and metallicity trends of d(V_phi)/dZ = 19 +/- 8 km/s/kpc and d[M/H]/dZ = -0.14 +/- 0.05 dex/kpc, respectively. These trends can be explained as a smooth transition between the different galactic components, although intrinsic gradients could not be excluded. In addition, a correlation d(V_phi)/d[M/H] = -45 +/- 12 km/s/dex between the orbital velocity and the metallicity of the thick disc is detected. This gradient is inconsistent with the SDSS photometric survey analysis, which did not detect any such trend, and challenges radial migration models of thick disc formation. Estimations of the scale heights and scale lengths for different metallicity bins of the thick disc result in consistent values, with hR~3.4 \pm 0.7 kpc, and hZ~694 \pm 45 pc, showing no evidence of relics of destroyed massive satellites.Comment: 19 pages, 15 figures, accepted for publication in A&

Chemical gradients in the Milky Way from the RAVE data

Aims. We aim at measuring the chemical gradients of the elements Mg, Al, Si, and Fe along the Galactic radius to provide new constraints on the chemical evolution models of the Galaxy and Galaxy models such as the Besancon model. Thanks to the large number of stars of our RAVE sample we can study how the gradients vary as function of the distance from the Galactic plane. Methods. We analysed three different samples selected from three independent datasets: a sample of 19 962 dwarf stars selected from the RAVE database, a sample of 10 616 dwarf stars selected from the Geneva-Copenhagen Survey (GCS) dataset, and a mock sample (equivalent to the RAVE sample) created by using the GALAXIA code, which is based on the Besancon model. The three samples were analysed by using the very same method for comparison purposes. We integrated the Galactic orbits and obtained the guiding radii (R-g) and the maximum distances from the Galactic plane reached by the stars along their orbits (Z(max)). We measured the chemical gradients as functions of R-g at different Z(max). Results. We found that the chemical gradients of the RAVE and GCS samples are negative and show consistent trends, although they are not equal: at Z(max) < 0.4 kpc and 4.5 < R-g(kpc) < 9.5, the iron gradient for the RAVE sample is d[Fe/H]/dR(g) = -0.065 dex kpc(-1), whereas for the GCS sample it is d[Fe/H]/dR(g) = -0.043 dex kpc(-1) with internal errors of +/-0.002 and +/-0.004 dex kpc(-1), respectively. The gradients of the RAVE and GCS samples become flatter at larger Z(max). Conversely, the mock sample has a positive iron gradient of d[Fe/H]/dR(g) = +0.053 +/- 0.003 dex kpc(-1) at Z(max) < 0.4 kpc and remains positive at any Z(max). These positive and unrealistic values originate from the lack of correlation between metallicity and tangential velocity in the Besancon model. In addition, the low metallicity and asymmetric drift of the thick disc causes a shift of the stars towards lower R-g and metallicity which, together with the thin-disc stars with a higher metallicity and R-g, generates a fictitious positive gradient of the full sample. The flatter gradient at larger Z(max) found in the RAVE and the GCS samples may therefore be due to the superposition of thin-and thick-disc stars, which mimicks a flatter or positive gradient. This does not exclude the possibility that the thick disc has no chemical gradient. The discrepancies between the observational samples and the mock sample can be reduced by i) decreasing the density; ii) decreasing the vertical velocity; and iii) increasing the metallicity of the thick disc in the Besancon model

Understanding Dwarf Galaxies in order to Understand Dark Matter

Much progress has been made in recent years by the galaxy simulation community in making realistic galaxies, mostly by more accurately capturing the effects of baryons on the structural evolution of dark matter halos at high resolutions. This progress has altered theoretical expectations for galaxy evolution within a Cold Dark Matter (CDM) model, reconciling many earlier discrepancies between theory and observations. Despite this reconciliation, CDM may not be an accurate model for our Universe. Much more work must be done to understand the predictions for galaxy formation within alternative dark matter models.Comment: Refereed contribution to the Proceedings of the Simons Symposium on Illuminating Dark Matter, to be published by Springe

<i>Gaia</i> Data Release 1. Summary of the astrometric, photometric, and survey properties

Context. At about 1000 days after the launch of Gaia we present the first Gaia data release, Gaia DR1, consisting of astrometry and photometry for over 1 billion sources brighter than magnitude 20.7. Aims. A summary of Gaia DR1 is presented along with illustrations of the scientific quality of the data, followed by a discussion of the limitations due to the preliminary nature of this release. Methods. The raw data collected by Gaia during the first 14 months of the mission have been processed by the Gaia Data Processing and Analysis Consortium (DPAC) and turned into an astrometric and photometric catalogue. Results. Gaia DR1 consists of three components: a primary astrometric data set which contains the positions, parallaxes, and mean proper motions for about 2 million of the brightest stars in common with the HIPPARCOS and Tycho-2 catalogues – a realisation of the Tycho-Gaia Astrometric Solution (TGAS) – and a secondary astrometric data set containing the positions for an additional 1.1 billion sources. The second component is the photometric data set, consisting of mean G-band magnitudes for all sources. The G-band light curves and the characteristics of ∼3000 Cepheid and RR-Lyrae stars, observed at high cadence around the south ecliptic pole, form the third component. For the primary astrometric data set the typical uncertainty is about 0.3 mas for the positions and parallaxes, and about 1 mas yr−1 for the proper motions. A systematic component of ∼0.3 mas should be added to the parallax uncertainties. For the subset of ∼94 000 HIPPARCOS stars in the primary data set, the proper motions are much more precise at about 0.06 mas yr−1. For the secondary astrometric data set, the typical uncertainty of the positions is ∼10 mas. The median uncertainties on the mean G-band magnitudes range from the mmag level to ∼0.03 mag over the magnitude range 5 to 20.7. Conclusions. Gaia DR1 is an important milestone ahead of the next Gaia data release, which will feature five-parameter astrometry for all sources. Extensive validation shows that Gaia DR1 represents a major advance in the mapping of the heavens and the availability of basic stellar data that underpin observational astrophysics. Nevertheless, the very preliminary nature of this first Gaia data release does lead to a number of important limitations to the data quality which should be carefully considered before drawing conclusions from the data

Thick disk kinematics from RAVE and the solar motion

Radial velocity surveys such as the Radial Velocity Experiment (RAVE) provide us with measurements of hundreds of thousands of nearby stars most of which belong to the Galactic thin, thick disk or halo. Ideally, to study the Galactic disks (both thin and thick) one should make use of the multi-dimensional phase-space and the whole pattern of chemical abundances of their stellar populations. In this paper, with the aid of the RAVE Survey, we study the thin and thick disks of the Milky Way, focusing on the latter. We present a technique to disentangle the stellar content of the two disks based on the kinematics and other stellar parameters such as the surface gravity of the stars. Using the Padova Galaxy Model, we checked the ability of our method to correctly isolate the thick disk component from the Galaxy mixture of stellar populations. We introduce selection criteria in order to clean the observed radial velocities from the Galactic differential rotation and to take into account the partial sky coverage of RAVE. We developed a numerical technique to statistically disentangle thin and thick disks from their mixture. We deduce the components of the solar motion relative to the Local Standard of Rest (LSR) in the radial and vertical direction, the rotational lag of the thick disk component relative to the LSR, and the square root of the absolute value of the velocity dispersion tensor for the thick disk alone. The analysis of the thin disk is presented in another paper. We find good agreement with previous independent parameter determinations. In our analysis we used photometrically determined distances. In the Appendix we show that similar values can be found for the thick disk alone as derived in the main sections of our paper even without the knowledge of photometric distances.Comment: accepted on A&A, please see companion paper "THIN disk kinem...