Large multi-object spectroscopic surveys require automated algorithms to optimise their observing strategy. One of the most ambitious upcoming spectroscopic surveys is the 4MOST survey. The 4MOST survey facility is a fibre-fed spectroscopic instrument on the VISTA telescope with a large enough field of view to survey a large fraction of the southern sky within a few years. Several Galactic and extragalactic surveys will be carried out simultaneously, so the combined target density will strongly vary. In this paper, we describe a new tiling algorithm that can naturally deal with the large target density variations on the sky and which automatically handles the different exposure times of targets. The tiling pattern is modelled as a marked point process, which is characterised by a probability density that integrates the requirements imposed by the 4MOST survey. The optimal tilling pattern with respect to the defined model is estimated by the tiles configuration that maximises the proposed probability density. In order to achieve this maximisation a simulated annealing algorithm is implemented. The algorithm automatically finds an optimal tiling pattern and assigns a tentative sky brightness condition and exposure time for each tile, while minimising the total execution time that is needed to observe the list of targets in the combined input catalogue of all surveys. Hence, the algorithm maximises the long-term observing efficiency and provides an optimal tiling solution for the survey. While designed for the 4MOST survey, the algorithm is flexible and can with simple modifications be applied to any other multi-object spectroscopic survey.

We study the late-time evolution of the central regions of two Milky Way-like simulations of galaxies formed in a cosmological context, one hosting a fast bar and the other a slow one. We find that bar length, R_b, measurements fluctuate on a dynamical timescale by up to 100%, depending on the spiral structure strength and measurement threshold. The bar amplitude oscillates by about 15%, correlating with R_b. The Tremaine-Weinberg-method estimates of the bars' instantaneous pattern speeds show variations around the mean of up to ~20%, typically anti-correlating with the bar length and strength. Through power spectrum analyses, we establish that these bar pulsations, with a period in the range ~60-200 Myr, result from its interaction with multiple spiral modes, which are coupled with the bar. Because of the presence of odd spiral modes, the two bar halves typically do not connect at exactly the same time to a spiral arm, and their individual lengths can be significantly offset. We estimated that in about 50% of bar measurements in Milky Way-mass external galaxies, the bar lengths of SBab type galaxies are overestimated by ~15% and those of SBbc types by ~55%. Consequently, bars longer than their corotation radius reported in the literature, dubbed "ultra-fast bars", may simply correspond to the largest biases. Given that the Scutum-Centaurus arm is likely connected to the near half of the Milky Way bar, recent direct measurements may be overestimating its length by 1-1.5 kpc, while its present pattern speed may be 5-10 km/s/kpc smaller than its time-averaged value.

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.

With backwards orbit integration we estimate birth locations of young stellar associations and moving groups identified in the solar neighborhood that are younger than 70 Myr. The birth locations of most of these stellar associations are at smaller galactocentric radius than the Sun, implying that their stars moved radially outwards after birth. Exceptions to this rule include are the Argus and Octans associations which formed outside the Sun's Galactocentric radius. Variations in birth heights of the stellar associations suggest that they were born in a corrugated disk of molecular clouds, similar to that inferred from the current filamentary molecular cloud distribution and dust extinction maps. Multiple spiral arm features with different but near corotation pattern speeds and at different heights could account for the stellar association birth sites. We find that the young stellar associations are located in between peaks in the UV stellar velocity distribution for stars in the solar neighborhood. This would be expected if they were born in a spiral arm which perturbs stellar orbits that cross it. In contrast, stellar associations seem to be located near peaks in the vertical phase space distribution, suggesting that the gas in which stellar associations are born moves vertically together with the low velocity dispersion disk stars.

Data-driven methods play an increasingly important role in the field of astrophysics. In the context of large spectroscopic surveys of stars, data-driven methods are key in deducing physical parameters for millions of spectra in a short time. Convolutional Neural Networks (CNNs) enable us to connect observables (e.g., spectra, stellar magnitudes) to physical properties (atmospheric parameters, chemical abundances, or labels in general). We trained a CNN, adopting stellar atmospheric parameters and chemical abundances from APOGEE DR16 (resolution R~22500) data as training set labels. As input we used parts of the intermediate-resolution RAVE DR6 spectra (R~7500) overlapping with the APOGEE DR16 data as well as broad-band ALL_WISE and 2MASS photometry together with Gaia DR2 photometry and parallaxes. We derived precise atmospheric parameters Teff, log(g), [M/H] and chemical abundances of [Fe/H], [alpha/M], [Mg/Fe], [Si/Fe], [Al/Fe], and [Ni/Fe] for 420165 RAVE spectra. The precision typically amounts to 60K in Teff, 0.06 in log(g), and 0.02-0.04 dex for individual chemical abundances. Incorporating photometry and astrometry as additional constraints substantially improves the results in terms of accuracy and precision of the derived labels, as long as we operate in parts of the parameter space that is well covered by the training sample. Scientific validation confirms the robustness of the CNN results. We provide a catalog of CNN-trained atmospheric parameters and abundances along with their uncertainties for 420165 stars in the RAVE survey. CNN-based methods provide a powerful way to combine spectroscopic, photometric, and astrometric data, without applying any priors in the form of stellar evolutionary models. The developed procedure can extend the scientific output of RAVE spectra beyond DR6 to on-going and planned surveys such as Gaia RVS and 4MOST.

We combine high-resolution spectroscopic data from APOGEE-2 Survey Data Release 16 (DR16) with broad-band photometric data from several sources, as well as parallaxes from {\it Gaia} Data Release 2 (DR2). Using the Bayesian isochrone-fitting code {\tt StarHorse}, we derive distances, extinctions and astrophysical parameters for around 388,815 APOGEE stars, achieving typical distance uncertainties of $\sim 6\%$ for APOGEE giants, $\sim 2\%$ for APOGEE dwarfs, as well as extinction uncertainties of $\sim 0.07$ mag when all photometric information is available, and $\sim 0.17$ mag if optical photometry is missing. {\tt StarHorse} uncertainties vary with the input spectroscopic catalogue, with the available photometry, and with the parallax uncertainties. To illustrate the impact of our results, we show that, thanks to {\it Gaia} DR2 and the now larger sky coverage of APOGEE-2 (including APOGEE-South), we obtain an extended map of the Galactic plane, providing unprecedented coverage of the disk close to the Galactic mid-plane ($|Z_{Gal}|<1$ kpc) from the Galactic Centre out to $R_{\rm Gal}\sim 20$ kpc. The improvements in statistics as well as distance and extinction uncertainties unveil the presence of the bar in stellar density, as well as the striking chemical duality in the innermost regions of the disk, now clearly extending to the inner bulge. We complement this paper with distances and extinctions for stars in other public released spectroscopic surveys: 324,999 in GALAH DR2, 4,928,715 in LAMOST DR5, 408,894 in RAVE DR6, and 6,095 in GES DR3

Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys.

The Milky Way provides an ideal laboratory to test our understanding of galaxy evolution, owing to our ability to observe our Galaxy over fine scales. However, connecting the Galaxy to the wider galaxy population remains difficult, due to the challenges posed by our internal perspective and to the different observational techniques employed. Here, we present a sample of galaxies identified as Milky Way Analogs (MWAs) on the basis of their stellar masses and bulge-to-total ratios, observed as part of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. We analyse the galaxies in terms of their stellar kinematics and populations as well as their ionised gas contents. We find our sample to contain generally young stellar populations in their outskirts. However, we find a wide range of stellar ages in their central regions, and we detect central AGN-like or composite-like activity in roughly half of the sample galaxies, with the other half consisting of galaxies with central star-forming emission or emission consistent with old stars. We measure gradients in gas metallicity and stellar metallicity that are generally flatter in physical units than those measured for the Milky Way; however, we find far better agreement with the Milky Way when scaling gradients by galaxies' disc scale lengths. From this, we argue much of the discrepancy in metallicity gradients to be due to the relative compactness of the Milky Way, with differences in observing perspective also likely to be a factor.

Simulating thin and extended galactic disks has long been a challenge in computational astrophysics. We introduce the NIHAO-UHD suite of cosmological hydrodynamical simulations of Milky Way mass galaxies and study stellar disk properties such as stellar mass, size and rotation velocity which agree well with observations of the Milky Way and local galaxies. In particular, the simulations reproduce the age-velocity dispersion relation and a multi-component stellar disk as observed for the Milky Way. Half of our galaxies show a double exponential vertical profile, while the others are well described by a single exponential model which we link to the disk merger history. In all cases, mono-age populations follow a single exponential whose scale height varies monotonically with stellar age and radius. The scale length decreases with stellar age while the scale height increases. The general structure of the stellar disks is already set at time of birth as a result of the inside-out and upside-down formation. Subsequent evolution modifies this structure by increasing both the scale length and height of all mono-age populations. Thus, our results put tight constraints on how much dynamical memory stellar disks can retain over cosmological timescales. Our simulations demonstrate that it is possible to form thin galactic disks in cosmological simulations provided there are no significant stellar mergers at low redshifts. Most of the stellar mass is formed in-situ with only a few percent ($\lesssim5\%$) brought in by merging satellites at early times. Redshift zero snapshots and halo catalogues are publicly available.

We present isochrone ages and initial bulk metallicities ($\rm [Fe/H]_{bulk}$, by accounting for diffusion) of 163,722 stars from the GALAH Data Release 2, mainly composed of main sequence turn-off stars and subgiants ($\rm 7000 K>T_{eff}>4000 K$ and $\rm log g>3$ dex). The local age-metallicity relationship (AMR) is nearly flat but with significant scatter at all ages; the scatter is even higher when considering the observed surface abundances. After correcting for selection effects, the AMR appear to have intrinsic structures indicative of two star formation events, which we speculate are connected to the thin and thick disks in the solar neighborhood. We also present abundance ratio trends for 16 elements as a function of age, across different $\rm [Fe/H]_{bulk}$ bins. In general, we find the trends in terms of [X/Fe] vs age from our far larger sample to be compatible with studies based on small ($\sim$ 100 stars) samples of solar twins but we now extend it to both sub- and super-solar metallicities. The $\alpha$-elements show differing behaviour: the hydrostatic $\alpha$-elements O and Mg show a steady decline with time for all metallicities while the explosive $\alpha$-elements Si, Ca and Ti are nearly constant during the thin disk epoch (ages $\lessapprox $ 12 Gyr). The s-process elements Y and Ba show increasing [X/Fe] with time while the r-process element Eu have the opposite trend, thus favouring a primary production from sources with a short time-delay such as core-collapse supernovae over long-delay events such as neutron star mergers.

We use Gaia DR2 astrometric and line-of-sight velocity information combined with two sets of distances obtained with a Bayesian inference method to study the 3D velocity distribution in the Milky Way disc. We search for variations in all Galactocentric cylindrical velocity components ($V_{\phi}$, $V_R$ and $V_z$) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent work showing that bulk vertical motions in the $R\text{-}z$ plane are consistent with a combination of breathing and bending modes. In the $x\text{-}y$ plane, we show that, although the amplitudes change, the structure produced by these modes is mostly invariant as a function of distance from the plane. Comparing to two different Galactic disc models, we demonstrate that the observed patterns can drastically change in short time intervals, showing the complexity of understanding the origin of vertical perturbations. A strong radial $V_R$ gradient was identified in the inner disc, transitioning smoothly from $16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $30^\circ<\phi<45^\circ$ ahead of the Sun-Galactic centre line, to $-16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $-45^\circ<\phi<-30^\circ$ lagging the solar azimuth. We use a simulation with no significant recent mergers to show that exactly the opposite trend is expected from a barred potential, but overestimated distances can flip this trend to match the data. Alternatively, using an $N$-body simulation of the Sagittarius dwarf-Milky Way interaction, we demonstrate that a major recent perturbation is necessary to reproduce the observations. Such an impact may have strongly perturbed the existing bar or even triggered its formation in the last $1\text{-}2$ Gyr.

