Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighboring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational potential to multiple stellar streams simultaneously by maximizing the clustering of the stream stars in action space. We apply this technique to members of the GD-1, Pal 5, Orphan and Helmi streams, exploiting both the individual and combined data sets. We describe the Galactic potential with a St\"ackel model, and vary up to five parameters simultaneously. We find that we can only constrain the enclosed mass, and that the strongest constraints come from the GD-1 and Pal 5 streams whose combined data set yields $M(< 20\ \mathrm{kpc}) = 3.47^{+0.95}_{-1.44} \times 10^{11} \ M_{\odot}$. When including the Orphan and Helmi stream in the data set, the mass uncertainty increases to $M(< 20\ \mathrm{kpc}) = 3.12^{+5.69}_{-1.07} \times 10^{11} \ M_{\odot}$, indicative of the hidden systematic errors in fits to individual streams. These systematics are likely due to insufficiencies in our St\"ackel model of the potential and to the limited phase space explored by the data set of 4 streams, so the larger uncertainty of the combined fit is the more robust measure of the actual uncertainty in the Milky Way potential.

The Gaia mission has provided the largest catalogue ever of sources with tangential velocity information. Using this catalogue for dynamical studies is, however, difficult because most of the stars lack line-of-sight velocity measurements. Recently, we presented a selection of $\sim 10^7$ halo stars with accurate distances that have been selected based on their photometry and proper motions. Using this sample, we model the tail of the velocity distribution in the stellar halo, locally and as a function of distance. Our goal is to measure the escape velocity, and herewith to constrain the mass of our Galaxy. We fit the tail of the velocity distribution with a power-law distribution, a commonly used approach first established by Leonard & Tremaine (1990). For the first time ever we use tangential velocities measured accurately for an unprecedented number of halo stars to estimate the escape velocity. For the solar neighbourhood we determine a very precise lower limit to the escape velocity: $497^{+40}_{-24}~{\rm km/s}$. Our best guess is that this value is underestimated by $\sim 10\%$. For a spherical NFW halo in a Milky Way potential, this translates into a lower limit on the mass of $M_{200} = 0.67^{+0.30}_{-0.15}\cdot10^{12}~{\rm M}_\odot$ and concentration parameter $c = 15.0^{2.6}_{-2.3}$. The behaviour of the escape velocity towards the inner Galaxy follows expectations from Milky Way dynamical models, but curiously we find that its value appears to increase with distance beyond the solar radius. This behaviour is due to a change in the shape of the velocity distribution, and could be related to the presence of velocity clumps. A tentative analysis of the escape velocity as a function of $(R,z)$ shows that it decreases too slowly compared to what is expected for a spherical halo using standard values for the characteristic parameters describing the galactic disc.

Circa 10 billion years ago the Milky Way merged with a massive satellite, Gaia-Enceladus. To gain insight into the properties of its debris we analyse in detail the suite of simulations from Villalobos & Helmi (2008), which includes an experiment that produces a good match to the kinematics of nearby halo stars inferred from Gaia data. We compare the kinematic distributions of stellar particles in the simulations and study the distribution of debris in orbital angular momentum, eccentricity and energy, and its relation to the mass-loss history of the simulated satellite. We confirm that Gaia-Enceladus probably fell in on a retrograde, 30$^\circ$ inclination orbit. We find that while 75% of the debris in our preferred simulation has large eccentricity ($> 0.8$), roughly 9% has eccentricity smaller than 0.6. Star particles lost early have large retrograde motions, and a subset of these have low eccentricity. Such stars would be expected to have lower metallicities as they stem from the outskirts of the satellite, and hence naively they could be confused with debris associated with a separate system. These considerations seem to apply to some of the stars from the postulated Sequoia galaxy. When a massive discy galaxy merges, it leaves behind debris with a complex phase-space structure, a large range of orbital properties, and a range of chemical abundances. Observationally, this results in substructures with very different properties, which can be misinterpreted as implying independent progeny. Detailed chemical abundances of large samples of stars and tailored hydrodynamical simulations are critical to resolving such conundrums.

We present part 2 of the 6th and final Data Release (DR6 or FDR) of the Radial Velocity Experiment (RAVE), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A) and span the complete time frame from the start of RAVE observations on 12 April 2003 to their completion on 4 April 2013. In the second of two publications, we present the data products derived from 518387 observations of 451783 unique stars using a suite of advanced reduction pipelines focussing on stellar atmospheric parameters, in particular purely spectroscopically derived stellar atmospheric parameters (Teff, log(g), and the overall metallicity), enhanced stellar atmospheric parameters inferred via a Bayesian pipeline using Gaia DR2 astrometric priors, and asteroseismically calibrated stellar atmospheric parameters for giant stars based on asteroseismic observations for 699 K2 stars. In addition, we provide abundances of the elements Fe, Al, and Ni, as well as an overall [alpha/Fe] ratio obtained using a new pipeline based on the GAUGUIN optimization method that is able to deal with variable signal-to-noise ratios. The RAVE DR6 catalogs are cross matched with relevant astrometric and photometric catalogs, and are complemented by orbital parameters and effective temperatures based on the infrared flux method. The data can be accessed via the RAVE Web site (http://rave-survey.org) or the Vizier database.

The Radial Velocity Experiment (RAVE) is a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A). The 6th and final data release (DR6 or FDR) is based on 518387 observations of 451783 unique stars. RAVE observations were taken between 12 April 2003 and 4 April 2013. Here we present the genesis, setup and data reduction of RAVE as well as wavelength-calibrated and flux-normalized spectra and error spectra for all observations in RAVE DR6. Furthermore, we present derived spectral classification and radial velocities for the RAVE targets, complemented by cross matches with Gaia DR2 and other relevant catalogs. A comparison between internal error estimates, variances derived from stars with more than one observing epoch and a comparison with radial velocities of Gaia DR2 reveals consistently that 68% of the objects have a velocity accuracy better than 1.4 km/s, while 95% of the objects have radial velocities better than 4.0 km/s. Stellar atmospheric parameters, abundances and distances are presented in subsequent publication. The data can be accessed via the RAVE Web (http://rave-survey.org) or the Vizier database.

We present the analysis of the FLAMES dataset targeting the central 25 arcmin region of the Sextans dSph. This dataset is the third major part of the high resolution spectroscopic section of the ESO large program 171.B-0588(A) obtained by the Dwarf galaxy Abundances and Radial-velocities Team (DART). Our sample is composed of red giant branch stars down to the level of the horizontal branch in Sextans. It allows to address questions related to both stellar nucleosynthesis and galaxy evolution. We provide metallicities for 81 stars, which cover the wide [Fe/H]=$-$3.2 to $-$1.5 dex range. The abundances of 10 other elements are derived: Mg, Ca, Ti, Sc, Cr, Mn, Co, Ni, Ba, and Eu. Despite its small mass, Sextans is a chemically evolved system, with evidence for the contribution of core-collapse and Type Ia supernovae as well as low metallicity AGBs. This new FLAMES sample offers a sufficiently large number of stars with chemical abundances derived at high accuracy to firmly establish the existence of a plateau in [$\alpha$/Fe] at $\sim 0.4$ dex, followed by a decrease above [Fe/H]$\sim-2$ dex. This is in stark similarity with the Fornax and Sculptor dSphs despite their very different masses and star formation histories. This suggests that these three galaxies had very similar star formation efficiencies in their early formation phases, probably driven by the early accretion of smaller galactic fragments, until the UV-background heating impacted them in different ways. The parallel between the Sculptor and Sextans dSph is also striking when considering Ba and Eu. Finally, as to the iron-peak elements, the decline of [Co/Fe] and [Ni/Fe] above [Fe/H]$\sim -2$ implies that the production yields of Ni and Co in SNeIa is lower than that of Fe. The decrease in [Ni/Fe] favours models of SNeIa based on the explosion of double degenerate sub-Chandrasekhar mass white dwarfs.

The Gaia mission has provided the largest ever astrometric chart of the Milky Way. Using it to map the Galactic halo is helpful for disentangling its merger history. The identification of halo stars in Gaia DR2 with reliable distance estimates requires special methods because such stars are typically farther away and scarce. We apply the reduced proper motion (RPM) method to identify halo main sequence stars on the basis of Gaia photometry and proper motions. Using the colour-absolute-magnitude relation for this type of stars, we calculate photometric distances. Our selection results in a set of $\sim10^7$ tentative main sequence halo stars with typical distance uncertainties of $7\%$ and with median velocity errors of 20 km/s. The median distance of our sample is $\sim 4.4$ kpc, with the faintest stars located at $\sim 16$ kpc. The spatial distribution of the stars in our sample is centrally concentrated. Visual inspection of the mean velocities of stars on the sky reveals large-scale patterns as well as clear imprints of the GD-1 stream and tentative hints of the Jhelum and Leiptr streams. Incompleteness and selection effects limit our ability to interpret the patterns reliably as well as to identify new substructures. We define a pseudo-velocity space by setting to zero the line-of-sight velocities of our sample stars. In this space, we recover several known structures such as the footprint of Gaia-Enceladus (i.e. the Gaia-Sausage) as well as the Helmi streams and some other retrograde substructures (Sequoia, Thamnos). We show that the two-point velocity correlation function reveals significant clustering on scales smaller than 100 km/s, of similar amplitude as found for the 6D Gaia halo sample. This clustering indicates the presence of nearby streams that are predominantly phase-mixed.

Recent studies using Gaia DR2 have identified a massive merger in the history of the Milky Way (MW) whose debris is markedly radial and counterrotating. This event, known as the Gaia-Enceladus/Gaia-Sausage (GE/GS), is also hypothesized to have built the majority of the inner stellar halo. We use the cosmological hydrodynamic simulation Illustris to place this merger in the context of galaxy assembly within $\Lambda$CDM. From $\sim$150 MW analogs, $\sim 80 \%$ have experienced at least one merger of similar mass and infall time as GE/GS. Within this sample, 37 have debris as radial as that of the GE/GS, which we dub the Ancient Radial Mergers (ARMs). Counterrotation is not rare among ARMs, with $43 \%$ having $>40 \%$ of their debris in counterrotating orbits. However, the compactness inferred for the GE/GS debris, given its large $\beta$ and its substantial contribution to the stellar halo, is more difficult to reproduce. The median radius of ARM debris is r$_{*,deb}\simeq 45$kpc, while GE/GS is thought to be mostly contained within $r\sim 30$ kpc. For most MW analogs, a few mergers are required to build the inner stellar halo, and ARM debris only accounts for $\sim 12 \%$ of inner accreted stars. Encouragingly, we find one ARM that is both compact and dominates the inner halo of its central, making it our best GE/GS analog. Interestingly, this merger deposits a significant number of stars (M$_*\simeq1.5 \times 10^9 M_\odot$) in the outer halo, suggesting that an undiscovered section of GE/GS may await detection.

The Sagittarius stream is one of the best tools that we currently have to estimate the mass and shape of our Galaxy. However, assigning membership and obtaining the phase-space distribution of the stars that form the tails is quite challenging. Our goal is to produce a catalogue of RR Lyrae stars of Sagittarius and obtain an empiric measurement of the trends along the stream in sky position, distance and tangential velocities. We generate two initial samples from the Gaia DR2 RR Lyrae catalogue: one, selecting only the stars within \pm20deg of the orbital plane of Sagittarius (Strip) and the other, the result of applying the Pole Count Map (nGC3) algorithm. We then use the model-independent, deterministic method developed in this work to remove most of the contamination by detecting and isolating the stream in distance and proper motions. The output is two empiric catalogues: the Strip sample (higher-completeness, lower-purity) which contains 11 677 stars, and the nGC3 sample (higher-purity, lower-completeness) with 6 608 stars. We characterise the changes along the stream in all the available dimensions, the 5 astrometric ones plus the metallicity, covering more than 2pi rad in the sky and obtain new estimates for the apocentres and the mean [Fe/H] of the RR Lyrae population. Also, we show the first map of the two components of the tangential velocity, thanks to the combination of distances and proper motions. Finally, we detect the bifurcation in the leading arm and report no significant difference between the two branches, either in metallicity, kinematics or distance. We provide the largest sample of RR Lyrae candidates of Sagittarius, which can be used as an input for a spectroscopic follow-up or as a reference for the new generation of models of the stream through the interpolators in distance and velocity that we have constructed.

The advent of Gaia's 2nd data release in combination with large spectroscopic surveys are revolutionizing our understanding of the Galaxy. Thanks to these and the knowledge accumulated thus far, a more mature picture of the evolution of the early Milky Way is emerging: * Two of the traditional Galactic components, i.e. the stellar halo and the thick disk, appear to be intimately linked: stars with halo-like kinematics originate in similar proportions, from a "heated" (thick) disk and from debris from a system named Gaia-Enceladus. Gaia-Enceladus was the last big merger event experienced by the Milky Way and probably completed around 10 Gyr ago. The puffed-up stars now present in the halo as a consequence of the merger have thus exposed the existence of a disk component at z ~ 1.8. * The Helmi streams, Sequoia, and Thamnos are amongst the newly uncovered or better characterized merger events. Knowledge of their progenitor's properties, star formation and chemical histories is still incomplete. * Debris' from different objects often overlap in phase-space. A task for the next years will be to use spectroscopic surveys for chemical labelling and to disentangle events from one another using dimensions other than only phase-space, metallicity or [alpha/Fe]. * These surveys will also provide line-of-sight velocities missing for faint stars and more accurate distance determinations for distant objects. The resulting samples of stars will cover a much wider volume of the Galaxy allowing, for example, linking kinematic substructures in the inner halo to spatial overdensities in the outer halo. * All the results obtained so far are in-line with expectations of current cosmological models. Yet, tailored hydrodynamical simulations as well as "constrained" cosmological simulations are needed to push our knowledge of the assembly of the Milky Way back to the earliest times. [abridged]

We aim to measure the proper motion along the Sagittarius stream that is the missing piece to determine its full 6D phase space coordinates. We conduct a blind search of over-densities in proper motion from Gaia DR2 in a broad region around the Sagittarius stream by applying wavelet transform techniques. We find that for most of the sky patches, the highest intensity peaks delineate the path of the Sagittarius stream. The 1500 peaks identified depict a continuous sequence spanning almost $2\pi$ in the sky, only obscured when the stream crosses the Galactic disk. Altogether, around $100\,000$ stars potentially belong to the stream as indicated by a coarse inspection of the colour-magnitude diagrams. From these stars, we determine the proper motion along the Sagittarius stream, making it the proper motion sequence with the largest span and continuity ever measured for a stream. A first comparison with existing N-body models of the stream reveals some discrepancies, especially near the pericentre of the trailing arm and an overestimation of the total proper motion for the leading arm. Our study can be the starting point for determining the variation of the population of stars along the stream, the distance to the stream with red clump stars, and the solar motion. It will also allow a much better measurement of the Milky Way potential.

Aims. We present the first three-dimensional internal motions for individual stars in the Draco dwarf spheroidal galaxy. Methods. By combining first-epoch $Hubble$ $Space$ $Telescope$ observations and second-epoch $Gaia$ Data Release 2 positions, we measured the proper motions of $149$ sources in the direction of Draco. We determined the line-of-sight velocities for a sub-sample of $81$ red giant branch stars using medium resolution spectra acquired with the DEIMOS spectrograph at the Keck II telescope. Altogether, this resulted in a final sample of $45$ Draco members with high-precision and accurate 3D motions, which we present as a table in this paper. Results. Based on this high-quality dataset, we determined the velocity dispersions at a projected distance of $\sim120$ pc from the centre of Draco to be $\sigma_{R} =11.0^{+2.1}_{-1.5}$ km/s, $\sigma_{T}=9.9^{+2.3}_{-3.1}$ km/s and $\sigma_{LOS}=9.0^{+1.1}_{-1.1}$ km/s in the projected radial, tangential, and line-of-sight directions. This results in a velocity anisotropy $\beta=0.25^{+0.47}_{-1.38}$ at $r \gtrsim120$ pc. Tighter constraints may be obtained using the spherical Jeans equations and assuming constant anisotropy and Navarro-Frenk-White (NFW) mass profiles, also based on the assumption that the 3D velocity dispersion should be lower than $\approx 1/3$ of the escape velocity of the system. In this case, we constrain the maximum circular velocity $V_{max}$ of Draco to be in the range of $10.2-17.0$ km/s. The corresponding mass range is in good agreement with previous estimates based on line-of-sight velocities only. Conclusions. Our Jeans modelling supports the case for a cuspy dark matter profile in this galaxy. Firmer conclusions may be drawn by applying more sophisticated models to this dataset and with new datasets from upcoming $Gaia$ releases.

Aims. Several kinematic and chemical substructures have been recently found amongst Milky Way halo stars with retrograde motions. It is currently unclear how these various structures are related to each other. This Letter aims to shed light on this issue. Methods. We explore the retrograde halo with an augmented version of the Gaia DR2 RVS sample, extended with data from three large spectroscopic surveys, namely RAVE, APOGEE and LAMOST. In this dataset, we identify several structures using the HDBSCAN clustering algorithm. We discuss their properties and possible links using all the available chemical and dynamical information. Results. In concordance with previous work, we find that stars with [Fe/H] $<-1$ have more retrograde motions than those with [Fe/H] $>-1$. The retrograde halo contains a mixture of debris from objects like Gaia-Enceladus, Sequoia, and even the chemically defined thick-disc. We find that the Sequoia has a smaller range in orbital energies than previously suggested and is confined to high-energy. Sequoia could be a small galaxy in itself, but since it overlaps both in integrals-of-motion space and chemical abundance space with the less bound debris of Gaia-Enceladus, its nature cannot be fully settled yet. In the low-energy part of the halo we find evidence for at least one more distinct structure: Thamnos. Stars in Thamnos are on low inclination, mildly eccentric retrograde orbits, moving at $v_{\phi}\approx-150$ km/s, and are chemically distinct from the other structures. Conclusions. Even with the excellent Gaia DR2 data it remains challenging to piece together all the fragments found in the retrograde halo. At this point, we are very much in need of large datasets with high-quality high-resolution spectra and tailored high-resolution hydrodynamical simulations of galaxy mergers.

$Context$. The assembly history experienced by the Milky Way is currently being unveiled thanks to the data provided by the $Gaia$ mission. It is likely that the globular cluster system of our Galaxy has followed a similarly intricate formation path. $Aims$. To constrain this formation path, we explore the link between the globular clusters and the known merging events that the Milky Way has experienced. $Methods$. To this end, we combined the kinematic information provided by $Gaia$ for almost all Galactic clusters, with the largest sample of cluster ages available after carefully correcting for systematic errors. To identify clusters with a common origin we analysed their dynamical properties, particularly in the space of integrals of motion. $Results$. We find that about 40% of the clusters likely formed in situ. A similarly large fraction, 35%, appear to be possibly associated to known merger events, in particular to $Gaia$-Enceladus (19%), the Sagittarius dwarf galaxy (5%), the progenitor of the Helmi streams (6%), and to the Sequoia galaxy (5%), although some uncertainty remains due to the degree of overlap in their dynamical characteristics. Of the remaining clusters, 16% are tentatively associated to a group with high binding energy, while the rest are all on loosely bound orbits and likely have a more heterogeneous origin. The resulting age-metallicity relations are remarkably tight and differ in their detailed properties depending on the progenitor, providing further confidence on the associations made. $Conclusions$. We provide a table listing the likely associations. Improved kinematic data by future Gaia data releases and especially a larger, systematic error-free sample of cluster ages would help to further solidify our conclusions.

We identify a simulated Milky Way analog in the EAGLE suite of cosmological hydrodynamical simulations. This galaxy not only shares similar global properties as the Milky Way, but was specifically selected because its merger history resembles that currently known for the Milky Way. In particular we find that this Milky Way analog has experienced its last significant merger (with a stellar mass ratio $\sim 0.2$) at $z\sim 1.2$. We show that this merger affected both the dynamical properties of the stars present at the time, contributing to the formation of a thick disk, and also leading to a significant increase in the star formation rate of the host. This object is thus particularly suitable for understanding the early evolutionary history of the Milky Way. It is also an ideal candidate for re-simulation with much higher resolution as this would allow addressing a plethora of interesting questions such as, for example, the specific distribution of dark matter near the Sun.

The velocity distribution of stars is a sensitive probe of the gravitational potential of the Galaxy, and hence of its dark matter distribution. In particular, the shape of the dark halo (e.g. spherical, oblate, or prolate) determines velocity correlations, and different halo geometries are expected to result in measurable differences. Here we explore and interpret the correlations in the $(v_R, v_z)$-velocity distribution as a function of position in the Milky Way. We selected a high-quality sample of stars from the Gaia DR2 catalogue and characterised the orientation of the velocity distribution or tilt angle over a radial distance range of $[4-13]~$kpc and up to $3.5~$kpc away from the Galactic plane while taking into account the effects of the measurement errors. We find that the tilt angles change from spherical alignment in the inner Galaxy ($R\sim4~$kpc) towards more cylindrical alignments in the outer Galaxy ($R\sim11~$kpc) when using distances that take a global zero-point offset in the parallax of $-29~\mu$as. However, if the amplitude of this offset is underestimated, then the inferred tilt angles in the outer Galaxy only appear shallower and are intrinsically more consistent with spherical alignment for an offset as large as $-54~\mu$as. We further find that the tilt angles do not seem to strongly vary with Galactic azimuth and that different stellar populations depict similar tilt angles. Therefore we introduce a simple analytic function that describes the trends found over the full radial range. Since the systematic parallax errors in Gaia DR2 depend on celestial position, magnitude, and colour in complex ways, it is not possible to fully correct for them. Therefore it will be particularly important for dynamical modelling of the Milky Way to thoroughly characterise the systematics in astrometry in future Gaia data releases.

We test the ability of Schwarzschild's orbit superposition method in measuring the mass content, scale radius and shape of a flattened dwarf spheroidal galaxy. Until now, most dynamical model efforts have assumed that dwarf spheroidal galaxies and their host halos are spherical. We use an Evans model (1993) to construct an isothermal mock galaxy whose properties somewhat resemble those of the Sculptor dwarf spheroidal galaxy. This mock galaxy contains flattened luminous and dark matter components, resulting in a logarithmic profile for the global potential. We have tested how well our Schwarzschild method could constrain the characteristic parameters of the system for different sample sizes, and also if the functional form of the potential was unknown. When assuming the true functional form of the potential, the Schwarzschild modelling technique is able to provide an accurate and precise measurement of the characteristic mass parameter of the system and reproduces well the light distribution and the stellar kinematics of our mock galaxy. When assuming a different functional form for the potential, such as a flattened NFW profile, we also constrain the mass and scale radius to their expected values. However in both cases, we find that the flattening parameter remains largely unconstrained. This is likely because the information content of the velocity dispersion on the geometric shape of the potential is too small, since $\sigma$ is constant across our mock dSph. Our results using Schwarzschild's method indicate that the mass enclosed can be derived reliably, even if the flattening parameter is unknown, and already for samples containing 2000 line-of-sight radial velocities, such as those currently available. Further applications of the method to more general distribution functions of flattened systems are needed to establish how well the flattening of dSph dark halos can be determined.

We present detailed chemical abundances for 99 red-giant branch stars in the centre of the Sculptor dwarf spheroidal galaxy, obtained from high-resolution VLT/FLAMES spectroscopy. This includes abundances of Li (detected in 1 out of 99 stars), Na, $\alpha$-elements (O, Mg, Si, Ca Ti), iron-peak elements (Sc, Cr, Fe, Co, Ni, Zn), r- and s-process elements (Ba, La, Nd, Eu). The sample, covering -2.3< [Fe/H]<-0.9, populates the whole metallicity distribution of the galaxy with the exception of the very low metallicity tail. There is a marked decrease in [$\alpha$/Fe] over our sample, from the Galactic halo plateau value at low [Fe/H] and then, after a `knee', a decrease to sub-solar [$\alpha$/Fe] at high [Fe/H]. This is consistent with products of core-collapse supernovae dominating at early times, followed by the onset of supernova type Ia as early as 12 Gyrs ago. The products from low-mass asymptotic giant branch stars, as traced by the s-process, also participate in the chemical evolution of Sculptor with a timescale comparable to that of supernovae type Ia. The r-process is consistent with having no time delay relative to core-collapse SN at the later stages of the chemical evolution in Sculptor. We derive empirical constraints on the relative importance of supernovae type II and Ia to the nucleosynthesis of individual iron-peak and $\alpha$-elements. The most important contribution of supernovae type Ia is to the iron-peak elements: Fe, Cr and Mn; but there is also a modest but non-negligible contribution to both the heavier $\alpha$-elements: S, Ca and Ti, and some of the iron-peak elements: Sc and Co. We see only very small or no contribution to O, Mg, Ni and Zn from supernovae type Ia in Sculptor. The observed chemical abundances in Sculptor show no evidence of a significantly different initial mass function, compared to that of the Milky Way.

We will study the formation history of the Milky Way, and the earliest phases of its chemical enrichment, with a sample of more than 1.5 million stars at high galactic latitude. Elemental abundances of up to 20 elements with a precision of better than 0.2 dex will be derived for these stars. The sample will include members of kinematically coherent substructures, which we will associate with their possible birthplaces by means of their abundance signatures and kinematics, allowing us to test models of galaxy formation. Our target catalogue is also expected to contain 30,000 stars at a metallicity of less than one hundredth that of the Sun. This sample will therefore be almost a factor of 100 larger than currently existing samples of metal-poor stars for which precise elemental abundances are available (determined from high-resolution spectroscopy), enabling us to study the early chemical evolution of the Milky Way in unprecedented detail.

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

The goal of this survey is to study the formation and evolution of the Milky Way halo to deduce its assembly history and the 3D distribution of mass in the Milky Way. The combination of multi-band photometry, Gaia proper motion and parallax data, and radial velocities and the metallicity and elemental abundances obtained from low-resolution spectra of halo giants with 4MOST, will yield an unprecedented characterisation of the Milky Way halo and its interface with the thick disc. The survey will produce a volume- and magnitude-limited complete sample of giant stars in the halo. It will cover at least 10,000 square degrees of high Galactic latitude, and measure line-of-sight velocities with a precision of 1-2 km/s as well as metallicities to within 0.2 dex.

A new moderately r-process-enhanced metal-poor star, RAVE J093730.5-062655, has been identified in the Milky Way halo as part of an ongoing survey by the R-Process Alliance. The temperature and surface gravity indicate that J0937-0626 is likely a horizontal branch star. At [Fe/H] = -1.86, J0937-0626 is found to have subsolar [X/Fe] ratios for nearly every light, alpha, and Fe-peak element. The low [alpha/Fe] ratios can be explained by an ~0.6 dex excess of Fe; J0937-0626 is therefore similar to the subclass of "iron-enhanced" metal-poor stars. A comparison with Milky Way field stars at [Fe/H] = -2.5 suggests that J0937-0626 was enriched in material from an event, possibly a Type Ia supernova, that created a significant amount of Cr, Mn, Fe, and Ni and smaller amounts of Ca, Sc, Ti, and Zn. The r-process enhancement of J0937-0626 is likely due to a separate event, which suggests that its birth environment was highly enriched in r-process elements. The kinematics of J0937-0626, based on Gaia DR2 data, indicate a retrograde orbit in the Milky Way halo; J0937-0626 was therefore likely accreted from a dwarf galaxy that had significant r-process enrichment.

We estimate the mass of the inner ($<20$ kpc) Milky Way and the axis ratio of its dark matter halo using globular clusters as tracers. At the same time, we constrain the phase-space distribution of the globular cluster system. We use the Gaia DR2 catalogue of 75 globular clusters' proper motions and recent measurements of the proper motions of another 20 distant clusters obtained with the Hubble Space Telescope. We describe the globular cluster system with a 2-component distribution function (DF), with a flat, rotating disc and a rounder, more extended halo. While fixing the Milky Way's disc and bulge, we let the mass and shape of the dark matter halo and we fit these two parameters, together with other six describing the DF, with a Bayesian method. We find the mass of the Galaxy within 20 kpc to be $M(<20{\,\rm kpc})=1.91^{+0.18}_{-0.17} \times 10^{11} M_\odot$, of which $M_{\rm DM}(<20{\,\rm kpc})=1.37^{+0.18}_{-0.17}\times 10^{11}M_\odot$ is in dark matter, and the density axis ratio of the dark matter halo to be $q=1.30 \pm 0.25$. This implies a virial mass $M_{\rm via} = 1.3 \pm 0.3 \times 10^{12} M_\odot$. Our analysis rules out oblate ($q<0.8$) and strongly prolate halos ($q>1.9$) with 99\% probability. Our preferred model reproduces well the observed phase-space distribution of globular clusters and has a disc component that closely resembles that of the Galactic thick disc. The halo component follows a power-law density profile $\rho \propto r^{-3.3}$, has a mean rotational velocity of $V_{\rm rot}\simeq -14\,\rm km\,s^{-1}$ at 20 kpc, and has a mildly radially biased velocity distribution ($\beta\simeq 0.2 \pm 0.07$, fairly constant outside 15 kpc). We also find that our distinction between disc and halo clusters resembles, although not fully, the observed distinction in metal-rich ([Fe/H]$>-0.8$) and metal-poor ([Fe/H]$\leq-0.8$) cluster populations.

The galactic halo likely grew over time in part by assembling smaller galaxies, the so-called building blocks. We investigate if the properties of these building blocks are reflected in the halo white dwarf (WD) population in the Solar neighborhood. Furthermore, we compute the halo WD luminosity functions (WDLFs) for four major building blocks of five cosmologically motivated stellar haloes. We couple the SeBa binary population synthesis model to the Munich-Groningen semi-analytic galaxy formation model, applied to the high-resolution Aquarius dark matter simulations. Although the semi-analytic model assumes an instantaneous recycling approximation, we model the evolution of zero-age main sequence stars to WDs, taking age and metallicity variations of the population into account. Although the majority of halo stars is old and metal-poor and therefore the WDs in the different building blocks have similar properties (including present-day luminosity), we find in our models that the WDs originating from building blocks that have young and/or metal-rich stars can be distinguished from WDs that were born in other building blocks. In practice however, it will be hard to prove that these WDs really originate from different building blocks, as the variations in the halo WD population due to binary WD mergers result in similar effects. The five joined stellar halo WD populations that we modelled result in WDLFs that are very similar to each other. We find that simple models with a Kroupa or Salpeter initial mass function (IMF) fit the observed luminosity function slightly better, since the Chabrier IMF is more top-heavy, although this result is dependent on our choice of the stellar halo mass density in the Solar neighborhood.

The halo of the Milky Way has long been hypothesized to harbour significant amounts of merger debris. This view has been supported over more than a decade by wide-field photometric surveys which have revealed the outer halo to be lumpy. The recent release of Gaia DR2 is allowing us to establish that mergers also have been important and possibly built up the majority of the inner halo. In this work we focus on the Helmi streams, a group of streams crossing the Solar vicinity and known for almost two decades. We characterize their properties and relevance for the build-up of the Milky Way's halo. We identify new members of the Helmi streams in an unprecedented dataset with full phase-space information combining Gaia DR2, and the APOGEE DR2, RAVE DR5 and LAMOST DR4 spectroscopic surveys. Based on the orbital properties of the stars, we find new stream members up to a distance of 5 kpc from the Sun, which we characterize using photometry and metallicity information. We also perform N-body experiments to constrain the time of accretion and properties of the progenitor of the streams. We find nearly 600 new members of the Helmi streams. Their HR diagram reveals a broad age range, from approximately 11 to 13 Gyr, while their metallicity distribution goes from $\sim$ 2.3 to $\sim$1.0, and peaks at [Fe/H] $\sim$1.5. These findings confirm that the streams originate in a dwarf galaxy. Furthermore, we find 7 globular clusters to be likely associated, and which follow a well-defined age-metallicity sequence whose properties suggest a relatively massive progenitor object. Our N-body simulations favour a system with a stellar mass of $\sim 10^8\,\mathrm{M}_\odot$ accreted $5 - 8$ Gyr ago. The debris from the Helmi streams is an important donor to the MilkyWay halo, contributing approximately 15\% of its mass in field stars and 10\% of its globular clusters.

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

This paper presents the detailed abundances and r-process classifications of 126 newly identified metal-poor stars as part of an ongoing collaboration, the R-Process Alliance. The stars were identified as metal-poor candidates from the RAdial Velocity Experiment (RAVE) and were followed-up at high spectral resolution (R~31,500) with the 3.5~m telescope at Apache Point Observatory. The atmospheric parameters were determined spectroscopically from Fe I lines, taking into account <3D> non-LTE corrections and using differential abundances with respect to a set of standards. Of the 126 new stars, 124 have [Fe/H]<-1.5, 105 have [Fe/H]<-2.0, and 4 have [Fe/H]<-3.0. Nine new carbon-enhanced metal-poor stars have been discovered, 3 of which are enhanced in r-process elements. Abundances of neutron-capture elements reveal 60 new r-I stars (with +0.3<=[Eu/Fe]<=+1.0 and [Ba/Eu]<0) and 4 new r-II stars (with [Eu/Fe]>+1.0). Nineteen stars are found to exhibit a `limited-r' signature ([Sr/Ba]>+0.5, [Ba/Eu]<0). For the r-II stars, the second- and third-peak main r-process patterns are consistent with the r-process signature in other metal-poor stars and the Sun. The abundances of the light, alpha, and Fe-peak elements match those of typical Milky Way halo stars, except for one r-I star which has high Na and low Mg, characteristic of globular cluster stars. Parallaxes and proper motions from the second Gaia data release yield UVW space velocities for these stars which are consistent with membership in the Milky Way halo. Intriguingly, all r-II and the majority of r-I stars have retrograde orbits, which may indicate an accretion origin.

We present mean absolute proper motion measurements for seven ultra-faint dwarf galaxies orbiting the Milky Way, namely Bo\"{o}tes III, Carina II, Grus II, Reticulum II, Sagittarius II, Segue 2 and Tucana IV. For four of these dwarfs our proper motion estimate is the first ever provided. The adopted astrometric data come from the second data release of the Gaia mission. We determine the mean proper motion for each galaxy starting from an initial guess of likely members, based either on radial velocity measurements or using stars on the Horizontal Branch identified in the Gaia ($G_{\rm BP}$-$G_{\rm RP}$,$G$) colour-magnitude diagram in the field of view towards the UFD. We then refine their membership iteratively using both astrometry and photometry. We take into account the full covariance matrix among the astrometric parameters when deriving the mean proper motions for these systems. Our procedure provides mean proper motions with typical uncertainties of $\sim0.1$ mas/yr, even for galaxies without prior spectroscopic information. In the case of Segue 2 we find that using radial velocity members only leads to biased results, presumably because of the small number of stars with measured radial velocities. Our procedure allows to maximize the number of member stars per galaxy regardless of the existence of prior spectroscopic information, and can therefore be applied on any faint or distant stellar system within reach of Gaia.

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

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

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.

We explore the phase-space structure of nearby halo stars identified kinematically from Gaia DR2 data. We focus on their distribution in velocity and in "integrals of motion" space as well as on their photometric properties. Our sample of stars selected to be moving at a relative velocity of at least 210 km/s with respect to the Local Standard of Rest, contains an important contribution from the low rotational velocity tail of the disk(s). The $V_R$-distribution of these stars depicts a small asymmetry similar to that seen for the faster rotating thin disk stars near the Sun. We also identify a prominent, slightly retrograde "blob", which traces the metal-poor halo main sequence reported by Gaia Collaboration et al. (2018d). We also find many small clumps especially noticeable in the tails of the velocity distribution of the stars in our sample. Their HR diagrams disclose narrow sequences characteristic of simple stellar populations. This stream-frosting confirms predictions from cosmological simulations, namely that substructure is most apparent amongst the fastest moving stars, typically reflecting more recent accretion events.

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.

There is long tradition extending more than a century on the identification of moving groups in the Solar neighborhood. However, with the advent of large kinematic surveys, and especially of the upcoming Gaia data releases there is a need for more sophisticated and automated substructure finders. We analyze the TGASxRAVE dataset to identify moving groups in the Galactic disk. These groups of stars may then be used to map dynamical and star formation processes in the vicinity of the Sun. We use the ROCKSTAR algorithm, a "friends-of-friends"-like substructure finder in 6D phase-space, and analyze the Hertzsprung-Russell diagrams of the groups identified. We find 125 moving groups within 300 pc of the Sun, containing on average 50 stars, and with 3D velocity dispersions smaller than 10 km/s. Most of these groups were previously unknown. Our photometric analysis allows us to isolate a subsample of 30 statistically significant groups likely composed of stars that were born together.

Stellar halos contain tracers of the assembly history of massive galaxies like our own. Exploiting the synergy between the TGAS and the spectroscopic RAVE surveys, Helmi et al. (2017) recently discovered several distinct substructures in the Solar neighbourhood, defined in integrals of motion space. Some of these substructures may be examples of the building blocks that built up the stellar halo. We analyse the chemical properties of stars in these substructures, with focus on their iron and $\alpha$-element abundances as provided by the RAVE survey chemical pipeline. We perform comparisons of the \feh and \mgfe distributions of the substructures to that of the entire halo sample defined in the TGAS$\times$RAVE dataset.We find that over half of the nine substructures have $\sigma_{\text{[Fe/H]}} \leq 0.3$~dex. Two of the substructures have $\sigma_{\text{[Fe/H]}} \leq 0.1$~dex, which makes them possible remnants of disrupted globular clusters. As expected most substructures and the vast majority of our stellar halo sample are $\alpha$-enhanced. Only one substructure shows a distinct [Mg/Fe] vs [Fe/H] abundance trend distinct from the rest of the halo stars in our sample.

Extended stellar haloes are a natural by-product of the hierarchical formation of massive galaxies. If merging is a non-negligible factor in the growth of our Galaxy, evidence of such events should be encoded in its stellar halo. Reliable identification of genuine halo stars is a challenging task however. The 1st Gaia data release contains the positions, parallaxes and proper motions for over 2 million stars, mostly in the Solar neighbourhood. Gaia DR2 will enlarge this sample to over 1.5 billion stars, the brightest ~5 million of which will have a full phase-space information. Our aim is to develop a machine learning model to reliably identify halo stars, even when their full phase-space information is not available. We use the Gradient Boosted Trees algorithm to build a supervised halo star classifier. The classifier is trained on a sample extracted from the Gaia Universe Model Snapshot, convolved with the errors of TGAS, as well as with the expected uncertainties of the upcoming Gaia DR2. We also trained our classifier on the cross-match between the TGAS and RAVE catalogues, where the halo stars are labelled in an entirely model independent way. We then use this model to identify halo stars in TGAS. When full phase- space information is available and for Gaia DR2-like uncertainties, our classifier is able to recover 90% of the halo stars with at most 30% distance errors, in a completely unseen test set, and with negligible levels of contamination. When line-of-sight velocity is not available, we recover ~60% of such halo stars, with less than 10% contamination. When applied to the TGAS data, our classifier detects 337 high confidence RGB halo stars. Although small, this number is consistent with the expectation from models given the data uncertainties. The large parallax errors are the biggest limitation to identify a larger number of halo stars in all the cases studied.

We combine parallaxes from the first Gaia data release with the spectrophotometric distance estimation framework for stars in the fifth RAVE survey data release. The combined distance estimates are more accurate than either determination in isolation - uncertainties are on average two times smaller than for RAVE-only distances (three times smaller for dwarfs), and 1.4 times smaller than TGAS parallax uncertainties (two times smaller for giants). We are also able to compare the estimates from spectrophotometry to those from Gaia, and use this to assess the reliability of both catalogues and improve our distance estimates. We find that the distances to the lowest log g stars are, on average, overestimated and caution that they may not be reliable. We also find that it is likely that the Gaia random uncertainties are smaller than the reported values. As a byproduct we derive ages for the RAVE stars, many with relative uncertainties less than 20 percent. These results for 219566 RAVE sources have been made publicly available, and we encourage their use for studies that combine the radial velocities provided by RAVE with the proper motions provided by Gaia. A sample that we believe to be reliable can be found by taking only the stars with the flag notification flag_any=0

We study the dynamical properties of halo stars located in the Solar Neighbourhood. Our goal is to explore how the properties of the halo depend on the selection criteria used for defining a sample of halo stars. Once this is understood we proceed to measure the shape and orientation of the halo's velocity ellipsoid and we use this information to put constraints on the gravitational potential of the Galaxy. We use the recently released Gaia DR1 catalogue cross-matched to the RAVE dataset for our analysis. We develop a dynamical criterion based on the distribution function of stars in various Galactic components, using action integrals to identify halo members, and compare this to metallicity and to kinematically selected samples. With this new method, we find 1156 stars in the Solar Neighbourhood to be likely members of the stellar halo. Our dynamically selected sample consists mainly of distant giants on elongated orbits. Their metallicity distribution is rather broad, with roughly half of the stars having [M/H] $\ge -1$ dex. The use of different selection criteria has an important impact on the characteristics of the velocity distributions obtained. Nonetheless, for our dynamically selected and for the metallicity selected samples, we find the local velocity ellipsoid to be aligned in spherical coordinates in a Galactocentric reference frame: this suggests that the total gravitational potential is rather spherical in the region spanned by the orbits of the halo stars in these samples.

It has recently been proposed that if the Galactic dark matter halo were triaxial it would induce lumpiness in the velocity distribution of halo stars in the Solar Neighbourhood through orbital resonances. These substructures could therefore provide a way of measuring its shape. We explore the robustness of this proposal by integrating numerically orbits starting from a realistic set of initial conditions in dark halo potentials of different shape. We have analysed the resulting velocity distributions in Solar neighbourhood-like volumes, and have performed statistical tests for the presence of kinematic substructures. Furthermore, we have characterized the particles' orbits using a Fourier analysis. The local velocity distributions obtained are relatively smooth, statistically consistent with being devoid of substructures even for a dark halo potential with significant but plausible triaxiality. Although resonances are indeed present and associated with specific regions of velocity space, the fraction of stars associated to these is relatively minor. The most significant imprint of the triaxiality of the dark halo is in fact, a variation in the shape of the velocity ellipsoid with spatial location.

We investigate the kinematics of red clump stars in the Solar Neighbourhood by combining data from TGAS and RAVE to constrain the local dark matter density. After calibrating the absolute magnitude of red clump stars, we characterize their velocity distribution over a radial distance range of $6$-$10$ kpc and up to $1.5$ kpc away from the Galactic plane. We then apply the axisymmetric Jeans equations on subsets representing the thin and thick disks to determine the (local) distribution of mass near the disk of our Galaxy. Our kinematic maps are well-behaved permitting a straightforward local determination of the vertical force, which we find to be $K_z^{\rm thin} = -2454 \pm 619$ and $K_z^{\rm thick} = -2141 \pm 774$ $(\mathrm{km/s})^2 \! /\mathrm{kpc}$ at $1.5$ kpc away from the Galactic plane for the thin and thick disk samples and for thin and thick disk scale heights of 0.28 kpc and 1.12 kpc respectively. These measurements can be translated into a local dark matter density $\rho_\textrm{DM} \sim 0.018 \pm 0.002 M_\odot/\textrm{pc}^3$. The systematic error on this estimate is much larger than the quoted statistical error, since even a 10% difference in the scale height of the thin disk leads to a 30% change in the value of $\rho_\textrm{DM}$, and a nearly equally good fit to the data.

We use data from the Radial Velocity Experiment (RAVE) and the Tycho-Gaia astrometric solution catalogue (TGAS) to compute the velocity fields yielded by the radial (VR), azimuthal (Vphi) and vertical (Vz) components of associated Galactocentric velocity. We search in particular for variation in all three velocity components with distance above and below the disc midplane, as well as how each component of Vz (line-of-sight and tangential velocity projections) modifies the obtained vertical structure. To study the dependence of velocity on proper motion and distance we use two main samples: a RAVE sample including proper motions from the Tycho-2, PPMXL and UCAC4 catalogues, and a RAVE-TGAS sample with inferred distances and proper motions from the TGAS and UCAC5 catalogues. In both samples, we identify asymmetries in VR and Vz. Below the plane we find the largest radial gradient to be dVR / dR = -7.01+- 0.61 km\s kpc, in agreement with recent studies. Above the plane we find a similar gradient with dVR / dR= -9.42+- 1.77 km\s kpc. By comparing our results with previous studies, we find that the structure in Vz is strongly dependent on the adopted proper motions. Using the Galaxia Milky Way model, we demonstrate that distance uncertainties can create artificial wave-like patterns. In contrast to previous suggestions of a breathing mode seen in RAVE data, our results support a combination of bending and breathing modes, likely generated by a combination of external or internal and external mechanisms.

The 3D motions of stars in small galaxies beyond our own are minute and yet they are crucial for our understanding of the nature of gravity and dark matter. Even for the dwarf galaxy Sculptor which is one of the best studied systems and inferred to be strongly dark matter dominated, there are conflicting reports on its mean motion around the Milky Way and the 3D internal motions of its stars have never been measured. Here we report, based on data from the Gaia space mission and the Hubble Space Telescope, a new precise measurement of Sculptor's mean proper motion. From this we deduce that Sculptor is currently at its closest approach to the Milky Way and moving on an elongated high-inclination orbit that takes it much farther away than previously thought. For the first time we are also able to measure the internal motions of stars in Sculptor. We find $\sigma_{R}=11.5 \pm 4.3$ km/s and $\sigma_{T}=8.5\pm3.2$ km/s along the projected radial and tangential directions, implying that the stars in our sample move preferentially on radial orbits as quantified by the anisotropy parameter, which we find to be $\beta\sim 0.86^{+0.12}_{-0.83}$ at a location beyond the core radius. Taken at face value such a high radial anisotropy requires abandoning conventional models for the mass distribution in Sculptor. Our sample is dominated by metal-rich stars and for these we find $\beta^{MR} \sim 0.95^{+0.04}_{-0.27}$, a value consistent with multi-component models where Sculptor is embedded in a cuspy dark halo as expected for cold dark matter.

4MOST is a new wide-field, high-multiplex spectroscopic survey facility for the VISTA telescope of ESO. Starting in 2022, 4MOST will deploy more than 2400 fibres in a 4.1 square degree field-of-view using a positioner based on the tilting spine principle. In this ontribution we give an outline of the major science goals we wish to achieve with 4MOST in the area of Galactic Archeology. The 4MOST Galactic Archeology surveys have been designed to address long-standing and far-reaching problems in Galactic science. They are focused on our major themes: 1) Near-field cosmology tests, 2) Chemo-dynamical characterisation of the major Milky Way stellar components, 3) The Galactic Halo and beyond, and 4) Discovery and characterisation of extremely metal-poor stars. In addition to a top-level description of the Galactic surveys we provide information about how the community will be able to join 4MOST via a call for Public Spectroscopic Surveys that ESO will launch.

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

