• ### Coma Berenices: first evidence for incomplete vertical phase-mixing in local velocity space with RAVE - confirmed with Gaia DR2(1804.07767)

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

Feb. 15, 2018 astro-ph.GA
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.
• ### PLATO as it is: a legacy mission for Galactic archaeology(1706.03778)

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

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

April 18, 2017 astro-ph.GA
We explore the correlations between velocity and metallicity and the possible distinct chemical signatures of the velocity over-densities of the local Galactic neighbourhood. We use the large spectroscopic survey RAVE and the Geneva Copenhagen Survey. We compare the metallicity distribution of regions in the velocity plane ($v_R,v_\phi$) with that of their symmetric counterparts ($-v_R,v_\phi$). We expect similar metallicity distributions if there are no tracers of a sub-population (e.g., a dispersed cluster, accreted stars), if the disk of the Galaxy is axisymmetric, and if the orbital effects of the spiral arms and the bar are weak. We find that the metallicity-velocity space of the solar neighbourhood is highly patterned. A large fraction of the velocity plane shows differences in the metallicity distribution when comparing symmetric $v_R$ regions. The typical differences in the median metallicity are of $0.05$ dex with a statistical significance of at least $95\%$, and with values up to $0.6$ dex. For low azimuthal velocity $v_\phi$, stars moving outwards in the Galaxy have on average higher metallicity than those moving inwards. These include stars in the Hercules and Hyades moving groups and other velocity branch-like structures. For higher $v_\phi$, the stars moving inwards have higher metallicity than those moving outwards. The most likely interpretation of the metallicity asymmetry is that it is due to the orbital effects of the bar and the radial metallicity gradient of the disk. We present a simulation that supports this idea. We have also discovered a positive gradient in $v_\phi$ with respect to metallicity at high metallicities, apart from the two known positive and negative gradients for the thick and thin disks, respectively.
• Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.
• ### Chromospherically Active Stars in the RAVE Survey. II. Young dwarfs in the Solar neighborhood(1612.02433)

Dec. 7, 2016 astro-ph.SR
A large sample of over 38,000 chromospherically active candidate solar-like stars and cooler dwarfs from the RAVE survey is addressed in this paper. An improved activity identification with respect to the previous study was introduced to build a catalog of field stars in the Solar neighborhood with an excess emission flux in the calcium infrared triplet wavelength region. The central result of this work is the calibration of the age--activity relation for the main sequence dwarfs in a range from a few $10 \; \mathrm{Myr}$ up to a few Gyr. It enabled an order of magnitude age estimation of the entire active sample. Almost 15,000 stars are shown to be younger than $1\;\mathrm{Gyr}$ and $\sim$2000 younger than $100\;\mathrm{Myr}$. The young age of the most active stars is confirmed by their position off the main sequence in the $J-K$ versus $N_{UV}-V$ diagram showing strong ultraviolet excess, mid-infrared excess in the $J-K$ versus $W_1-W_2$ diagram and very cool temperatures ($J-K>0.7$). They overlap with the reference pre-main sequence RAVE stars often displaying X-ray emission. The activity level increasing with the color reveals their different nature from the solar-like stars and probably represents an underlying dynamo generating magnetic fields in cool stars. 50\% of the RAVE objects from DR5 are found in the TGAS catalog and supplemented with accurate parallaxes and proper motions by Gaia. This makes the database of a large number of young stars in a combination with RAVE's radial velocities directly useful as a tracer of the very recent large-scale star formation history in the Solar neighborhood. The data are available online in the Vizier database.
• ### Identification of Globular Cluster Stars in RAVE data II: Extended tidal debris around NGC 3201(1601.02282)

Jan. 10, 2016 astro-ph.GA
We report the identification of extended tidal debris potentially associated with the globular cluster NGC 3201, using the RAVE catalogue. We find the debris stars are located at a distance range of 1-7 kpc based on the forthcoming RAVE distance estimates. The derived space velocities and integrals of motion show interesting connections to NGC 3201, modulo uncertainties in the proper motions. Three stars, which are among the 4 most likely candidates for NGC 3201 tidal debris, are separated by 80 degrees on the sky yet are well matched by the 12 Gyr, [Fe/H] = -1.5 isochrone appropriate for the cluster. This is the first time tidal debris around this cluster has been reported over such a large spatial extent, with implications for the cluster$'$s origin and dynamical evolution.
• ### RAVE stars tidally stripped/ejected from $\omega$ Centauri globular cluster(1508.02839)

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

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