Galactic astrophysics is now in the process of building a multi-dimensional map of the Galaxy. For such a map, stellar ages are the essential ingredient. Ages are however measured only indirectly by comparing observational data with models. It is often difficult to provide a single age value for a given star, as several non-overlapping solutions are possible. We aim at recovering the underlying log(age) distribution from the measured log(age) probability density function for an arbitrary set of stars. We build an age inversion method, namely, we represent the measured log(age) probability density function as a weighted sum of probability density functions of mono-age populations. Weights in that sum give the underlying log(age) distribution. Mono-age populations are simulated so that the distribution of stars on the log g-[Fe/H] plane is close to that of the observed sample. We tested the age inversion method on simulated data, demonstrating that it is capable of properly recovering the true log(age) distribution for a large (N > 103) sample of stars. The method was further applied to large public spectroscopic surveys. For RAVE-on, LAMOST and APOGEE we also applied age inversion to mono-metallicity samples, successfully recovering age-metallicity trends present in higher-precision APOGEE data and chemical evolution models. We conclude that applying an age inversion method as presented in this work is necessary to recover the underlying age distribution of a large (N > 103 ) set of stars. These age distributions can be used to explore for instance age-metallicity relations.

The next decade affords tremendous opportunity to achieve the goals of Galactic archaeology. That is, to reconstruct the evolutionary narrative of the Milky Way, based on the empirical data that describes its current morphological, dynamical, temporal and chemical structures. Here, we describe a path to achieving this goal. The critical observational objective is a Galaxy-scale, contiguous, comprehensive mapping of the disk's phase space, tracing where the majority of the stellar mass resides. An ensemble of recent, ongoing, and imminent surveys are working to deliver such a transformative stellar map. Once this empirical description of the dust-obscured disk is assembled, we will no longer be operationally limited by the observational data. The primary and significant challenge within stellar astronomy and Galactic archaeology will then be in fully utilizing these data. We outline the next-decade framework for obtaining and then realizing the potential of the data to chart the Galactic disk via its stars. One way to support the investment in the massive data assemblage will be to establish a Galactic Archaeology Consortium across the ensemble of stellar missions. This would reflect a long-term commitment to build and support a network of personnel in a dedicated effort to aggregate, engineer, and transform stellar measurements into a comprehensive perspective of our Galaxy.

Combining the precise parallaxes and optical photometry delivered by Gaia's second data release (Gaia DR2) with the photometric catalogues of PanSTARRS-1, 2MASS, and AllWISE, we derive Bayesian stellar parameters, distances, and extinctions for 265 million stars brighter than G=18. Because of the wide wavelength range used, our results substantially improve the accuracy and precision of previous extinction and effective temperature estimates. After cleaning our results for both unreliable input and output data, we retain 137 million stars, for which we achieve a median precision of 5% in distance, 0.20 mag in V-band extinction, and 245 K in effective temperature for G<14, degrading towards fainter magnitudes (12%, 0.20 mag, and 245 K at G=16; 16%, 0.23 mag, and 260 K at G=17, respectively). We find a very good agreement with the asteroseismic surface gravities and distances of 7000 stars in the Kepler, the K2-C3, and the K2-C6 fields, with stellar parameters from the APOGEE survey, as well as with distances to star clusters. Our results are available through the ADQL query interface of the Gaia mirror at the Leibniz-Institut f\"{u}r Astrophysik Potsdam (gaia.aip.de) and as binary tables at data.aip.de. As a first application, in this paper we provide distance- and extinction-corrected colour-magnitude diagrams, extinction maps as a function of distance, and extensive density maps, demonstrating the potential of our value-added dataset for mapping the three-dimensional structure of our Galaxy. In particular, we see a clear manifestation of the Galactic bar in the stellar density distributions, an observation that can almost be considered a direct imaging of the Galactic bar.

Galactic disc chemical evolution models generally ignore azimuthal surface density variation that can introduce chemical abundance azimuthal gradients. Recent observations, however, have revealed chemical abundance changes with azimuth in the gas and stellar components of both the Milky Way and external galaxies. To quantify the effects of spiral arm density fluctuations on the azimuthal variations of the oxygen and iron abundances in disc galaxies. We develop a new 2D galactic disc chemical evolution model, capable of following not just radial but also azimuthal inhomogeneities. The density fluctuations resulting from a Milky Way-like N-body disc formation simulation produce azimuthal variations in the oxygen abundance gradients of the order of 0.1 dex. Moreover, in agreement with the most recent observations in external galaxies, the azimuthal variations are more evident in the outer galactic regions. Using a simple analytical model, we show that the largest fluctuations with azimuth result near the spiral structure corotation resonance, where the relative speed between spiral and gaseous disc is the slowest. In conclusion we provided a new 2D chemical evolution model capable of following azimuthal density variations. Density fluctuations extracted from a Milky Way-like dynamical model lead to a scatter in the azimuthal variations of the oxygen abundance gradient in agreement with observations in external galaxies. We interpret the presence of azimuthal scatter at all radii by the presence of multiple spiral modes moving at different pattern speeds, as found in both observations and numerical simulations.

Simpson's paradox, or Yule-Simpson effect, arises when a trend appears in different subsets of data but disappears or reverses when these subsets are combined. We describe here seven cases of this phenomenon for chemo-kinematical relations believed to constrain the Milky Way disk formation and evolution. We show that interpreting trends in relations, such as the radial and vertical chemical abundance gradients, the age-metallicity relation, and the metallicity-rotational velocity relation (MVR), can lead to conflicting conclusions about the Galaxy past if analyses marginalize over stellar age and/or birth radius. It is demonstrated that the MVR in RAVE giants is consistent with being always strongly negative, when narrow bins of [Mg/Fe] are considered. This is directly related to the negative radial metallicity gradients of stars grouped by common age (mono-age populations) due to the inside out disk formation. The effect of the asymmetric drift can then give rise to a positive MVR trend in high-[alpha/Fe] stars, with a slope dependent on a given survey's selection function and observational uncertainties. We also study the variation of lithium abundance, A(Li), with [Fe/H] of AMBRE:HARPS dwarfs. A strong reversal in the positive A(Li)-[Fe/H] trend of the total sample is found for mono-age populations, flattening for younger groups of stars. Dissecting by birth radius shows strengthening in the positive A(Li)-[Fe/H] trend, shifting to higher [Fe/H] with decreasing birth radius; these observational results suggest new constraints on chemical evolution models. This work highlights the necessity for precise age estimates for large stellar samples covering wide spatial regions.

The mechanisms of the formation and evolution of the Milky Way are encoded in the orbits, chemistry and ages of its stars. With the 4MOST MIlky way Disk And BuLgE Low-Resolution Survey (4MIDABLE-LR) we aim to study kinematic and chemical substructures in the Milky Way disc and bulge region with samples of unprecedented size out to larger distances and greater precision than conceivable with Gaia alone or any other ongoing or planned survey. Gaia gives us the unique opportunity for target selection based almost entirely on parallax and magnitude range, hence increasing the efficiency in sampling larger Milky Way volumes with well-defined and effective selection functions. Our main goal is to provide a detailed chrono-chemo-kinematical extended map of our Galaxy and the largest Gaia follow-up down to $G = 19$ magnitudes (Vega). The complex nature of the disc components (for example, large target densities and highly structured extinction distribution in the Milky Way bulge and disc area), prompted us to develop a survey strategy with five main sub-surveys that are tailored to answer the still open questions about the assembly and evolution of our Galaxy, while taking full advantage of the Gaia data.

The signatures of the formation and evolution of a galaxy are imprinted in its stars. Their velocities, ages, and chemical compositions present major constraints on models of galaxy formation, and on various processes such as the gas inflows and outflows, the accretion of cold gas, radial migration, and the variability of star formation activity. Understanding the evolution of the Milky Way requires large observational datasets of stars via which these quantities can be determined accurately. This is the science driver of the 4MOST MIlky way Disc And BuLgE High-Resolution (4MIDABLE-HR) survey: to obtain high-resolution spectra at $R \sim 20\,000$ and to provide detailed elemental abundances for large samples of stars in the Galactic disc and bulge. High data quality will allow us to provide accurate spectroscopic diagnostics of two million stellar spectra: precise radial velocities; rotation; abundances of many elements, including those that are currently only accessible in the optical, such as Li, s-, and r-process; and multi-epoch spectra for a sub-sample of stars. Synergies with complementary missions like Gaia and TESS will provide masses, stellar ages and multiplicity, forming a multi-dimensional dataset that will allow us to explore and constrain the origin and structure of the Milky Way.

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.

We propose here that the lithium decrease at super-solar metallicities observed in high resolution spectroscopic surveys can be explained by the interplay of mixed populations, coming from the inner regions of the Milky Way disc. The lower lithium content of these stars is a consequence of inside-out disc formation, plus radial migration. In this framework, local stars with super-solar metallicities would have migrated to the solar vicinity and depleted their original lithium during their travel time. To arrive to such a result, we took advantage of the AMBRE catalog of lithium abundances combined with chemical evolution models which take into account the contribution to the lithium enrichment by different nucleosynthetic sources. A large proportion of migrated stars can explain the observed lower lithium abundance at super-solar metallicities. We stress that nowadays, there is no stellar model able to predict Li-depletion for such super-solar metallicity stars, and the Solar Li-depletion has to be assumed. In addition, it currently exists no solid quantitative estimate of the proportion of migrated stars in the Solar neighborhood and their travel time. Our results illustrate how important it is to properly include radial migration when comparing chemical evolution models to observations, and that in this case, the lithium decrease at larger metallicities does not necessarily imply that stellar yields have to be modified, contrary to previous claims in literature.

We analyse an N-body simulation of the interaction of the Milky Way (MW) with a Sagittarius-like dSph (Sgr), looking for signatures which may be attributed to its orbital history in the phase space volume around the Sun in light of $Gaia$ DR2 discoveries. The repeated impacts of Sgr excite coupled vertical and radial oscillations in the disc which qualitatively, and to a large degree quantitatively are able to reproduce many features in the 6D $Gaia$ DR2 samples, from the median $V_{R},V_{\phi}, V_{z} $ velocity maps to the local $\delta\rho(v_{z},z)$ phase-space spiral which is a manifestation of the global disc response to coupled oscillations within a given volume. The patterns in the large-scale velocity field are well described by tightly wound spirals and vertical corrugations excited from Sgr's impacts. We show that the last pericentric passage of Sgr resets the formation of the local present-day $\delta\rho(v_z,z)$ spiral and situate its formation around 500-800 Myr. As expected $\delta\rho(vz,z)$ grows in size and decreases in woundedness as a function of radius in both the $Gaia$ DR2 data and simulations. This is the first N-body model able to explain so many of the features in the data on different scales. We demonstrate how to use the full extent of the Galactic disc to date perturbations dating from Myr to Gyr, probe the underlying potential and constrain the mass-loss history of Sgr. $\delta\rho(vz,z)$ looks the same in all stellar populations age bins down to the youngest ages which rules out a bar buckling origin.

We present and discuss the stellar kinematics and populations of the S0 galaxy FCC 170 (NGC 1381) in the Fornax cluster, using deep MUSE data from the Fornax 3D survey. We show the maps of the first four moments of the stellar line-of-sight velocity distribution and of the mass-weighted mean stellar age, metallicity and [Mg/Fe] abundance ratio. The high-quality MUSE stellar kinematic measurements unveil the structure of this massive galaxy: a nuclear disk, a bar seen as a boxy bulge with a clear higher-velocity-dispersion X shape, a fast-rotating and flaring thin disk and a slower rotating thick disk. Whereas their overall old age makes it difficult to discuss differences in the formation epoch between these components, we find a clear-cut distinction between metal-rich and less [Mg/Fe]-enhanced populations in the thin-disk, boxy-bulge and nuclear disk, and more metal-poor and [Mg/Fe]-enhanced stars in the thick disk. Located in the densest region of the Fornax cluster, where signs of tidal stripping have been recently found, the evolution of FCC 170 might have been seriously affected by its environment. We discuss the possibility of its "pre-processing" in a subgroup before falling into the present-day cluster, which would have shaped this galaxy a long time ago. The thick disk displays a composite star formation history, as a significant fraction of younger stars co-exist with the main older thick-disk population. The former sub-population is characterized by even lower-metallicity and higher-[Mg/Fe] values, suggesting that these stars formed later and faster in a less chemically evolved satellite, which was subsequently accreted. Finally, we discuss evidence that metal-rich and less [Mg/Fe]-enhanced stars were brought in the outer parts of the thick disk by the flaring of the thin disk.