We present a novel analysis of the metal-poor star sample in the complete Radial Velocity Experiment (RAVE) Data Release 5 catalog with the goal of identifying and characterizing all very metal-poor stars observed by the survey. Using a three-stage method, we first identified the candidate stars using only their spectra as input information. We employed an algorithm called t-SNE to construct a low-dimensional projection of the spectrum space and isolate the region containing metal-poor stars. Following this step, we measured the equivalent widths of the near-infrared CaII triplet lines with a method based on flexible Gaussian processes to model the correlated noise present in the spectra. In the last step, we constructed a calibration relation that converts the measured equivalent widths and the color information coming from the 2MASS and WISE surveys into metallicity and temperature estimates. We identified 877 stars with at least a 50% probability of being very metal-poor $(\rm [Fe/H] < -2\,\rm dex)$, out of which 43 are likely extremely metal-poor $(\rm [Fe/H] < -3\,\rm dex )$. The comparison of the derived values to a small subsample of stars with literature metallicity values shows that our method works reliably and correctly estimates the uncertainties, which typically have values $\sigma_{\rm [Fe/H]} \approx 0.2\,\mathrm{dex}$. In addition, when compared to the metallicity results derived using the RAVE DR5 pipeline, it is evident that we achieve better accuracy than the pipeline and therefore more reliably evaluate the very metal-poor subsample. Based on the repeated observations of the same stars, our method gives very consistent results. The method used in this work can also easily be extended to other large-scale data sets, including to the data from the Gaia mission and the upcoming 4MOST survey.

We explore the correlations between velocity and metallicity and the possible distinct chemical signatures of the velocity over-densities of the local Galactic neighbourhood. We use the large spectroscopic survey RAVE and the Geneva Copenhagen Survey. We compare the metallicity distribution of regions in the velocity plane ($v_R,v_\phi$) with that of their symmetric counterparts ($-v_R,v_\phi$). We expect similar metallicity distributions if there are no tracers of a sub-population (e.g., a dispersed cluster, accreted stars), if the disk of the Galaxy is axisymmetric, and if the orbital effects of the spiral arms and the bar are weak. We find that the metallicity-velocity space of the solar neighbourhood is highly patterned. A large fraction of the velocity plane shows differences in the metallicity distribution when comparing symmetric $v_R$ regions. The typical differences in the median metallicity are of $0.05$ dex with a statistical significance of at least $95\%$, and with values up to $0.6$ dex. For low azimuthal velocity $v_\phi$, stars moving outwards in the Galaxy have on average higher metallicity than those moving inwards. These include stars in the Hercules and Hyades moving groups and other velocity branch-like structures. For higher $v_\phi$, the stars moving inwards have higher metallicity than those moving outwards. The most likely interpretation of the metallicity asymmetry is that it is due to the orbital effects of the bar and the radial metallicity gradient of the disk. We present a simulation that supports this idea. We have also discovered a positive gradient in $v_\phi$ with respect to metallicity at high metallicities, apart from the two known positive and negative gradients for the thick and thin disks, respectively.

We characterize the selection function of RAVE using 2MASS as our underlying population, which we assume represents all stars which could have potentially been observed. We evaluate the completeness fraction as a function of position, magnitude, and color in two ways: first, on a field-by-field basis, and second, in equal-size areas on the sky. Then, we consider the effect of the RAVE stellar parameter pipeline on the final resulting catalogue, which in principle limits the parameter space over which our selection function is valid. Our final selection function is the product of the completeness fraction and the selection function of the pipeline. We then test if the application of the selection function introduces biases in the derived parameters. To do this, we compare a parent mock catalogue generated using Galaxia with a mock-RAVE catalogue where the selection function of RAVE has been applied. We conclude that for stars brighter than I = 12, between $4000 \rm K < T_{\rm eff} < 8000 \rm K$ and $0.5 < \rm{log}\,g < 5.0$, RAVE is kinematically and chemically unbiased with respect to expectations from Galaxia.

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.

Aims: We present the first measurement of the proper motion and orbit of the very distant and intriguing globular cluster NCG 2419. Methods: We have combined data from HST and Gaia DR1 to derive the relative proper motions of stars in the direction to the cluster. To tie to an absolute reference frame we have used a background galaxy located in the field. Results: We find the absolute proper motion of NGC 2419 to be $(\mu_{\alpha}\cos(\delta)$, $\mu_{\delta}$)=($-0.17\pm0.26,-0.49\pm0.17$) mas/yr. We have integrated the orbit of the cluster in a Galactic potential and found it to oscillate between $\sim$53 kpc and $\sim$98 kpc on a nearly polar orbit. This makes it very likely that NGC 2419 is a former cluster of the Sagittarius dwarf spheroidal galaxy, also because it shares the same sense of rotation around the Milky Way.

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.

The hierarchical structure formation model predicts that stellar halos should form, at least partly, via mergers. If this was a predominant formation channel for the Milky Way's halo, imprints of this merger history in the form of moving groups or streams should exist also in the vicinity of the Sun. Here we study the kinematics of halo stars in the Solar neighbourhood using the very recent first data release from the Gaia mission, and in particular the TGAS dataset, in combination with data from the RAVE survey. Our aim is to determine the amount of substructure present in the phase-space distribution of halo stars that could be linked to merger debris. To characterise kinematic substructure, we measure the velocity correlation function in our sample of halo (low metallicity) stars. We also study the distribution of these stars in the space of energy and two components of the angular momentum, in what we call "Integrals of Motion" space. The velocity correlation function reveals substructure in the form of an excess of pairs of stars with similar velocities, well above that expected for a smooth distribution. Comparison to cosmological simulations of the formation of stellar halos indicate that the levels found are consistent with the Galactic halo having been built fully via accretion. Similarly, the distribution of stars in the space of "Integrals of motion" is highly complex. A strikingly high fraction (between 58% and upto 73%) of the stars that are somewhat less bound than the Sun are on (highly) retrograde orbits. A simple comparison to Milky Way-mass galaxies in cosmological hydrodynamical simulations suggests that less than 1% have such prominently retrograde outer halos. We also identify several other statistically significant structures in "Integrals of Motion" space that could potentially be related to merger events.