April 29, 2015 astro-ph.GA
We present the identification of potential members of nearby Galactic globular clusters using radial velocities from the RAdial Velocity Experiment Data Release 4 (RAVE-DR4) survey database. Our identifications are based on three globular clusters -- NGC 3201, NGC 5139 ($\omega$ Cen) and NGC 362 -- all of which are shown to have |RV|>100 km/s. The identification of globular cluster stars in RAVE DR4 data offers a unique opportunity to test the precision and accuracy of the stellar parameters determined with the currently available Stellar Parameter Pipelines (SPPs) used in the survey, as globular clusters are ideal testbeds for the validation of stellar atmospheric parameters, abundances, distances and ages. For both NGC 3201 and $\omega$ Cen, there is compelling evidence for numerous members (> 10) in the RAVE database; in the case of NGC 362 the evidence is more ambiguous, and there may be significant foreground and/or background contamination in our kinematically-selected sample. A comparison of the RAVE-derived stellar parameters and abundances with published values for each cluster and with BASTI isochrones for ages and metallicities from the literature reveals overall good agreement, with the exception of the apparent underestimation of surface gravities for giants, in particular for the most metal-poor stars. Moreover, if the selected members are part of the main body of each cluster our results would also suggest that the distances from Binney et al. 2013, where only isochrones more metal-rich than -0.9 dex were used, are typically underestimated by ~ 40% with respect to the published distances for the clusters, while the distances from Zwitter et al. 2010 show stars ranging from 1 to ~ 6.5 kpc -- with indications of a trend toward higher distances at lower metallicities -- for the three clusters analysed in this study.
• ### The imprints of the Galactic bar on the thick disk with RAVE(1503.02072)

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

Dec. 19, 2014 astro-ph.GA
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.
• ### Characterizing the High-Velocity Stars of RAVE: The Discovery of a Metal-Rich Halo Star Born in the Galactic Disk(1412.1484)

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

Oct. 2, 2014 astro-ph.CO, astro-ph.GA
We determine the Galactic potential in the solar neigbourhood from RAVE observations. We select red clump stars for which accurate distances, radial velocities, and metallicities have been measured. Combined with data from the 2MASS and UCAC catalogues, we build a sample of 4600 red clump stars within a cylinder of 500 pc radius oriented in the direction of the South Galactic Pole, in the range of 200 pc to 2000 pc distances. We deduce the vertical force and the total mass density distribution up to 2 kpc away from the Galactic plane by fitting a distribution function depending explicitly on three isolating integrals of the motion in a separable potential locally representing the Galactic one with four free parameters. Because of the deep extension of our sample, we can determine nearly independently the dark matter mass density and the baryonic disc surface mass density. We find (i) at 1kpc Kz/(2piG) = 68.5 pm 1.0 Msun/pc2, and (ii) at 2 kpc Kz/(2piG) = 96.9 pm 2.2 Msun/pc2. Assuming the solar Galactic radius at R0 = 8.5 kpc, we deduce the local dark matter density rhoDM (z=0) = 0.0143 pm 0.0011Msun pc3 = 0.542 pm 0.042 Gev/cm3 and the baryonic surface mass density Sigma = 44.4 pm 4.1 Msun/pc2 . Our results are in agreement with previously published Kz determinations up to 1 kpc, while the extension to 2 kpc shows some evidence for an unexpectedly large amount of dark matter. A flattening of the dark halo of order 0.8 can produce such a high local density in combination with a circular velocity of 240 km/s . Another explanation, allowing for a lower circular velocity, could be the presence of a secondary dark component, a very thick disc resulting either from the deposit of dark matter from the accretion of multiple small dwarf galaxies, or from the presence of an effective phantom thick disc in the context of effective galactic-scale modifications of gravity.
• ### Constraining the Galaxy's dark halo with RAVE stars(1406.4130)

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

Aug. 26, 2014 astro-ph.GA, astro-ph.SR
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.
• ### APASS Landolt-Sloan BVgri photometry of RAVE stars. I. Data, effective temperatures and reddenings(1408.5476)

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

June 27, 2014 astro-ph.GA
We provide new constraints on the chemo-dynamical models of the Milky Way by measuring the radial and vertical chemical gradients for the elements Mg, Al, Si, Ti, and Fe in the Galactic disc and the gradient variations as a function of the distance from the Galactic plane ($Z$). We selected a sample of giant stars from the RAVE database using the gravity criterium 1.7$<$log g$<$2.8. We created a RAVE mock sample with the Galaxia code based on the Besan\c con model and selected a corresponding mock sample to compare the model with the observed data. We measured the radial gradients and the vertical gradients as a function of the distance from the Galactic plane $Z$ to study their variation across the Galactic disc. The RAVE sample exhibits a negative radial gradient of $d[Fe/H]/dR=-0.054$ dex kpc$^{-1}$ close to the Galactic plane ($|Z|<0.4$ kpc) that becomes flatter for larger $|Z|$. Other elements follow the same trend although with some variations from element to element. The mock sample has radial gradients in fair agreement with the observed data. The variation of the gradients with $Z$ shows that the Fe radial gradient of the RAVE sample has little change in the range $|Z|\lesssim0.6$ kpc and then flattens. The iron vertical gradient of the RAVE sample is slightly negative close to the Galactic plane and steepens with $|Z|$. The mock sample exhibits an iron vertical gradient that is always steeper than the RAVE sample. The mock sample also shows an excess of metal-poor stars in the [Fe/H] distributions with respect to the observed data. These discrepancies can be reduced by decreasing the number of thick disc stars and increasing their average metallicity in the Besan\c con model.
• ### Kinematic modelling of the Milky Way using the RAVE and GCS stellar surveys(1405.7435)

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

Jan. 28, 2014 astro-ph.GA
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.
• ### A new stellar chemo-kinematic relation reveals the merger history of the Milky Way disc(1310.5145)

Dec. 25, 2013 astro-ph.CO, astro-ph.GA
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.
• ### The RAVE survey: the Galactic escape speed and the mass of the Milky Way(1309.4293)

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