We present [C/N]-[Fe/H] abundance trends from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey, Data Release 14 (DR14), for red giant branch stars across the Milky Way Galaxy (MW, 3 kpc $<$ R $<$ 15 kpc). The carbon-to-nitrogen ratio (often expressed as [C/N]) can indicate the mass of a red giant star, from which an age can be inferred. Using masses and ages derived by Martig et al., we demonstrate that we are able to interpret the DR14 [C/N]-[Fe/H] abundance distributions as trends in age-[Fe/H] space. Our results show that an anti-correlation between age and metallicity, which is predicted by simple chemical evolution models, is not present at any Galactic zone. Stars far from the plane ($|$Z$|$ $>$ 1 kpc) exhibit a radial gradient in [C/N] ($\sim$ $-$0.04 dex/kpc). The [C/N] dispersion increases toward the plane ($\sigma_{[C/N]}$ = 0.13 at $|$Z$|$ $>$ 1 kpc to $\sigma_{[C/N]}$ = 0.18 dex at $|$Z$|$ $<$ 0.5 kpc). We measure a disk metallicity gradient for the youngest stars (age $<$ 2.5 Gyr) of $-$0.060 dex/kpc from 6 kpc to 12 kpc, which is in agreement with the gradient found using young CoRoGEE stars by Anders et al. Older stars exhibit a flatter gradient ($-$0.016 dex/kpc), which is predicted by simulations in which stars migrate from their birth radii. We also find that radial migration is a plausible explanation for the observed upturn of the [C/N]-[Fe/H] abundance trends in the outer Galaxy, where the metal-rich stars are relatively enhanced in [C/N].