We present a set of 87 RAVE stars with detected solar like oscillations, observed during Campaign 1 of the K2 mission (RAVE K2-C1 sample). This dataset provides a useful benchmark for testing the gravities provided in RAVE Data Release 4 (DR4), and is key for the calibration of the RAVE Data Release 5 (DR5). In the present work, we use two different pipelines, GAUFRE (Valentini et al. 2013) and Sp_Ace (Boeche & Grebel 2016), to determine atmospheric parameters and abundances by fixing log(g) to the seismic one. Our strategy ensures highly consistent values among all stellar parameters, leading to more accurate chemical abundances. A comparison of the chemical abundances obtained here with and without the use of seismic log(g) information has shown that an underestimated (overestimated) gravity leads to an underestimated (overestimated) elemental abundance (e.g. [Mg/H] is underestimated by ~0.25 dex when the gravity is underestimated by 0.5 dex). We then perform a comparison between the seismic gravities and the spectroscopic gravities presented in the RAVE DR4 catalogue, extracting a calibration for log(g) of RAVE giants in the colour interval 0.50<(J - Ks)<0.85. Finally, we show a comparison of the distances, temperatures, extinctions (and ages) derived here for our RAVE K2-C1 sample with those derived in RAVE DR4 and DR5.DR5 performs better than DR4 thanks to the seismic calibration, although discrepancies can still be important for objects for which the difference between DR4/DR5 and seismic gravities differ by more than ~0.5 dex. The method illustrated in this work will be used for analysing RAVE targets present in the other K2 campaigns, in the framework of Galactic Archaeology investigations.

Data Release 5 (DR5) of the Radial Velocity Experiment (RAVE) is the fifth data release from a magnitude-limited (9< I < 12) survey of stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra ($R\sim7500$) covering the Ca-triplet region (8410-8795\AA) span the complete time frame from the start of RAVE observations in 2003 to their completion in 2013. Radial velocities from 520,781 spectra of 457,588 unique stars are presented, of which 255,922 stellar observations have parallaxes and proper motions from the Tycho-Gaia astrometric solution (TGAS) in Gaia DR1. For our main DR5 catalog, stellar parameters (effective temperature, surface gravity, and overall metallicity) are computed using the RAVE DR4 stellar pipeline, but calibrated using recent K2 Campaign 1 seismic gravities and Gaia benchmark stars, as well as results obtained from high-resolution studies. Also included are temperatures from the Infrared Flux Method, and we provide a catalogue of red giant stars in the dereddened color $(J-Ks)_0$ interval (0.50,0.85) for which the gravities were calibrated based only on seismology. Further data products for sub-samples of the RAVE stars include individual abundances for Mg, Al, Si, Ca, Ti, Fe, and Ni, and distances found using isochrones. Each RAVE spectrum is complemented by an error spectrum, which has been used to determine uncertainties on the parameters. The data can be accessed via the RAVE Web site or the Vizier database.

In the $\Lambda$CDM model of structure formation, a stellar spheroid grows by the assembly of smaller galaxies, the so-called building blocks. Combining the Munich-Groningen semi-analytical model of galaxy formation with the high resolution Aquarius simulations of dark matter haloes, we study the assembly history of the stellar spheroids of six Milky Way-mass galaxies, focussing on building block properties such as mass, age and metallicity. These properties are compared to those of the surviving satellites in the same models. We find that the building blocks have higher star formation rates on average, and this is especially the case for the more massive objects. At high redshift these dominate in star formation over the satellites, whose star formation timescales are longer on average. These differences ought to result in a larger $\alpha$-element enhancement from Type II supernovae in the building blocks (compared to the satellites) by the time Type Ia supernovae would start to enrich them in iron, explaining the observational trends. Interestingly, there are some variations in the star formation timescales of the building blocks amongst the simulated haloes, indicating that [$\alpha$/Fe] abundances in spheroids of other galaxies might differ from those in our own Milky Way.

In the current concordance cosmology small halos are expected to be completely dark and can significantly perturb low-mass galaxies during minor merger interactions. These interactions may well contribute to the diversity of the dwarf galaxy population. Dwarf galaxies in the field are often observed to have peculiarities in their structure, morphology, and kinematics, as well as strong bursts of star formation without apparent cause. We aim to characterize the signatures of minor mergers of dwarf galaxies with dark satellites to aid their observational identification. We explore and quantify a variety of structural, morphological, and kinematic indicators of merging dwarf galaxies and their remnants using a suite of hydrodynamical simulations. The most sensitive indicators of mergers with dark satellites are large asymmetries in the gaseous and stellar distributions, enhanced central surface brightness and starbursts, and velocity offsets and misalignments between the cold gas and stellar components. In general, merging systems span a wide range of values of the most commonly used indicators, while isolated objects tend to have more confined values. Interestingly, we find in our simulations that a significantly off-centered burst of star formation can pinpoint the location of the dark satellite. Observational systems with such characteristics are perhaps the most promising for unveiling the presence of the hitherto, missing satellites.

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

Stellar streams result from the tidal disruption of satellites and star clusters as they orbit a host galaxy, and can be very sensitive probes of the gravitational potential of the host system. We select and study narrow stellar streams formed in a Milky-Way-like dark matter halo of the Aquarius suite of cosmological simulations, to determine if these streams can be used to constrain the present day characteristic parameters of the halo's gravitational potential. We find that orbits integrated in static spherical and triaxial NFW potentials both reproduce the locations and kinematics of the various streams reasonably well. To quantify this further, we determine the best-fit potential parameters by maximizing the amount of clustering of the stream stars in the space of their actions. We show that using our set of Aquarius streams, we recover a mass profile that is consistent with the spherically-averaged dark matter profile of the host halo, although we ignored both triaxiality and time evolution in the fit. This gives us confidence that such methods can be applied to the many streams that will be discovered by the Gaia mission to determine the gravitational potential of our Galaxy.

We model the time evolution of gaps in tidal streams caused by the impact of a dark matter subhalo, while both orbit a spherical gravitational potential. To this end, we make use of the simple behaviour of orbits in action-angle space. A gap effectively results from the divergence of two nearby orbits whose initial phase-space separation is, for very cold thin streams, largely given by the impulse induced by the subhalo. We find that in a spherical potential the size of a gap increases linearly with time for sufficiently long timescales. We have derived an analytic expression that shows how the growth rate depends on the mass of the perturbing subhalo, its scale and its relative velocity with respect to the stream. We have verified these scalings using N-body simulations and find excellent agreement. For example, a subhalo of mass 10^8 Msun directly impacting a very cold thin stream on an inclined orbit can induce a gap that may reach a size of several tens of kpc after a few Gyr. The gap size fluctuates importantly with phase on the orbit, and it is largest close to pericentre. This indicates that it may not be fully straightforward to invert the spectrum of gaps present in a stream to recover the mass spectrum of the subhalos.

We perform the first self-consistent measurement of the rate of interactions between stellar tidal streams created by disrupting satellites and dark subhalos in a cosmological simulation of a Milky-Way-mass galaxy. Using a retagged version of the Aquarius A dark-matter-only simulation, we selected 18 streams of tagged star particles that appear thin at the present day and followed them from the point their progenitors accrete onto the main halo, recording in each snapshot the characteristics of all dark-matter subhalos passing within several distance thresholds of any tagged star particle in each stream. We considered distance thresholds corresponding to constant impact parameters (1, 2, and 5 kpc), as well as those proportional to the region of influence of each subhalo (one and two times its half-mass radius $r_{1/2}$). We then measured the age and present-day, phase-unwrapped length of each stream in order to compute the interaction rate in different mass bins and for different thresholds, and compared these to analytic predictions from the literature. We measure a median rate of $1.5^{+3.0}_{-1.1}\ (9.1^{+17.5}_{-7.1},\ 61.8^{+211}_{-40.6})$ interactions within 1 (2, 5) kpc of the stream per 10 kpc of stream length per 10 Gyr. Resolution effects (both time and particle number) affect these estimated rates by lowering them.

We have studied the behaviour of stellar streams in the Aquarius fully cosmological N-body simulations of the formation of Milky Way halos. In particular, we have characterised the streams in angle/frequency spaces derived using an approximate but generally well-fitting spherical potential. We have also run several test-particle simulations to understand and guide our interpretation of the different features we see in the Aquarius streams. Our goal is both to establish which deviations of the expected action-angle behaviour of streams exist because of the approximations made on the potential, but also to derive to what degree we can use these coordinates to model streams reliably. We have found that many of the Aquarius streams wrap in angle space along relatively straight lines, and also in frequency space. On the other hand, from our controlled simulations we have been able to establish that deviations from spherical symmetry, the use of incorrect potentials and the inclusion of self-gravity lead to streams in angle space to still be along relatively straight lines but also to depict wiggly behaviour whose amplitude increases as the approximation to the true potential becomes worse. In frequency space streams typically become thicker and somewhat distorted. Therefore, our analysis explains most of the features seen in the approximate angle and frequency spaces for the Aquarius streams with the exception of their somewhat `noisy' and `patchy' morphologies. These are likely due to the interactions with the large number of dark matter subhalos present in the cosmological simulations. Since the measured angle-frequency misalignments of the Aquarius streams can largely be attributed to using the wrong (spherical) potential, determining the mass growth history of these halos will only be feasible once the true potential has been determined robustly.

Stellar halos and globular cluster (GC) systems contain valuable information regarding the assembly history of their host galaxies. Motivated by the detection of a significant rotation signal in the outer halo GC system of M31, we investigate the likelihood of detecting such a rotation signal in projection, using cosmological simulations. To this end we select subsets of tagged particles in the halos of the Aquarius simulations to represent mock GC systems, and analyse their kinematics. We find that GC systems can exhibit a non-negligible rotation signal provided the associated stellar halo also has a net angular momentum. The ability to detect this rotation signal is highly dependent on the viewing perspective, and the probability of seeing a signal larger than that measured in M31 ranges from 10% to 90% for the different halos in the Aquarius suite. High values are found from a perspective such that the projected angular momentum of the GC system is within 40 deg of the rotation axis determined via the projected positions and line-of-sight velocities of the GCs. Furthermore, the true 3D angular momentum of the outer stellar halo is relatively well aligned, within 35 deg, with that of the mock GC systems. We argue that the net angular momentum in the mock GC systems arises naturally when the majority of the material is accreted from a preferred direction, namely along the dominant dark matter filament of the large-scale structure that the halos are embedded in. This, together with the favourable edge-on view of M31's disk suggests that it is not a coincidence that a large rotation signal has been measured for its outer halo GC system.

We use the high-resolution Aquarius cosmological dark matter simulations coupled to the semi-analytic model by Starkenburg et al. (2013) to study the HI content and velocity width properties of field galaxies at the low mass end in the context of $\Lambda$CDM. We compare our predictions to the observed ALFALFA survey HI mass and velocity width functions, and find very good agreement without fine-tuning, when considering central galaxies. Furthermore, the properties of the dark matter halos hosting galaxies, characterised by their peak velocity and circular velocity at 2 radial disk scalelengths overlap perfectly with the inferred values from observations. This suggests that our galaxies are placed in the right dark matter halos, and consequently at face value, we do not find any discrepancy with the predictions from the $\Lambda$CDM model. Our analysis indicates that previous tensions, apparent when using abundance matching models, arise because this technique cannot be straightforwardly applied for objects with masses $M_{vir} < 10^{10} M_{\odot}$.

We present the high resolution spectroscopic study of five -3.9<=[Fe/H]<=-2.5 stars in the Local Group dwarf spheroidal, Sculptor, thereby doubling the number of stars with comparable observations in this metallicity range. We carry out a detailed analysis of the chemical abundances of alpha, iron peak, light and heavy elements, and draw comparisons with the Milky Way halo and the ultra faint dwarf stellar populations. We show that the bulk of the Sculptor metal-poor stars follows the same trends in abundance ratios versus metallicity as the Milky Way stars. This suggests similar early conditions of star formation and a high degree of homogeneity of the interstellar medium. We find an outlier to this main regime, which seems to miss the products of the most massive of the TypeII supernovae. In addition to its value to help refining galaxy formation models, this star provides clues to the production of cobalt and zinc. Two of our sample stars have low odd-to-even barium isotope abundance ratios, suggestive of a fair proportion of s-process; we discuss the implication for the nucleosynthetic origin of the neutron capture elements.

It has been proposed that mergers induce starbursts and lead to important morphological changes in galaxies. Most studies so far have focused on large galaxies, but dwarfs might also experience such events, since the halo mass function is scale-free in the concordance cosmological model. Notably, because of their low mass, most of their interactions will be with dark satellites. In this paper we follow the evolution of gas-rich disky dwarf galaxies as they experience a minor merger with a dark satellite. We aim to characterize the effects of such an interaction on the dwarf's star formation, morphology, and kinematical properties. We performed a suite of carefully set-up hydrodynamical simulations of dwarf galaxies that include dark matter, gas, and stars merging with a satellite consisting solely of dark matter. For the host system we vary the gas fraction, disk size and thickness, halo mass, and concentration, while we explore different masses, concentrations, and orbits for the satellite. We find that the interactions cause strong starbursts of both short and long duration in the dwarfs. Their star formation rates increase by factors of a few to 10 or more. They are strongest for systems with extended gas disks and high gas fractions merging with a high-concentration satellite on a planar, radial orbit. In contrast to analogous simulations of Milky Way-mass galaxies, many of the systems experience strong morphological changes and become spheroidal even in the presence of significant amounts of gas. The simulated systems compare remarkably well with the observational properties of a large selection of irregular dwarf galaxies and blue compact dwarfs. This implies that mergers with dark satellites might well be happening but not be fully evident, and may thus play a role in the diversity of the dwarf galaxy population.

We study the evolution of streams in a time-dependent spherical gravitational potential. Our goal is to establish what are the imprints of this time evolution on the properties of streams as well as their observability. To this end, we have performed a suite of numerical experiments for a host system that doubles its mass during the integration time and for a variety of initial conditions. In these experiments we found that the most striking imprint is a misalignment of 10 degrees in the angular location of the apocentres of the streams compared to the static case (and to the orbit of the centre of mass), which only becomes apparent for sufficiently long streams. We have also developed an analytic model using action-angle variables which allows us to explain this behaviour and to identify the most important signature of time evolution, namely a difference in the slope defined by the distribution of particles along a stream in frequency and in angle space. Although a difference in slope can arise when the present-day potential is not correctly modelled, this shortcoming can be by-passed because in this case, streams are no longer straight lines in angle space, but depict a wiggly appearance and an implausible energy gradient. The difference in slope due to time-evolution is small, typically $\sim10^{-2}$ and its amplitude depends on the growth rate of the potential, but nonetheless we find that it could be observable if accurate full-space information for nearby long streams is available.

We compare the dark matter halos' structural parameters derived for four Milky Way dwarf spheroidal galaxies to those of subhalos found in cosmological $N$-body simulations. We confirm that estimates of the mass at a single fixed radius are fully consistent with the observations. However, when a second structural parameter such as the logarithmic slope of the dark halo density profile measured close to the half-light radius is included in the comparison, we find little to no overlap between the satellites and the subhalos. Typically the right mass subhalos have steeper profiles at these radii than measurements of the dSph suggest. Using energy arguments we explore if it is possible to solve this discrepancy by invoking baryonic effects. Assuming that feedback from supernovae can lead to a reshaping of the halos, we compute the required efficiency and find entirely plausible values for a significant fraction of the subhalos and even as low as 0.1%. This implies that care must be taken not to exaggerate the effect of supernovae feedback as this could make the halos too shallow. These results could be used to calibrate and possibly constrain feedback recipes in hydrodynamical simulations.

We study the kinematics of a local sample of stars, located within a cylinder of 500 pc radius centered on the Sun, in the RAVE dataset. We find clear asymmetries in the $v_R$-$v_\phi$ velocity distributions of thin and thick disk stars: here are more stars moving radially outwards for low azimuthal velocities and more radially inwards for high azimuthal velocities. Such asymmetries have been previously reported for the thin disk as being due to the Galactic bar, but this is the first time that the same type of structures are seen in the thick disk. Our findings imply that the velocities of thick disk stars should no longer be described by Schwarzschild's, multivariate Gaussian or purely axisymmetric distributions. Furthermore, the nature of previously reported substructures in the thick disk needs to be revisited as these could be associated with dynamical resonances rather than to accretion events. It is clear that dynamical models of the Galaxy must fit the 3D velocity distributions of the disks, rather than the projected 1D, if we are to understand the Galaxy fully.

We present a proper-motion study on models of the dwarf spheroidal galaxy Sculptor, based on the predicted proper-motion accuracy of Gaia measurements. Gaia will measure proper motions of several hundreds of stars for a Sculptor-like system. Even with an uncertainty on the proper motion of order 1.5 times the size of an individual proper-motion value of ~10 mas/century, we find that it is possible to recover Sculptor's systemic proper motion at its distance of 79 kpc.

We present a new method for constraining the Milky Way halo gravitational potential by simultaneously fitting multiple tidal streams. This method requires full three-dimensional positions and velocities for all stars to be fit, but does not require identification of any specific stream or determination of stream membership for any star. We exploit the principle that the action distribution of stream stars is most clustered when the potential used to calculate the actions is closest to the true potential. Clustering is quantified with the Kullback-Leibler Divergence (KLD), which also provides conditional uncertainties for our parameter estimates. We show, for toy Gaia-like data in a spherical isochrone potential, that maximizing the KLD of the action distribution relative to a smoother distribution recovers the true values of the potential parameters. The precision depends on the observational errors and the number of streams in the sample; using KIII giants as tracers, we measure the enclosed mass at the average radius of the sample stars accurate to 3% and precise to 20-40%. Recovery of the scale radius is precise to 25%, and is biased 50% high by the small galactocentric distance range of stars in our mock sample (1-25 kpc, or about three scale radii, with mean 6.5 kpc). About 15 streams, with at least 100 stars per stream, are needed to obtain upper and lower bounds on the enclosed mass and scale radius when observational errors are taken into account; 20-25 streams are required to stabilize the size of the confidence interval. If radial velocities are provided for stars out to 100 kpc (10 scale radii), all parameters can be determined with 10% accuracy and 20% precision (1.3% accuracy in the case of the enclosed mass), underlining the need for ground-based spectroscopic follow-up to complete the radial velocity catalog for faint halo stars observed by Gaia.

Using the RAdial Velocity Experiment fourth data release (RAVE DR4), and a new metallicity calibration that will be also taken into account in the future RAVE DR5, we investigate the existence and the properties of super-solar metallicity stars ([M/H] > +0.1 dex) in the sample, and in particular in the Solar neighbourhood. We find that RAVE is rich in super-solar metallicity stars, and that the local metallicity distribution function declines remarkably slowly up to +0.4 dex. Our results show that the kinematics and height distributions of the super-solar metallicity stars are identical to those of the [M/H] < 0 thin-disc giants that we presume were locally manufactured. The eccentricities of the super-solar metallicity stars indicate that half of them are on a roughly circular orbit (e < 0.15), so under the assumption that the metallicity of the interstellar medium at a given radius never decreases with time, they must have increased their angular momenta by scattering at corotation resonances of spiral arms from regions far inside the Solar annulus. The likelihood that a star will migrate radially does not seem to decrease significantly with increasing amplitude of vertical oscillations within range of oscillation amplitudes encountered in the disc.

We aim to characterize high-velocity (HiVel) stars in the solar vicinity both chemically and kinematically using the fourth data release of the RAdial Velocity Experiment (RAVE). We used a sample of 57 HiVel stars with Galactic rest-frame velocities larger than 275 km s$^{-1}$. With 6D position and velocity information, we integrated the orbits of the HiVel stars and found that, on average, they reach out to 13 kpc from the Galactic plane and have relatively eccentric orbits consistent with the Galactic halo. Using the stellar parameters and [$\alpha$/Fe] estimates from RAVE, we found the metallicity distribution of the HiVel stars peak at [M/H] = -1.2 dex and is chemically consistent with the inner halo. There are a few notable exceptions that include a hypervelocity star (HVS) candidate, an extremely high-velocity bound halo star, and one star that is kinematically consistent with the halo but chemically consistent with the disk. High-resolution spectra were obtained for the metal-rich HiVel star candidate and the second highest velocity star in the sample. Using these high-resolution data, we report the discovery of a metal-rich halo star that has likely been dynamically ejected into the halo from the Galactic thick disk. This discovery could aid in explaining the assembly of the most metal-rich component of the Galactic halo.

In the context of the current $\Lambda$CDM cosmological model small dark matter haloes are abundant and satellites of dwarf galaxies are expected to be predominantly dark. Since low mass galaxies have smaller baryon fractions interactions with these satellites may leave particularly dramatic imprints. We uncover the influence of the most massive of these dark satellites on disky dwarf galaxies and the possible dynamical and morphological transformations that result from these interactions. We use a suite of carefully set-up, controlled simulations of isolated dwarf galaxies. The primary dwarf galaxies have solely a stellar disk in the dark matter halo and the secundaries are completely devoid of baryons. We vary the disk mass, halo concentration, initial disk thickness and inclination of the satellite orbit. The disky dwarf galaxies are heated and disrupted due to the minor merger event, more extremely for higher satellite over disk mass ratios, and the morphology and kinematics are significantly altered. Moreover, for less concentrated haloes the minor merger can completely destroy the disk leaving a low-luminosity spheroidal-like galaxy instead. We conclude that dwarf galaxies are very much susceptible to being disturbed by dark galaxies and that even a minor merger event can significantly disrupt and alter the structure and kinematics of a dwarf galaxy. This process may be seen as a new channel for the formation of dwarf spheroidal galaxies.

We use the kinematics of $\sim200\,000$ giant stars that lie within $\sim 1.5$ kpc of the plane to measure the vertical profile of mass density near the Sun. We find that the dark mass contained within the isodensity surface of the dark halo that passes through the Sun ($(6\pm0.9)\times10^{10}\,\mathrm{M_\odot}$), and the surface density within $0.9$ kpc of the plane ($(69\pm10)\,\mathrm{M_\odot\,pc^{-2}}$) are almost independent of the (oblate) halo's axis ratio $q$. If the halo is spherical, 46 per cent of the radial force on the Sun is provided by baryons, and only 4.3 per cent of the Galaxy's mass is baryonic. If the halo is flattened, the baryons contribute even less strongly to the local radial force and to the Galaxy's mass. The dark-matter density at the location of the Sun is $0.0126\,q^{-0.89}\,\mathrm{M_\odot\,pc^{-3}}=0.48\,q^{-0.89}\,\mathrm{GeV\,cm^{-3}}$. When combined with other literature results we find hints for a mildly oblate dark halo with $q \simeq 0.8$. Our value for the dark mass within the solar radius is larger than that predicted by cosmological dark-matter-only simulations but in good agreement with simulations once the effects of baryonic infall are taken into account. Our mass models consist of three double-exponential discs, an oblate bulge and a Navarro-Frenk-White dark-matter halo, and we model the dynamics of the RAVE stars in the corresponding gravitational fields by finding distribution functions $f(\mathbf{J})$ that depend on three action integrals. Statistical errors are completely swamped by systematic uncertainties, the most important of which are the distance to the stars in the photometric and spectroscopic samples and the solar distance to the Galactic centre. Systematics other than the flattening of the dark halo yield overall uncertainties $\sim 15$ per cent.

We develop a new method to account for the finite lifetimes of stars and trace individual abundances within a semi-analytic model of galaxy formation. At variance with previous methods, based on the storage of the (binned) past star formation history of model galaxies, our method projects the information about the metals produced by each simple stellar population (SSP) in the future. Using this approach, an accurate accounting of the timings and properties of the individual SSPs composing model galaxies is possible. We analyse the dependence of our chemical model on various ingredients, and apply it to six simulated haloes of roughly Milky Way mass and with no massive close neighbour at z=0. For all models considered, the [Fe/H] distributions of the stars in the disc component are in good agreement with Milky Way data, while for the spheroid component (whose formation we model only through mergers) these are offset low with respect to observational measurements for the Milky Way bulge. This is a consequence of narrow star formation histories, with relatively low rates of star formation. The slow recycling of gas and energy from supernovae in our chemical model has important consequences on the predicted star formation rates, which are systematically lower than the corresponding rates in the same physical model but with an instantaneous recycling approximation. The halo that resembles most our Galaxy in terms of its global properties also reproduces the observed relation between the average metallicity and luminosity of the Milky Way satellites, albeit with a slightly steeper slope.

Fornax is one of the most massive dwarf spheroidal galaxies in the Local Group. The Fornax field star population is dominated by intermediate age stars but star formation was going on over almost its entire history. It has been proposed that Fornax experienced a minor merger event. Despite recent progress, only the high metallicity end of Fornax field stars ([Fe/H]>-1.2 dex) has been sampled in larger number via high resolution spectroscopy. We want to better understand the full chemical evolution of this galaxy by better sampling the whole metallicity range, including more metal poor stars. We use the VLT-FLAMES multi-fibre spectrograph in high-resolution mode to determine the abundances of several alpha, iron-peak and neutron-capture elements in a sample of 47 individual Red Giant Branch stars in the Fornax dwarf spheroidal galaxy. We combine these abundances with accurate age estimates derived from the age probability distribution from the colour-magnitude diagram of Fornax. Similar to other dwarf spheroidal galaxies, the old, metal-poor stars of Fornax are typically alpha-rich while the young metal-rich stars are alpha-poor. In the classical scenario of the time delay between SNe II and SNe Ia, we confirm that SNe Ia started to contribute to the chemical enrichment at [Fe/H] between -2.0 and -1.8 dex. We find that the onset of SNe Ia took place between 12-10 Gyrs ago. The high values of [Ba/Fe], [La/Fe] reflect the influence of SNe Ia and AGB stars in the abundance pattern of the younger stellar population of Fornax. Our findings of low [alpha/Fe] and enhanced [Eu/Mg] are compatible with an initial mass function that lacks the most massive stars and with star formation that kept going on throughout the whole history of Fornax. We find that massive stars kept enriching the interstellar medium in alpha-elements, although they were not the main contributor to the iron enrichment.

In previous work we have presented Schwarzschild models of the Sculptor dSph, demonstrating that this system could be embedded in dark matter halos that are either cusped or cored. Here we show that the non-parametric distribution function recovered through Schwarschild's method is bimodal in energy and angular momentum space for all best fitting mass models explored. We demonstrate that this bimodality is directly related to the two components known to be present in Sculptor through stellar populations analysis, although our method is purely dynamical in nature and does not use this prior information. It therefore constitutes independent confirmation of the existence of two physically distinct dynamical components in Sculptor and suggests a rather complex assembly history for this dwarf galaxy.

Aims: We study single and binary white dwarfs in the inner halo of the Milky Way in order to learn more about the conditions under which the population of halo stars was born, such as the initial mass function (IMF), the star formation history, or the binary fraction. Methods: We simulate the evolution of low-metallicity halo stars at distances up to ~ 3 kpc using the binary population synthesis code SeBa. We use two different white dwarf cooling models to predict the present-day luminosities of halo white dwarfs. We determine the white dwarf luminosity functions (WDLFs) for eight different halo models and compare these with the observed halo WDLF of white dwarfs in the SuperCOSMOS Sky Survey. Furthermore, we predict the properties of binary white dwarfs in the halo and determine the number of halo white dwarfs that is expected to be observed with the Gaia satellite. Results: By comparing the WDLFs, we find that a standard IMF matches the observations more accurately than a top-heavy one, but the difference with a bottom-heavy IMF is small. A burst of star formation 13 Gyr ago fits slightly better than a star formation burst 10 Gyr ago and also slightly better than continuous star formation $10-13$ Gyrs ago. Gaia will be the first instrument to constrain the bright end of the field halo WDLF, where contributions from binary WDs are considerable. Many of these will have He cores, of which a handful have atypical surface gravities ($\log g < 6$) and reach luminosities $\log(L/L_\odot) > 0$ in our standard model for WD cooling. These so called pre-WDs, if observed, can help us to constrain white dwarf cooling models and might teach us something about the fraction of halo stars that reside in binaries.

Stellar population studies of globular clusters have suggested that the brightest clusters in the Galaxy might actually be the remnant nuclei of dwarf spheroidal galaxies. If the present Galactic globular clusters formed within larger stellar systems, they are likely surrounded by extra-tidal halos and/or tails made up of stars that were tidally stripped from their parent systems. The stellar surroundings around globular clusters are therefore one of the best places to look for the remnants of an ancient dwarf galaxy. Here an attempt is made to search for tidal debris around the supernovae enriched globular clusters M22 and NGC 1851 as well as the kinematically unique cluster NGC 3201. The stellar parameters from the Radial Velocity Experiment (RAVE) are used to identify stars with RAVE metallicities, radial velocities and elemental-abundances consistent with the abundance patterns and properties of the stars in M22, NGC 1851 and NGC 3201. The discovery of RAVE stars that may be associated with M22 and NGC 1851 are reported, some of which are at projected distances of ~10 degrees away from the core of these clusters. Numerous RAVE stars associated with NGC 3201 suggest that either the tidal radius of this cluster is underestimated, or that there are some unbound stars extending a few arc minutes from the edge of the cluster's radius. No further extra-tidal stars associated with NGC 3201 could be identified. The bright magnitudes of the RAVE stars make them easy targets for high resolution follow-up observations, allowing an eventual further chemical tagging to solidify (or exclude) stars outside the tidal radius of the cluster as tidal debris. In both our radial velocity histograms of the regions surrounding NGC 1851 and NGC 3201, a peak of stars at 230 km/s is seen, consistent with extended tidal debris from omega Centauri.

We provide APASS photometry in the Landolt BV and Sloan g'r'i' bands for all the 425,743 stars included in the latest 4th RAVE Data Release. The internal accuracy of the APASS photometry of RAVE stars, expressed as error of the mean of data obtained and separately calibrated over a median of 4 distinct observing epochs and distributed between 2009 and 2013, is 0.013, 0.012, 0.012, 0.014 and 0.021 mag for B, V, g', r' and i' band, respectively. The equally high external accuracy of APASS photometry has been verified on secondary Landolt and Sloan photometric standard stars not involved in the APASS calibration process, and on a large body of literature data on field and cluster stars, confirming the absence of offsets and trends. Compared with the Carlsberg Meridian Catalog (CMC-15), APASS astrometry of RAVE stars is accurate to a median value of 0.098 arcsec. Brightness distribution functions for the RAVE stars have been derived in all bands. APASS photometry of RAVE stars, augmented by 2MASS JHK infrared data, has been chi2 fitted to a densely populated synthetic photometric library designed to widely explore in temperature, surface gravity, metallicity and reddening. Resulting Teff and E(B-V), computed over a range of options, are provided and discussed, and will be kept updated in response to future APASS and RAVE data releases. In the process it is found that the reddening caused by an homogeneous slab of dust, extending for 140 pc on either side of the Galactic plane and responsible for E(B-V,poles)=0.036 +/- 0.002 at the galactic poles, is a suitable approximation of the actual reddening encountered at Galactic latitudes |b|>=25 deg.

The diffuse interstellar bands (DIBs) are absorption lines observed in visual and near infrared spectra of stars. Understanding their origin in the interstellar medium (ISM) is one of the oldest problems in astronomical spectroscopy, as DIBs ahave been known since 1922. In a completely new approach to understanding DIBs, we combined information from nearly 500,000 stellar spectra obtained by the massive spectroscopic survey RAVE (Radial Velocity Experiment) to produce the first pseudo three-dimensional map of the strength of the DIB at 8620~{\AA}ngstroms covering the nearest 3 kiloparsecs from the Sun, and show that it follows our independently constructed spatial distribution of extinction by interstellar dust along the Galactic plane. Despite having a similar distribution in the Galactic plane, the DIB 8620 carrier has a significantly larger vertical scale height than the dust. Even if one DIB may not represent the general DIB population, our observations outline the future direction of DIB research.

Aims: We investigate whether the cylindrical (galactocentric) radial velocity gradient of ~ -3 km/s/kpc, directed radially from the Galactic center and recently observed in the stars of the solar neighborhood with the RAVE survey, can be explained by the resonant effects of the bar near the solar neighborhood. Methods: We compared the results of test particle simulations of the Milky Way with a potential that includes a rotating bar with observations from the RAVE survey. To this end we applied the RAVE selection function to the simulations and convolved these with the characteristic RAVE errors. We explored different "solar neighborhoods" in the simulations, as well as different bar models Results: We find that the bar induces a negative radial velocity gradient at every height from the Galactic plane, outside the outer Lindblad resonance and for angles from the long axis of the bar compatible with the current estimates. The selection function and errors do not wash away the gradient, but often make it steeper, especially near the Galactic plane, because this is where the RAVE survey is less radially extended. No gradient in the vertical velocity ispresent in our simulations, from which we may conclude that this cannot be induced by the bar.

We present a study of the chemical abundances of Solar neighbourhood stars associated to dynamical structures in the Milky Way's (thick) disk. These stars were identified as overdensity in the eccentricity range 0.3< ecc < 0.5 in the Copenhagen-Geneva Survey by Helmi et al. (2006). We find that the stars with these dynamical characteristics do not constitute a homogeneous population. A relatively sharp transition in dynamical and chemical properties appears to occur at a metallicity of [Fe/H] ~ -0.4. Stars with [Fe/H] > -0.4 have mostly lower eccentricities, smaller vertical velocity dispersions, are alpha-enhanced and define a rather narrow sequence in [alpha/Fe] vs [Fe/H], clearly distinct from that of the thin disk. Stars with [Fe/H] < -0.4 have a range of eccentricities, are hotter vertically, and depict a larger spread in [alpha/Fe]. We have also found tentative evidence of substructure possibly associated to the disruption of a metal-rich star cluster. The differences between these populations of stars is also present in e.g. [Zn/Fe], [Ni/Fe] and [SmII/Fe], suggesting a real physical distinction.