One of the fundamental assumptions of chemical evolution models (CEMs) of the Milky Way (MW) and other spirals is that higher gas accretion rates are expected in the past, and in the inner regions of the Galaxy. This leads to the so-called `inside-out disc formation scenario'. Yet, these are probably the most unconstrained inputs of such models. In the present paper, we aim at investigating these main assumptions by studying how gas is accreted in four simulated MW-like galaxies assembled within the $\Lambda$CDM scenario. The galaxies were obtained using two different simulation techniques, cosmological setups and initial conditions. Two of them are MW candidates corresponding to the chemodynamical model of Minchev et al. (2013, 2014) (known as MCM) and the Local Group cosmological simulation of Nuza et al. (2014). We investigate vertical and radial gas accretion on to galaxy discs as a function of cosmic time and disc radius. We find that accretion in the MW-like galaxies seem to happen in two distinct phases, namely: an early, more violent period; followed by a subsequent, slowly declining phase. Our simulations seem to give support to the assumption that the amount of gas incorporated into the MW disc exponentially decreases with time, leading to current net accretion rates of $0.6-1\,$M$_\odot\,$yr$^{-1}$. In particular, accretion timescales on to the simulated thin-disc-like structures are within $\sim5-7\,$Gyr, consistent with expectations from CEMs. Moreover, our simulated MW discs are assembled from the inside-out with gas in the inner disc regions accreted in shorter timescales than in external ones, in qualitative agreement with CEMs of the Galaxy. However, this type of growth is not general to all galaxies and it is intimately linked to their particular merger and gas accretion history.

HST/ACS observations along the major axis of M33 show that the mean age of its stars decreases with increasing distance from the galaxy center. Such a behavior is consistent with an inside-out growth of the disc. However, in the outermost observed field, at r$\simeq$11.6 kpc, a reversal of this gradient is detected, with old stars found in high percentages beyond this radius. In this work we investigate the origin of such a reversal in stellar age gradient, by using a simulated M33 analogue from the Constrained Local UniversE Simulations (CLUES). The simulated M33 is similar to the observed one in terms of mass, rotation velocity, surface brightness and, similar to what has been reported in observations, shows a stellar age turnaround at large radii. We demonstrate that this reversal is mostly a result of stellar accretion from old satellite galaxies and, to a lesser extent, of stellar migration of in-situ stars. The old accreted stars, with formation times tf < 4 Gyrs, are kinematically hot and can be differentiated from the in-situ stars by their high velocity dispersion and the fact that they do not have rotationally-supported orbits. In the future, obtaining kinematic information of the stars in the outskirt of M33 will help to verify this scenario.

We present a semi-empirical, largely model-independent approach for estimating Galactic birth radii, r_birth, for Milky Way disk stars. The technique relies on the justifiable assumption that a negative radial metallicity gradient in the interstellar medium (ISM) existed for most of the disk lifetime. Stars are projected back to their birth positions according to the observationally derived age and [Fe/H] with no kinematical information required. Applying our approach to the AMBRE:HARPS and HARPS-GTO local samples, we show that we can constrain the ISM metallicity evolution with Galactic radius and cosmic time, [Fe/H]_ISM(r, t), by requiring a physically meaningful r_birth distribution. We find that the data are consistent with an ISM radial metallicity gradient that flattens with time from ~-0.15 dex/kpc at the beginning of disk formation, to its measured present-day value (-0.07 dex/kpc). We present several chemo-kinematical relations in terms of mono-r_birth populations. One remarkable result is that the kinematically hottest stars would have been born locally or in the outer disk, consistent with thick disk formation from the nested flares of mono-age populations and predictions from cosmological simulations. This phenomenon can be also seen in the observed age-velocity dispersion relation, in that its upper boundary is dominated by stars born at larger radii. We also find that the flatness of the local age-metallicity relation (AMR) is the result of the superposition of the AMRs of mono-r_birth populations, each with a well-defined negative slope. The solar birth radius is estimated to be 7.3+-0.6 kpc, for a current Galactocentric radius of 8 kpc.

The solar neighborhood contains disc stars that have recently crossed spiral arms in the Galaxy. We propose that boundaries in local velocity distributions separate stars that have recently crossed or been perturbed by a particular arm from those that haven't. Ridges in the stellar velocity distributions constructed from the second Gaia data release trace orbits that could have touched nearby spiral arms at apocentre or pericentre. The multiple ridges and arcs seen in local velocity distributions are consistent with the presence of multiple spiral features and different pattern speeds and imply that the outer Galaxy is flocculent rather than grand design.

Before the publication of the Gaia DR2 we confirmed with RAVE and TGAS an observation recently made with the GALAH survey by Quillen ey al. concerning the Coma Berenices moving group in the Solar neighbourhood, namely that it is only present at negative Galactic latitudes. This allowed us to show that it is coherent in vertical velocity, providing a first evidence for incomplete vertical phase-mixing. We estimated for the first time from dynamical arguments that the moving group must have formed at most ~ 1.5 Gyr ago, and related this to a pericentric passage of the Sagittarius dwarf satellite galaxy. The present note is a rewritten version of the original arXiv post on this result now also including a confirmation of our finding with Gaia DR2.

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

In chemodynamical evolution models it is usually assumed that the Milky Way galaxy forms from the inside-out implying that gas inflows onto the disk decrease with galactocentric distance. Similarly, to reproduce differences between chemical abundances of the thick disk and bulge with respect to those of the thin disk, higher accretion fluxes at early times are postulated. By using a suite of Milky Way-like galaxies extracted from cosmological simulations, we investigate the accretion of gas on the simulated stellar disks during their whole evolution. In general, we find that the picture outlined above holds, although the detailed behavior depends on the assembly history of the Galaxy and the complexities inherent to the physics of galaxy formation.

We explore the connections between stellar age, chemistry, and kinematics across a Galactocentric distance of $7.5 < R\,(\mathrm{kpc}) < 9.0$, using a sample of $\sim 12\,000$ intermediate-mass (FGK) turnoff stars observed with the RAdial Velocity Experiment (RAVE) survey. The kinematics of this sample are determined using radial velocity measurements from RAVE, and parallax and proper motion measurements from the Tycho-Gaia Astrometric Solution (TGAS). In addition, ages for RAVE stars are determined using a Bayesian method, taking TGAS parallaxes as a prior. We divide our sample into young ($0 < \tau < 3$ Gyr) and old ($8 < \tau < 13$ Gyr) populations, and then consider different metallicity bins for each of these age groups. We find significant differences in kinematic trends of young and old, metal-poor and metal-rich, stellar populations. In particular, we find a strong metallicity dependence in the mean Galactocentric radial velocity as a function of radius ($\partial {V_{\rm R}}/\partial R$) for young stars, with metal-rich stars having a much steeper gradient than metal-poor stars. For $\partial {V_{\phi}}/\partial R$, young, metal-rich stars significantly lag the LSR with a slightly positive gradient, while metal-poor stars show a negative gradient above the LSR. We interpret these findings as correlations between metallicity and the relative contributions of the non-axisymmetries in the Galactic gravitational potential (the spiral arms and the bar) to perturb stellar orbits.

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.

Using tracer particles embedded in self-gravitating shearing sheet N-body simulations, we investigate the distance in guiding centre radius that stars or star clusters can migrate in a few orbital periods. The standard deviations of guiding centre distributions and maximum migration distances depend on the Toomre or critical wavelength and the contrast in mass surface density caused by spiral structure. Comparison between our simulations and estimated guiding radii for a few young super-solar metallicity open clusters, including NGC 6583, suggests that the contrast in mass surface density in the solar neighbourhood has standard deviation (in the surface density distribution) divided by mean of about 1/4 and larger than measured using COBE data by Drimmel and Spergel. Our estimate is consistent with a standard deviation of $\sim$0.07 dex in the metallicities measured from high-quality spectroscopic data for 38 young open clusters (<1 Gyr) with mean galactocentric radius 7-9 kpc.

The aim of Galactic Archaeology is to recover the evolutionary history of the Milky Way from its present day kinematical and chemical state. Because stars move away from their birth sites, the current dynamical information alone is not sufficient for this task. The chemical composition of stellar atmospheres, on the other hand, is largely preserved over the stellar lifetime and, together with accurate ages, can be used to recover the birthplaces of stars currently found at the same Galactic radius. In addition to the availability of large stellar samples with accurate 6D kinematics and chemical abundance measurements, this requires detailed modeling with both dynamical and chemical evolution taken into account. An important first step is to understand the variety of dynamical processes that can take place in the Milky Way, including the perturbative effects of both internal (bar and spiral structure) and external (infalling satellites) agents. We discuss here (1) how to constrain the Galactic bar, spiral structure, and merging satellites by their effect on the local and global disc phase-space, (2) the effect of multiple patterns on the disc dynamics, and (3) the importance of radial migration and merger perturbations for the formation of the Galactic thick disc. Finally, we discuss the construction of Milky Way chemo-dynamical models and relate to observations.

We live in an age where an enormous amount of astrometric, photometric, asteroseismic, and spectroscopic data of Milky Way stars are being acquired, many orders of magnitude larger than about a decade ago. Thanks to the Gaia astrometric mission and followup ground-based spectroscopic surveys in the next 5-10 years about 10-20 Million stars will have accurate 6D kinematics and chemical composition measurements. KEPLER-2, PLATO, and TESS will provide asteroseismic ages for a good fraction of those. In this article we outline some outstanding problems concerning the formation and evolution of the Milky Way and argue that, due to the complexity of physical processes involved in the formation of disk galaxies, numerical simulations in the cosmological context are needed for the interpretation of Milky Way observations. We also discuss in some detail the formation of the Milky Way thick disk, chemodynamical models, and the effects of radial migration.

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

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.

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

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

Correlations between stellar chemistry and kinematics have long been used to gain insight into the evolution of the Milky Way Galaxy. Orbital angular momentum is a key physical parameter and it is often estimated from three-dimensional space motions. We here demonstrate the lower uncertainties that can be achieved in the estimation of one component of velocity through selection of stars in key directions and use of line-of-sight velocity alone (i.e. without incorporation of proper motion data). In this first paper we apply our technique to stars observed in the direction of Galactic rotation in the APOGEE survey. We first derive the distribution of azimuthal velocities, $v_\phi$, then from these and observed radial coordinates, estimate the stellar guiding centre radii, $R_g$, within $6.9\leq R \leq 10$ kpc with uncertainties smaller than (or of the order of) 1kpc. We show that there is no simple way to select a clean stellar sample based on low errors on proper motions and distances to obtain high-quality 3D velocities and hence one should pay particular attention when trying to identify kinematically peculiar stars based on velocities derived using the proper motions. Using our azimuthal velocity estimations, we investigate the joint distribution of elemental abundances and rotational kinematics free from the blurring effects of epicyclic motions, and we derive the $\partial v_\phi / \partial R_g$ trends for the thin and thick discs as a function of radius. Our analysis provides further evidence for radial migration within the thin disc and hints against radial migration playing a significant role in the evolution of the thick disc.

Using combined asteroseismic and spectroscopic observations of 418 red-giant stars close to the Galactic disc plane (6 kpc $<R_{\rm Gal}\lesssim13$ kpc, $|Z_{\rm Gal}|<0.3$ kpc), we measure the age dependence of the radial metallicity distribution in the Milky Way's thin disc over cosmic time. The slope of the radial iron gradient of the young red-giant population ($-0.058\pm0.008$ [stat.] $\pm0.003$ [syst.] dex/kpc) is consistent with recent Cepheid measurements. For stellar populations with ages of $1-4$ Gyr the gradient is slightly steeper, at a value of $-0.066\pm0.007\pm0.002$ dex/kpc, and then flattens again to reach a value of $\sim-0.03$ dex/kpc for stars with ages between 6 and 10 Gyr. Our results are in good agreement with a state-of-the-art chemo-dynamical Milky-Way model in which the evolution of the abundance gradient and its scatter can be entirely explained by a non-varying negative metallicity gradient in the interstellar medium, together with stellar radial heating and migration. We also offer an explanation for why intermediate-age open clusters in the Solar Neighbourhood can be more metal-rich, and why their radial metallicity gradient seems to be much steeper than that of the youngest clusters. Already within 2 Gyr, radial mixing can bring metal-rich clusters from the innermost regions of the disc to Galactocentric radii of 5 to 8 kpc. We suggest that these outward-migrating clusters may be less prone to tidal disruption and therefore steepen the local intermediate-age cluster metallicity gradient. Our scenario also explains why the strong steepening of the local iron gradient with age is not seen in field stars. In the near future, asteroseismic data from the K2 mission will allow for improved statistics and a better coverage of the inner-disc regions, thereby providing tighter constraints on the evolution of the central parts of the Milky Way.

Studying the Milky Way disk structure using stars in narrow bins of [Fe/H] and [alpha/Fe] has recently been proposed as a powerful method to understand the Galactic thick and thin disk formation. It has been assumed so far that these mono-abundance populations (MAPs) are also coeval, or mono-age, populations. Here we study this relationship for a Milky Way chemo-dynamical model and show that equivalence between MAPs and mono-age populations exists only for the high-[alpha/Fe] tail, where the chemical evolution curves of different Galactic radii are far apart. At lower [alpha/Fe]-values a MAP is composed of stars with a range in ages, even for small observational uncertainties and a small MAP bin size. Due to the disk inside-out formation, for these MAPs younger stars are typically located at larger radii, which results in negative radial age gradients that can be as large as 2 Gyr/kpc. Positive radial age gradients can result for MAPs at the lowest [alpha/Fe] and highest [Fe/H] end. Such variations with age prevent the simple interpretation of observations for which accurate ages are not available. Studying the variation with radius of the stellar surface density and scale-height in our model, we find good agreement to recent analyses of the APOGEE red-clump (RC) sample when 1-4 Gyr old stars dominate (as expected for the RC). Our results suggest that the APOGEE data are consistent with a Milky Way model for which mono-age populations flare for all ages. We propose observational tests for the validity of our predictions and argue that using accurate age measurements, such as from asteroseismology, is crucial for putting constraints on the Galactic formation and evolution.

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.

Understanding the nature of spiral structure in disk galaxies is one of the main, and still unsolved questions in galactic astronomy. However, theoretical works are proposing new testable predictions whose detection is becoming feasible with recent development in instrumentation. In particular, streaming motions along spiral arms are expected to induce azimuthal variations in the chemical composition of a galaxy at a given galactic radius. In this letter we analyse the gas content in NGC 6754 with VLT/MUSE data to characterise its 2D chemical composition and H$\alpha$ line-of-sight velocity distribution. We find that the trailing (leading) edge of the NGC 6754 spiral arms show signatures of tangentially-slower, radially-outward (tangentially-faster, radially-inward) streaming motions of metal-rich (poor) gas over a large range of radii. These results show direct evidence of gas radial migration for the first time. We compare our results with the gas behaviour in a $N$-body disk simulation showing spiral morphological features rotating with a similar speed as the gas at every radius, in good agreement with the observed trend. This indicates that the spiral arm features in NGC 6754 may be transient and rotate similarly as the gas does at a large range of radii.

In the Milky Way, the thick disk can be defined using individual stellar abundances, kinematics, or age; or geometrically, as stars high above the mid-plane. In nearby galaxies, where only a geometric definition can be used, thick disks appear to have large radial scale-lengths, and their red colors suggest that they are uniformly old. The Milky Way's geometrically thick disk is also radially extended, but it is far from chemically uniform: alpha-enhanced stars are confined within the inner Galaxy. In simulated galaxies, where old stars are centrally concentrated, geometrically thick disks are radially extended, too. Younger stellar populations flare in the simulated disks' outer regions, bringing those stars high above the mid-plane. The resulting geometrically thick disks therefore show a radial age gradient, from old in their central regions to younger in their outskirts. Based on our age estimates for a large sample of giant stars in the APOGEE survey, we can now test this scenario for the Milky Way. We find that the geometrically-defined thick disk in the Milky Way has indeed a strong radial age gradient: the median age for red clump stars goes from ~9 Gyr in the inner disk to 5 Gyr in the outer disk. We propose that at least some nearby galaxies could also have thick disks that are not uniformly old, and that geometrically thick disks might be complex structures resulting from different formation mechanisms in their inner and outer parts.

With the advent of the space missions CoRoT and Kepler, it has become feasible to determine precise asteroseismic masses and ages for large samples of red-giant stars. In this paper, we present the CoRoGEE dataset -- obtained from CoRoT lightcurves for 606 red giant stars in two fields of the Galactic disc which have been co-observed for an ancillary project of APOGEE. We have used the Bayesian parameter estimation code PARAM to calculate distances, extinctions, masses, and ages for these stars in a homogeneous analysis, resulting in relative statistical uncertainties of $\sim2\%$ in distance, $\sim4\%$ in radius, $\sim9\%$ in mass and $\sim25\%$ in age. We also assess systematic age uncertainties due to different input physics and mass loss. We discuss the correlation between ages and chemical abundance patterns of field stars over a large radial range of the Milky Way's disc (5 kpc $<R_{\rm Gal}<$ 14 kpc), focussing on the [$\alpha$/Fe]-[Fe/H]-age plane in five radial bins of the Galactic disc. We find an overall agreement with the expectations of chemical-evolution models computed before the present data were available, especially for the outer regions. However, our data also indicate that a significant fraction of stars now observed near and beyond the Solar Neighbourhood migrated from inner regions. Mock CoRoGEE observations of a chemo-dynamical Milky Way disc model show that the number of high-metallicity stars in the outer disc is too high to be accounted for even by the strong radial mixing present in the model. The mock observations also reveal that the age distribution of the [$\alpha$/Fe]-enhanced sequence in the CoRoGEE inner-disc field is much broader than expected from a combination of radial mixing and observational errors. We suggest that a thick disc/bulge component that formed stars for more than 3 Gyr may account for these discrepancies.

Recent improvements in the age dating of stellar populations and single stars allow us to study the ages and abundance of stars and galaxies with unprecedented accuracy. We here compare the relation between age and \alpha-element abundances for stars in the solar neighborhood to that of local, early-type galaxies. We find both relations to be very similar. Both fall into two regimes with a flat slope for ages younger than ~9 Gyr and a steeper slope for ages older than that value. This quantitative similarity seems surprising, given the different types of galaxies and scales involved. For the sample of early-type galaxies we also show that the data are inconsistent with literature delay time distributions of either single or double Gaussian shape. The data are consistent with a power law delay time distribution. We thus confirm that the delay time distribution inferred for the Milky Way from chemical evolution arguments also must apply to massive early-type galaxies. We also offer a tentative explanation for the seeming universality of the age-[\alpha/Fe] relation as the manifestation of averaging of different stellar populations with varying chemical evolution histories.

During the past decade, several studies reported a correlation between chemical abundances of stars and condensation temperature (also known as Tc trend). However, the real astrophysical nature of this correlation is still debated. The main goal of this work is to explore the possible dependence of the Tc trend on stellar Galactocentric distances, Rmean. We used high-quality spectra of about 40 stars observed with the HARPS and UVES spectrographs to derive precise stellar parameters, chemical abundances, and stellar ages. A differential line-by-line analysis was applied to achieve the highest possible precision in the chemical abundances. We confirm previous results that [X/Fe] abundance ratios depend on stellar age and that for a given age, some elements also show a dependence on Rmean. When using the whole sample of stars, we observe a weak hint that the Tc trend depends on Rmean. The observed dependence is very complex and disappears when only stars with similar ages are considered. To conclude on the possible dependence of the Tc trend on the formation place of stars, a larger sample of stars with very similar atmospheric parameters and stellar ages observed at different Galactocentric distances is needed

We use a high resolution cosmological zoom simulation of a Milky Way-sized halo to study the observable features in velocity and metallicity space associated with the dynamical influence of spiral arms. For the first time, we demonstrate that spiral arms, that form in a disc in a fully cosmological environment with realistic galaxy formation physics, drive large-scale systematic streaming motions. In particular, on the trailing edge of the spiral arms the peculiar galacto-centric radial and azimuthal velocity field is directed radially outward and azimuthally backward, whereas it is radially inward and azimuthally forward on the leading edge. Owing to the negative radial metallicity gradient, this systematic motion drives, at a given radius, an azimuthal variation in the residual metallicity that is characterised by a metal rich trailing edge and a metal poor leading edge. We show that these signatures are theoretically observable in external galaxies with Integral Field Unit instruments such as VLT/MUSE, and if detected, would provide evidence for large-scale systematic radial migration driven by spiral arms.

In a companion paper, we have presented the combined asteroseismic-spectroscopic dataset obtained from CoRoT lightcurves and APOGEE infra-red spectra for 678 solar-like oscillating red giants in two fields of the Galactic disc (CoRoGEE). We have measured chemical abundance patterns, distances, and ages of these field stars which are spread over a large radial range of the Milky Way's disc. Here we show how to simulate this dataset using a chemodynamical Galaxy model. We also demonstrate how the observation procedure influences the accuracy of our estimated ages.

The main goal of Galactic Archaeology is understanding the formation and evolution of the basic Galactic components. This requires sophisticated chemo-dynamical modeling, where disc asymmetries (e.g., perturbations from the central bar, spirals arms, and infalling satellites) and non-equilibrium processes are taken into account self-consistently. Here we discuss the current status of Galactic chemo-dynamical modeling and focus on a recent hybrid technique, which helps circumvent traditional problems with chemical enrichment and star formation encountered in fully self-consistent cosmological simulations. We show that this model can account for a number of chemo-kinematic relations in the Milky Way. In addition, we demonstrate that (1) our model matches well the observed age-[alpha/Fe] and age-[Fe/H] relations and (2) that the scatter in the age-[Fe/H] relation cannot be simply explained by blurring (stars on apo- and pericenters visiting the solar vicinity) but significant radial migration (stars born elsewhere but ending up at the solar vicinity today because of a change in guiding radius) is needed. We emphasize the importance of accurate stellar ages, such as those obtained through asteroseismology by the CoRoT and Kepler missions, for breaking the degeneracy among different Galactic evolution scenarios.

Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends. Methods. We analyse five simulated spiral galaxies with masses comparable to the Milky Way. The simulations span a range of star formation and energy feedback strengths and prescriptions, particle- and grid-based hydrodynamical implementations, as well as initial conditions/assembly history. Results. Consistently, we find that the steeper, negative, radial metallicity gradients seen in the mid-plane flatten with increasing height away from the plane. In simulations with stronger (and/or more spatially-extended) feedback, the negative radial gradients invert, becoming positive for heights in excess of 1 kpc. Such behaviour is consistent with that inferred from recent observations. Our measurements of the vertical metallicity gradients show no clear correlation with galactocentric radius, and are in good agreement with those observed in the Milky Way's thick disc (locally). Conclusions. Simulations employing stronger/more extended feedback prescriptions possess radial and vertical metallicity and kinematic gradients more in line with recent observations. The inverted, positive, radial metallicity gradients seen in the simulated thick stellar discs originate from a population of younger, more metal-rich, stars formed in-situ, superimposed upon a background population of older migrators from the inner disc; the contrast provided by the former increases radially, due to the inside-out growth of the disc. A similar behaviour may be responsible for the same flattening seen in the radial gradients with scaleheight in the Milky Way.

Aims: We study the effects of the cosmological assembly history on the chemical and dynamical properties of the discs of spiral galaxies as a function of radius. Methods: We make use of the simulated Milky-Way mass, fully-cosmological discs, from {\tt RaDES} (Ramses Disc Environment Study). We analyse their assembly history by examining the proximity of satellites to the galactic disc, instead of their merger trees, to better gauge which satellites impact the disc. We present stellar age and metallicity profiles, Age-Metallicity Relation (AMR), Age-Velocity dispersion Relation (AVR), and Stellar Age Distribution (SAD) in several radial bins for the simulated galaxies. Results: Assembly histories can be divided into three different stages: i) a merger dominated phase, when a large number of mergers with mass ratios of $\sim$1:1 take place (lasting $\sim$3.2$\pm$0.4 Gyr on average); ii) a quieter phase, when $\sim$1:10 mergers take place (lasting $\sim$4.4$\pm$2.0 Gyr) - these two phases are able to kinematically heat the disc and produce a disc that is chemically mixed over its entire radial extension; iii) a "secular" phase where the few mergers that take place have mass ratios below 1:100, and which do not affect the disc properties (lasting $\sim$5.5$\pm$2.0 Gyr). Phase ii ends with a final merger event (at time $t_\mathrm{jump}$) marking the onset of important radial differences in the AMR, AVR, and SAD. Conclusions: Inverted AMR trends in the outer parts of discs, for stars younger than $t_\mathrm{jump}$, are found as the combined effect of radial motions and star formation in satellites temporarily located in these outer parts. "U-shaped" stellar age profiles change to an old plateau ($\sim$10 Gyr) in the outer discs for the entire {\tt RaDES} sample. This shape is a consequence of inside-out growth of the disc, radial motions of disc stars ... [abridged]

We study the individual evolution histories of three nearby low-mass edge-on galaxies (IC 5052, NGC4244, and NGC5023). Using resolved stellar populations, we constructed star count density maps for populations of different ages and analyzed the change of structural parameters with stellar age within each galaxy. We do not detect a separate thick disk in any of the three galaxies, even though our observations cover a wider range in equivalent surface brightness than any integrated light study. While scale heights increase with age, each population can be well described by a single disk. Two of the galaxies contain a very weak additional component, which we identify as the faint halo. The mass of these faint halos is lower than 1% of the mass of the disk. The three galaxies show low vertical heating rates, which are much lower than the heating rate of the Milky Way. This indicates that heating agents, such as giant molecular clouds and spiral structure, are weak in low-mass galaxies. All populations in the three galaxies exhibit no or only little flaring. While this finding is consistent with previous integrated light studies, it poses strong constraints on galaxy simulations, where strong flaring is often found as a result of interactions or radial migration.

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

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

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

We estimate the size and distribution of the parent populations for the 6 largest (at least 20 stars in the Solar neighborhood) chemical groups identified in the Chemical Tagging experiment by Mitschang et al.~2014. Stars in the abundance groups tend to lie near a boundary in angular momentum versus eccentricity space where the probability is highest for a star to be found in the Solar neighborhood and where orbits have apocenter approximately equal to the Sun's galactocentric radius. Assuming that the parent populations are uniformly distributed at all azimuthal angles in the Galaxy, we estimate that the parent populations of these abundance groups contain at least 200,000 members. The spread in angular momentum of the groups implies that the assumption of a uniform azimuthal distribution only fails for the two youngest groups and only for the highest angular momentum stars in them. The parent populations of three thin disk groups have narrow angular momentum distributions, but tails in the eccentricity and angular momentum distributions suggest that only a small fraction of stars have migrated and increased in eccentricity. In contrast, the parent populations of the thick disk groups exhibit both wide angular momentum and eccentricity distributions implying that both heating and radial migration has taken place.

Recent spectroscopic observations in the Milky Way suggest that the chemically defined thick disk (stars with high [alpha/Fe] ratios and thus old) has a significantly smaller scale-length than the thin disk. This is in apparent contradiction with observations of external edge-on galaxies, where the thin and thick components have comparable scale-lengths. Moreover, while observed disks do not flare (scale-height does not increase with radius), numerical simulations suggest that disk flaring is unavoidable, resulting from both environmental effects and secular evolution. Here we address these problems by studying two different suites of simulated galactic disks formed in the cosmological context. We show that the scale-heights of coeval populations always increase with radius. However, the total population can be decomposed morphologically into thin and thick disks, which do not flare. We relate this to the disk inside-out formation, where younger populations have increasingly larger scale-lengths and flare at progressively larger radii. In this new picture, thick disks are composed of the imbedded flares of mono-age stellar populations. Assuming that disks form inside out, we predict that morphologically defined thick disks must show a decrease in age (or [alpha/Fe] ratios) with radius and that coeval populations should always flare. This also explains the observed inversion in the metallicity and [alpha/Fe] gradients for stars away from the disk midplane in the Milky Way. The results of this work are directly linked to, and can be seen as evidence of, inside-out disk growth.

We report the discovery of a group of apparently young CoRoT red-giant stars exhibiting enhanced [alpha/Fe] abundance ratios (as determined from APOGEE spectra) with respect to Solar values. Their existence is not explained by standard chemical evolution models of the Milky Way, and shows that the chemical-enrichment history of the Galactic disc is more complex. We find similar stars in previously published samples for which isochrone-ages could be robustly obtained, although in smaller relative numbers, which could explain why these stars have not received prior attention. The young [alpha/Fe]-rich stars are much more numerous in the CoRoT-APOGEE (CoRoGEE) inner-field sample than in any other high-resolution sample available at present, as only CoRoGEE can explore the inner-disc regions and provide ages for its field stars. The kinematic properties of the young [$\alpha$/Fe]-rich stars are not clearly thick-disc like, despite their rather large distances from the Galactic mid-plane. Our tentative interpretation of these and previous intriguing observations in the Milky Way is that these stars were formed close to the end of the Galactic bar, near corotation -- a region where gas can be kept inert for longer times, compared to other regions shocked more frequently by the passage of spiral arms. Moreover, that is where the mass return from older inner-disc stellar generations should be maximal (according to an inside-out disc-formation scenario), further diluting the in-situ gas. Other possibilities to explain these observations (e.g., a recent gas-accretion event) are also discussed.

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

In the last few decades many efforts have been made to understand the effect of spiral arms on the gas and stellar dynamics in the Milky Way disc. One of the fundamental parameters of the spiral structure is its angular velocity, or pattern speed $\Omega_p$, which determines the location of resonances in the disc and the spirals' radial extent. The most direct method for estimating the pattern speed relies on backward integration techniques, trying to locate the stellar birthplace of open clusters. Here we propose a new method based on the interaction between the spiral arms and the stars in the disc. Using a sample of around 500 open clusters from the {\it New Catalogue of Optically Visible Open Clusters and Candidates}, and a sample of 500 giant stars observed by APOGEE, we find $\Omega_p = 23.0\pm0.5$ km s$^{-1}$ kpc$^{-1}$, for a local standard of rest rotation $V_0=220$~km s$^{-1}$ and solar radius $R_0=8.0$~kpc. Exploring a range in $V_0$ and $R_0$ within the acceptable values, 200-240 km s$^{-1}$ and 7.5-8.5 kpc, respectively, results only in a small change in our estimate of $\Omega_p$, that is within the error. Our result is in close agreement with a number of studies which suggest values in the range 20-25 km s$^{-1}$ kpc$^{-1}$. An advantage of our method is that we do not need knowledge of the stellar age, unlike in the case of the birthplace method, which allows us to use data from large Galactic surveys. The precision of our method will be improved once larger samples of disk stars with spectroscopic information will become available thanks to future surveys such as 4MOST.

We have observed an N-body/Smoothed Particle Hydrodynamics simulation of a Milky Way like barred spiral galaxy. We present a simple method that samples N-body model particles into mock Gaia stellar observations and takes into account stellar populations, dust extinction and Gaia's science performance estimates. We examine the kinematics around a nearby spiral arm at a similar position to the Perseus arm at three lines of sight in the disc plane; (l,b)=(90,0), (120,0) and (150,0) degrees. We find that the structure of the peculiar kinematics around the co-rotating spiral arm, which is found in Kawata et al. (2014b), is still visible in the observational data expected to be produced by Gaia despite the dust extinction and expected observational errors of Gaia. These observable kinematic signatures will enable testing whether the Perseus arm of the Milky Way is similar to the co-rotating spiral arms commonly seen in N-body simulations.

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.

[Abridge] In the first paper of this series (paper I) we presented a new approach for studying the chemo-odynamical evolution in disk galaxies, focusing on the Milky Way. Here we extend these results to different distances from the Galactic center, looking for variations of observables that can be related to on-going and future spectroscopic surveys. By separating the effects of kinematic heating and radial migration, we show that migration is much more important, even for the oldest and hottest stellar population. The distributions of stellar birth guiding radii and final guiding radii (signifying contamination from migration and heating, respectively) widen with increasing distance from the Galactic center. As a result, the slope in the age-metallicity relation flattens significantly at Galactic radii larger than solar. The radial metallicity and [Mg/Fe] gradients in our model show significant variations with height above the plane due to changes in the mixture of stellar ages. An inversion in the radial metallicity gradient is found from negative to weakly positive (at r<10 kpc), and from positive to negative for the [Mg/Fe] gradient, with increasing distance from the disk plane. We relate this to the combined effect of (i) the predominance of young stars close to the disk plane and old stars away from it, (ii) the more concentrated older stellar component, and (iii) the flaring of mono-age disk populations. We also investigate the effect of recycled gas flows on the mean [Fe/H] and find that in the region 4<r<12 kpc the introduced errors are less than 0.05-0.1 dex, related to the fact that inward and outward flows mostly cancel in that radial range. We show that radial migration cannot compete with the inside-out formation of the disk, exposed by the more centrally concentrated older disk populations, and consistent with recent observations.

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

Galactic Archaeology, i.e. the use of chemo-dynamical information for stellar samples covering large portions of the Milky Way to infer the dominant processes involved in its formation and evolution, is now a powerful method thanks to the large recently completed and ongoing spectroscopic surveys. It is now important to ask the right questions when analyzing and interpreting the information contained in these rich datasets. To this aim, we have developed a chemodynamical model for the Milky Way that provides quantitative predictions to be compared with the chemo-kinematical properties extracted from the stellar spectra. Three key parameters are needed to make the comparison between data and model predictions useful in order to advance in the field, namely: precise proper-motions, distances and ages. The uncertainties involved in the estimate of ages and distances for field stars are currently the main obstacles in the Galactic Archaeology method. Two important developments might change this situation in the near future: asteroseismology and the now launched Gaia. When combined with the large datasets from surveys like RAVE, SEGUE, LAMOST, Gaia-ESO, APOGEE , HERMES and the future 4MOST we will have the basic ingredients for the reconstruction of the MW history in hands. In the light of these observational advances, the development of detailed chemo-dynamical models tailored to the Milky Way is urgently needed in the field. Here we show the steps we have taken, both in terms of data analysis and modelling. The examples shown here illustrate how powerful can the Galactic Archaeology method become once ages and distances are known with better precision than what is currently feasible.

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 study seven simulated disc galaxies, three with a quiescent merger history, and four with mergers in their last 9 Gyr of evolution. We compare their structure at z=0 by decomposing them into "mono-age populations" (MAPs) of stars within 500 Myr age bins. All studied galaxies undergo a phase of merging activity at high redshift, so that stars older than 9 Gyr are found in a centrally concentrated component, while younger stars are mostly found in discs. We find that most MAPs have simple exponential radial and vertical density profiles, with a scale-height that typically increases with age. Because a large range of merger histories can create populations with simple structures, this suggests that the simplicity of the structure of mono-abundance populations observed in the Milky Way by Bovy et al. (2012b,c) is not necessarily a direct indicator of a quiescent history for the Milky Way. Similarly, the anti-correlation between scale-length and scale-height does not necessarily imply a merger-free history. However, mergers produce discontinuities between thin and thick disc components, and jumps in the age-velocity relation. The absence of a structural discontinuity between thin and thick disc observed in the Milky Way would seem to be a good indicator that no merger with a mass ratio larger than 1:15-1:10 occurred in the last 9 Gyr. Mergers at higher redshift might nevertheless be necessary to produce the thickest, hottest components of the Milky Way's disc.

We study the relation between stellar ages and vertical velocity dispersion (the age-velocity relation, or AVR) in a sample of seven simulated disc galaxies. In our simulations, the shape of the AVR for stars younger than 9 Gyr depends strongly on the merger history at low redshift, with even 1:10 - 1:15 mergers being able to create jumps in the AVR (although these jumps might not be detectable if the errors on stellar ages are on the order of 30%). For galaxies with a quiescent history at low redshift, we find that the vertical velocity dispersion rises smoothly for ages up to 8-9 Gyr, following a power law with a slope of ~0.5, similar to what is observed in the solar neighbourhood by the Geneva-Copenhagen Survey. For these galaxies, we show that the slope of the AVR is not imprinted at birth, but is the result of subsequent heating. By contrast, in all our simulations, the oldest stars form a significantly different population, with a high velocity dispersion. These stars are usually born kinematically hot in a turbulent phase of intense mergers at high redshift, and also include some stars accreted from satellites. This maximum in velocity dispersion is strongly decreased when age errors are included, suggesting that observations can easily miss such a jump with the current accuracy of age measurements.

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

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.

Despite the recent advancements in the field of galaxy formation and evolution, fully self-consistent simulations are still unable to make the detailed predictions necessary for the planned and ongoing large spectroscopic and photometry surveys of the Milky Way disc. These difficulties arise from the very uncertain nature of sub-grid physical energy feedback within models, affecting both star formation rates and chemical enrichment. To avoid these problems, we have introduced a new approach which consists of fusing disc chemical evolution models with compatible numerical simulations. We demonstrate the power of this method by showing that a range of observational results can be explained by our new model. We show that due to radial migration from mergers at high redshift and the central bar at later times, a sizable fraction of old metal-poor, high-[alpha/Fe] stars can reach the solar vicinity. This naturally accounts for a number of recent observations related to both the thin and thick discs, despite the fact that we use thin-disc chemistry only. Within the framework of our model, the MW thick disc has emerged naturally from (i) stars born with high velocity dispersions at high redshift, (ii) stars migrating from the inner disc very early on due to strong merger activity, and (iii) further radial migration driven by the bar and spirals at later times. A significant fraction of old stars with thick-disc characteristics could have been born near the solar radius.

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.

[Abridged] In this first paper of this series, we present a new approach for studying the chemo-dynamical evolution in disk galaxies, which consists of fusing disk chemical evolution models with compatible numerical simulations of galactic disks. This method avoids known star formation and chemical enrichment problems encountered in simulations. Here we focus on the Milky Way, by using a detailed thin-disk chemical evolution model and a simulation in the cosmological context, with dynamical properties close to those of our Galaxy. We show that, due to radial migration from mergers at high redshift and the central bar at later times, a sizable fraction of old metal-poor high-[alpha/Fe] stars reaches the solar vicinity. This naturally accounts for a number of recent observations related to both the thin and thick disks, despite the fact that we use thin-disk chemistry only. Although significant radial mixing is present, the slope in the age-metallicity relation is only weakly affected, with a scatter compatible with recent observational work. While we find a smooth density distribution in the [O/Fe]-[Fe/H] plane, we can recover the observed discontinuity by selecting particles according to kinematic criteria used in high-resolution samples to define the thin and thick disks. We outline a new method for estimating the birth place of the Sun and predict that the most likely radius lies in the range 4.4 < r < 7.7 kpc (for a location at r = 8 kpc). A new, unifying model for the Milky Way thick disk is offered, where both mergers and radial migration play a role at different stages of the disk evolution. We show that in the absence of early-on massive mergers the vertical velocity dispersion of the oldest stars is underestimated by a factor of ~2 compared with observations. We can, therefore, argue that the Milky Way thick disk is unlikely to have been formed through a quiescent disk evolution.

We explore a second order Hamiltonian vertical resonance model for X-shaped or peanut-shaped galactic bulges. The X-shape is caused by the 2:1 vertical Lindblad resonance with the bar, with two vertical oscillation periods per orbital period in the bar frame. We examine N-body simulations and find that due to the bar slowing down and disk thickening during bar buckling, the resonance and associated peanut-shape moves outward. The peanut-shape is consistent with the location of the vertical resonance, independent of whether the bar buckled or not. We estimate the resonance width from the potential m=4 Fourier component and find that the resonance is narrow, affecting orbits over a narrow range in the angular momentum distribution, dL/L ~ 0.05. As the resonance moves outward, stars originally in the mid plane are forced out of the mid plane into orbits just within the resonance separatrix. The height of the separatrix orbits, estimated from the Hamiltonian model, is approximately consistent with the peanut-shape height. The X-shape is comprised of stars in the vicinity of the resonance separatrix. The velocity distributions from the simulations illustrate that low inclination orbits are depleted within resonance. Within resonance, the vertical velocity distribution is broad, consistent with resonant heating caused by the passage of the resonance through the disk. In the Milky Way bulge we relate the azimuthally averaged mid-plane mass density near the vertical resonance to the rotation curve and bar pattern speed. At an estimated vertical resonance galactocentric radius of ~1.3 kpc, we confirm a mid-plane density of ~5x10^8 Msol/kpc^3, consistent with recently estimated mass distributions. We find that the rotation curve, bar pattern speed, 2:1 vertical resonance location, X-shape tips, and mid-plane mass density, are all self-consistent in the Milky Way galaxy bulge.

By selecting in the RAVE-DR4 survey the stars located between 1 and 2 kpc above the Galactic plane, we question the consistency of the simplest three-component model (thin disc, thick disc, halo) for the Milky Way. We confirm that the metallicity and azimuthal velocity distribution functions of the thick disc are not Gaussian. In particular, we find that the thick disc has an extended metallicity tail going at least down to [M/H]=-2 dex, contributing roughly at 3% of the entire thick disc population and having a shorter scale-length compared to the canonical thick disc. The mean azimuthal velocity of these metal-poor stars allows us to estimate the correlation between the metallicity and the orbital velocity, which is an important constraint on the formation mechanisms of the Galactic thick disc. Given our simple approach, we find dVphi/d[M/H] ~ 50 km/s/dex, which is in a very good agreement with previous literature values. We complete the study with a brief discussion on the implications of the formation scenarios for the thick disc, and suggest that given the above mentioned characteristics, a thick disc mainly formed by radial migration mechanisms seems unlikely.

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.

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 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.

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.

Recently, Widrow and collaborators announced the discovery of vertical density waves in the Milky Way disk. Here we investigate a scenario where these waves were induced by the Sagittarius dwarf galaxy as it plunged through the Galaxy. Using numerical simulations, we find that the Sagittarius impact produces North-South asymmetries and vertical wave-like behavior that qualitatively agrees with what is observed. The extent to which vertical modes can radially penetrate into the disc, as well as their amplitudes, depend on the mass of the perturbing satellite. We show that the mean height of the disc is expected to vary more rapidly in the radial than in the azimuthal direction. If the observed vertical density asymmetry is indeed caused by vertical oscillations, we predict radial and azimuthal variations of the mean vertical velocity, correlating with the spatial structure. These variations can have amplitudes as large as 8 km/s.

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.

Non-axisymmetric components, such as spirals and central bars, play a major role in shaping galactic discs. An important aspect of the disc secular evolution driven by these perturbers is the radial migration of stars. It has been suggested recently that migration can populate a thick-disc component from inner-disc stars with high vertical energies. Since this has never been demonstrated in simulations, we study in detail the effect of radial migration on the disc velocity dispersion and disc thickness, by separating simulated stars into migrators and non-migrators. We apply this method to three isolated barred Tree-SPH N-body galaxies with strong radial migration. Contrary to expectations, we find that as stellar samples migrate, on the average, their velocity dispersion change (by as much as 50%) in such a way as to approximately match the non-migrating population at the radius at which they arrive. We show that, in fact, migrators suppress heating in parts of the disc. To confirm the validity of our findings, we also apply our technique to three cosmological re-simulations, which use a completely different simulation scheme and, remarkably, find very similar results. We believe the inability of migration to thicken discs is a fundamental property of internal disc evolution, irrespective of the migration mechanism at work. We explain this with the approximate conservation of the (average) vertical and radial actions rather than the energy. This "action mixing" can be used to constrain the migration rate in the Milky Way: estimates of the average vertical action in observations for different populations of stars should reveal flattening with radius for older groups of stars.

We present the metallicity distribution function (MDF) for 24,270 G and 16,847 K dwarfs at distances from 0.2 to 2.3 kpc from the Galactic plane, based on spectroscopy from the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey. This stellar sample is significantly larger in both number and volume than previous spectroscopic analyses, which were limited to the solar vicinity, making it ideal for comparison with local volume-limited samples and Galactic models. For the first time, we have corrected the MDF for the various observational biases introduced by the SEGUE target selection strategy. The SEGUE sample is particularly notable for K dwarfs, which are too faint to examine spectroscopically far from the solar neighborhood. The MDF of both spectral types becomes more metal-poor with increasing |Z|, which reflects the transition from a sample with small [alpha/Fe] values at small heights to one with enhanced [alpha/Fe] above 1 kpc. Comparison of our SEGUE distributions to those of two different Milky Way models reveals that both are more metal-rich than our observed distributions at all heights above the plane. Our unbiased observations of G and K dwarfs provide valuable constraints over the |Z|-height range of the Milky Way disk for chemical and dynamical Galaxy evolution models, previously only calibrated to the solar neighborhood, with particular utility for thin- and thick-disk formation models.

We investigate the evolution of galactic disks in N-body Tree-SPH simulations. We find that disks, initially truncated at three scale-lengths, can triple their radial extent, solely driven by secular evolution. Both Type I (single exponential) and Type II (down-turning) observed disk surface-brightness profiles can be explained by our findings. We relate these results to the strong angular momentum outward transfer, resulting from torques and radial migration associated with multiple patterns, such as central bars and spiral waves of different multiplicity. We show that even for stars ending up on cold orbits, the changes in angular momentum exhibit complex structure as a function of radius, unlike the expected effect of transient spirals alone. Focussing on one of our models, we find evidence for non-linear coupling among m=1, 2, 3 and 4 density waves, where m is the pattern multiplicity. We suggest that the naturally occurring larger resonance widths at galactic radii beyond four scale-lengths may have profound consequences on the formation and location of breaks in disk density profiles, provided spirals are present at such large distances. We also consider the effect of gas inflow and show that when in-plane smooth gas accretion of ~5 M_sun/yr is included, the outer disks become more unstable, leading to a strong increase in the stellar velocity dispersion. This, in turn, causes the formation of a Type III (up-turning) profile in the old stellar population. We propose that observations of Type III surface brightness profiles, combined with an up-turn in the stellar velocity dispersions beyond the disk break, could be a signature of ongoing gas-accretion. The results of this study suggest that disk outskirts comprised of stars migrated from the inner disk would have relatively large radial velocity dispersions, and significant thickness when seen edge-on. [Abridged]

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.

It is now known that minor mergers are capable of creating structure in the phase-space distribution of their host galaxy's disc. In order to search for such imprints in the Milky Way, we analyse the SEGUE F/G-dwarf and the Schuster et al. (2006) stellar samples. We find similar features in these two completely independent stellar samples, consistent with the predictions of a Milky Way minor-merger event. We next apply the same analyses to high-resolution, idealised N-body simulations of the interaction between the Sagittarius dwarf galaxy and the Milky Way. The energy distributions of stellar particle samples in small spatial regions in the host disc reveal strong variations of structure with position. We find good matches to the observations for models with a mass of Sagittarius' dark matter halo progenitor $\lessapprox 10^{11}$ M$_{\odot}$. Thus, we show that this kind of analysis could be used to provide unprecedentedly tight constraints on Sagittarius' orbital parameters, as well as place a lower limit on its mass.

NGC 6791 is a unique stellar system among Galactic open clusters being at the same time one of the oldest open clusters and the most metal rich. Combination of its properties is puzzling and poses question of its origin. One possible scenario is that the cluster formed close to the Galactic Center and later migrated outwards to its current location. In this work we study the cluster's orbit and investigate the possible migration processes which might have displaced NGC 6791 to its present-day position, under the assumption that it actually formed in the inner disk. To this aim we performed integrations of NGC 6791's orbit in a potential consistent with the main Milky Way parameters. In addition to analytical expressions for halo, bulge and disk, we also consider the effect of bar and spiral arm perturbations, which are expected to be very important for the disk dynamical evolution, especially inside the solar circle. Starting from state-of-the art initial conditions for NGC 6791, we calculate 1000 orbits back in time for about 1 Gyr turning on and off different non-axisymmetric components of the global potential. We then compare statistical estimates of the cluster's recent orbital parameters with the orbital parameters of 10^4 test-particles originating close to the Galactic Center (having initial galocentric radii in the range of 3-5 kpc) and undergoing radial migration during 8 Gyr of forward integration. We find that a model which incorporates a strong bar and spiral arm perturbations can indeed be responsible for the migration of NGC 6791 from the inner disk (galocentric radii of 3-5 kpc) to its present-day location. Such a model can provide orbital parameters which are close enough to the observed ones. However, the probability of this scenario as it results from our investigations is very low.

By means of N-body simulations we study the response of a galactic disc to a minor merger event. We find that non-self-gravitating, spiral-like features are induced in the thick disc. As we have shown in a previous work, this "ringing" also leaves an imprint in velocity space (the u-v plane) in small spatial regions, such as the solar neighbourhood. As the disc relaxes after the event, clumps in the u-v plane get closer with time, allowing us to estimate the time of impact. In addition to confirming the possibility of this diagnostic, here we show that in a more realistic scenario, the in-fall trajectory of the perturber gives rise to an azimuthal dependence of the structure in phase-space. We also find that the space defined by the energy and angular momentum of stars is a better choice than velocity space, as clumps remain visible even in large local volumes. This makes their observational detection much easier since one need not be restricted to a small spatial volume. We show that information about the time of impact, the mass of the perturber, and its trajectory is stored in the kinematics of disc stars.

Disc non-axisymmetrc components, such as spirals and central bars, are nowadays known to play an important role in shaping galactic discs. Here we use Tree-SPH N-body simulations to examine the effect of these perturbers on two aspects: the occurrence of multiple patterns in discs and the effects of radial migration on disc thickening. We find that, in addition to a central bar, multiple spiral patterns and lopsided modes develop in all models. Interaction among these asymmetric features results in a large scale stellar migration. However, we show that, despite the strong radial mixing, discs cannot be thickened sufficiently to match observed thick discs. We relate this to the adiabatic cooling as stars migrate radially outwards. We also find that the bulge contribution to a thick-disc component for an Sa-type galaxy at ~2.5 disc scale-lengths is less than 1% and zero in the case of a Milky Way-like, Sb-type. Our findings cast doubt on the plausibility of thick disc formation via stellar radial migration.

Bars and spirals are among the main drivers of the secular evolution of galactic disks, and it is, thus, of prime importance to better understand their exact nature and dynamics. In this respect, the Milky Way provides a unique laboratory in which to study a snapshot of their detailed dynamical influence in six-dimensional phase-space. With the advent of present and future astrometric and spectroscopic surveys, such six-dimensional phase-space data can be obtained for stars in an increasingly large volume around the Sun. Here, we review the signatures of - and constraints on - disk non-axisymmetries obtained from recent stellar kinematical data sets, including (i) the detection of resonant moving groups in the solar neighbourhood, (ii) the non-zero value of the Oort constants C and K, and (iii) the detection of a radial velocity gradient of 4 km/s/kpc in the extended solar neighbourhood (d<2 kpc).

An N-body hybrid simulation, integrating both massive and tracer particles, of a Galactic disk is used to study the stellar phase space distribution or velocity distributions in different local neighborhoods. Pattern speeds identified in Fourier spectrograms suggest that two-armed and three-armed spiral density waves, a bar and a lopsided motion are coupled in this simulation, with resonances of one pattern lying near resonances of other patterns. We construct radial and tangential (uv) velocity distributions from particles in different local neighborhoods. More than one clump is common in these local velocity distributions regardless of the position in the disk. Features in the velocity distribution observed at one galactic radius are also seen in nearby neighborhoods (at larger and smaller radii) but with shifted mean v values. This is expected if the v velocity component of a clump sets the mean orbital galactic radius of its stars. We find that gaps in the velocity distribution are associated with the radii of kinks or discontinuities in the spiral arms. These gaps also seem to be associated with Lindblad resonances with spiral density waves and so denote boundaries between different dominant patterns in the disk. We discuss implications for interpretations of the Milky Way disk based on local velocity distributions. Velocity distributions created from regions just outside the bar's Outer Lindblad resonance and with the bar oriented at 45 degrees from the Sun-Galactic center line more closely resemble that seen in the solar neighborhood (triangular in shape at lower uv and with a Hercules like stream) when there is a strong nearby spiral arm, consistent with the observed Centaurus Arm tangent, just interior to the solar neighborhood.

The Hyades stream has long been thought to be a dispersed vestige of the Hyades cluster. However, recent analyses of the parallax distribution, of the mass function, and of the action-space distribution of stream stars have shown it to be rather composed of orbits trapped at a resonance of a density disturbance. This resonant scenario should leave a clearly different signature in the element abundances of stream stars than the dispersed cluster scenario, since the Hyades cluster is chemically homogeneous. Here, we study the metalllicty as well as the element abundances of Li, Na, Mg, Fe, Zr, Ba, La, Ce, Nd, and Eu for a random sample of stars belonging to the Hyades stream, and compare them with those of stars from the Hyades cluster. From this analysis: (i) we independently confirm that the Hyades stream cannot be solely composed of stars originating in the Hyades cluster; (ii) we show that some stars from the Hyades stream nevertheless have abundances compatible with an origin in the cluster; (iii) we emphasize that the use of Li as a chemical tag of the cluster origin of main-sequence stars is very efficient in the range 5500 - 6200 K, since the Li sequence in the Hyades cluster is very tight; (iv) we show that, while this evaporated population has a metallicity excess of ~0.2 dex w.r.t. the local thin disk population, identical to that of the Hyades cluster, the remainder of the Hyades stream population has still a metallicity excess of ~0.06 to 0.15 dex, consistent with an origin in the inner Galaxy; (v) we show that the Hyades stream can be interpreted as an inner 4:1 resonance of the spiral pattern: this then also reproduces an orbital family compatible with the Sirius stream, and places the origin of the Hyades stream up to 1 kpc inwards from the solar radius, which might explain the observed metallicity excess of the stream population.

We have recently identified a new radial migration mechanism resulting from the overlap of spiral and bar resonances in galactic disks. Here we confirm the efficiency of this mechanism in fully self-consistent, Tree-SPH simulations, as well as high-resolution pure N-body simulations. In all barred cases we clearly identify the effect of spiral-bar resonance overlap by measuring a bimodality in the changes of angular momentum in the disk, dL, whose maxima are near the bar's corotation and outer Lindblad resonance. This contrasts with the smooth distribution of dL for a simulation with no stable bar present, where strong radial migration is induced by multiple spirals. The presence of a disk gaseous component appears to increase the rate of angular momentum exchange by about 20%. The efficiency of this mechanism is such that galactic stellar disks can extend to over 10 scale-lengths within 1-3 Gyr in both Milky Way size and low-mass galaxies (circular velocity ~100 km/s). We also show that metallicity gradients can flatten in less than 1 Gyr rendering mixing in barred galaxies an order of magnitude more efficient than previously thought.

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.

We have recently identified a previously unknown radial migration mechanism resulting from the overlap of spiral and bar resonances in galactic discs (Minchev & Famaey 2010, Minchev et al. 2010). This new mechanism is much more efficient than mixing by transient spirals and its presence is unavoidable in all barred galaxies, such as our own Milky Way. The consequences of this are a strong flattening in the metallicity gradient in the disc, an extended disc profile, and the formation of a thick disc component, all taking place in only a couple of Gyr. This timescale is drastically shorter than previously expected and thus can put strong constraints on the longevity, strength and pattern speeds of the Galactic bar and Spiral Structure.

While it has long been known that a large number of short-lived transient spirals can cause stellar migration, here we report that another mechanism is also effective at mixing disks of barred galaxies. The resonance overlap of the bar and spiral structure induces a nonlinear response leading to a strong redistribution of angular momentum in the disk. We find that, depending on the amplitudes of the perturbers, the changes in angular momentum, dL, could occur up to an order of magnitude faster than in the case of recurrent spirals. The signature of this mechanism is a bimodality in dL with maxima near the bar's corotation and its outer Lindblad resonance; this is independent of the properties of the spiral structure. For parameters consistent with the Milky Way the disk mixes in about 3 Gyr and the stellar velocity dispersion increases with time in a manner roughly consistent with observations. This new mechanism could account for both the observed age-velocity relation and the absence of age-metallicity relation in the solar neighborhood. Spiral-bar interaction could also explain observations showing that strongly barred galaxies have weaker metallicity gradients than weakly barred or non-barred galaxies.

We find that a steady state bar induces transient features at low velocities in the solar neighborhood velocity distribution due to the initial response of the disc, following the formation of the bar. We associate these velocity streams with two quasi-periodic orbital families, librating around the stable x_1(1) and x_1(2) orbits near the bar's outer Lindblad resonance (OLR). In a reference frame moving with the bar, these otherwise stationary orbits precess on a timescale dependent on the strength of the bar, consistent with predictions from a simple Hamiltonian model for the resonance. This behavior allows the two orbital families to reach the solar neighborhood and manifest themselves as clumps in the u-v plane moving away from (x_1(2)), and toward (x_1(1)) the Galactic center. Depending on the bar parameters and time since its formation, this model is consistent with the Pleiades and Coma Berenices, or Pleiades and Sirius moving groups seen in the Hipparcos stellar velocity distribution, if the Milky Way bar angle is 30<phi_0<45 [deg] and its pattern speed is Omega_b/Omega_0=1.82\pm 0.07, where Omega_0 is the angular velocity of the local standard of rest (LSR). Since the process is recurrent, we can achieve a good match about every six LSR rotations. However, to be consistent with the fraction of stars in the Pleiades, we estimate that the Milky Way bar formed ~2 Gyr ago. This model argues against a common dynamical origin for the Hyades and Pleiades moving groups.

Previous work has related the Galactic bar to structure in the local stellar velocity distribution. Here we show that the bar also influences the spatial gradients of the velocity vector via the Oort constants. By numerical integration of test-particles we simulate measurements of the Oort C value in a gravitational potential including the Galactic bar. We account for the observed trend that C is increasingly negative for stars with higher velocity dispersion. By comparing measurements of C with our simulations we improve on previous models of the bar, estimating that the bar pattern speed is Omega_b/Omega_0=1.87\pm0.04, where Omega_0 is the local circular frequency, and the bar angle lies within 20[deg]<phi_0<45[deg]. We find that the Galactic bar affects measurements of the Oort constants A and B less than ~2 km/s/kpc for the hot stars.

Using test particle simulations we examine the structure of the outer Galactic disk as it is perturbed by a satellite in a tight eccentric orbit about the Galaxy. A satellite of mass a few times 10^9 Msol can heat the outer Galactic disk, excite spiral structure and a warp and induce streams in the velocity distribution. We examine particle eccentricity versus the change in mean radius between initial and current orbits. Correlations between these quantities are reduced after a few satellite pericenter passages. Stars born in the outer galaxy can be moved in radius from their birth positions and be placed in low eccentricity orbits inside their birth radii. We propose that mergers and perturbations from satellite galaxies and subhalos can induce radial mixing in the stellar metallicity distribution.

We perform numerical simulations of a stellar galactic disk with initial conditions chosen to represent an unrelaxed population which might have been left following a merger. Stars are unevenly distributed in radial action angle, though the disk is axisymmetric. The velocity distribution in the simulated Solar neighborhood exhibits waves traveling in the direction of positive v, where u,v are the radial and tangential velocity components. As the system relaxes and structure wraps in phase space, the features seen in the uv-plane move closer together. We show that these results can be obtained also by a semi-analytical method. We propose that this model could provide an explanation for the high velocity streams seen in the Solar neighborhood at approximate v in km/s, of -60 (HR 1614), -80 (Arifyanto and Fuchs 2006), -100 (Arcturus), and -160 (Klement et al. 2008). In addition, we predict four new features at v ~ -140, -120, 40 and 60 km/s. By matching the number and positions of the observed streams, we estimate that the Milky Way disk was strongly perturbed ~1.9 Gyr ago. This event could have been associated with Galactic bar formation.

The morphology of the outer rings of early-type spiral galaxies is compared to integrations of massless collisionless particles initially in nearly circular orbits. Particles are perturbed by a quadrupolar gravitational potential corresponding to a growing and secularly evolving bar. We find that outer rings with R1R2 morphology and pseudorings are exhibited by the simulations even though they lack gaseous dissipation. Simulations with stronger bars form pseudorings earlier and more quickly than those with weaker bars. We find that the R1 ring, perpendicular to the bar, is fragile and dissolves after a few bar rotation periods if the bar pattern speed increases by more than ~ 8%, bar strength increases (by >~ 140%) after bar growth, or the bar is too strong (Q_T>0.3). If the bar slows down after formation, pseudoring morphology persists and the R2 ring perpendicular to the bar is populated due to resonance capture. The R2 ring remains misaligned with the bar and increases in ellipticity as the bar slows down. The R2 ring becomes scalloped and does not resemble any ringed galaxies if the bar slows down more than 3.5% suggesting that bars decrease in strength before they slow down this much. We compare the morphology of our simulations to B-band images of 9 ringed galaxies from the Ohio State University Bright Spiral Galaxy Survey, and we find a reasonable match in morphologies to R1R2' pseudorings seen within a few bar rotation periods of bar formation. Some of the features previously interpreted in terms of dissipative models may be due to transient structure associated with recent bar growth and evolution.

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.

HD 179821 is an evolved star of unknown progenitor mass range (either post-Asymptotic Giant Branch or post-Red Supergiant) exhibiting a double peaked spectral energy distribution (SED) with a sharp rise from $\sim8-20$ $\mu$m. Such features have been associated with ejected dust shells or inwardly truncated circumstellar discs. In order to compare SEDs from both systems, we employ a spherically symmetric radiative transfer code and compare it to a radiative, inwardly truncated disc code. As a case study, we model the broad-band SED of HD 179821 using both codes. Shortward of 40 $\mu$m, we find that both models produce equivalent fits to the data. However, longward of 40 $\mu$m, the radial density distribution and corresponding broad range of disc temperatures produce excess emission above our spherically symmetric solutions and the observations. For HD 179821, our best fit consists of a $T_{eff}=7000$ K central source characterized by $\tau_V\sim1.95$ and surrounded by a radiatively driven, spherically symmetric dust shell. The extinction of the central source reddens the broad-band colours so that they resemble a $T_{eff}=5750$ K photosphere. We believe that HD 179821 contains a hotter central star than previously thought. Our results provide an initial step towards a technique to distinguish geometric differences from spectral modeling.

Previous work has related the Galactic Bar to structure in the local stellar velocity distribution. Here we show that the Bar also influences the spatial gradients of the velocity vector via the Oort constants. By numerical integration of test-particles we simulate measurements of the Oort C value in a gravitational potential including the Galactic Bar. We account for the observed trend that C is increasingly negative for stars with higher velocity dispersion. By comparing measurements of C with our simulations we improve on previous models of the Bar, estimating that the Bar pattern speed is Omega_b/Omega_0=1.87\pm0.04, where Omega_0 is the local circular frequency, and the Bar angle lies within 20<phi_0<45 deg. We find that the Galactic Bar affects measurements of the Oort constants A and B less than ~2 km/s/kpc for the hot stars.

We investigate the effect of spiral structure on the Galactic disk as viewed by pencil beams centered on the Sun, relevant to upcoming surveys such as ARGOS, SEGUE, and GAIA. We create synthetic Galactic maps which we call Pencil Beam Maps (PBMs) of the following observables: line-of-sight velocities, the corresponding velocity dispersion, and the stellar number density that are functions of distance from the observer. We show that such maps can be used to infer spiral structure parameters, such as pattern speed, solar phase angle, and number of arms. The mean line-of-sight velocity and velocity dispersion are affected by up to ~35 km/s which is well within the detectable limit for forthcoming radial velocity surveys. One can measure the pattern speed by searching for imprints of resonances. In the case of a two-armed spiral structure it can be inferred from the radius of a high velocity dispersion ring situated at the 2:1 ILR. This information, however, must be combined with information related to the velocities and stellar number density in order to distinguish from a four-armed structure. If the pattern speed is such that the 2:1 ILR is hidden inside the Galactic bulge the 2:1 OLR will be present in the outer Galaxy and thus can equivalently be used to estimate the pattern speed. Once the pattern speed is known the solar angle can be estimated from the line-of-sight velocities and the number density PBMs. Forthcoming radial velocity surveys are likely to provide powerful constraints of the structure of the Milky Way disk.

We perform test-particle simulations in a 2D, differentially rotating stellar disk, subjected to a two-armed steady state spiral density wave perturbation in order to estimate the influence of spiral structure on the local velocity field. By using Levenberg-Marquardt least-squares fit we decompose the local velocity field (as a result of our simulations) into Fourier components to fourth order. Thus we obtain simulated measurements of the Oort constants, A, B, C, and K. We get relations between the Fourier coefficients and some galactic parameters, such as the phase angle of the Solar neighborhood and the spiral pattern speed. We show that systematic errors due to the presence of spiral structure are likely to affect the measurements of the Oort constants. Moderate strength spiral structure causes errors of order 5 km/s/kpc in A and B. We find variations of the Fourier coefficients with velocity dispersion, pattern speed, and sample depth. For a sample at an average heliocentric distance of 0.8 kpc we can summarize our findings as follows:(i) if our location in the Galaxy is near corotation then we expect a vanishing value for C for all phase angles;(ii) for a hot disk, spiral structure induced errors for all Oort constants vanish at, and just inward of the corotation radius; (iii) as one approaches the 4:1 indblad resonances |C| increases and so does its variation with galactic azimuth; (iv) for all simulations |C|, on average, is larger for lower stellar velocity dispersions, contrary to recent measurements.

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 consider a differentially rotating, 2D stellar disk perturbed by two steady state spiral density waves moving at different patterns speeds. Our investigation is based on direct numerical integration of initially circular test-particle orbits. We examine a range of spiral strengths and spiral speeds and show that stars in this time dependent gravitational field can be heated (their random motions increased).This is particularly noticeable in the simultaneous propagation of a 2-armed spiral density wave near the corotation resonance (CR), and a weak 4-armed one near the inner and outer 4:1 Lindblad resonances. In simulations with 2 spiral waves moving at different pattern speeds we find: (1) the variance of the radial velocity, sigma_R^2, exceeds the sum of the variances measured from simulations with each individual pattern; (2) sigma_R^2 can grow with time throughout the entire simulation; (3) sigma_R^2 is increased over a wider range of radii compared to that seen with one spiral pattern; (4) particles diffuse radially in real space whereas they don't when only one spiral density wave is present. Near the CR with the stronger, 2-armed pattern, test particles are observed to migrate radially. These effects take place at or near resonances of both spirals so we interpret them as the result of stochastic motions. This provides a possible new mechanism for increasing the stellar velocity dispersion in galactic disks. If multiple spiral patterns are present in the Galaxy we predict that there should be large variations in the stellar velocity dispersion as a function of radius.

With 2D hydrodynamical simulations of disks perturbed externally by stars, brown dwarfs or planets we investigate possible scenarios that can account for the spiral structure in circumstellar disks. We consider two scenarios, spiral structure driven by an external bound planet or low mass star and that excited by a previous stellar close encounter or flyby. We find that both scenarios produce morphology similar to that observed in the outer disks of HD 141569A and HD 100546; moderately open 2-armed outer spiral structure. The outer two-armed spiral structure observed in the disk of HD 141569A is qualitatively reproduced with tidal perturbations from its companion binary HD 141569B,C on a prograde orbit near periapse. Our simulation accounts for the outer spiral arms, but is less successful than the secular model of Augereau and Papaloizou at matching the lopsidedness or asymmetry of the disk edge at 300AU. The disk has been previously truncated by the tidal force from the binary. A bound object (stellar or planetary) is unlikely to explain the spiral structure in HD 100546. A co-eval planet or brown dwarf in the disk of sufficient mass to account for the amplitude of the spiral structure would be detectable in NICMOS and STIS images, however existing images reveal no such object. A previous encounter could explain the observed structure, provided that the encounter occurred less than a few thousand year ago. The object responsible for causing the spiral structure must then be within a few arcminutes of the star. However, the USNO-B proper motion survey reveals no candidate object. Moreover, the probability that a field star encountered HD 100546 in the past few thousand years is very low.

Clumps in the solar neighborhood's stellar velocity distribution could be caused by spiral density waves. In the solar neighborhood, stellar velocities corresponding to orbits that are nearly closed in the frame rotating with a spiral pattern represent likely regions for stellar concentrations. Via particle integration, we show that orbits can intersect the solar neighborhood when they are excited by Lindblad resonances with a spiral pattern. We find that a two-armed spiral density wave with pattern speed placing the Sun near the 4:1 Inner Lindblad Resonance (ILR) can cause two families of nearly closed orbits in the solar neighborhood. One family corresponds to square shaped orbits aligned so their peaks lie on top of, and support, the two dominant stellar arms. The second family correspond to orbits 45 degrees out of phase with the other family. Such a spiral density pattern could account for two major clumps in the solar neighborhood's velocity distribution. The Pleiades/Hyades moving group corresponds to the first family of orbits and the Coma Berenices moving group corresponds to the second family. This model requires a spiral pattern speed of approximately 0.66 +- 0.03 times the angular rotation rate of the Sun or 18.1 +- 0.8 km/s/kpc.