We investigate the kinematic parameters of the Milky Way disc using the RAVE and GCS stellar surveys. We do this by fitting a kinematic model to the data taking the selection function of the data into account. For stars in the GCS we use all phase-space coordinates, but for RAVE stars we use only $(l,b,v_{\rm los})$. Using MCMC technique, we investigate the full posterior distributions of the parameters given the data. We investigate the `age-velocity dispersion' relation for the three kinematic components ($\sigma_R,\sigma_{\phi},\sigma_z$), the radial dependence of the velocity dispersions, the Solar peculiar motion ($U_{\odot},V_{\odot}, W_{\odot} $), the circular speed $\Theta_0$ at the Sun and the fall of mean azimuthal motion with height above the mid-plane. We confirm that the Besan\c{c}on-style Gaussian model accurately fits the GCS data, but fails to match the details of the more spatially extended RAVE survey. In particular, the Shu distribution function (DF) handles non-circular orbits more accurately and provides a better fit to the kinematic data. The Gaussian distribution function not only fits the data poorly but systematically underestimates the fall of velocity dispersion with radius. We find that correlations exist between a number of parameters, which highlights the importance of doing joint fits. The large size of the RAVE survey, allows us to get precise values for most parameters. However, large systematic uncertainties remain, especially in $V_{\odot}$ and $\Theta_0$. We find that, for an extended sample of stars, $\Theta_0$ is underestimated by as much as $10\%$ if the vertical dependence of the mean azimuthal motion is neglected. Using a simple model for vertical dependence of kinematics, we find that it is possible to match the Sgr A* proper motion without any need for $V_{\odot}$ being larger than that estimated locally by surveys like GCS.

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

Probability density functions are determined from new stellar parameters for the distance moduli of stars for which the RAdial Velocity Experiment (RAVE) has obtained spectra with S/N>=10. Single-Gaussian fits to the pdf in distance modulus suffice for roughly half the stars, with most of the other half having satisfactory two-Gaussian representations. As expected, early-type stars rarely require more than one Gaussian. The expectation value of distance is larger than the distance implied by the expectation of distance modulus; the latter is itself larger than the distance implied by the expectation value of the parallax. Our parallaxes of Hipparcos stars agree well with the values measured by Hipparcos, so the expectation of parallax is the most reliable distance indicator. The latter are improved by taking extinction into account. The effective temperature absolute-magnitude diagram of our stars is significantly improved when these pdfs are used to make the diagram. We use the method of kinematic corrections devised by Schoenrich, Binney & Asplund to check for systematic errors for general stars and confirm that the most reliable distance indicator is the expectation of parallax. For cool dwarfs and low-gravity giants <pi> tends to be larger than the true distance by up to 30 percent. The most satisfactory distances are for dwarfs hotter than 5500 K. We compare our distances to stars in 13 open clusters with cluster distances from the literature and find excellent agreement for the dwarfs and indications that we are over-estimating distances to giants, especially in young clusters.

We analyse the kinematics of ~400000 RAVE stars. We split the sample into hot and cold dwarfs, red-clump and non-clump giants. The kinematics of the clump giants are consistent with being identical with those of non-clump giants. We fit Gaussian velocity ellipsoids to the meridional-plane components of velocity of each star class and give formulae from which the shape and orientation of the velocity ellipsoid can be determined at any location. The data are consistent with the giants and the cool dwarfs sharing the same velocity ellipsoids; sigma_z rises from 21 kms in the plane to sim 55 kms at |z|=2 kpc, while sigma_r rises from 37 kms to 82 kms. At (R,z) the longest axis of one of these velocity ellipsoids is inclined to the Galactic plane by an angle ~0.8 arctan(z/R). We use a novel formula to obtain precise fits to the highly non-Gaussian distributions of v_phi components. We compare the observed velocity distributions with the predictions of a dynamical model fitted to the velocities of stars that lie within ~150 pc of the Sun and star counts towards the Galactic pole. The model accurately reproduces the non-Gaussian nature of the v_r and v_z distributions and provides excellent fits to the data for v_z at all locations. The model v_phi distributions for the cool dwarfs fit the data extremely well, while those for the hot dwarfs have displacements to low v_phi that grow with |z| from very small values near the plane. At |z|>0.5 kpc, the theoretical v_phi distributions for giants show a deficit of stars with large v_phi and the model v_r distributions are too narrow. Systematically over-estimating distances by 20 per cent introduces asymmetry into the model v_r and v_z distributions near the plane and but significantly improves the fits to the data at |z|>0.5 kpc. The quality of the fits lends credence to the assumed, disc-dominated, gravitational potential.

We present a new simple functional form to model the evolution of a spherical mass distribution in a cosmological context. Two parameters control the growth of the system and this is modelled using a redshift dependent exponential for the scale mass and scale radius. In this new model, systems form inside out and the mass of a given shell can be made to never decrease, as generally expected. This feature makes it more suitable for studying the smooth growth of galactic potentials or cosmological halos than other parametrizations often used in the literature. This is further confirmed through a comparison to the growth of dark matter halos in the Aquarius simulations.

We analyze the Aquarius simulations to characterize the shape of dark matter halos with peak circular velocity in the range 8<Vmax<200 km/s, and perform a convergence study using the various Aquarius resolution levels. For the converged objects, we determine the principal axis (a<b<c) of the normalized inertia tensor as a function of radius. We find that the triaxiality of field halos is an increasing function of halo mass, so that the smallest halos in our sample are ~40-50% rounder than Milky Way-like objects at the radius where the circular velocity peaks, rmax. We find that the distribution of subhalo axis ratios is consistent with that of field halos of comparable Vmax. Inner and outer contours within each object are well aligned, with the major axis preferentially pointing in the radial direction for subhalos closest to the center of their host halo. We also analyze the dynamical structure of subhalos likely to host luminous satellites comparable to the classical dwarf spheroidals in the Local Group. These halos have axis ratios that increase with radius, and which are mildly triaxial with <b/a>~0.75 and <c/a>~0.60 at r~1 kpc. Their velocity ellipsoid become strongly tangentially biased in the outskirts as a consequence of tidal stripping.

Non-axisymmetries in the Galactic potential (spiral arms and bar) induce kinematic groups such as the Hercules stream. Assuming that Hercules is caused by the effects of the Outer Lindblad Resonance of the Galactic bar, we model analytically its properties as a function of position in the Galaxy and its dependence on the bar's pattern speed and orientation. Using data from the RAVE survey we find that the azimuthal velocity of the Hercules structure decreases as a function of Galactocentric radius, in a manner consistent with our analytical model. This allows us to obtain new estimates of the parameters of the Milky Way's bar. The combined likelihood function of the bar's pattern speed and angle has its maximum for a pattern speed of $\Omega_b=(1.89 \pm 0.08) \times \Omega_0$ where $\Omega_0$ is the local circular frequency. Assuming a Solar radius of 8.05 kpc and a local circular velocity of 238 km/s, this corresponds to $\Omega_b=56\pm 2$ km/s/kpc. On the other hand, the bar's orientation $\phi_b$ cannot be constrained with the available data. In fact, the likelihood function shows that a tight correlation exists between the pattern speed and the orientation, implying that a better description of our best fit results is given by the linear relation $\Omega_b/\Omega_0=1.91 + 0.0044[\phi_b(\deg)-48]$, with standard deviation of 0.02. For example, for an angle of $\phi_b=30\deg$ the pattern speed is $54.0\pm0.5$ km/s/kpc. These results are not very sensitive to the other Galactic parameters such as the circular velocity curve or the peculiar motion of the Sun, and are robust to biases in distance.

The velocity dispersions of stars near the Sun are known to increase with stellar age, but age can be difficult to determine so a proxy like the abundance of alpha elements (e.g., Mg) with respect to iron, [alpha/Fe], is used. Here we report an unexpected behavior found in the velocity dispersion of a sample of giant stars from the RAdial Velocity Experiment (RAVE) survey with high quality chemical and kinematical information, in that it decreases strongly for stars with [Mg/Fe] > 0.4 dex (i.e., those that formed in the first Gyr of the Galaxy's life). These findings can be explained by perturbations from massive mergers in the early Universe, which have affected more strongly the outer parts of the disc, and the subsequent radial migration of stars with cooler kinematics from the inner disc. Similar reversed trends in velocity dispersion are also found for different metallicity subpopulations. Our results suggest that the Milky Way disc merger history can be recovered by relating the observed chemo-kinematic relations to the properties of past merger events.

Aims: To study the imprints of a rotating bar on the kinematics of stars in the thin and thick disks throughout the Galaxy. Methods: We perform test particle numerical simulations of the thin and thick disks in a 3D Galactic potential that includes a halo, a bulge, thin and thick disks, and a Ferrers bar. We analyze the resulting velocity distributions of populations corresponding to both disks, for different positions in the Galaxy and for different structural parameters of the bar. Results: We find that the velocity distributions of the disks are affected by the bar, and that strong transient effects are present for approximately 10 bar rotations after this is introduced adiabatically. On long (more realistic) timescales, the effects of the bar are strong on the kinematics of thin disk stars, and weaker on those in the thick disk, but in any case significant. Furthermore, we find that it is possible to trace the imprints of the bar also vertically and at least up to z ~ 1 kpc for the thin disk and z ~ 2 kpc for the thick disk.

We construct new estimates on the Galactic escape speed at various Galactocentric radii using the latest data release of the Radial Velocity Experiment (RAVE DR4). Compared to previous studies we have a database larger by a factor of 10 as well as reliable distance estimates for almost all stars. Our analysis is based on the statistical analysis of a rigorously selected sample of 90 high-velocity halo stars from RAVE and a previously published data set. We calibrate and extensively test our method using a suite of cosmological simulations of the formation of Milky Way-sized galaxies. Our best estimate of the local Galactic escape speed, which we define as the minimum speed required to reach three virial radii $R_{340}$, is $533^{+54}_{-41}$ km/s (90% confidence) with an additional 5% systematic uncertainty, where $R_{340}$ is the Galactocentric radius encompassing a mean over-density of 340 times the critical density for closure in the Universe. From the escape speed we further derive estimates of the mass of the Galaxy using a simple mass model with two options for the mass profile of the dark matter halo: an unaltered and an adiabatically contracted Navarro, Frenk & White (NFW) sphere. If we fix the local circular velocity the latter profile yields a significantly higher mass than the un-contracted halo, but if we instead use the statistics on halo concentration parameters in large cosmological simulations as a constraint we find very similar masses for both models. Our best estimate for $M_{340}$, the mass interior to $R_{340}$ (dark matter and baryons), is $1.3^{+0.4}_{-0.3} \times 10^{12}$ M$_\odot$ (corresponding to $M_{200} = 1.6^{+0.5}_{-0.4} \times 10^{12}$ M$_\odot$). This estimate is in good agreement with recently published independent mass estimates based on the kinematics of more distant halo stars and the satellite galaxy Leo I.

We use the very high resolution, fully cosmological simulations from the Aquarius project, coupled to a semi-analytical model of galaxy formation, to study the phase-space distribution of halo stars in "solar neighbourhood"-like volumes. We find that this distribution is very rich in substructure in the form of stellar streams for all five stellar haloes we have analysed. These streams can be easily identified in velocity space, as well as in spaces of pseudo-conserved quantities such as E vs. Lz. In our best-resolved local volumes, the number of identified streams ranges from ~ 300 to 600, in very good agreement with previous analytical predictions, even in the presence of chaotic mixing. The fraction of particles linked to (massive) stellar streams in these volumes can be as large as 84%. The number of identified streams is found to decrease as a power-law with galactocentric radius. We show that the strongest limitation to the quantification of substructure in our poorest-resolved local volumes is particle resolution rather than strong diffusion due to chaotic mixing.

The RAVE survey, combined with proper motions and distance estimates, can be used to study in detail stellar kinematics in the extended solar neighbourhood (solar suburb). Using the red clump, we examine the mean velocity components in 3D between an R of 6 and 10 kpc and a Z of -2 to 2 kpc, concentrating on North-South differences. Simple parametric fits to the R, Z trends for VPHI and the velocity dispersions are presented. We confirm the recently discovered gradient in mean Galactocentric radial velocity, VR, finding that the gradient is more marked below the plane, with a Z gradient also present. The vertical velocity, VZ, also shows clear structure, with indications of a rarefaction-compression pattern, suggestive of wave-like behaviour. We perform a rigorous error analysis, tracing sources of both systematic and random errors. We confirm the North-South differences in VR and VZ along the line-of-sight, with the VR estimated independent of the proper motions. The complex three-dimensional structure of velocity space presents challenges for future modelling of the Galactic disk, with the Galactic bar, spiral arms and excitation of wave-like structures all probably playing a role.

We present the stellar atmospheric parameters (effective temperature, surface gravity, overall metallicity), radial velocities, individual abundances and distances determined for 425 561 stars, which constitute the fourth public data release of the RAdial Velocity Experiment (RAVE). The stellar atmospheric parameters are computed using a new pipeline, based on the algorithms of MATISSE and DEGAS. The spectral degeneracies and the 2MASS photometric information are now better taken into consideration, improving the parameter determination compared to the previous RAVE data releases. The individual abundances for six elements (magnesium, aluminum, silicon, titanium, iron and nickel) are also given, based on a special-purpose pipeline which is also improved compared to that available for the RAVE DR3 and Chemical DR1 data releases. Together with photometric information and proper motions, these data can be retrieved from the RAVE collaboration website and the Vizier database.

We analyse the kinematics of disc stars observed by the RAVE survey in and beyond the Solar neighbourhood.We detect significant overdensities in the velocity distributions using a technique based on the wavelet transform.We find that the main local kinematic groups are large scale features, surviving at least up to ~1 kpc from the Sun in the direction of anti-rotation, and also at ~700 pc below the Galactic plane.We also find that for regions located at different radii than the Sun, the known groups appear shifted in the velocity plane. For example, the Hercules group has a larger azimuthal velocity for regions inside the Solar circle and a lower value outside. We have also discovered a new group at (U, V) = (92,-22) km/s in the Solar neighbourhood and confirmed the significance of other previously found groups. Some of these trends detected for the first time are consistent with dynamical models of the effects of the bar and the spiral arms. More modelling is required to definitively characterise the non-axisymmetric components of our Galaxy using these groups.

We demonstrate the existence of distribution functions that can be used to represent spherical massless cored stellar systems embedded in cuspy dark matter halos with constant mildly tangential velocity anisotropy. In particular, we derive analytically the functional form of the distribution function for a Plummer stellar sphere in a Hernquist dark halo, for \beta_0 = -0.5 and for different degrees of embedding. This particular example satisfies the condition that the central logarithmic slope of the light profile \gamma_0 > 2 \beta_0. Our models have velocity dispersion profiles similar to those observed in nearby dwarf spheroidal galaxies. Hence they can be used to generate initial conditions for a variety of problems, including N-body simulations that may represent dwarf galaxies in the Local Group.

We have developed spherically symmetric dynamical models of dwarf spheroidal galaxies using Schwarzschild's orbit superposition 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 (e.g. velocity anisotropy). This method is thus less prone to biases introduced by assumptions in comparison to the more commonly used Jeans modelling, and it allows us to derive the dark matter content in a robust way. Here we present our results for the Sculptor dwarf spheroidal galaxy, after testing our methods on mock data sets. We fit both the second and fourth velocity moment profile to break the mass-anisotropy degeneracy. For an NFW dark matter halo profile, we find that the mass of Sculptor within 1 kpc is M_1kpc = 1.03 \pm 0.07 x 10^8 Msol, and that its velocity anisotropy profile is tangentially biased and nearly constant with radius. The preferred concentration (c ~ 15) is low for its dark matter mass but consistent within the scatter found in N-body cosmological simulations. When we let the value of the central logarithmic slope \alpha{} vary, we find that the best-fit model has \alpha{} = 0, although an NFW cusp or shallower is consistent at 1 \sigma{} confidence level. On the other hand, very cuspy density profiles with logarithmic central slopes \alpha{} < -1.5 are strongly disfavoured for Sculptor.

We propose a new model for the dark matter halo of the Milky Way that fits the properties of the stellar stream associated with the Sagittarius dwarf galaxy. Our dark halo is oblate with q_z = 0.9 for r < 10 kpc, and can be made to follow the Law & Majewski model at larger radii. However, we find that the dynamical perturbations induced by the Large Magellanic Cloud on the orbit of Sgr cannot be neglected when modeling its streams. When taken into account, this leads us to constrain the Galaxy's outer halo shape to have minor-to-major axis ratio (c/a)_\Phi = 0.8 and intermediate-to-major axis ratio (b/a)_\Phi = 0.9, in good agreement with cosmological expectations.

The class of tidal features around galaxies known variously as "shells" or "umbrellas" comprises debris that has arisen from high-mass-ratio mergers with low impact parameter; the nearly radial orbits of the debris give rise to a unique morphology, a universal density profile, and a tight correlation between positions and velocities of the material. As such they are accessible to analytical treatment, and can provide a relatively clean system for probing the gravitational potential of the host galaxy. In this work we present a simple analytical model that describes the density profile, phase-space distribution, and geometry of a shell, and whose parameters are directly related to physical characteristics of the interacting galaxies. The model makes three assumptions: that their orbit is radial, that the potential of the host is spherical at the shell radii, and that the satellite galaxy had a Maxwellian velocity distribution. We quantify the error introduced by the first two assumptions and show that selecting shells by their appearance on the sky is a sufficient basis to assume that these simplifications are valid. We further demonstrate that (1) given only an image of a shell, the radial gravitational force at the shell edge and the phase-space density of the satellite are jointly constrained, (2) that combining the image with measurements of either point line-of-sight velocities or integrated spectra will yield an independent estimate of the gravitational force at a shell, and (3) that an independent measurement of this force is obtained for each shell observed around a given galaxy, potentially enabling a determination of the galactic mass distribution.

We review our current understanding of the internal dynamical properties of the dwarf spheroidal galaxies surrounding the Milky Way. These are the most dark matter dominated galaxies, and as such may be considered ideal laboratories to test the current concordance cosmological model, and in particular provide constraints on the nature of the dominant form of dark matter. We discuss the latest observations of the kinematics of stars in these systems, and how these may be used to derive their mass distribution. We tour through the various dynamical techniques used, with emphasis on the complementarity and limitations, and discuss what the results imply also in the context of cosmological models. Finally we provide an outlook on exciting developments in this field.

We use orbit based dynamical models to fit the 2nd and 4th moments of the line of sight velocity distributions of the Fornax, Sculptor, Carina and Sextans dwarf spheroidal galaxies. Our goal is to compare dark matter profile models of these four systems using Bayesian evidence. We consider NFW, Einasto and several cored profiles for their dark halos and present the probability distribution functions of the model parameters. When considering each system separately, we find there is no preference for one of these specific parametric density profiles. However, the combined evidence shows that it is unlikely that all galaxies are embedded in the same type of cored profiles of the form \rho_{DM} \propto 1/(1 + r^2)^{\beta/2}, where \beta = 3, 4. For each galaxy, we also obtain an almost model independent, and therefore accurate, measurement of the logarithmic slope of the dark matter density distribution at a radius ~ r_{-3}, i.e. where the logarithmic slope of the stellar density profile is -3. This slope ranges from ~ -1.4 for Fornax to ~ -1.1 for Sextans, both at ~ 1 kpc. All our best fit models essentially have the same mass distribution over a large range in radius (from just below r_{-3} to the last measured data point). This remarkable finding likely implies much stronger constraints on the characteristics that subhalos extracted from cosmological simulations should have in order to host the dSph galaxies around the Milky Way.

Aims: We study the relations between stellar kinematics and chemical abundances of a large sample of RAVE giants in search for selection criteria needed for disentangling different Galactic stellar populations. Methods: We select a sample of 2167 giant stars with signal-to-noise per spectral measurements above 75 from the RAVE chemical catalogue and follow the analysis performed by Gratton and colleagues on 150 subdwarf stars spectroscopically observed at high-resolution. We then use a larger sample of 9131 giants (with signal-to-noise above 60) to investigate the chemo-kinematical characteristics of our stars by grouping them into nine subsamples with common eccentricity ($e$) and maximum distance achieved above the Galactic plane ($Z_max$). Results: The RAVE kinematical and chemical data proved to be reliable by reproducing the results by Gratton et al. obtained with high-resolution spectroscopic data. Our analysis, based on the $e$-$Z_max$ plane combined with additional orbital parameters and chemical information, provides an alternative way of identifying different populations of stars. In addition to extracting canonical thick- and thin-disc samples, we find a group of stars in the Galactic plane ($Z_max<1$ kpc and 0.4 $< e < $0.6), which show homogeneous kinematics but differ in their chemical properties. We interpret this as a clear sign that some of these stars have experienced the effects of heating and/or radial migration, which have modified their original orbits. The accretion origin of such stars cannot be excluded.

We combine the six high-resolution Aquarius dark matter simulations with a semi-analytic galaxy formation model to investigate the properties of the satellites of Milky Way-like galaxies. We find good correspondence with the observed luminosity function, luminosity-metallicity relation and radial distribution of the Milky Way satellites. The star formation histories of the dwarf galaxies in our model vary widely, in accordance with what is seen observationally. Ram-pressure stripping of hot gas from the satellites leaves a clear imprint of the environment on the characteristics of a dwarf galaxy. We find that the fraction of satellites dominated by old populations of stars matches observations well. However, the internal metallicity distributions of the model satellites appear to be narrower than observed. This may indicate limitations in our treatment of chemical enrichment, which is based on the instantaneous recycling approximation. Our model works best if the dark matter halo of the Milky Way has a mass of ~8 x 10^11 Msun, in agreement with the lower estimates from observations. The galaxy that resembles the Milky Way the most also has the best matching satellite luminosity function, although it does not contain an object as bright as the SMC or LMC. Compared to other semi-analytic models and abundance matching relations we find that central galaxies reside in less massive haloes, but the halo mass-stellar mass relation in our model is consistent both with hydrodynamical simulations and with recent observations.

In preparation for future, large-scale, multi-object, high-resolution spectroscopic surveys of the Galaxy, we present a series of tests of the precision in radial velocity and chemical abundances that any such project can achieve at a 4m class telescope. We briefly discuss a number of science cases that aim at studying the chemo-dynamical history of the major Galactic components (bulge, thin and thick disks, and halo) - either as a follow-up to the Gaia mission or on their own merits. Based on a large grid of synthetic spectra that cover the full range in stellar parameters of typical survey targets, we devise an optimal wavelength range and argue for a moderately high-resolution spectrograph. As a result, the kinematic precision is not limited by any of these factors, but will practically only suffer from systematic effects, easily reaching uncertainties <1 km/s. Under realistic survey conditions (namely, considering stars brighter than r=16 mag with reasonable exposure times) we prefer an ideal resolving power of R~20000 on average, for an overall wavelength range (with a common two-arm spectrograph design) of [395;456.5] nm and [587;673] nm. We show for the first time on a general basis that it is possible to measure chemical abundance ratios to better than 0.1 dex for many species (Fe, Mg, Si, Ca, Ti, Na, Al, V, Cr, Mn, Co, Ni, Y, Ba, Nd, Eu) and to an accuracy of about 0.2 dex for other species such as Zr, La, and Sr. While our feasibility study was explicitly carried out for the 4MOST facility, the results can be readily applied to and used for any other conceptual design study for high-resolution spectrographs.

We present a new analysis of the Aquarius simulations done in combination with a semi-analytic galaxy formation model. Our goal is to establish whether the subhalos present in LCDM simulations of Milky Way-like systems could host the dwarf spheroidal (dSph) satellites of our Galaxy. Our analysis shows that, contrary to what has been assumed in most previous work, the mass profiles of subhalos are generally not well fit by NFW models but that Einasto profiles are preferred. We find that for shape parameters alpha = 0.2 - 0.5 and Vmax = 10 - 30 km/s there is very good correspondence with the observational constraints obtained for the nine brightest dSph of the Milky Way. However, to explain the internal dynamics of these systems as well as the number of objects of a given circular velocity the total mass of the Milky Way should be ~ 8x10^11 Msun, a value that is in agreement with many recent determinations, and at the low mass end of the range explored by the Aquarius simulations. Our simulations show important scatter in the number of bright satellites, even when the Aquarius Milky Way-like hosts are scaled to a common mass, and we find no evidence for a missing population of massive subhalos in the Galaxy. This conclusion is also supported when we examine the dynamics of the satellites of M31.

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.

The 4MOST consortium is currently halfway through a Conceptual Design study for ESO with the aim to develop a wide-field (>3 square degree, goal >5 square degree), high-multiplex (>1500 fibres, goal 3000 fibres) spectroscopic survey facility for an ESO 4m-class telescope (VISTA). 4MOST will run permanently on the telescope to perform a 5 year public survey yielding more than 20 million spectra at resolution R~5000 ({\lambda}=390-1000 nm) and more than 2 million spectra at R~20,000 (395-456.5 nm & 587-673 nm). The 4MOST design is especially intended to complement three key all-sky, space-based observatories of prime European interest: Gaia, eROSITA and Euclid. Initial design and performance estimates for the wide-field corrector concepts are presented. We consider two fibre positioner concepts, a well-known Phi-Theta system and a new R-Theta concept with a large patrol area. The spectrographs are fixed configuration two-arm spectrographs, with dedicated spectrographs for the high- and low-resolution. A full facility simulator is being developed to guide trade-off decisions regarding the optimal field-of-view, number of fibres needed, and the relative fraction of high-to-low resolution fibres. Mock catalogues with template spectra from seven Design Reference Surveys are simulated to verify the science requirements of 4MOST. The 4MOST consortium aims to deliver the full 4MOST facility by the end of 2018 and start delivering high-level data products for both consortium and ESO community targets a year later with yearly increments.

One of the strongest predictions of the LambdaCDM cosmological model is the presence of dark satellites orbiting all types of galaxies. We focus here on the dynamical effects of such satellites on disky dwarf galaxies, and demonstrate that these encounters can be dramatic. Although mergers with M_sat > M_d are not very common, because of the lower baryonic content they occur much more frequently on the dwarf scale than for L_*-galaxies. As an example, we present a numerical simulation of a 20% (virial) mass ratio merger between a dark satellite and a disky dwarf (akin to the Fornax dwarf galaxy in luminosity) that shows that the merger remnant has a spheroidal morphology. We conclude that perturbations by dark satellites provide a plausible path for the formation of dSph systems and also could trigger starbursts in gas rich dwarf galaxies. Therefore the transition from disky to the often amorphous, irregular, or spheroidal morphologies of dwarfs could be a natural consequence of the dynamical heating of hitherto unobservable dark satellites.

The RAdial Velocity Experiment (RAVE) is a medium resolution R~7500 spectroscopic survey of the Milky Way which already obtained over half a million stellar spectra. They present a randomly selected magnitude-limited sample, so it is important to use a reliable and automated classification scheme which identifies normal single stars and discovers different types of peculiar stars. To this end we present a morphological classification of 350,000 RAVE survey stellar spectra using locally linear embedding, a dimensionality reduction method which enables representing the complex spectral morphology in a low dimensional projected space while still preserving the properties of the local neighborhoods of spectra. We find that the majority of all spectra in the database ~90-95% belong to normal single stars, but there is also a significant population of several types of peculiars. Among them the most populated groups are those of various types of spectroscopic binary and chromospherically active stars. Both of them include several thousands of spectra. Particularly the latter group offers significant further investigation opportunities since activity of stars is a known proxy of stellar ages. Applying the same classification procedure to the sample of normal single stars alone shows that the shape of the projected manifold in two dimensional space correlates with stellar temperature, surface gravity and metallicity.

The detailed abundances of 23 elements in nine bright RGB stars in the Carina dSph are presented based on high resolution spectra gathered at the VLT and Magellan telescopes. A spherical model atmospheres analysis is applied using standard methods to spectra ranging from 380 to 680 nm. The stars in this analysis range from -2.9 < [Fe/H] < -1.3, and adopting the ages determined by Lemasle et al. (2012), we are able to examine the chemical evolution of Carina's old and intermediate-aged populations. One of the main results from this work is the evidence for inhomogeneous mixing in Carina; a large dispersion in [Mg/Fe] indicates poor mixing in the old population, an offset in the [alpha/Fe] ratios between the old and intermediate-aged populations (when examined with previously published results) suggests that the second star formation event occurred in alpha-enriched gas, and one star, Car-612, seems to have formed in a pocket enhanced in SN Ia/II products. This latter star provides the first direct link between the formation of stars with enhanced SN Ia/II ratios in dwarf galaxies to those found in the outer Galactic halo (Ivans et al. 2003). Another important result is the potential evidence for SN II driven winds. We show that the very metal-poor stars in Carina have not been enhanced in AGB or SN Ia products, and therefore their very low ratios of [Sr/Ba] suggests the loss of contributions from the early SNe II. Low ratios of [Na/Fe], [Mn/Fe], and [Cr/Fe] in two of these stars support this scenario, with additional evidence from the low [Zn/Fe] upper limit for one star. It is interesting that the chemistry of the metal-poor stars in Carina is not similar to those in the Galaxy, most of the other dSphs, or the UFDs, and suggests that Carina may be at the critical mass where some chemical enrichment events are lost through SN II driven winds.

We characterize substructure in the simulated stellar halos of Cooper et al. (2010) which were formed by the disruption of satellite galaxies within the cosmological N-body simulations of galactic halos of the Aquarius Project. These stellar halos exhibit a wealth of tidal features: broad overdensities and very narrow faint streams akin to those observed around the Milky Way. The substructures are distributed anisotropically on the sky, a characteristic that should become apparent in the next generation of photometric surveys. The normalized RMS of the density of stars on the sky appears to be systematically larger for our halos compared to the value estimated for the Milky Way from main sequence turn-off stars in the Sloan Digital Sky Survey. We show that this is likely to be due in part to contamination by faint QSOs and redder main sequence stars, and might suggest that ~10% of the Milky Way halo stars have formed in-situ.

With the advent of the Gaia data, the unprecedented kinematic census of great part of the Milky Way disc will allow us to characterise the local kinematic groups and new groups in different disc neighbourhoods. First, we show here that the models predict a stellar kinematic response to the spiral arms and bar strongly dependent on disc position. For example, we find that the kinematic groups induced by the spiral arm models change significantly if one moves only ~ 0.6 kpc in galactocentric radius, but ~ 2 kpc in azimuth. There are more and stronger groups as one approaches the spiral arms. Depending on the spiral pattern speed, the kinematic imprints are more intense in nearby vicinities or far from the Sun. Secondly, we present a preliminary study of the kinematic groups observed by RAVE. This sample will allow us, for the first time, to study the dependence on Galactic position of the (thin and thick) disc moving groups. In the solar neighbourhood, we find the same kinematics groups as detected in previous surveys, but now with better statistics and over a larger spatial volume around the Sun. This indicates that these structures are indeed large scale kinematic features.

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.

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.

We present a method to constrain the potential of the non-axisymmetric components of the Galaxy using the kinematics of stars in the solar neighborhood. The basic premise is that dynamical substructures in phase-space (i.e. due to the bar and/or spiral arms) are associated with families of periodic or irregular orbits, which may be easily identified in orbital frequency space. We use the "observed" positions and velocities of stars as initial conditions for orbital integrations in a variety of gravitational potentials. We then compute their characteristic frequencies, and study the structure present in the frequency maps. We find that the distribution of dynamical substructures in velocity- and frequency-space is best preserved when the integrations are performed in the "true" gravitational potential.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We present chemical elemental abundances for $36,561$ stars observed by the RAdial Velocity Experiment (RAVE), an ambitious spectroscopic survey of our Galaxy at Galactic latitudes $|$b$|>25^{\circ}$ and with magnitudes in the range 9$<I_{DENIS}<$13. RAVE spectra cover the Ca-triplet region at 8410--8795\AA\ with resolving power R$\sim$7500. This first data release of the RAVE chemical catalogue is complementary to the third RAVE data release of radial velocities and stellar parameters, and it contains chemical abundances for the elements Mg, Al, Si, Ca, Ti, Fe and Ni, with a mean error of $\sim$0.2 dex, as judged from accuracy tests performed on synthetic and real spectra. Abundances are estimated through a dedicated processing pipeline in which the curve of growth of individual lines is obtained from a library of absorption-line equivalent widths to construct a model spectrum that is then matched to the observed spectrum via a $\chi^2$-minimization technique. We plan to extend this pipeline to include estimates for other elements, such as oxygen and sulfur, in future data releases.

We describe a correlation function statistic that quantifies the amount of spatial and kinematic substructure in the stellar halo. We test this statistic using model stellar halo realizations constructed from the Aquarius suite of six high-resolution N-body simulations in combination with the Galform semi-analytic galaxy formation model. These simulations show considerable scatter in the properties of stellar haloes. We find that our statistic can distinguish between these plausible alternatives for the global structure of the Milky Way stellar halo. We compare with observational data and show that pencil beam surveys of ~100 tracer stars (such as the Spaghetti Survey) are not sufficient to constrain the degree of structure in the Milky Way halo with this statistic. Larger area surveys with >1000 tracer stars (such as BHB stars in the Sloan Digital Sky Survey) provide much tighter constraints on comparisons between models and data. In our simulations, we find examples of haloes with spatial and kinematic substructure consistent with the available Milky Way data.

We use cosmological simulations from the Aquarius Project to study the orbital history of the Large Magellanic Cloud (LMC) and its potential association with other satellites of the Milky Way (MW). We search for dynamical analogs to the LMC and find a subhalo that matches the LMC position and velocity at either of its two most recent pericentric passages. This suggests that the LMC is not necessarily on its first approach to the MW, provided that the virial mass of the Milky Way is M_200 ~1.8e12 M_sun. The simulation results yield specific predictions for the position and velocity of systems associated with the LMC prior to infall. If on first approach, most should lie close to the LMC because the Galactic tidal field has not yet had enough time to disperse them. If on second approach, the list of potential associates increases substantially. Interestingly, our analysis rules out an LMC association for Draco and Ursa Minor, two of the dwarf spheroidals suggested by Lynden-Bell & Lynden-Bell to form part of the "Magellanic Ghostly Stream". Our results also indicate that the direction of the orbital angular momentum is a powerful test of LMC association. This test, however, requires precise proper motions, which are unavailable for most MW satellites. Of the 4 satellites with published proper motions, only the Small Magellanic Cloud is clearly associated with the LMC. Taken at face value, the proper motions of Carina, Fornax and Sculptor rule them out as potential associates, but this conclusion should be revisited when better data become available. The dearth of satellites clearly associated with the Clouds might be solved by wide-field imaging surveys that target its surroundings, a region that may prove a fertile hunting ground for faint, previously unnoticed MW satellites.

We use the high resolution cosmological N-body simulations from the Aquarius project to investigate in detail the mechanisms that determine the shape of Milky Way-type dark matter haloes. We find that, when measured at the instantaneous virial radius, the shape of individual haloes changes with time, evolving from a typically prolate configuration at early stages to a more triaxial/oblate geometry at the present day. This evolution in halo shape correlates well with the distribution of the infalling material: prolate configurations arise when haloes are fed through narrow filaments, which characterizes the early epochs of halo assembly, whereas triaxial/oblate configurations result as the accretion turns more isotropic at later times. Interestingly, at redshift z=0, clear imprints of the past history of each halo are recorded in their shapes at different radii, which also exhibit a variation from prolate in the inner regions to triaxial/oblate in the outskirts. Provided that the Aquarius haloes are fair representatives of Milky Way-like 10^12 Msun objects, we conclude that the shape of such dark matter haloes is a complex, time-dependent property, with each radial shell retaining memory of the conditions at the time of collapse.

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.

We apply the method of Burnett & Binney (2010) for the determination of stellar distances and parameters to the internal catalogue of the Radial Velocity Experiment (Steinmetz et al. 2006). Subsamples of stars that either have Hipparcos parallaxes or belong to well-studied clusters, inspire confidence in the formal errors. Distances to dwarfs cooler than ~6000 K appear to be unbiased, but those to hotter dwarfs tend to be too small by ~10% of the formal errors. Distances to giants tend to be too large by about the same amount. The median distance error in the whole sample of 216,000 stars is 28% and the error distribution is similar for both giants and dwarfs. Roughly half the stars in the RAVE survey are giants. The giant fraction is largest at low latitudes and in directions towards the Galactic Centre. Near the plane the metallicity distribution is remarkably narrow and centred on [M/H]-0.04 dex; with increasing |z| it broadens out and its median moves to [M/H] ~ -0.5. Mean age as a function of distance from the Galactic centre and distance |z| from the Galactic plane shows the anticipated increase in mean age with |z|.

We have undertaken the study of the elemental abundances and kinematic properties of a metal-poor sample of candidate thick-disk stars selected from the RAVE spectroscopic survey of bright stars to differentiate among the present scenarios of the formation of the thick disk. In this paper, we report on a sample of 214 red giant branch, 31 red clump/horizontal branch, and 74 main-sequence/sub-giant branch metal-poor stars, which serves to augment our previous sample of only giant stars. We find that the thick disk [alpha/Fe] ratios are enhanced, and have little variation (<0.1 dex), in agreement with our previous study. The augmented sample further allows, for the first time, investigation of the gradients in the metal-poor thick disk. For stars with [Fe/H] < -1.2, the thick disk shows very small gradients, <0.03 +/- 0.02 dex/kpc, in alpha-enhancement, while we find a +0.01 +/- 0.04 dex/kpc radial gradient and a -0.09 +/- 0.05 dex/kpc vertical gradient in iron abundance. In addition, we show that the peak of the distribution of orbital eccentricities for our sample agrees better with models in which the stars that comprise the thick disk were formed primarily in the Galaxy, with direct accretion of stars contributing little. Our results thus disfavor direct accretion of stars from dwarf galaxies into the thick disk as a major contributor to the thick disk population, but cannot discriminate between alternative models for the thick disk, such as those that invoke high-redshift (gas-rich) mergers, heating of a pre-existing thin stellar disk by a minor merger, or efficient radial migration of stars.

Repeated spectroscopic observations of stars in the Radial Velocity Experiment (RAVE) database are used to identify and examine single-lined binary (SB1) candidates. The RAVE latest internal database (VDR3) includes radial velocities, atmospheric and other parameters for approximately quarter million of different stars with little less than 300,000 observations. In the sample of ~20,000 stars observed more than once, 1333 stars with variable radial velocities were identified. Most of them are believed to be SB1 candidates. The fraction of SB1 candidates among stars with several observations is between 10% and 15% which is the lower limit for binarity among RAVE stars. Due to the distribution of time spans between the re-observation that is biased towards relatively short timescales (days to weeks), the periods of the identified SB1 candidates are most likely in the same range. Because of the RAVE's narrow magnitude range most of the dwarf candidates belong to the thin Galactic disk while the giants are part of the thick disk with distances extending to up to a few kpc. The comparison of the list of SB1 candidates to the VSX catalog of variable stars yielded several pulsating variables among the giant population with the radial velocity variations of up to few tens of km/s. There are 26 matches between the catalog of spectroscopic binary orbits (SB9) and the whole RAVE sample for which the given periastron time and the time of RAVE observation were close enough to yield a reliable comparison. RAVE measurements of radial velocities of known spectroscopic binaries are consistent with their published radial velocity curves.

We present the third data release of the RAdial Velocity Experiment (RAVE) which is the first milestone of the RAVE project, releasing the full pilot survey. The catalog contains 83,072 radial velocity measurements for 77,461 stars in the southern celestial hemisphere, as well as stellar parameters for 39,833 stars. This paper describes the content of the new release, the new processing pipeline, as well as an updated calibration for the metallicity based upon the observation of additional standard stars. Spectra will be made available in a future release. The data release can be accessed via the RAVE webpage: http://www.rave-survey.org.

We present the results of a dedicated search for extremely metal-poor stars in the Fornax, Sculptor and Sextans dSphs. Five stars were selected from two earlier VLT/Giraffe and HET/HRS surveys and subsequently followed up at high spectroscopic resolution with VLT/UVES. All of them turned out to have [Fe/H] <= -3 and three stars are below [Fe/H] -3.5. This constitutes the first evidence that the classical dSphs Fornax and Sextans join Sculptor in containing extremely metal-poor stars and suggests that all of the classical dSphs contain extremely metal-poor stars. One giant in Sculptor at [Fe/H]=-3.96 +- 0.10 is the most metal-poor star ever observed in an external galaxy. We carried out a detailed analysis of the chemical abundances of the alpha, iron peak, and the heavy elements, and we performed a comparison with the Milky Way halo and the ultra faint dwarf stellar populations. Carbon, barium and strontium show distinct features characterized by the early stages of galaxy formation and can constrain the origin of their nucleosynthesis.

We identify a new, nearby (0.5 < d < 10 kpc) stream in data from the RAdial Velocity Experiment (RAVE). As the majority of stars in the stream lie in the constellation of Aquarius we name it the Aquarius Stream. We identify 15 members of the stream lying between 30 < l < 75 and -70< b <-50, with heliocentric line-of-sight velocities V_los~-200 km/s. The members are outliers in the radial velocity distribution, and the overdensity is statistically significant when compared to mock samples created with both the Besan\c{c}on Galaxy model and newly-developed code Galaxia. The metallicity distribution function and isochrone fit in the log g - T_eff plane suggest the stream consists of a 10 Gyr old population with [m/H]~-1.0. We explore relations to other streams and substructures, finding the stream cannot be identified with known structures: it is a new, nearby substructure in the Galaxy's halo. Using a simple dynamical model of a dissolving satellite galaxy we account for the localization of the stream. We find that the stream is dynamically young and therefore likely the debris of a recently disrupted dwarf galaxy or globular cluster. The Aquarius stream is thus a specimen of ongoing hierarchical Galaxy formation, rare for being right in the solar suburb.

Using a sample of 213,713 stars from the Radial Velocity Experiment (RAVE) survey, limited to a distance of 2 kpc from the Sun and to |z|<1 kpc, we report the detection of a velocity gradient of disc stars in the fourth quadrant, directed radially from the Galactic centre. In the direction of the Galactic centre, we apply a simple method independent of stellar proper motions and of Galactic parameters to assess the existence of this gradient in the RAVE data. This velocity gradient corresponds to |K+C| < 3 km/s/kpc, where K and C are the Oort constants measuring the local divergence and radial shear of the velocity field, respectively. In order to illustrate the effect, assuming a zero radial velocity of the Local Standard of Rest we then reconstruct the two-dimensional Galactocentric velocity maps using two different sets of proper motions and photometric distances based either on isochrone fitting or on K-band magnitudes, and considering two sets of values for the Galactocentric radius of the Sun and local circular speed. Further observational confirmation of our finding with line-of-sight velocities of stars at low latitudes, together with further modelling, should help constrain the non-axisymmetric components of the Galactic potential, including the bar, the spiral arms and possibly the ellipticity of the dark halo.

Very metal-poor stars are of obvious importance for many problems in chemical evolution, star formation, and galaxy evolution. Finding complete samples of such stars which are also bright enough to allow high-precision individual analyses is of considerable interest. We demonstrate here that stars with iron abundances [Fe/H] < -2 dex, and down to below -4 dex, can be efficiently identified within the Radial Velocity Experiment (RAVE) survey of bright stars, without requiring additional confirmatory observations. We determine a calibration of the equivalent width of the Calcium triplet lines measured from the RAVE spectra onto true [Fe/H], using high spectral resolution data for a subset of the stars. These RAVE iron abundances are accurate enough to obviate the need for confirmatory higher-resolution spectroscopy. Our initial study has identified 631 stars with [Fe/H] <= -2, from a RAVE database containing approximately 200,000 stars. This RAVE-based sample is complete for stars with [Fe/H] < -2.5, allowing statistical sample analysis. We identify three stars with [Fe/H] <= -4. Of these, one was already known to be `ultra metal-poor', one is a known carbon-enhanced metal-poor star, but we obtain [Fe/H]= -4.0, rather than the published [Fe/H]=-3.3, and derive [C/Fe] = +0.9, and [N/Fe] = +3.2, and the third is at the limit of our S/N. RAVE observations are on-going and should prove to be a rich source of bright, easily studied, very metal-poor stars.

The RAdial Velocity Experiment (RAVE) is a spectroscopic survey of the Milky Way. We use the subsample of spectra with spectroscopically determined values of stellar parameters to determine the distances to these stars. The list currently contains 235,064 high quality spectra which show no peculiarities and belong to 210,872 different stars. The numbers will grow as the RAVE survey progresses. The public version of the catalog will be made available through the CDS services along with the ongoing RAVE public data releases. The distances are determined with a method based on the work by Breddels et al.~(2010). Here we assume that the star undergoes a standard stellar evolution and that its spectrum shows no peculiarities. The refinements include: the use of either of the three isochrone sets, a better account of the stellar ages and masses, use of more realistic errors of stellar parameter values, and application to a larger dataset. The derived distances of both dwarfs and giants match within ~21% to the astrometric distances of Hipparcos stars and to the distances of observed members of open and globular clusters. Multiple observations of a fraction of RAVE stars show that repeatability of the derived distances is even better, with half of the objects showing a distance scatter of \simlt 11%. RAVE dwarfs are ~300 pc from the Sun, and giants are at distances of 1 to 2 kpc, and up to 10 kpc. This places the RAVE dataset between the more local Geneva-Copenhagen survey and the more distant and fainter SDSS sample. As such it is ideal to address some of the fundamental questions of Galactic structure and evolution in the pre-Gaia era. Individual applications are left to separate papers, here we show that the full 6-dimensional information on position and velocity is accurate enough to discuss the vertical structure and kinematic properties of the thin and thick disks.

We use VLT/FLAMES intermediate resolution (R~6500) spectra of individual red giant branch stars in the near-infrared CaII triplet (CaT) region to investigate the wide-area metallicity properties and internal kinematics of the Sextans dwarf spheroidal galaxy (dSph). Our final sample consists of 174 probable members of Sextans with accurate line-of-sight velocities (+- 2 km/s) and CaT [Fe/H] measurements (+- 0.2 dex). We use the MgI line at 8806.8 \AA\, as an empirical discriminator for distinguishing between probable members of the dSph (giant stars) and probable Galactic contaminants (dwarf stars). Sextans shows a similar chemo-dynamical behaviour to other Milky Way dSphs, with its central regions being more metal rich than the outer parts and with the more metal-rich stars displaying colder kinematics than the more metal-poor stars. Hints of a velocity gradient are found along the projected major axis and along an axis at P.A.=191 deg, however a larger and more spatially extended sample may be necessary to pin down the amplitude and direction of this gradient. We detect a cold kinematic substructure at the centre of Sextans, consistent with being the remnant of a disrupted very metal poor stellar cluster. We derive the most extended line-of-sight velocity dispersion profile for Sextans, out to a projected radius of 1.6 deg. From Jeans modelling of the observed line-of-sight velocity dispersion profile we find that this is consistent with both a cored dark matter halo with large core radius and cuspy halo with low concentration. The mass within the last measured point is in the range 2-4 x 10^8 M_sun, giving very large mass-to-light ratios, from 460 to 920 (M/L)_(V,sun).

We study the eccentricity distribution of a thick disc sample of stars observed in the Radial Velocity Experiment (RAVE) and compare it to that expected in four simulations of thick disc formation in the literature (accretion of satellites, heating of a primordial thin disc during a merger, radial migration, and gas-rich mergers), as compiled by Sales et al. (2009). We find that the distribution of our sample is peaked at low eccentricities and falls off smoothly and rather steeply to high eccentricities. This distribution is fairly robust to changes in distances, thin disc contamination, and the particular thick disc sample used. Our results are inconsistent with what is expected for the pure accretion simulation, since we find that the dynamics of local thick disc stars implies that the majority must have formed "in situ". Of the remaining models explored, the eccentricity distribution of our stars appears to be most consistent with the gas-rich merger case.

We report on the discovery of new members of nearby young moving groups, exploiting the full power of combining the RAVE survey with several stellar age diagnostic methods and follow-up high-resolution optical spectroscopy. The results include the identification of one new and five likely members of the beta Pictoris moving group, ranging from spectral types F9 to M4 with the majority being M dwarfs, one K7 likely member of the epsilon Cha group and two stars in the Tuc-Hor association. Based on the positive identifications we foreshadow a great potential of the RAVE database in progressing toward a full census of young moving groups in the solar neighbourhood.

We use the Aquarius simulation series to study the imprint of assembly history on the structure of Galaxy-mass cold dark matter halos. Our results confirm earlier work regarding the influence of mergers on the mass density profile and the inside-out growth of halos. The inner regions that contain the visible galaxies are stable since early times and are significantly affected only by major mergers. Particles accreted diffusely or in minor mergers are found predominantly in the outskirts of halos. Our analysis reveals trends that run counter to current perceptions of hierarchical halo assembly. For example, major mergers (i.e. those with progenitor mass ratios greater than 1:10) contribute little to the total mass growth of a halo, on average less than 20 per cent for our six Aquarius halos. The bulk is contributed roughly equally by minor mergers and by "diffuse" material which is not resolved into individual objects. This is consistent with modeling based on excursion-set theory which suggests that about half of this diffuse material should not be part of a halo of any scale. Interestingly, the simulations themselves suggest that a significantly fraction is not truly diffuse, since it was ejected from earlier halos by mergers prior to their joining the main system. The Aquarius simulations resolve halos to much lower mass scales than are expected to retain gas or form stars. These results thus confirm that most of the baryons from which visible galaxies form are accreted diffusely, rather than through mergers, and they suggest that only relatively rare major mergers will affect galaxy structure at later times.

Theories of thick disk formation can be differentiated by measurements of stellar elemental abundances. We have undertaken a study of metal-poor stars selected from the RAVE spectroscopic survey of bright stars to establish whether or not there is a significant population of metal-poor thick-disk stars ([Fe/H] <~ -1.0) and to measure their elemental abundances. In this paper, we present abundances of four alpha-elements (Mg, Si, Ca, Ti) and iron for a subsample of 212 RGB and 31 RC/HB stars from this study. We find that the [alpha/Fe] ratios are enhanced implying that enrichment proceeded by purely core-collapse supernovae. This requires that star formation in each star forming region had a short duration. The relative lack of scatter in the [alpha/Fe] ratios implies good mixing in the ISM prior to star formation. In addition, the ratios resemble that of the halo, indicating that the halo and thick disk share a similar massive star IMF. We conclude that the alpha enhancement of the metal-poor thick disk implies that direct accretion of stars from dwarf galaxies similar to surviving dwarf galaxies today did not play a major role in the formation of the thick disk.

Wide-field multi-object spectroscopy is a high priority for European astronomy over the next decade. Most 8-10m telescopes have a small field of view, making 4-m class telescopes a particularly attractive option for wide-field instruments. We present a science case and design drivers for a wide-field multi-object spectrograph (MOS) with integral field units for the 4.2-m William Herschel Telescope (WHT) on La Palma. The instrument intends to take advantage of a future prime-focus corrector and atmospheric-dispersion corrector that will deliver a field of view 2 deg in diameter, with good throughput from 370 to 1,000 nm. The science programs cluster into three groups needing three different resolving powers R: (1) high-precision radial-velocities for Gaia-related Milky Way dynamics, cosmological redshift surveys, and galaxy evolution studies (R = 5,000), (2) galaxy disk velocity dispersions (R = 10,000) and (3) high-precision stellar element abundances for Milky Way archaeology (R = 20,000). The multiplex requirements of the different science cases range from a few hundred to a few thousand, and a range of fibre-positioner technologies are considered. Several options for the spectrograph are discussed, building in part on published design studies for E-ELT spectrographs. Indeed, a WHT MOS will not only efficiently deliver data for exploitation of important imaging surveys planned for the coming decade, but will also serve as a test-bed to optimize the design of MOS instruments for the future E-ELT.

For the first time we show the detailed late-stage chemical evolution history of small nearby dwarf spheroidal galaxy in the Local Group. We present the results of a high resolution (R$\sim$20000) FLAMES/GIRAFFE abundance study at ESO/VLT of 81 photometrically selected red giant branch stars in the central 25$'$ of the Fornax dwarf spheroidal galaxy. We present abundances of \alfe\ (Mg, Si, Ca and Ti), iron-peak elements (Fe, Ni and Cr) and heavy elements (Y, Ba, La, Nd and Eu). Our sample was randomly selected, and is clearly dominated by the younger and more metal rich component of Fornax which represents the major fraction of stars in the central region. This means that the majority of our stars are 1$-$4 Gyr old, and thus represent the end phase of chemical evolution in this system. Our sample of stars has unusually low [$\alpha$/Fe], [Ni/Fe] and [Na/Fe] compared to the Milky Way stellar populations at the same [Fe/H]. The particularly important role of stellar winds from low metallicity AGB stars in the creation of s-process elements is clearly seen from the high [Ba/Y]. Furthermore, we present evidence for an s-process contribution to Eu.

We devise a new method for the detection of double-lined binary stars in a sample of the Radial Velocity Experiment (RAVE) survey spectra. The method is both tested against extensive simulations based on synthetic spectra, and compared to direct visual inspection of all RAVE spectra. It is based on the properties and shape of the cross-correlation function, and is able to recover ~80% of all binaries with an orbital period of order 1 day. Systems with periods up to 1 year are still within the detection reach. We have applied the method to 25,850 spectra of the RAVE second data release and found 123 double-lined binary candidates, only eight of which are already marked as binaries in the SIMBAD database. Among the candidates, there are seven that show spectral features consistent with the RS CVn type (solar type with active chromosphere) and seven that might be of W UMa type (over-contact binaries). One star, HD 101167, seems to be a triple system composed of three nearly identical G-type dwarfs. The tested classification method could also be applicable to the data of the upcoming Gaia mission.

We perform collisionless N-body simulations to investigate the evolution of the structural and kinematical properties of simulated thick disks induced by the growth of an embedded thin disk. The thick disks used in the present study originate from cosmologically-common 5:1 encounters between initially-thin primary disk galaxies and infalling satellites. The growing thin disks are modeled as static gravitational potentials and we explore a variety of growing-disk parameters that are likely to influence the response of thick disks. We find that the final thick-disk properties depend strongly on the total mass and radial scale-length of the growing thin disk, and much less sensitively on its growth timescale and vertical scale-height as well as the initial sense of thick-disk rotation. Overall, the growth of an embedded thin disk can cause a substantial contraction in both the radial and vertical direction, resulting in a significant decrease in the scale-lengths and scale-heights of thick disks. Kinematically, a growing thin disk can induce a notable increase in the mean rotation and velocity dispersions of thick-disk stars. We conclude that the reformation of a thin disk via gas accretion may play a significant role in setting the structure and kinematics of thick disks, and thus it is an important ingredient in models of thick-disk formation.

The NIR Ca II triplet absorption lines have proven to be an important tool for quantitative spectroscopy of individual red giant branch stars in the Local Group, providing a better understanding of metallicities of stars in the Milky Way and dwarf galaxies and thereby an opportunity to constrain their chemical evolution processes. An interesting puzzle in this field is the significant lack of extremely metal-poor stars, below [Fe/H]=-3, found in classical dwarf galaxies around the Milky Way using this technique. The question arises whether these stars are really absent, or if the empirical Ca II triplet method used to study these systems is biased in the low-metallicity regime. Here we present results of synthetic spectral analysis of the Ca II triplet, that is focused on a better understanding of spectroscopic measurements of low-metallicity giant stars. Our results start to deviate strongly from the widely-used and linear empirical calibrations at [Fe/H]<-2. We provide a new calibration for Ca II triplet studies which is valid for -0.5<[Fe/H]<-4. We subsequently apply this new calibration to current data sets and suggest that the classical dwarf galaxies are not so devoid of extremely low-metallicity stars as was previously thought.

We model the formation of the Galactic stellar halo via the accretion of satellite galaxies onto a time-dependent semi-cosmological galactic potential. Our goal is to characterize the substructure left by these accretion events in a close manner to what may be possible with the {\it Gaia} mission. We have created a synthetic {\it Gaia} Solar Neighbourhood catalogue by convolving the 6D phase-space coordinates of stellar particles from our disrupted satellites with the latest estimates of the {\it Gaia} measurement errors, and included realistic background contamination due to the Galactic disc(s) and bulge. We find that, even after accounting for the expected observational errors, the resulting phase-space is full of substructure. We are able to successfully isolate roughly 50% of the different satellites contributing to the `Solar Neighbourhood' by applying the Mean-Shift clustering algorithm in energy-angular momentum space. Furthermore, a Fourier analysis of the space of orbital frequencies allows us to obtain accurate estimates of time since accretion for approximately 30% of the recovered satellites.

We present six simulations of galactic stellar haloes formed by the tidal disruption of accreted dwarf galaxies in a fully cosmological setting. Our model is based on the Aquarius project, a suite of high resolution N-body simulations of individual dark matter haloes. We tag subsets of particles in these simulations with stellar populations predicted by the Galform semi-analytic model. Our method self-consistently tracks the dynamical evolution and disruption of satellites from high redshift. The luminosity function and structural properties of surviving satellites, which agree well with observations, suggest that this technique is appropriate. We find that accreted stellar haloes are assembled between 1<z<7 from less than 5 significant progenitors. These progenitors are old, metal-rich satellites with stellar masses similar to the brightest Milky Way dwarf spheroidals. In contrast to previous stellar halo simulations, we find that several of these major contributors survive as self-bound systems to the present day. Both the number of these significant progenitors and their infall times are inherently stochastic. This results in great diversity among our stellar haloes, which amplifies small differences between the formation histories of their dark halo hosts. The masses and density/surface-brightness profiles of the stellar haloes are consistent with expectations from the Milky Way and M31. Each halo has a complex structure, consisting of well-mixed components, tidal streams, shells and other subcomponents. This structure is not adequately described by smooth models. We find one example of an accreted thick disk. Contrasts in age and metallicity between halo stars and those in surviving satellites are in broad agreement with recent observations. [Abridged]

(Abridged) Aims:We develop a method for deriving distances from spectroscopic data and obtaining full 6D phase-space coordinates for the RAVE survey's second data release. Methods: We used stellar models combined with atmospheric properties from RAVE (Teff, logg and [Fe/H]) and (J-Ks) photometry from archival sources to derive absolute magnitudes. We are able to derive the full 6D phase-space coordinates for a large sample of RAVE stars. This method is tested with artificial data, Hipparcos trigonometric parallaxes and observations of the open cluster M67. Results: When we applied our method to a set of 16 146 stars, we found that 25% (4 037) of the stars have relative (statistical) distance errors of < 35%, while 50% (8 073) and 75% (12 110) have relative (statistical) errors smaller than 45% and 50%, respectively. Our various tests show that we can reliably estimate distances for main-sequence stars, but there is an indication of potential systematic problems with giant stars. For the main-sequence star sample (defined as those with log(g) > 4), 25% (1 744) have relative distance errors < 31%, while 50% (3 488) and 75% (5 231) have relative errors smaller than 36% and 42%, respectively. Our full dataset shows the expected decrease in the metallicity of stars as a function of distance from the Galactic plane. The known kinematic substructures in the U and V velocity components of nearby dwarf stars are apparent in our dataset, confirming the accuracy of our data and the reliability of our technique. We provide independent measurements of the orientation of the UV velocity ellipsoid and of the solar motion, and they are in very good agreement with previous work. Conclusions: The distance catalogue for the RAVE second data release is available at http://www.astro.rug.nl/~rave

We combine a series of high-resolution simulations with semi-analytic galaxy formation models to follow the evolution of a system resembling the Milky Way and its satellites. The semi-analytic model is based on that developed for the Millennium Simulation, and successfully reproduces the properties of galaxies on large scales, as well as those of the Milky Way. In this model, we are able to reproduce the luminosity function of the satellites around the Milky Way by preventing cooling in haloes with Vvir < 16.7 km/s (i.e. the atomic hydrogen cooling limit) and including the impact of the reionization of the Universe. The physical properties of our model satellites (e.g. mean metallicities, ages, half-light radii and mass-to-light ratios) are in good agreement with the latest observational measurements. We do not find a strong dependence upon the particular implementation of supernova feedback, but a scheme which is more efficient in galaxies embedded in smaller haloes, i.e. shallower potential wells, gives better agreement with the properties of the ultra-faint satellites. Our model predicts that the brightest satellites are associated with the most massive subhaloes, are accreted later (z $\lta$ 1), and have extended star formation histories, with only 1 per cent of their stars made by the end of the reionization. On the other hand, the faintest satellites were accreted early, are dominated by stars with age > 10 Gyr, and a few of them formed most of their stars before the reionization was complete. Objects with luminosities comparable to those of the classical MW satellites are associated with dark matter subhaloes with a peak circular velocity $\gta$ 10 km/s, in agreement with the latest constraints.

We study the evolution of satellite debris to establish the most suitable space to identify past merger events. We confirm that the space of orbital frequencies is very promising in this respect. In frequency space individual streams can be easily identified, and their separation provides a direct measurement of the time of accretion. We are able to show for a few idealised gravitational potentials that these features are preserved also in systems that have evolved strongly in time. Furthermore, this time evolution is imprinted in the distribution of streams in frequency space. We have also tested the power of the orbital frequencies in a fully self-consistent (live) N-body simulation of the merger between a disk galaxy and a massive satellite. Even in this case streams can be easily identified and the time of accretion of the satellite can be accurately estimated.

We use N-body simulations to study the effects of tides on the kinematical structure of satellite galaxies orbiting a Milky Way-like potential. Our work is motivated by observations of dwarf spheroidal galaxies in the Local Group, for which often a distinction is possible between a cold centrally concentrated metal rich and a hot, extended metal poor population. We find that an important attenuation of the initial differences in the distribution of the two stellar components occurs for orbits with small pericentric radii (r_per < 20 kpc). This is mainly due to: i) the loss of the gravitational support provided by the dark matter component after tidal stripping takes place, which forces a re-configuration of the luminous components, and ii) tides preferentially affect the more extended stellar component, leading to a net decrease in its velocity dispersion as a response for the mass loss, which thus shrinks the kinematical gap. We apply these ideas to the Sculptor and Carina dwarf spheroidals. Differences in their orbits might help to explain, under the assumption of similar initial configurations, why in the former a clear kinematical separation between metal poor and metal rich stars is apparent, while in Carina this segregation is significantly more subtle.

We study the orbital properties of stars in four (published) simulations of thick disks formed by: i) accretion from disrupted satellites, ii) heating of a pre-existing thin disk by a minor merger, iii) radial migration and iv) gas rich mergers. We find that the distribution of orbital eccentricities are predicted to be different for each model: a prominent peak at low eccentricity is expected for the heating, migration and gas-rich merging scenarios, while the eccentricity distribution is broader and shifted towards higher values for the accretion model. These differences can be traced back to whether the bulk of the stars in each case is formed 'in-situ' or is 'accreted', and are robust to the peculiarities of each model. A simple test based on the eccentricity distribution of nearby thick disk stars may thus help elucidate the dominant formation mechanism of the Galactic thick disk.

Galaxy formation models typically assume that the size and rotation speed of galaxy disks are largely dictated by the mass, concentration, and spin of their surrounding dark matter haloes. Equally important, however, are the fraction of baryons in the halo that collect into the central galaxy, as well as the net angular momentum that they are able to retain during its assembly process. We explore the latter using a set of four large cosmological N-body/gasdynamical simulations drawn from the OWLS (OverWhelmingly Large Simulations) project. These runs differ only in their implementation of feedback from supernovae. We find that, when expressed as fractions of their virial values, galaxy mass and net angular momentum are tightly correlated. Galaxy mass fractions, m_d=M_gal/M_vir, depend strongly on feedback, but only weakly on halo mass or spin over the halo mass range explored here (M_vir>1e11 h^{-1}M_sun). The angular momentum of a galaxy, j_d=J_gal/J_vir, correlates with m_d in a manner that is insensitive to feedback and that deviates strongly from the simple j_d = m_d assumption often adopted in semi-analytic models of galaxy formation. The m_d-j_d correlation implies that, in a given halo, galaxy disk size is maximal when the central galaxy makes up a substantial fraction (~20%-30%) of all baryons within the virial radius. At z=2, such systems may host gaseous disks with radial scale lengths as large as those reported for star-forming disks by the SINS survey, even in moderately massive haloes of average spin. Extended disks at z=2 may thus signal the presence of systems where galaxy formation has been particularly efficient, rather than the existence of haloes with unusually high spin parameter.

CONTEXT: The RAVE spectroscopic survey for galactic structure and evolution obtains 8400-8800 Ang spectra at 7500 resolving power at the UK Schmidt Telescope using the 6dF multi-fiber positioner. More than 300,000 9<I<12 and |b|>25 deg southern stars have been observed to date. AIMS: This paper presents the first intrinsic examination of stellar spectra from the RAVE survey, aimed at evaluating their diagnostic potential for peculiar stars and at contributing to the general understanding of Luminous Blue Variables (LBVs). METHODS: We used the multi-epoch spectra for all seven LBVs observed, between 2005 and 2008, in the Large Magellanic Cloud (LMC) by the RAVE survey. RESULTS: We demonstrate that RAVE spectra possess significant diagnostic potential when applied to peculiar stars and, in particular, LBVs. The behaviour of the radial velocities for both emission and absorption lines, and the spectral changes between outburst and quiescence states are described and found to agree with evidence gathered at more conventional wavelengths. The wind outflow signatures and their variability are investigated, with multi-components detected in S Doradus. Photoionisation modelling of the rich emission line spectrum of R 127 shows evidence of a massive detached ionised shell that was ejected during the 1982-2000 outburst. Surface inhomogeneities in the nuclear-processed material, brought to the surface by heavy mass loss, could have been observed in S Doradus, even if alternative explanations are possible. We also detect the transition from quiescence to outburst state in R 71. Finally, our spectrum of R 84 offers one of the clearest views of its cool companion.

We have measured the amount of kinematic substructure in the Galactic halo using the final data set from the Spaghetti project, a pencil-beam high latitude sky survey. Our sample contains 101 photometrically selected and spectroscopically confirmed giants with accurate distance, radial velocity and metallicity information. We have developed a new clustering estimator: the "4distance" measure, which when applied to our data set leads to the identification of 1 group and 7 pairs of clumped stars. The group, with 6 members, can confidently be matched to tidal debris of the Sagittarius dwarf galaxy. Two pairs match the properties of known Virgo structures. Using models of the disruption of Sagittarius in Galactic potentials with different degrees of dark halo flattening, we show that this favors a spherical or prolate halo shape, as demonstrated by Newberg et al. (2007) using SDSS data. One additional pair can be linked to older Sagittarius debris. We find that 20% of the stars in the Spaghetti data set are in substructures. From comparison with random data sets we derive a very conservative lower limit of 10% to the amount of substructure in the halo. However, comparison to numerical simulations shows that our results are also consistent with a halo entirely built up from disrupted satellites, provided the dominating features are relatively broad due to early merging or relatively heavy progenitor satellites.

We use a hybrid approach that combines high-resolution simulations of the formation of a Milky Way-like halo with a semi-analytic model of galaxy formation to study the mass content of dwarf galaxies in the concordance $\Lambda$CDM cosmology. We find that the mass within 600 pc of dark matter haloes hosting luminous satellites has a median value of $\sim 3.2\times 10^7$ Msun with very little object-to-object scatter. In contrast, the present day total luminosities of the model satellites span nearly five orders of magnitude. These findings are in very good agreement with the results recently reported in the literature for the dwarf spheroidal galaxies of the Milky Way. In our model, dwarf irregular galaxies like the Small Magellanic Cloud, are predicted to have similar or slightly larger dark matter mass within 600 pc.

We study the mass, velocity dispersion, and anisotropy profiles of $\Lambda$CDM halos using a suite of N-body simulations of unprecedented numerical resolution (the {\it Aquarius Project}). Our analysis confirms a number of results claimed by earlier work, and clarifies a few issues where conflicting claims may be found in the recent literature. The spherically-averaged density profile becomes progressively shallower inwards and, at the innermost resolved radius, the logarithmic slope is $\gamma \equiv -$d$\ln\rho/$d$\ln r \simlt 1$. Asymptotic inner slopes as steep as the recently claimed $\rho \propto r^{-1.2}$ are clearly ruled out. The radial dependence of $\gamma$ is well approximated by a power-law, $\gamma \propto r^{\alpha}$ (the Einasto profile). The shape parameter, $\alpha$, varies slightly but significantly from halo to halo, implying that the mass profiles of $\Lambda$CDM halos are not strictly universal: different halos cannot, in general, be rescaled to look identical. Departures from similarity are also seen in velocity dispersion profiles and correlate with those in density profiles so as to preserve a power-law form for the spherically averaged pseudo-phase-space density, $\rho/\sigma^3\propto r^{-1.875}$. Our conclusions are reliable down to radii below 0.4% of the virial radius, providing well-defined predictions for halo structure when baryonic effects are neglected, and thus an instructive theoretical template against which the modifications induced by the baryonic components of real galaxies can be judged.

We present a large sample of fully self-consistent hydrodynamical Nbody/Tree-SPH simulations of isolated dwarf spheroidal galaxies (dSphs). It has enabled us to identify the key physical parameters and mechanisms at the origin of the observed variety in the Local Group dSph properties. The initial total mass (gas + dark matter) of these galaxies is the main driver of their evolution. Star formation (SF) occurs in series of short bursts. In massive systems, the very short intervals between the SF peaks mimic a continuous star formation rate, while less massive systems exhibit well separated SF bursts, as identified observationally. The delay between the SF events is controlled by the gas cooling time dependence on galaxy mass. The observed global scaling relations, luminosity-mass and luminosity-metallicity, are reproduced with low scatter. We take advantage of the unprecedentedly large sample size and data homogeneity of the ESO Large Programme DART, and add to it a few independent studies, to constrain the star formation history of five Milky Way dSphs, Sextans, LeoII, Carina, Sculptor and Fornax. For the first time, [Mg/Fe] vs [Fe/H] diagrams derived from high-resolution spectroscopy of hundreds of individual stars are confronted with model predictions. We find that the diversity in dSph properties may well result from intrinsic evolution. We note, however, that the presence of gas in the final state of our simulations, of the order of what is observed in dwarf irregulars, calls for removal by external processes.

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

The fact that dark matter (DM), thus far, has revealed itself only on scales of galaxies and larger, again thrusts onto astrophysics the opportunity and the responsibility to confront the age old mystery "What is the nature of matter?" By deriving basic data on the nature of DM - e.g., mass of its particle(s), present mean temperature, distribution in galaxies and other structures in the universe, and capacity for dissipational collapse - we will be uncovering the properties of the dominant species of matter in the universe and significantly extending the standard models of particle physics. Determining the mass of the DM particle to an order of magnitude would help to sort out the particle family to which it (or they) belongs. Beyond mass, there are issues of stability. The DM particle may be unstable with a measurable half-life, or it may become unstable after absorbing a certain amount of energy from collisions. In both cases it would contribute to the present hot dark matter component. Some key parameters of DM can most accurately be measured in the very nearby universe because DM dominates the mass in the outer Milky Way (MW), in other galaxies in the Local Group, and in the Local Group in its entirety. The presence and amount of DM can be quantified by study of dynamical processes observable in fine detail within these entities. Precise measurements of 3-D velocities for stars, coherent star streams, and stars in satellite stellar systems out to the edge of the Galaxy can reveal "what is the shape, orientation, density law, and lumpiness of the dark matter halo" as well as "what is the total mass of the Galaxy?"

The concordance Lambda Cold Dark Matter (Lambda-CDM) model for the formation of structure in the Universe, while remarkably successful at describing observations of structure on large scales, continues to be challenged by observations on galactic scales. Fortunately, CDM models and their various proposed alternatives make a rich variety of testable predictions that make the Local Group and our own Milky Way Galaxy key laboratories for exploring dark matter (DM) in this regime. However, some of the most definitive tests of local DM require microarcsecond astrometry of faint sources, an astrometric regime that is a unique niche of SIM Lite. This chapter explores the important and distinct contributions that can be made by SIM Lite in the exploration of galaxy dynamics and DM on galaxy scales and that have cosmological consequences. Key areas of potential SIM Lite exploration include (1) measuring the shape, orientation, density law, and lumpiness of the dark halo of the Milky Way and other nearby galaxies, (2) determining the orbits of Galactic satellites, which may be representatives of late infall from the hierarchical formation of the Milky Way, (3) ascertaining the distribution of angular momentum and orbital anisotropy of stars and globular clusters to the outer reaches of the Galactic halo, dynamical properties that hold clues to the early hierarchical formation of the Galaxy, (4) measuring the physical nature of DM by placing strong constraints on the phase space density in the cores of nearby dSph galaxies, and (5) reconstructing the dynamical history of the Local Group through the determination of orbits and masses of its constituent galaxies.

We analyse the phase-space structure of simulated thick discs that are the result of a significant merger between a disc galaxy and a satellite. Our main goal is to establish what would be the characteristic imprints of a merger origin for the Galactic thick disc. We find that the spatial distribution predicted for thick disc stars is asymmetric, seemingly in agreement with recent observations of the Milky Way thick disc. Near the Sun, the accreted stars are expected to rotate more slowly, to have broad velocity distributions, and to occupy preferentially the wings of the line-of-sight velocity distributions. The majority of the stars in our model thick discs have low eccentricity orbits (in clear reference to the pre-existing heated disc) which gives rise to a characteristic (sinusoidal) pattern for their line of sight velocities as function of galactic longitude. The z-component of the angular momentum of thick disc stars provides a clear discriminant between stars from the pre-existing disc and those from the satellite, particularly at large radii. These results are robust against the particular choices of initial conditions made in our simulations, and thus provide clean tests of the disc heating via a minor merger scenario for the formation of thick discs.

Using a sample of 248 metal-poor stars (RR Lyraes, red giants and RHB stars) which is remarkable for the accuracy of its 6-D kinematical data, we find a new component for the local halo which has an axial ratio c/a ~ 0.2, a similar flattening to the thick disk. It has a small prograde rotation but is supported by velocity anisotropy, and contains more intermediate-metallicity stars (with -1.5 < [Fe/H] < -1.0) than the rest of our sample. We suggest that this component was formed quite late, during or after the formation of the disk. It formed either from the gas that was accreted by the last major mergers experienced by the Galaxy, or by dynamical friction of massive infalling satellite(s) with the halo and possibly the stellar disk or thick disk. The remainder of the stars in our sample exhibit a clumpy distribution in energy and angular momentum, suggesting that the early, chaotic conditions under which the inner halo formed were not violent enough to erase the record of their origins. The clumpy structure suggests that a relatively small number of progenitors were responsible for building up the inner halo, in line with theoretical expectations. We find a difference in mean binding energy between the RR Lyrae variables and the red giants in our sample, suggesting that more of the RR Lyraes in the sample belong to the outer halo, and that the outer halo may be somewhat younger, as first suggested by Searle and Zinn (1978). We also find that the RR Lyrae mean rotation is more negative than the red giants, which is consistent with the recent result of Carollo et al.(2007) that the outer halo has a retrograde rotation and with the difference in kinematics seen between RR Lyraes and BHB stars by Kinman et al.(2007).

We study predictions for dark matter phase-space structure near the Sun based on high-resolution simulations of six galaxy halos taken from the Aquarius Project. The local DM density distribution is predicted to be remarkably smooth; the density at the Sun differs from the mean over a best-fit ellipsoidal equidensity contour by less than 15% at the 99.9% confidence level. The local velocity distribution is also very smooth, but it differs systematically from a (multivariate) Gaussian distribution. This is not due to the presence of individual clumps or streams, but to broad features in the velocity modulus and energy distributions that are stable both in space and time and reflect the detailed assembly history of each halo. These features have a significant impact on the signals predicted for WIMP and axion searches. For example, WIMP recoil rates can deviate by ~10% from those expected from the best-fit multivariate Gaussian models. The axion spectra in our simulations typically peak at lower frequencies than in the case of multivariate Gaussian velocity distributions. Also in this case, the spectra show significant imprints of the formation of the halo. This implies that once direct DM detection has become routine, features in the detector signal will allow us to study the dark matter assembly history of the Milky Way. A new field, "dark matter astronomy", will then emerge.

(Abridged) We use cosmological N-body simulations to study the properties of substructure halos in galaxy-sized dark matter halos. We extend prior work on the subject by considering the whole population of subhalos physically associated with the main system. These are defined as subhalos that have at some time in the past been within the virial radius of the halo's main progenitor and that have survived as self-bound entities to $z=0$. We find that this population extends beyond {\it three times} the virial radius, and contains objects on extreme orbits. We trace the origin of these unorthodox orbits to the tidal dissociation of bound groups of subhalos, which results in the ejection of some subhalos along tidal streams. Ejected subhalos are primarily low-mass systems, leading to mass-dependent biases in their spatial distribution and kinematics: the lower the subhalo mass at accretion time, the less centrally concentrated and kinematically hotter their descendant population. The bias is strongest amongst the most massive subhalos, but disappears at the low-mass end. Our findings imply that subhalos identified within the virial radius represent an incomplete census of the substructure physically related to a halo: only about {\it one half} of all associated subhalos are found today within the virial radius of a halo. These results may explain the age dependence of the clustering of low-mass halos, and has implications for (i) the interpretation of the structural parameters and assembly histories of halos neighboring massive systems; (ii) the existence of low-mass dynamical outliers in the Local Group; and (iii) the presence of evidence for evolutionary effects well outside the traditional virial boundary of a galaxy system.

The Arcturus moving group is a well-populated example of phase space substructure within the thick disk of our Galaxy. Traditionally regarded as the remains of a dissolved open cluster, it has recently been suggested to be a remnant of a satellite accreted by our Galaxy. To investigate the origin of the group we undertook a high-resolution spectroscopic abundance study of Arcturus group members and candidates. Examining abundance of Fe, Mg, Ca, Ti, Cr, Ni, Zn, Ce, Nd, Sm and Gd for 134 stars we found that the group is chemically similar to disk stars and does not exhibit a clear chemical homogeneity. Furthermore, we confirm the existence of the Arcturus group using the Nordstroem et al. (2004), Schuster et al. (2006) and RAdial Velocity Experiment (RAVE) surveys (Steinmetz et al. 2006). The origin of the group still remains unresolved: the chemical results are consistent with a dynamical origin but do not entirely rule out a merger one. Certainly, the Arcturus group provides a challenge to our understanding of the nature and origin of the Galaxy's thick disk.

We present simulations of the formation of thick disks via the accretion of two-component satellites onto a pre-existing thin disk. Our goal is to establish the detailed characteristics of the thick disks obtained in this way, as well as their dependence on the initial orbital and internal properties of the accreted objects. We find that mergers with 10-20% mass of the mass of the host lead to the formation of thick disks whose characteristics are similar, both in morphology as in kinematics, to those observed. Despite the relatively large mass ratios, the host disks are not fully destroyed by the infalling satellites: a remaining kinematically cold and thin component containing ~15-25% of the mass can be identified, which is embedded in a hotter and thicker disk. This may for example, explain the existence of a very old thin disk stars in the Milky Way. The final scale-heights of the disks depend both on the initial inclination and properties of the merger, but the fraction of satellite stellar particles at ~4 scale-heights directly measures the mass ratio between the satellite and host galaxy. Our thick disks typically show boxy isophotes at very low surface brightness levels (>6 magnitudes below their peak value). Kinematically, the velocity ellipsoids of the simulated thick disks are similar to that of the Galactic thick disk at the solar radius. The trend of sigma_Z/sigma_R with radius is found to be a very good discriminant of the initial inclination of the accreted satellite. In the Milky Way, the possible existence of a vertical gradient in the rotational velocity of the thick disk as well as the observed value of sigma_Z/sigma_R at the solar vicinity appear to favour the formation of the thick disk by a merger with either low or intermediate orbital inclination.

We have performed the largest ever particle simulation of a Milky Way-sized dark matter halo, and present the most comprehensive convergence study for an individual dark matter halo carried out thus far. We have also simulated a sample of 6 ultra-highly resolved Milky-way sized halos, allowing us to estimate the halo-to-halo scatter in substructure statistics. In our largest simulation, we resolve nearly 300,000 gravitationally bound subhalos within the virialized region of the halo. Simulations of the same object differing in mass resolution by factors up to 1800 accurately reproduce the largest subhalos with the same mass, maximum circular velocity and position, and yield good convergence for the abundance and internal properties of dark matter substructures. We detect up to four generations of subhalos within subhalos, but contrary to recent claims, we find less substructure in subhalos than in the main halo when regions of equal mean overdensity are compared. The overall substructure mass fraction is much lower in subhalos than in the main halo. Extrapolating the main halo's subhalo mass spectrum down to an Earth mass, we predict the mass fraction in substructure to be well below 3% within 100 kpc, and to be below 0.1% within the Solar Circle. The inner density profiles of subhalos show no sign of converging to a fixed asymptotic slope and are well fit by gently curving profiles of Einasto form. The mean concentrations of isolated halos are accurately described by the fitting formula of Neto et al. down to maximum circular velocities of 1.5 km/s, an extrapolation over some 5 orders of magnitude in mass. However, at equal maximum circular velocity, subhalos are more concentrated than field halos, with a characteristic density that is typically ~2.6 times larger and increases towards the halo centre.

Dark matter is the dominant form of matter in the universe, but its nature is unknown. It is plausibly an elementary particle, perhaps the lightest supersymmetric partner of known particle species. In this case, annihilation of dark matter in the halo of the Milky Way should produce gamma-rays at a level which may soon be observable. Previous work has argued that the annihilation signal will be dominated by emission from very small clumps (perhaps smaller even than the Earth) which would be most easily detected where they cluster together in the dark matter halos of dwarf satellite galaxies. Here we show, using the largest ever simulation of the formation of a galactic halo, that such small-scale structure will, in fact, have a negligible impact on dark matter detectability. Rather, the dominant and likely most easily detectable signal will be produced by diffuse dark matter in the main halo of the Milky Way. If the main halo is strongly detected, then small dark matter clumps should also be visible, but may well contain no stars, thereby confirming a key prediction of the Cold Dark Matter (CDM) model.

We present a measure of the inclination of the velocity ellipsoid at 1 kpc below the Galactic plane using a sample of red clump giants from the RAVE DR2 release. We find that the velocity ellipsoid is tilted towards the Galactic plane with an inclination of 7.3 +/-1.8 degree. We compare this value to computed inclinations for two mass models of the Milky Way. We find that our measurement is consistent with a short scale length of the stellar disc (Rd ~2 kpc) if the dark halo is oblate or with a long scale length (Rd~3 kpc) if the dark halo is prolate. Once combined with independent constraints on the flattening of the halo, our measurement suggests that the scale length is approximately halfway between these two extreme values, with a preferred range [2.5-2.7] kpc for a nearly spherical halo. Nevertheless, no model can be clearly ruled out. With the continuation of the RAVE survey, it will be possible to provide a strong constraint on the mass distribution of the Milky Way using refined measurements of the orientation of the velocity ellipsoid.

We use a series of high-resolution simulations of a `Milky-Way' halo coupled to semi-analytic methods to study the formation of our own Galaxy and of its stellar halo. The physical properties of our model Milky Way, as well as the age and metallicity distribution of stars in the different components, are in relatively good agreement with observational measurements. Assuming that the stellar halo builds up from the cores of the satellite galaxies that merged with the Milky Way over its life-time, we are able to study the physical and structural properties of this component. In agreement with previous work, we find that the largest contribution to the stellar halo should come from a few relatively massive (10^8 - 10^10 Msun) satellites, accreted at early times. Our "stellar halo" does not exhibit any clear metallicity gradient, but higher metallicity stars are more centrally concentrated than stars of lower abundance. This indicates that the probability of observing low-metallicity halo stars increases with distance from the Galactic centre. We find that the proposed "dual" nature of the Galactic stellar halo can be explained in our model as a result of a mass-metallicity relation imprinted in the building blocks of this component.

Stellar halos may hold some of the best preserved fossils of the formation history of galaxies. They are a natural product of the merging processes that probably take place during the assembly of a galaxy, and hence may well be the most ubiquitous component of galaxies, independently of their Hubble type. This review focuses on our current understanding of the spatial structure, the kinematics and chemistry of halo stars in the Milky Way. In recent years, we have experienced a change in paradigm thanks to the discovery of large amounts of substructure, especially in the outer halo. I discuss the implications of the currently available observational constraints and fold them into several possible formation scenarios. Unraveling the formation of the Galactic halo will be possible in the near future through a combination of large wide field photometric and spectroscopic surveys, and especially in the era of Gaia.

We have used spectra of hot stars from the RAVE Survey in order to investigate the visibility and properties of five diffuse interstellar bands previously reported in the literature. The RAVE spectroscopic survey for Galactic structure and kinematics records CCD spectra covering the 8400-8800 Ang wavelength region at 7500 resolving power. The spectra are obtained with the UK Schmidt at the AAO, equipped with the 6dF multi-fiber positioner. The DIB at 8620.4 Ang is by far the strongest and cleanest of all DIBs occurring within the RAVE wavelength range, with no interference by underlying absorption stellar lines in hot stars. It correlates so tightly with reddening that it turns out to be a reliable tool to measure it, following the relation E(B-V) = 2.72 (+/- 0.03) x E.W.(Ang), valid throughout the general interstellar medium of our Galaxy. The presence of a DIB at 8648 Ang is confirmed. Its intensity appears unrelated to reddening, in agreement with scanty and preliminary reports available in the literature, and its measurability is strongly compromised by severe blending with underlying stellar HeI doublet at 8649 Ang. The two weak DIBS at 8531 and 8572 Ang do not appear real and should actually be blends of underlying stellar lines. The very weak DIB at 8439 Ang cannot be resolved within the profile of the much stronger underlying hydrogen Paschen 18 stellar line.

We use N-body simulations to explore the origin and a plausible orbit for the Orphan Stream, one of the faintest substructures discovered so far in the outer halo of our Galaxy. We are able to reproduce its position, velocity and distance measurements by appealing to a single wrap of a double-component satellite galaxy. We find that the progenitor of the Orphan Stream could have been an object similar to today's Milky Way dwarfs, such as Carina, Draco, Leo II or Sculptor; and unlikely to be connected to Complex A or Ursa Major II. Our models suggest that such progenitors, if accreted on orbits with apocenters smaller than ~35 kpc, are likely to give rise to very low surface brightness streams, which may be hiding in the outer halo and remain largely undetected with current techniques. The systematic discovery of these ghostly substructures may well require wide field spectroscopic surveys of the Milky Way's outer stellar halo.

We report the discovery that substructures/subhaloes of a galaxy-size halo tend to fall in together in groups in cosmological simulations, something that may explain the oddity of the MW satellite distribution. The original clustering at the time of infall is still discernible in the angular momenta of the subhaloes even for events which took place up to eight Gyrs ago, $z \sim 1$. This phenomenon appears to be rather common since at least 1/3 of the present-day subhaloes have fallen in groups in our simulations. Hence, this may well explain the Lynden-Bell & Lynden-Bell ghostly streams. We have also found that the probability of building up a flattened distribution similar to the MW satellites is as high as $\sim 80%$ if the MW satellites were from only one group and $\sim 20%$ when five groups are involved. Therefore, we conclude that the `peculiar' distribution of satellites around the MW can be expected with the CDM structure formation theory. This non-random assignment of satellites to subhaloes implies an environmental dependence on whether these low-mass objects are able to form stars, possibly related to the nature of reionization in the early Universe.

We present a new determination of the mass content of the Sculptor dwarf spheroidal galaxy, based on a novel approach which takes into account the two distinct stellar populations present in this galaxy. This method helps to partially break the well-known mass-anisotropy degeneracy present in the modelling of pressure-supported stellar systems.

We use a series of high-resolution N-body simulations of a `Milky-Way' halo, coupled to semi-analytic techniques, to study the formation of our own Galaxy and of its stellar halo. Our model Milky Way galaxy is a relatively young system whose physical properties are in quite good agreement with observational determinations. In our model, the stellar halo is mainly formed from a few massive satellites accreted early on during the galaxy's lifetime. The stars in the halo do not exhibit any metallicity gradient, but higher metallicity stars are more centrally concentrated than stars with lower abundances. This is due to the fact that the most massive satellites contributing to the stellar halo are also more metal rich, and dynamical friction drags them closer to the inner regions of the host halo.

We present the second data release of the Radial Velocity Experiment (RAVE), an ambitious spectroscopic survey to measure radial velocities (RVs) and stellar atmosphere parameters of up to one million stars using the 6dF multi-object spectrograph on the 1.2-m UK Schmidt Telescope of the Anglo-Australian Observatory (AAO). It is obtaining medium resolution spectra (median R=7,500) in the Ca-triplet region (8,410--8,795 \AA) for southern hemisphere stars in the magnitude range 9<I<12. Following the first data release (Steinmetz et al. 2006) the current release doubles the sample of published RVs, now containing 51,829 RVs for 49,327 individual stars observed on 141 nights between April 11 2003 and March 31 2005. Comparison with external data sets shows that the new data collected since April 3 2004 show a standard deviation of 1.3 km/s, about twice better than for the first data release. For the first time this data release contains values of stellar parameters from 22,407 spectra of 21,121 individual stars. They were derived by a penalized \chi^2 method using an extensive grid of synthetic spectra calculated from the latest version of Kurucz models. From comparison with external data sets, our conservative estimates of errors of the stellar parameters (for a spectrum with S/N=40) are 400 K in temperature, 0.5 dex in gravity, and 0.2 dex in metallicity. We note however that the internal errors estimated from repeat RAVE observations of 822 stars are at least a factor 2 smaller. We demonstrate that the results show no systematic offsets if compared to values derived from photometry or complementary spectroscopic analyses. The data release includes proper motion and photometric measurements. It can be accessed via the RAVE webpage: http://www.rave-survey.org and through CDS.

We present VLT/FLAMES spectroscopic observations (R~6500) in the CaII triplet region for 470 probable kinematic members of the Sculptor (Scl) dwarf spheroidal galaxy. The accurate velocities (+- 2 km/s) and large area coverage of Scl allow us to measure a velocity gradient of 7.6_{-2.2}^{+3.0} km/s deg^{-1} along the Scl projected major axis, likely a signature of intrinsic rotation. We also use our kinematic data to measure the mass distribution within this system. By considering independently the kinematics of the distinct stellar components known to be present in Scl, we are able to relieve known degeneracies, and find that the observed velocity dispersion profiles are best fitted by a cored dark matter halo with core radius r_c= 0.5 kpc and mass enclosed within the last measured point M(< 1.8 kpc)=3.4 +- 0.7 x 10^8 M_sun, assuming an increasingly radially anisotropic velocity ellipsoid. This results in a mass-to-light ratio of 158+-33 (M/L)_sun inside 1.8 kpc. An NFW profile with concentration c=20 and mass M(< 1.8 kpc) = 2.2_{-0.7}^{+1.0} x 10^8 M_sun is also statistically consistent with the observations, but it tends to yield poorer fits for the metal rich stars.

We analyse the dynamical properties of substructures in a high-resolution dark matter simulation of the formation of a Milky Way-like halo in a $\Lambda$CDM cosmology. Our goal is to shed light on the dynamical peculiarities of the Milky Way satellites. Our simulations show that about 1/3 of the subhalos have been accreted in groups. We quantify this clustering by measuring the alignment of the angular momentum of subhalos in a group. We find that this signal is visible even for objects accreted up to $z \sim 1$, i.e. 8 Gyr ago, and long after the spatial coherence of the groups has been lost due the host tidal field. This group infall may well explain the ghostly streams proposed by Lynden-Bell & Lynden-Bell to orbit the Milky Way. Our analyses also show that if most satellites originate in a few groups, the disk-like distribution of the Milky Way satellites would be almost inevitable. This non-random assignment of satellites to subhalos implies an environmental dependence on whether these low-mass objects are able to form stars, possibly related to the nature of reionization in the early Universe. With this picture, both the ``ghostly streams'' and the ``disk-like configuration'' are manifestations of the same phenomenon: the hierarchical growth of structure down to the smallest scales.

We analyze the distribution of G and K type stars towards the Galactic poles using RAVE and ELODIE radial velocities, 2MASS photometric star counts, and UCAC2 proper motions. The combination of photometric and 3D kinematic data allows us to disentangle and describe the vertical distribution of dwarfs, sub-giants and giants and their kinematics. We identify discontinuities within the kinematics and magnitude counts that separate the thin disk, thick disk and a hotter component. The respective scale heights of the thin disk and thick disk are 225$\pm$10 pc and 1048$\pm$36 pc. We also constrain the luminosity function and the kinematic distribution function. The existence of a kinematic gap between the thin and thick disks is incompatible with the thick disk having formed from the thin disk by a continuous process, such as scattering of stars by spiral arms or molecular clouds. Other mechanisms of formation of the thick disk such as `created on the spot' or smoothly `accreted' remain compatible with our findings.

We present a new and completely general technique for calculating the fine-grained phase-space structure of dark matter throughout the Galactic halo. Our goal is to understand this structure on the scales relevant for direct and indirect detection experiments. Our method is based on evaluating the geodesic deviation equation along the trajectories of individual DM particles. It requires no assumptions about the symmetry or stationarity of the halo formation process. In this paper we study general static potentials which exhibit more complex behaviour than the separable potentials studied previously. For ellipsoidal logarithmic potentials with a core, phase mixing is sensitive to the resonance structure, as indicated by the number of independent orbital frequencies. Regions of chaotic mixing can be identified by the very rapid decrease in the real space density of the associated dark matter streams. We also study the evolution of stream density in ellipsoidal NFW halos with radially varying isopotential shape, showing that if such a model is applied to the Galactic halo, at least $10^5$ streams are expected near the Sun. The most novel aspect of our approach is that general non-static systems can be studied through implementation in a cosmological N-body code. Such an implementation allows a robust and accurate evaluation of the enhancements in annihilation radiation due to fine-scale structure such as caustics. We embed the scheme in the current state-of-the-art code GADGET-3 and present tests which demonstrate that N-body discreteness effects can be kept under control in realistic configurations.

We demonstrate that low resolution CaII triplet (CaT) spectroscopic estimates of the overall metallicity ([Fe/H]) of individual Red Giant Branch (RGB) stars in two nearby dwarf spheroidal galaxies (dSphs) agree to +-0.1-0.2 dex with detailed high resolution spectroscopic determinations for the same stars over the range -2.5 < [Fe/H] < -0.5. For this study we used a sample of 129 stars observed in low and high resolution mode with VLT/FLAMES in the Sculptor and Fornax dSphs. We also present the data reduction steps we used in our low resolution analysis and show that the typical accuracy of our velocity and CaT [Fe/H] measurement is ~2 km/s and 0.1 dex respectively. We conclude that CaT-[Fe/H] relations calibrated on globular clusters can be applied with confidence to RGB stars in composite stellar populations over the range -2.5 < [Fe/H] < -0.5.

We study the phase-space behaviour of nearby trajectories in integrable potentials. We show that the separation of nearby orbits initially diverges very fast, mimicking a nearly exponential behaviour, while at late times it grows linearly. This initial exponential phase, known as Miller's instability, is commonly found in N-body simulations, and has been attributed to short-term (microscopic) N-body chaos. However we show here analytically that the initial divergence is simply due to the shape of an orbit in phase-space. This result confirms previous suspicions that this transient phenomenon is not related to an instability in the sense of non-integrable behaviour in the dynamics of N-body systems.

We have searched for in-falling stellar streams on to the local Milky Way disc in the CORAVEL and RAVE surveys. The CORAVEL survey consists of local dwarf stars (Nordstrom et al. Geneva-Copenhagen survey) and local Famaey et al. giant stars. We select RAVE stars with radial velocities that are sensitive to the Galactic vertical space velocity (Galactic latitude b < -45 deg). Kuiper statistics have been employed to test the symmetry of the Galactic vertical velocity distribution functions in these samples for evidence of a net vertical flow that could be associated with a (tidal?) stream of stars with vertically coherent kinematics. In contrast to the `Field of Streams' found in the outer halo, we find that the local volumes of the solar neighbourhood sampled by the CORAVEL dwarfs (complete within ~3 x 10^-4 kpc^3), CORAVEL giants (complete within ~5 x 10^-2 kpc^3) and RAVE (5-15% complete within ~8 kpc^3) are devoid of any vertically coherent streams containing hundreds of stars. This is sufficiently sensitive to allow our RAVE sample to rule out the passing of the tidal stream of the disrupting Sagittarius (Sgr) dwarf galaxy through the solar neighbourhood. This agrees with the most recent determination of its orbit and dissociates it from the Helmi et al. halo stream. Our constraints on the absence of the Sgr stream near the Sun could prove a useful tool for discriminating between Galactic potential models. The lack of a net vertical flow through the solar neighbourhood in the CORAVEL giants and RAVE samples argues against the Virgo overdensity crossing the disc near the Sun. There are no vertical streams in the CORAVEL giants and RAVE samples with stellar densities >1.6 x 10^4 and 1.5 x 10^3 stars kpc^-3 respectively and therefore no evidence for locally enhanced dark matter.

We have assembled a sample of halo stars in the solar neighborhood to look for halo substructure in velocity and angular momentum space. Our sample includes red giants, RR Lyrae, and red horizontal branch stars within 2.5 kpc of the Sun with [Fe/H] less than -1.0. It was chosen to include stars with accurate distances, space velocities, and metallicities as well as well-quantified errors. We confirm the existence of the streams found by Helmi and coworkers, which we refer to as the H99 streams. These streams have a double-peaked velocity distribution in the z direction. We use the results of modeling of the H99 streams by Helmi and collaborators to test how one might use v_z velocity information and radial velocity information to detect kinematic substructure in the halo. We find that detecting the H99 streams with radial velocities alone would require a large sample. We use the velocity distribution of the H99 streams to estimate their age. From our model of the progenitor of the H99 streams, we determine that it was accreted between 6 and 9 Gyr ago. The H99 streams have [alpha/Fe] abundances similar to other halo stars in the solar neighborhood, suggesting that the gas that formed these stars were enriched mostly by Type II SNe. We have also discovered in angular momentum space two other possible substructures, which we refer to as the retrograde and prograde outliers. The retrograde outliers are likely to be halo substructure, but the prograde outliers are most likely part of the smooth halo. The retrograde outliers have significant structure in the v_phi direction and show a range of [alpha/Fe]. The methods presented in this paper can be used to exploit the kinematic information present in future large databases like RAVE, SDSSII/SEGUE, and Gaia.

We report new constraints on the local escape speed of our Galaxy. Our analysis is based on a sample of high velocity stars from the RAVE survey and two previously published datasets. We use cosmological simulations of disk galaxy formation to motivate our assumptions on the shape of the velocity distribution, allowing for a significantly more precise measurement of the escape velocity compared to previous studies. We find that the escape velocity lies within the range $498\kms < \ve < 608 \kms$ (90 per cent confidence), with a median likelihood of $544\kms$. The fact that $\ve^2$ is significantly greater than $2\vc^2$ (where $\vc=220\kms$ is the local circular velocity) implies that there must be a significant amount of mass exterior to the Solar circle, i.e. this convincingly demonstrates the presence of a dark halo in the Galaxy. For a simple isothermal halo, one can calculate that the minimum radial extent is $\sim58$ kpc. We use our constraints on $\ve$ to determine the mass of the Milky Way halo for three halo profiles. For example, an adiabatically contracted NFW halo model results in a virial mass of $1.42^{+1.14}_{-0.54}\times10^{12}M_\odot$ and virial radius of $305^{+66}_{-45}$ kpc (90 per cent confidence). For this model the circular velocity at the virial radius is $142^{+31}_{-21}\kms$. Although our halo masses are model dependent, we find that they are in good agreement with each other.

Recent observations show that the number of stars with very low metallicities in the dwarf spheroidal satellites of the Milky Way is low, despite the low average metallicities of stars in these systems. We undertake numerical simulations of star formation and metal enrichment of dwarf galaxies in order to verify whether this result can be reproduced with "standard" assumptions. The answer is likely to be negative, unless some selection bias against very low metallicity stars is present in the observations.

As part of the Dwarf galaxies Abundances and Radial-velocities Team (DART) program, we have measured the metallicities of a large sample of stars in four nearby dwarf spheroidal galaxies (dSph's): Sculptor, Sextans, Fornax, and Carina. The low mean metal abundances and the presence of very old stellar populations in these galaxies have supported the view that they are fossils from the early universe. However, contrary to naive expectations, we find a significant lack of stars with metallicities below [Fe/H] ~ -3 dex in all four systems. This suggests that the gas that made up the stars in these systems had been uniformly enriched prior to their formation. Furthermore, the metal-poor tail of the dSph metallicity distribution is significantly different from that of the Galactic halo. These findings show that the progenitors of nearby dSph's appear to have been fundamentally different from the building blocks of the Milky Way, even at the earliest epochs.

As part of the DART project we have used the ESO/2.2m Wide Field Imager in conjunction with the VLT/FLAMES GIRAFFE spectrograph to study the detailed properties of the resolved stellar population of the Fornax dwarf spheroidal galaxy out to and beyond its tidal radius. We re-derived the structural parameters of the Fornax dwarf spheroidal using our wide field imaging covering the galaxy out to its tidal radius, and analysed the spatial distribution of the Fornax stars of different ages as selected from Colour-Magnitude Diagram analysis. We have obtained accurate velocities and metallicities from spectra in the CaII triplet wavelength region for 562 Red Giant Branch stars which have velocities consistent with membership in Fornax dwarf spheroidal. We have found evidence for the presence of at least three distinct stellar components: a young population (few 100 Myr old) concentrated in the centre of the galaxy, visible as a Main Sequence in the Colour-Magnitude Diagram; an intermediate age population (2-8 Gyr old); and an ancient population (> 10Gyr),which are distinguishable from each other kinematically, from the metallicity distribution and in the spatial distribution of stars found in the Colour-Magnitude Diagram. From our spectroscopic analysis we find that the ``metal rich'' stars ([Fe/H]> -1.3) show a less extended and more concentrated spatial distribution, and display a colder kinematics than the ``metal poor'' stars ([Fe/H<-1.3). There is tentative evidence that the ancient stellar population in the centre of Fornax does not exhibit equilibrium kinematics. This could be a sign of a relatively recent accretion of external material.

We present the first data release of the Radial Velocity Experiment (RAVE), an ambitious spectroscopic survey to measure radial velocities and stellar atmosphere parameters (temperature, metallicity, surface gravity) of up to one million stars using the 6dF multi-object spectrograph on the 1.2-m UK Schmidt Telescope of the Anglo-Australian Observatory (AAO). The RAVE program started in 2003, obtaining medium resolution spectra (median R=7,500) in the Ca-triplet region ($\lambda\lambda$ 8,410--8,795 \AA) for southern hemisphere stars drawn from the Tycho-2 and SuperCOSMOS catalogs, in the magnitude range 9<I<12. The first data release is described in this paper and contains radial velocities for 24,748 individual stars (25,274 measurements when including re-observations). Those data were obtained on 67 nights between 11 April 2003 to 03 April 2004. The total sky coverage within this data release is $\sim$4,760 square degrees. The average signal to noise ratio of the observed spectra is 29.5, and 80% of the radial velocities have uncertainties better than 3.4 km/s. Combining internal errors and zero-point errors, the mode is found to be 2 km/s. Repeat observations are used to assess the stability of our radial velocity solution, resulting in a variance of 2.8 km/s. We demonstrate that the radial velocities derived for the first data set do not show any systematic trend with color or signal to noise. The RAVE radial velocities are complemented in the data release with proper motions from Starnet 2.0, Tycho-2 and SuperCOSMOS, in addition to photometric data from the major optical and infrared catalogs (Tycho-2, USNO-B, DENIS and 2MASS). The data release can be accessed via the RAVE webpage: http://www.rave-survey.org.

We search for signatures of past accretion events in the Milky Way in the recently published catalogue by Nordstrom et al.(2004), containing accurate spatial and kinematic information as well as metallicities for 13240 nearby stars. To optimize our strategy, we use numerical simulations to characterize the properties of the debris from disrupted satellites. We find that stars with a common progenitor should show distinct correlations between their orbital parameters; in particular, between the apocentre A and pericentre P, as well as their z-angular momentum (L_z). In the APL-space, such stars are expected to cluster around regions of roughly constant eccentricity. The APL space for the Nordstrom catalogue exhibits a wealth of substructure, much of which can be linked to dynamical perturbations induced by spiral arms and the Galactic bar. However, our analysis also reveals a statistically significant excess of stars on orbits of common (moderate) eccentricity, analogous to the pattern expected for merger debris. Besides being dynamically peculiar, the 274 stars in these substructures have very distinct metallicity and age distributions, providing further evidence of their extra-Galactic provenance. It is possible to identify among these stars, three coherent Groups with characteristic ages and metallicities, that, in all likelihood, correspond to the remains of disrupted satellites. [abridged]

We have compiled a new sample of 240 halo objects with accurate distance and radial velocity measurements, including globular clusters, satellite galaxies, field blue horizontal branch stars and red giant stars from the Spaghetti survey. The new data lead to a significant increase in the number of known objects for Galactocentric radii beyond 50 kpc, which allows a reliable determination of the radial velocity dispersion profile out to very large distances. The radial velocity dispersion shows an almost constant value of 120 km/s out to 30 kpc and then continuously declines down to 50 km/s at about 120 kpc. This fall-off puts important constraints on the density profile and total mass of the dark matter halo of the Milky Way. For a constant velocity anisotropy, the isothermal profile is ruled out, while both a dark halo following a truncated flat model of mass $1.2^{+1.8}_{-0.5}\times 10^{12}$ M_sun and an NFW profile of mass $0.8^{+1.2}_{-0.2}\times 10^{12}$ M_sun and c=18 are consistent with the data. The significant increase in the number of tracers combined with the large extent of the region probed by these has allowed a more precise determination of the Milky Way mass in comparison to previous works. We also show how different assumptions for the velocity anisotropy affect the performance of the mass models.

The Gaia mission is designed as a Galaxy explorer, and will measure simultaneously, in a survey mode, the five or six phase space parameters of all stars brighter than 20th magnitude, as well as providing a description of their astrophysical characteristics. These measurements are obtained by combining an astrometric instrument with micro-arcsecond capabilities, a photometric system giving the magnitudes and colours in 15 bands and a medium resolution spectrograph named the Radial Velocity Spectrometer (RVS). The latter instrument will produce spectra in the 848 to 874 nm wavelength range, with a resolving power R = 11 500, from which radial velocities, rotational velocities, atmospheric parameters and abundances can be derived. A companion paper (Katz et al. 2004) presents the characteristics of the RVS and its performance. This paper details the outstanding scientific impact of this important part of the Gaia satellite on some key open questions in present day astrophysics. The unbiased and simultaneous acquisition of multi-epoch radial velocities and individual abundances of key elements in parallel with the astrometric parameters is essential for the determination of the dynamical state and formation history of our Galaxy. Moreover, for stars brighter than V=15, the resolving power of the RVS will give information about most of the effects which influence the position of a star in the Hertzsprung-Russell diagram, placing unprecedented constraints on the age, internal structure and evolution of stars of all types. Finally, the RVS multi-epoch observations are ideally suited to the identification, classification and characterisation of the many types of double, multiple and variable stars.

We have analyzed the catalog of 2106 non-kinematically selected metal poor stars in the solar neighborhood published by Beers et al. (2000), with the goal of quantifying the amount of substructure in the motions of the fastest halo stars. We have computed the two-point velocity correlation function for a subsample of halo stars within 1-2 kpc of the Sun, and found statistical evidence of substructure, with a similar amplitude to that predicted by high resolution CDM simulations. The signal is due to a small kinematic group whose dynamical properties are compared to the stellar "stream", previously discovered by Helmi et al. (1999). If real, this high velocity moving group would provide further support to the idea that substructures remain as fossils from the formation of the Galaxy as expected in the CDM scenario.

We have found evidence for the presence of two distinct ancient stellar components (both geq 10 Gyr old) in the Sculptor dwarf spheroidal galaxy. We used the ESO Wide Field Imager (WFI) in conjunction with the VLT/FLAMES spectrograph to study the properties of the resolved stellar population of Sculptor out to and beyond the tidal radius. We find that two components are discernible in the spatial distribution of Horizontal Branch stars in our imaging, and in the [Fe/H] and v_hel distributions for our large sample of spectroscopic measurements. They can be generally described as a ``metal-poor'' component ([Fe/H]< -1.7) and a ``metal-rich'' component ([Fe/H]> -1.7). The metal-poor stars are more spatially extended than the metal-rich stars, and they also appear to be kinematically distinct. These results provide an important insight into the formation processes of small systems in the early universe and the conditions found there. Even this simplest of galaxies appears to have had a surprisingly complex early evolution.

The definition and optimisation studies for the Gaia satellite spectrograph, the Radial Velocity Spectrometer (RVS), converged in late 2002 with the adoption of the instrument baseline. This paper reviews the characteristics of the selected configuration and presents its expected performance. The RVS is a 2.0 by 1.6 degree integral field spectrograph, dispersing the light of all sources entering its field of view with a resolving power R=11 500 over the wavelength range [848, 874] nm. The RVS will continuously and repeatedly scan the sky during the 5 years of the Gaia mission. On average, each source will be observed 102 times over this period. The RVS will collect the spectra of about 100-150 million stars up to magnitude V~17-18. At the end of the mission, the RVS will provide radial velocities with precisions of ~2 km/s at V=15 and \~15-20 km/s at V=17, for a solar metallicity G5 dwarf. The RVS will also provide rotational velocities, with precisions (at the end of the mission) for late type stars of sigma_vsini ~5 km/s at V~15 as well as atmospheric parameters up to V~14-15. The individual abundances of elements such as Silicon and Magnesium, vital for the understanding of Galactic evolution, will be obtained up to V~12-13. Finally, the presence of the 862.0 nm Diffuse Interstellar Band (DIB) in the RVS wavelength range will make it possible to derive the three dimensional structure of the interstellar reddening.

Recently, radial velocities have been measured for a large sample of M giants from the 2MASS catalog, selected to be part of the Sgr dwarf leading and trailing streams. Here we present a comparison of their kinematics to models of the Sgr dwarf debris orbiting Galactic potentials, with halo components of varying degrees of flattening and elongation. This comparison shows that the portion of the trailing stream mapped so far is dynamically young and hence does not provide very stringent constraints on the shape of the Galactic dark-matter halo. The leading stream, however, contains slightly older debris, and its kinematics provide for the first time direct evidence that the dark-matter halo of the Galaxy has a prolate shape with an average density axis ratio within the orbit of Sgr close to 5/3.

It has been recently claimed that the confined structure of the debris from the Sagittarius dwarf implies that the dark-matter halo of our Galaxy should be nearly spherical, in strong contrast with predictions from cold dark-matter (CDM) simulations, where dark halos are found to have typical density axis ratios of 0.6 - 0.8. In this letter, numerical simulations are used to show that the Sagittarius streams discovered thus far are too young dynamically to be sensitive to the shape of the dark halo of the Milky Way. The data presently available are entirely consistent with a Galactic dark-matter halo that could either be oblate or prolate, with minor-to-major density axis ratios c/a as low as 0.6 within the region probed by the orbit of the Sagittarius dwarf.

We reanalyze the group of stars associated kinematically with Arcturus (Eggen 1971), and confirm that they constitute a peculiar grouping of metal-poor stars with similar apocentric radius, common angular momentum, and distinct metal abundance patterns. These properties are consistent with those expected for a group of stars originating from the debris of a disrupted satellite. The Arcturus group stands out clearly in the catalog of Beers et al. (2000) as well as in other compilations of metal-poor stars in the solar neighborhood. Its angular momentum appears too low to arise from dynamical perturbations induced by the Galactic bar, and coincides with that of the kinematically-peculiar population of stars identified above and below the Galactic plane by Gilmore, Wyse & Norris (2002): the Arcturus group is in all likelihood the solar neighborhood extension of such population. Further analysis is needed to confirm our proposal that Arcturus, together with some of the brightest stars in the night sky, may have formed beyond the confines of the Galaxy. However, validating the extragalactic origin of the Arcturus group would lend support to the prevalence of accretion events envisioned in hierarchical models of galaxy formation, and would suggest that early mergers may have contributed a substantial fraction of the old and metal-deficient stars in the Galactic disk.

We examine the tidal disruption of satellite galaxies in a cosmological simulation of the formation of a disk galaxy in the $\Lambda$CDM scenario. We find that the disruption of satellite galaxies in orbits roughly coplanar with the disk leads naturally to the formation of ring-like stellar structures similar to that recently discovered in the outer disk of the Milky Way. Two interpretations appear plausible in this context. One is that the ring is a transitory, localized radial density enhancement reflecting the apocenter of particles stripped from a satellite during a recent pericentric passage (a ``tidal arc'' reminiscent of the tidal arms seen in disk galaxy mergers). In the second scenario, the ring is analogous to the ``shells'' found around some elliptical galaxies, and would result from a minor merger that took place several Gyr ago. The two interpretations differ in several ways. Tidal arcs are expected to span a limited longitude range; may carry a substantial fraction of the original mass of the satellite; should exhibit a significant velocity gradient with Galactic longitude; and are in all likelihood asymmetric above and below the plane of the disk. Shells, on the other hand, may comprise at most a small fraction of the original mass of the satellite and, due to their more relaxed state, ought to be symmetric above and below the plane, with no discernible velocity gradients across the structure. If confirmed as a tidal feature, the ring discovered by SDSS in the outer Galactic disk would strengthen the view--supported by numerical simulations--that minor mergers have played a significant role building up not only the stellar halo, but also the disk components of the Galaxy.

We present a method for photometric selection of metal-poor halo giants from the imaging data of the Sloan Digital Sky Survey (SDSS). These stars are offset from the stellar locus in the (g-r) vs. (u-g) color-color diagram. Based on a sample of 29 candidates for which spectra were taken, we derive a selection efficiency of the order of 50%, for stars brighter than $r \sim 17^m$. The candidates selected in 400 deg$^2$ of sky from the SDSS Early Data Release trace the known halo structures (tidal streams from the Sagittarius dwarf galaxy, the Draco dwarf spheroidal galaxy), indicating that such a color-selected sample can be used to study the halo structure even without spectroscopic information. This method, and supplemental techniques for selecting halo stars, such as RR Lyrae stars and other blue horizontal branch stars, can produce an unprecedented three-dimensional map of the Galactic halo based on the SDSS imaging survey.

We study the formation of the Milky Way's halo in a $\Lambda$CDM cosmology by scaling down a high resolution simulation of the formation of a cluster of galaxies. We determine how much phase-space substructure is left over from the objects that merge to build up the present galaxy. We study the debris streams originating from such objects and find that their evolution can be explained simply in terms of the conservation of phase-space density. Analysing the mass growth history of our halo we find that its inner regions have been in place for more than 10 Gyr, but that the growth of the halo as a whole is more gradual, in agreement with other high resolution simulations of dark-matter halos. Recent accretion contributes to the inner 10 kpc of the halo only at the 10$^{-4}$ level. Finally we determine the number of dark-matter streams as a function of distance from the centre of the halo. In the equivalent of the ``Solar vicinity'', we find that the dark-matter is smoothly distributed in space, and that the velocity ellipsoid is formed by hundreds of thousands of streams, most of which have velocity dispersions of the order of 1 km/s or less.

We study the phase-space structure of a dark-matter halo formed in a high resolution simulation of a Lambda CDM cosmology. Our goal is to quantify how much substructure is left over from the inhomogeneous growth of the halo, and how it may affect the signal in experiments aimed at detecting the dark matter particles directly. If we focus on the equivalent of ``Solar vicinity'', we find that the dark-matter is smoothly distributed in space. The probability of detecting particles bound within dense lumps of individual mass less than 10^7 M_\sun h^{-1} is small, less than 10^{-2}. The velocity ellipsoid in the Solar neighbourhood deviates only slightly from a multivariate Gaussian, and can be thought of as a superposition of thousands of kinematically cold streams. The motions of the most energetic particles are, however, strongly clumped and highly anisotropic. We conclude that experiments may safely assume a smooth multivariate Gaussian distribution to represent the kinematics of dark-matter particles in the Solar neighbourhood. Experiments sensitive to the direction of motion of the incident particles could exploit the expected anisotropy to learn about the recent merging history of our Galaxy.

The last 10-20 years has seen a profound shift in views of how the Galaxy's halo formed. The idea of a monolithic early collapse of a single system (Eggen, Lynden-Bell and Sandage 1962) has been challenged by observations at high redshift and by cosmological models of structure formation. These findings imply that we should see clear evidence of hierarchical formation processes in nearby galaxies. Recent studies of our Galaxy, made possible by large-scale CCD surveys such as the Sloan Digital Sky Survey (SDSS), have begun to reveal tantalizing evidence of substructure in the outer halo. We review evidence for tidal streams associated with known Milky Way satellites and for star streams whose progenitors are still unknown. This includes results from the SDSS and our own ongoing pencil-beam halo survey, the Spaghetti survey.

We analyse a high resolution simulation of the formation of a cluster's dark-matter halo in a $\Lambda$CDM cosmology (Springel et al. 2001). The resolution achieved allows us to map the phase-space structure in detail, and characterize its evolution and degree of lumpiness. Scaling down the cluster halo to a Milky-Way size halo, we probe the substructure expected in the solar neighbourhood. Here we specifically address the relevance of such substructure for direct detection experiments aimed at determining the nature of dark-matter.

As part of the Spaghetti Project Survey (SPS) we have detected a concentration of giant stars well above expectations for a smooth halo model. The position (l~350, b~50) and distance (~50 kpc) of this concentration match those of the Northern over-density detected by SDSS (Yanny et al. 2000, Ivezic et al. 2000). We find additional evidence for structure at ~80 kpc in the same direction. We present radial velocities for many of these stars, including the first published results from the 6.5m Magellan telescope. The radial velocities for stars in these structures are in excellent agreement with models of the dynamical evolution of the Sgr dwarf tidal debris, whose center is 60 degrees away. The metallicity of stars in these streams is lower than that of the main body of the Sgr dwarf, which may indicate a radial metallicity gradient prior to disruption.

We present two simple dynamical models for Sagittarius based on N-body simulations of the progressive disruption of a satellite galaxy orbiting for 12.5 Gyr within a realistic Galactic potential. In both models the satellite initially has observable properties similar to those of current outlying dwarfs; in one case it is purely stellar while in the other it is embedded in an extended massive halo. The purely stellar progenitor is a King model with a total velocity dispersion of 18.9 km/s, a core radius of 0.44 kpc and a tidal radius of 3 kpc. The initial stellar distribution in the other case follows a King profile with the same core radius, a slightly larger total velocity dispersion and similar extent. Both these models are consistent with all published data on the current Sagittarius system, they match not only the observed properties of the main body of Sagittarius, but also those reported for unbound debris at larger distances.

In this contribution we discuss the expected properties of our halo if our Galaxy had been built from mergers of smaller systems. Using both analytic arguments and high-resolution simulations, we find that the Galactic halo in the Solar neighbourhood should have a smooth spatial distribution and that its stellar velocity ellipsoid should consist of several hundred kinematically cold streams, relicts from the different merging events. We also discuss observational evidence which supports this hierarchical picture: two stellar streams originating in the same system that probably fell onto the Milky Way about ten Gyr ago. Finally we address what future astrometric missions like FAME or GAIA may reveal by observing the motions of nearby halo stars.

We run numerical simulations of the disruption of satellite galaxies in a Galactic potential to build up the entire stellar halo, in order to investigate what the next generation of astrometric satellites will reveal by observing the halo of the Milky Way. We generate artificial DIVA, FAME and GAIA halo catalogues, in which we look for the signatures left by the accreted satellites. We develop a method based on the standard Friends-of-Friends algorithm applied to the space of integrals of motion. We find this simple method can recover about 50% of the different accretion events, when the observational uncertainties expected for GAIA are taken into account, even when the exact form of the Galactic potential is unknown. The recovery rate for DIVA and FAME is much smaller, but these missions, like GAIA, should be able to test the hierarchical formation paradigm on our Galaxy by measuring the amount of halo substructure in the form of nearby kinematically cold streams with for example, a two-point correlation function in velocity space.

It is now generally believed that galaxies were built up through gravitational amplification of primordial fluctuations and the subsequent merging of smaller precursor structures. The stars of the structures that assembled to form the Milky Way should now make up much or all of its bulge and halo, in which case one hopes to find "fossil" evidence for those precursor structures in the present distribution of halo stars. Confirmation that this process is continuing came with the discovery of the Sagittarius dwarf galaxy, which is being disrupted by the Milky Way, but direct evidence that this process provided the bulk of the Milky Way's population of old stars has so far been lacking. Here we show that about ten per cent of the metal--poor stars in the halo of the Milky Way, outside the radius of the Sun's orbit, come from a single coherent structure that was disrupted during or soon after the Galaxy's formation. This object had a highly inclined orbit about the Milky Way at a maximum distance of $\sim$ 16 kpc, and it probably resembled the Fornax and Sagittarius dwarf spheroidal galaxies.

We study numerical simulations of satellite galaxy disruption in a potential resembling that of the Milky Way. Our goal is to assess whether a merger origin for the stellar halo would leave observable fossil structure in the phase-space distribution of nearby stars. We show how mixing of disrupted satellites can be quantified using a coarse-grained entropy. Although after 10 Gyr few obvious asymmetries remain in the distribution of particles in configuration space, strong correlations are still present in velocity space. We give a simple analytic description of these effects, based on a linearised treatment in action-angle variables, which shows how the kinematic and density structure of the debris stream changes with time. By applying this description we find that a single satellite of current luminosity $10^8 L_\sun$ disrupted 10 Gyr ago from an orbit circulating in the inner halo (mean apocentre $\sim 12$ kpc) would contribute about $\sim 30$ kinematically cold streams with internal velocity dispersions below 5 km/s to the local stellar halo. If the whole stellar halo were built by disrupted satellites, it should consist locally of 300 - 500 such streams. Clear detection of all these structures would require a sample of a few thousand stars with 3-D velocities accurate to better than 5 km/s. Even with velocity errors several times worse than this, the expected clumpiness should be quite evident. We apply our formalism to a group of stars detected near the North Galactic Pole, and derive an order of magnitude estimate for the initial properties of the progenitor system.

We study numerical simulations of satellite galaxy disruption in a potential resembling that of the Milky Way. Our goal is to assess whether a merger origin for the stellar halo would leave observable fossil structure in the phase-space distribution of nearby stars. We show how mixing of disrupted satellites can be quantified using a coarse-grained entropy. Although after 10 billion years few obvious asymmetries remain in the distribution of particles in configuration space, strong correlations are still present in velocity space. We briefly describe how we can understand these effects based on the conservation of fine-grained phase-space density in an action-angle formalism. We also discuss the implications of our results on the known properties of the stellar halo.

We investigate what the proposed ESA astrometric satellite GAIA will reveal by observing the halo of the Milky Way. Specifically, we look for halo streams which are the signatures left by the merging/accretion events experienced by a typical galaxy like the Milky Way. We run numerical simulations of the disruption of satellite galaxies in a Galactic potential to generate artificial GAIA halo catalogs. We recover the streams by searching for peaks in angular momentum space.