$\texttt{galkin}$: a new compilation of the Milky Way rotation curve data
Miguel Pato, Fabio Iocco

TL;DR
galkin is a comprehensive, user-friendly compilation of Milky Way rotation curve data from various sources, designed to aid dynamical studies and dark matter research.
Contribution
It provides a new, curated dataset of 2780 kinematic measurements and a tool for easy data handling, covering 3 to 20 kpc in the Galaxy.
Findings
Enables detailed dynamical analysis of the Milky Way
Facilitates constraints on dark matter distribution
Supports tests of modified gravity theories
Abstract
We present , a novel compilation of kinematic measurements tracing the rotation curve of our Galaxy, together with a tool to treat the data. The compilation is optimised to Galactocentric radii between 3 and 20 kpc and includes the kinematics of gas, stars and masers in a total of 2780 measurements carefully collected from almost four decades of literature. A simple, user-friendly tool is provided to select, treat and retrieve the data of all source references considered. This tool is especially designed to facilitate the use of kinematic data in dynamical studies of the Milky Way with various applications ranging from dark matter constraints to tests of modified gravity.
| tracer type | [kpc] | quadrants | tracers | |
| gas kinematics | HI terminal velocities | |||
| Fich+ ’89 Fich et al. (1989) | 2.1 – 8.0 | 1,4 | 149/149 | |
| Malhotra ’95 Malhotra (1995) | 2.1 – 7.5 | 1,4 | 110/110 | |
| McClure-Griffiths & Dickey ’07 McClure-Griffiths and Dickey (2007) | 2.8 – 7.6 | 4 | 701/761 | |
| HI thickness method | ||||
| Honma & Sofue ’97 Honma and Sofue (1997) | 6.8 – 20.2 | – | 13/13 | |
| CO terminal velocities | ||||
| Burton & Gordon ’78 Burton and Gordon (1978) | 1.4 – 7.9 | 1 | 284/284 | |
| Clemens ’85 Clemens (1985) | 1.9 – 8.0 | 1 | 143/143 | |
| Knapp+ ’85 Knapp et al. (1985) | 0.6 – 7.8 | 1 | 37/37 | |
| Luna+ ’06 Luna et al. (2006) | 2.0 – 8.0 | 4 | 272/457 | |
| HII regions | ||||
| Blitz ’79 Blitz (1979) | 8.7 – 11.0 | 2,3 | 3/3 | |
| Fich+ ’89 Fich et al. (1989) | 9.4 – 12.5 | 3 | 5/104 | |
| Turbide & Moffat ’93 Turbide and Moffat (1993) | 11.8 – 14.7 | 3 | 5/8 | |
| Brand & Blitz ’93 Brand and Blitz (1993) | 5.2 – 16.5 | 1,2,3,4 | 148/206 | |
| Hou+ ’09 Hou et al. (2009) | 3.5 – 15.5 | 1,2,3,4 | 274/815 | |
| giant molecular clouds | ||||
| Hou+ ’09 Hou et al. (2009) | 6.0 – 13.7 | 1,2,3,4 | 30/963 | |
| star kinematics | open clusters | |||
| Frinchaboy & Majewski ’08 Frinchaboy and Majewski (2008) | 4.6 – 10.7 | 1,2,3,4 | 60/71 | |
| planetary nebulae | ||||
| Durand+ ’98 Durand et al. (1998) | 3.6 – 12.6 | 1,2,3,4 | 79/867 | |
| classical cepheids | ||||
| Pont+ ’94 Pont et al. (1994) | 5.1 – 14.4 | 1,2,3,4 | 245/278 | |
| Pont+ ’97 Pont et al. (1997) | 10.2 – 18.5 | 2,3,4 | 32/48 | |
| carbon stars | ||||
| Demers & Battinelli ’07 Demers and Battinelli (2007) | 9.3 – 22.2 | 1,2,3 | 55/103 | |
| Battinelli+ ’13 Battinelli et al. (2013) | 12.1 – 24.8 | 1,2 | 35/36 | |
| masers | masers | |||
| Reid+ ’14 Reid et al. (2014) | 4.0 – 15.6 | 1,2,3,4 | 80/103 | |
| Honma+ ’12 Honma et al. (2012) | 7.7 – 9.9 | 1,2,3,4 | 11/52 | |
| Stepanishchev & Bobylev ’11 Stepanishchev and Bobylev (2011) | 8.3 | 3 | 1/1 | |
| Xu+ ’13 Xu et al. (2013) | 7.9 | 4 | 1/30 | |
| Bobylev & Bajkova ’13 Bobylev and Bajkova (2013) | 4.7 – 9.4 | 1,2,4 | 7/31 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
galkin: a new compilation of the Milky Way rotation curve data
Miguel Pato
The Oskar Klein Centre for Cosmoparticle Physics, Department of Physics, Stockholm University, AlbaNova, SE-106 91 Stockholm, Sweden
Fabio Iocco
ICTP South American Institute for Fundamental Research, and Instituto de Física Teórica - Universidade Estadual Paulista (UNESP), Rua Dr. Bento Teobaldo Ferraz 271, 01140-070 São Paulo, SP Brazil
Abstract
We present galkin, a novel compilation of kinematic measurements tracing the rotation curve of our Galaxy, together with a tool to treat the data. The compilation is optimised to Galactocentric radii between and kpc and includes the kinematics of gas, stars and masers in a total of 2780 measurements carefully collected from almost four decades of literature. A simple, user-friendly tool is provided to select, treat and retrieve the data of all source references considered. This tool is especially designed to facilitate the use of kinematic data in dynamical studies of the Milky Way with various applications ranging from dark matter constraints to tests of modified gravity.
rotation curve, Galactic dynamics, Milky Way
I Introduction
The rotation curve of a spiral galaxy provides far-reaching insight into its properties, as noticeably explored for decades now (see e.g. Refs. Freeman (1970); Roberts and Rots (1973); Bosma (1978); Rubin et al. (1978); Persic and Salucci (1995)). Data on the rotation curve of our Galaxy, a spiral itself, have also been available for several decades Blitz (1979); Burton and Gordon (1978); Clemens (1985); Knapp et al. (1985); Fich et al. (1989); Pont et al. (1994). However, the data are rather disperse throughout the literature and groups of references are often neglected. We therefore set out to assemble a comprehensive compilation of the decades-long observational effort to pinpoint the rotation curve of the Milky Way. The compilation, named galkin111To download your copy of galkin, please refer to our GitHub page github.com/galkintool/galkin or contact us at [email protected]., improves upon existing ones (e.g. Refs. Sofue et al. (2009a); Bhattacharjee et al. (2014)) on several aspects, including most notably: (i) an enlarged database of observations appropriately treated for unified use, and (ii) the release of a simple out-of-the-box tool to retrieve the data. This compilation has been presented in Ref. Iocco et al. (2015) and later adopted in other works in the literature. Without venturing into any analysis of the Galactic structure or dynamics (as done in galpy Bovy (2015)), here we provide instead a thorough description of the data sets (Sec. II) as well as the features of an out-of-the-box tool (Sec. III) to access the database and output the desired data for independent analyses. The open source code provided is simple, flexible and can be easily modified to include new datasets or other types of measurements. The latter feature is particularly relevant on the eve of the precision era soon to be introduced by the Gaia satellite de Bruijne (2012) and an array of optical and near-infrared ground-based surveys such as APOGEE-2 apo ; Majewski et al. (2015), GALAH De Silva et al. (2015), WEAVE wea and 4MOST de Jong et al. (2012). Our compilation can be regarded as a step forward in unifying the current state of the art, yet it is certainly susceptible of further inclusions – please see our own extensive caveats and notes throughout the manuscript. We encourage the community to adopt galkin and participate in its extension as new datasets arise.
II Compilation
The rationale behind our literature survey is to be as exhaustive as possible, but still selective enough to put together a clean and reliable sample of kinematic tracers. For that reason, we have decided to exclude tracers with kinematic distance determination only, important asymmetric drift or large random motions. Whereas we have taken utter care in making such compilation exhaustive, we cannot exclude that some data sets may be missing; the reader compelled to add extra data sets is welcome to do so by modifying the open source code described in Sec. III. Our compilation is focussed on the range of Galactocentric radii kpc and is not intended for use far outside this range. In particular, compilations dedicated to the very internal regions of the Milky Way are available elsewhere in the literature Sofue (2013). However, we warn the reader that the use of kinematic tracers at kpc may be problematic for certain applications Chemin et al. (2015).
Note as well that kinematic tracers beyond the ones in our compilation exist and are discussed at length in the literature, including tracers in the outer halo Dehnen et al. (2006); Xue et al. (2008); Kafle et al. (2012); Bhattacharjee et al. (2014); Kafle et al. (2014), disc Bovy et al. (2012); Bovy and Rix (2013); Loebman et al. (2014); Russeil (2003) and off the Galactic plane Kuijken and Gilmore (1991); Holmberg and Flynn (2004); Moni Bidin et al. (2012); Bovy and Tremaine (2012); Moni Bidin et al. (2015). For careful analyses of the implications of these and other tracers on the mass distribution of the Milky Way, we refer the reader to the works cited above and also to Refs. Garbari et al. (2011); Smith et al. (2012); Garbari et al. (2012); Zhang et al. (2013); Read (2014); Dehnen and Binney (1998); Sofue et al. (2009b); Catena and Ullio (2010); Weber and de Boer (2010); Salucci et al. (2010); Iocco et al. (2011); McMillan (2011); Nesti and Salucci (2013); Sofue (2015), where the authors have often times built their own tracer compilations. Partial compilations are easily reproducible with the galkin tool given the possibility to (de-)select individual references. Again, we remind the interested reader that our code galkin is modular and easily modifiable by the user for the addition of new data sets.
We divide up the data in three main categories:
gas kinematics, including
- •
HI terminal velocities Fich et al. (1989); Malhotra (1995); McClure-Griffiths and Dickey (2007),
- •
HI thickness Honma and Sofue (1997),
- •
CO terminal velocities Burton and Gordon (1978); Clemens (1985); Knapp et al. (1985); Luna et al. (2006),
- •
HII regions Blitz (1979); Fich et al. (1989); Turbide and Moffat (1993); Brand and Blitz (1993); Hou et al. (2009), and
- •
giant molecular clouds Hou et al. (2009); 2. 2.
star kinematics, including
- •
open clusters Frinchaboy and Majewski (2008),
- •
planetary nebulae Durand et al. (1998),
- •
classical cepheids Pont et al. (1994, 1997), and
- •
carbon stars Demers and Battinelli (2007); Battinelli et al. (2013); and 3. 3.
masers Reid et al. (2014); Honma et al. (2012); Stepanishchev and Bobylev (2011); Xu et al. (2013); Bobylev and Bajkova (2013).
Tab. 1 recaps the key features of each data set. We refer to the original source references and in particular to the supplementary information of Ref. Iocco et al. (2015) for in-depth descriptions of the different object types. In order to obtain a clean sample, we have imposed various data selection cuts on the available sources following closely the recommendations of each original reference (notice, however, that the interested user can easily override these software cuts hacking the source code directly). The Appendix at the end of the present paper gives a full account of our data selection and treatment for each reference listed in Tab. 1.
Our compilation consists of 2780 individual tracers distributed in Galactocentric radius , Galactic longitude and height above Galactic plane as shown in Fig. 1. Each object is specified by its coordinates , heliocentric distance and heliocentric line-of-sight velocity . The uncertainties on and are largely unimportant and hence neglected, whereas the uncertainties on and are taken from the original references (cf. details in the Appendix). In radio observations, it is customary to report measurements of in terms of the line-of-sight velocity in the local standard of rest (LSR) for a fixed peculiar solar motion . In these cases, we infer by subtracting the peculiar solar motion used in the reference off the reported (cf. the Appendix). Once is obtained, this is summed to the adopted peculiar solar motion to get the final LSR line-of-sight velocity . Each object has an associated measurement . The corresponding Galactocentric radius follows from simple geometry as
[TABLE]
where is the distance of the Sun to the Galactic centre. Under the assumption of circular orbits, the angular circular velocity of the object is found by inverting
[TABLE]
where is the local circular velocity. The uncertainties on and are propagated to and , respectively. We shall also provide the familiar circular velocity and corresponding uncertainties, but note that the errors of and are strongly positively correlated, while those of and are independent. All uncertainties currently implemented in galkin are symmetric following the information available in each reference; future data might provide the full distribution of observables, which would then be treated in upcoming versions of galkin and would be of great value for Bayesian studies. The procedure described above is common to all object types in Tab. 1, with some modifications in two cases. For terminal velocities, we set and (or, equivalently, ) in Eqs. (1) and (2), and each measurement reads . For the HI thickness method, the measured quantity is instead of , so each data point is defined by , cf. Refs. Merrifield (1992); Honma and Sofue (1997). We also process the proper motions when available, as is the case of open clusters and masers. All other details of data treatment are exhaustively documented in the Appendix.
Note that Eqs. (1) and (2) can be used to infer only upon fixing of , and , which are quantities still affected by sizeable uncertainties – see Refs. Gillessen et al. (2009); Ando et al. (2011); Malkin (2012); Reid et al. (2014) for , Refs. Reid and Brunthaler (2004); Reid et al. (2009); Bovy et al. (2009); McMillan and Binney (2010); Bovy et al. (2012); Reid et al. (2014) for and Refs. Dehnen and Binney (1998); Schönrich et al. (2010); Bovy et al. (2012); Reid et al. (2014) for . Making the reasonable choice of Galactic parameters kpc, km/s and km/s Schönrich et al. (2010), we show in Fig. 2 the positions and circular velocities of all tracers implemented. This is the main output of galkin.
III Tool description
We now turn to the description of the galkin tool, whose function is to allow the user to access the data described Sec. II in a customisable way. The tool is written in Python and the code package is available through our GitHub page, github.com/galkintool/galkin.
The distribution has a very simple structure. The parent directory contains the setup file setup.py together with the usual README file. The galkin package is provided under galkin/, where galkin/data/ contains the original data as extracted from the corresponding references222Note that the main point of this tool is to handle in a unified way data from references using different Galactic parameters, different definitions and often overlapping in observed sources. The data files under data/ contain the original published data of each reference; please do not use these files directly unless you are fully aware of all details of each reference., and example scripts can be found under bin/. The main dependencies and installation instructions are described in the README file, whereas here we simply summarise the usage of the tool assuming it is properly installed and running. Note that galkin adopts astropy Astropy Collaboration et al. (2013) for coordinate conversion if required by the user, but the code can also be used without installing this package, which is the default option.
The goal of galkin is to provide ready-to-use data files containing all the necessary kinematic tracer information for constraining the rotation curve of the Galaxy. The user may choose the values for the Galactic parameters and , as well as (de)select either entire classes of tracers (gas, stars, masers) or single references in Tab. 1. It is also possible to add a given systematic uncertainty to the line-of-sight velocity of each tracer. This is done with the help of the script bin/galkin_data.py either through a graphic interface (see Fig. 3) by launching the tool from bin/ with the command
python galkin_data.py
or through a customisable input file by typing
python galkin_data.py inputpars.txt
The code then processes the original data sets of the selected references converting each data set consistently for the chosen values of the Galactic parameters. The output is stored in three data files:
- bin/output/posdata.dat with the position information of each tracer, namely ;
- bin/output/vcdata.dat with the rotation curve measurements, namely ; and
- bin/output/rawdata.dat with the raw data measurements, namely .
In the first two files, the values of the Galactic parameters chosen by the user are reported in the first line and the source reference for each tracer is indicated in the last column. For testing purposes, we provide under bin/output/ the sample files corresponding to our baseline choice of Galactic parameters used to produce Fig. 2 (for the entire database and single classes objects, separately) as well as for a variation of the Galactic parameters (for the entire database only).
The tool also includes the script bin/galkin_plotter.py to read and visualise the output described above. This can be launched from bin/ with the command
python galkin_plotter.py output/vcdata.dat output/posdata.dat
which produces a set of demonstrative plots including the spatial distribution of the tracers and the inferred rotation curve.
Finally, we provide the script bin/galkin_data_fast.py in order to illustrate how to use galkin inside another code without dealing with input nor output files. This is a faster version of the data processing pipeline that is specifically designed for applications that need to use galkin repeateadly, e.g. in scans over Galactic parameters.
The structure of galkin is purposely minimal and modular. The code can be easily adapted to replace or modify existing data sets or single tracers and to add new data sets as they become available. The idea behind galkin is to provide an up-to-date and user-friendly compilation of tracers of the kinematics of the Milky Way over the coming years.
Acknowledgements. M. P. acknowledges the support from Wenner-Gren Stiftelserna in Stockholm; F. I. from the Simons Foundation and FAPESP process 2014/11070-2.
Appendix A Data treatment
Tracers of the rotation curve of the Milky Way usually adopted in the literature can be roughly divided into three categories: gas kinematics, kinematics of stellar objects, and masers. For each of these classes of objects, different methodologies are used to infer positions, distances and kinematics. Without attempting here a review of the properties of each tracer type, we point to the supplementary information of Ref. Iocco et al. (2015), where this collection of data has first been presented. The reader will find there an in-depth description of all tracer types as well as appropriate references to the original literature. In this Appendix, we recap the original source references and fully document our data treatment. Let us notice that halo stars are currently not included in this version of galkin, due to the additional assumptions needed. Future versions of the code will properly include halo stars as an extra tracer type.
Tab. 1 displays the key features of all data sets used in galkin. In the following we provide a detailed description of the treatment applied to each data set. Equation, figure and table numbering refers to the original source references.
A.1 Gas kinematics
A.1.1 HI terminal velocities
Fich+ ’89 Fich et al. (1989)
From Tab. 2, we take the terminal velocities measured at different longitudes and assume an overall velocity uncertainty km/s following Sec. II.b.i. The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation333 Notice that the original reference fails to indicate explicitely the adopted peculiar solar motion. Following Ref. Reid et al. (2009), we assume the adopted value is the old standard solar motion as defined in the Appendix of Ref. Reid et al. (2009).. We further correct for the peculiar LSR motion of km/s in the radial direction, cf. Sec. II.b.i and Eq. (4)444 Notice that this correction depends on the assumed peculiar solar motion, so this procedure is strictly valid only when the user selects the same peculiar solar motion as in the original reference; in practice, however, this is a small correction and we ignore the slight inconsistency..
Malhotra ’95 Malhotra (1995)
From Fig. 7, we take the terminal velocities measured at different Galactocentric radii and assume an overall velocity uncertainty km/s in line with the dispersions computed in Sec. 3.4 for both the first and fourth quadrants. The Galatocentric radii are converted into longitudes through depending on the quadrant (first quadrant: circles and triangles in top panel of Fig. 7; fourth quandrant: squares in bottom panel of Fig. 7). The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation.
McClure-Griffiths & Dickey ’07 McClure-Griffiths and
Dickey (2007)
From Tab. 1 (online version), we take the terminal velocities measured at different longitudes . According to Secs. 3.3 and 4.3.1 and the caption of Fig. 8, the velocity uncertainty amounts to km/s for , km/s for and km/s for ; we conservatively assume km/s for the whole range . The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We exclude from the full data set regions with discrete HI clouds at (cf. Sec. 4; the width of the intervals is fixed a posteriori in order to eliminate the spikes in velocity) and also the region at because there (cf. Sec. 4.2; this last cut is actually already performed in the online version of Tab. 1).
A.1.2 HI thickness method
Honma & Sofue ’97 Honma and Sofue (1997)
From Tab. 1, we take the Galactocentric radii fitted through the thickness method for different velocities and assume an overall uncertainty km/s following Sec. 2.4. The authors of this reference find that the method of Merrifield ’92 Merrifield (1992) (i.e. method 1 in Tab. 1) covers the largest range of and is the most accurate, so we use the results of that method for our compilation.
A.1.3 CO terminal velocities
Burton & Gordon ’78 Burton and Gordon (1978)
From Fig. 2, we take the terminal velocities measured at different longitudes , where already includes the line width correction, cf. Sec. 3 and Eq. (4a). The resolution of the original CO data presented in this reference is km/s, while it amounts to km/s for the other data sets Bania (1977); Liszt et al. (1977) plotted in Fig. 2 (cf. Sec. 2); therefore, we conservatively assume km/s all along. The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation.
Clemens ’85 Clemens (1985)
From Tab. 2, we take the terminal velocities measured at different longitudes and the corresponding uncertainties , which actually represent 0.67 of the standard deviation, cf. footnote c in Tab. 2. The terminal velocities in Tab. 2 are uncorrected for the line width, so the true terminal velocities read km/s (cf. Sec. II.c) and the corresponding uncertainties are , where km/s is the typical error in the line width correction, cf. Sec. II.c. The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We further correct for the peculiar LSR motion of km/s in the azimuthal direction, cf. Sec. II.d.
Knapp+ ’85 Knapp et al. (1985)
From Fig. 5, we take the terminal velocities (top panel) and the velocity dispersions (central panel) measured at different Galactocentric radii . The statistical uncertainty of the terminal velocities is summed in quadrature to the dispersion, i.e. . The Galatocentric radii are converted into longitudes through considering that the data refers to the first quadrant only. The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation.
Luna+ ’06 Luna et al. (2006)
From Fig. 1 (top panel), we take the terminal velocities measured at different longitudes and assume an overall velocity uncertainty km/s following Sec. 3.1. The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We exclude from the full data set the region at due to the influence of the Carina arm and Centaurus, cf. Sec. III.a.iii in Ref. Alvarez et al. (1990).
A.1.4 HII regions
Blitz ’79 Blitz (1979)
From Tab. 2, we take the distances and line-of-sight velocities corresponding to objects towards different Galactic coordinates . The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We exclude from the full data set all objects in (or coincident to objects in) Fich+ ’89, Turbide & Moffat ’93, Brand & Blitz ’93 and Hou+ ’09; only three objects remain: S148, S306 and Mon R2.
Fich+ ’89 Fich et al. (1989)
From Tab. 1, we take the distances and line-of-sight velocities corresponding to objects towards different Galactic coordinates and add a random motion of km/s in quadrature to following Sec. II.i.b. The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We further correct for the peculiar LSR motion of km/s in the radial direction, cf. Sec. II.b.i and Eq. (4). We exclude from the full data set all objects in (or coincident to objects in) Turbide & Moffat ’93 and Brand & Blitz ’93.
Turbide & Moffat ’93 Turbide and Moffat (1993)
From Tab. 5 (for ), we take the distances and line-of-sight velocities corresponding to objects towards different Galactic coordinates . Note that the reported velocities correspond to Eq. (11) with km/s and a peculiar solar motion km/s, which is consistently subtracted off in our compilation. For the objects in Tab. 5 based on Refs. [1,3] therein, we further correct for the peculiar LSR motion of km/s in the radial direction, cf. Sec. 4 and Eq. (11). We exclude from the full data set all objects in (or coincident to objects in) Brand & Blitz ’93.
Brand & Blitz ’93 Brand and Blitz (1993)
From Tab. 1, we take the distances and line-of-sight velocities corresponding to objects towards different Galactic coordinates and add a random motion of km/s in quadrature to following Sec. 3.1. The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We exclude from the full data set objects near to the Galactic centre or anti Galactic centre (cf. Sec. 3.3) and nearby objects at kpc (cf. Sec. 3.3).
Hou+ ’09 Hou et al. (2009)
From Tab. A1, we take the distances and line-of-sight velocities corresponding to objects towards different Galactic coordinates and assume an overall velocity uncertainty km/s. The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We exclude from the full data set objects without stellar distances nor at the tangent points, near to the Galactic centre or anti Galactic centre and objects in (or coincident to objects in) Brand & Blitz ’93.
A.1.5 Giant molecular clouds
Hou+ ’09 Hou et al. (2009)
From Tab. A2, we take the distances and line-of-sight velocities corresponding to objects towards different Galactic coordinates and adopt an overall velocity uncertainty km/s and an overall relative distance uncertainty . The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We exclude from the full data set objects without stellar distances and near to the Galactic centre or anti Galactic centre .
A.2 Star kinematics
A.2.1 Open clusters
Frinchaboy & Majewski ’08 Frinchaboy and
Majewski (2008)
From Tab. 1, we take the distances corresponding to objects towards different equatorial coordinates (then converted to Galactic coordinates ) based on Ref. Dias et al. (2002) and assume an overall relative distance uncertainty . From Tab. 12, we take the heliocentric line-of-sight velocities , where whenever possible we use the bulk kinematics derived with the three-dimensional and radial velocity membership criterion (3D+RV), and also the proper motions , . We exclude from the full data set objects near to the anti Galactic centre ; NGC 1513, NGC 7654 (only one member identified, cf. Sec. 6.1); Collinder 258, Lynga 1, NGC 6250 (probably wrong memberships, cf. Sec. 6.2.3); and NGC 6416 (high residual velocity).
A.2.2 Planetary nebulae
Durand+ ’98 Durand et al. (1998)
From Tab. 2 (online version), we take the heliocentric line-of-sight velocities corresponding to objects towards different equatorial coordinates (then converted to Galactic coordinates ). The distances are found by cross-matching Tab. 2 to Tabs. 1 and 3 of Ref. Zhang (1995), which report individually determined and statistical distances, respectively. When available, we prefer the individually determined distances (i.e. Tab. 1 of Ref. Zhang (1995)), based on the methods presented in Ref. Zhang (1993) with a relative uncertainty ranging from 15 to 40% (cf. Sec. 3.3 in Ref. Zhang (1993)) so an overall relative distance uncertainty is adopted; when not possible, we use the statistical distance scale (i.e. Tab. 3 of Ref. Zhang (1995)), for which the relative uncertainty is (cf. Sec. 5.2 in Ref. Zhang (1995)). We further correct for the -term of km/s (cf. Eq. (4) and Tab. 3) and model the asymmetric drift with respect to the peculiar solar motion in use (i.e. in the radial direction and in the azimuthal direction; cf. Eq. (4) and Tab. 3). The error of , and (cf. Tab. 3) are propagated to the line-of-sight velocity . We exclude from the full data set objects near to the Galactic centre and off the Galactic plane pc (cf. Sec. 4.2); NGC 6565, M2-29, M1-46, Sa 1-8, NGC 6741, NGC 7293, Vy 2-2, NGC 6702, NGC 7094, NGC 2392, NGC 4071, NGC 6026, NGC 6302, M2-7, Th 3-14 (deficient statistical distances, cf. Sec. 5.2 in Ref. Zhang (1995); notice typo in this reference where NGC 6320 is mentioned instead of NGC 6302); BoBn 1 (atypical motion; cf. Sec. 4.1) and NGC 6567, A 8, Pu 1, M1-5 (high residual velocity).
A.2.3 Classical cepheids
Pont+ ’94 Pont et al. (1994)
From Tab. 3, we take the distance moduli (FW) and heliocentric line-of-sight velocities corresponding to objects towards different Galactic coordinates and assume an overall distance modulus uncertainty mag following Sec. 11.3. The distance moduli are converted to distances through . The velocity uncertainty reads , where km/s depending on the method used to calculate the radial velocity (cf. Sec. 11.3) and km/s is the contribution of the velocity ellipsoid (cf. Sec. 11.3). We further correct for the -term of km/s, cf. Sec. 11.6. We use the reduced sample of 266 stars (cf. Sec. 11.4 and Tab. 4) excluding objects near to the anti Galactic centre and two further objects due to high residual velocity.
Pont+ ’97 Pont et al. (1997)
From Tab. 1, we take the heliocentric line-of-sight velocities , period and colours and corresponding to objects towards different Galactic coordinates . The distance moduli are found with the period-luminosity-colour relation described in Ref. Pont et al. (1994) (FW) and the zero point given in Sec. 3.1; an overall distance modulus uncertainty mag is assumed following Sec. 3.3. The distance moduli are converted to distances through . The velocity uncertainty reads , where km/s for velocity measurements (cf. Tab. 1) and km/s is the contribution of the velocity ellipsoid (cf. Sec. 11.3 of Ref. Pont et al. (1994)). We further add a km/s systematic uncertainty (cf. Sec. 6; see also Sec. 5.3) to the derived circular velocity. We exclude from the full data set objects near to the anti Galactic centre and with no measured radial velocity or no .
A.2.4 Carbon stars
Demers & Battinelli ’07 Demers and
Battinelli (2007)
From Tab. 4 (online version), we take the heliocentric line-of-sight velocities corresponding to objects towards different Galactic coordinates and add a random motion of km/s in quadrature to following Sec. 5.1. From Tab. 1 (online version), we take the distances and assume an overall relative distance uncertainty following Sec. 5.1. We exclude from the full data set objects near to the anti Galactic centre ; stars no. 20 and 23 (nearby fast stars; cf. Sec. 5.2) and stars no. 17, 18, 27, 28, 35, 42, 52, 56, 58 and 60 (possibly belonging to Canis Major; cf. Sec. 5.2).
Battinelli+ ’13 Battinelli et al. (2013)
From Tab. 1, we take the Galactocentric radii (computed with kpc) and heliocentric line-of-sight velocities corresponding to objects towards different Galactic coordinates and add a random motion of km/s in quadrature to following Sec. 5.1 of Ref. Demers and Battinelli (2007). The Galactocentric radii are converted to heliocentric distances by inverting Eq. (1) using kpc. We exclude from the full data set star no. 712 (probably belongs to the Sagittarius dwarf galaxy; cf. Sec. 4).
A.3 Masers
Reid+ ’14 Reid et al. (2014)
From Tab. 1, we take the parallaxes , line-of-sight velocities and equatorial proper motions , corresponding to objects towards different equatorial coordinates (then converted to Galactic coordinates ) and add a virial motion of km/s in quadrature to following Sec. 3 in Ref. Reid et al. (2009). The parallaxes are converted to distances through , while the equatorial proper motions are converted to Galactic proper motions , . The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We further corrected for the mean peculiar motion of the masers, i.e. km/s (cf. Tab. 4, A5) and km/s (cf. Sec. 4.4). We exclude from the full data set eight masers with kpc (cf. footnote 3) and 15 outlier masers (cf. footnote 4).
Honma+ ’12 Honma et al. (2012)
From Tab. 1, we take the parallaxes , line-of-sight velocities and equatorial proper motions , corresponding to objects towards different Galactic coordinates and add a virial motion of km/s in quadrature to following Sec. 3 in Ref. Reid et al. (2009). The parallaxes are converted to distances through , while the equatorial proper motions are converted to Galactic proper motions , . The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We further corrected for the mean peculiar motion of the masers, i.e. km/s (cf. Tab. 4, A5 in Ref. Reid et al. (2014)) and km/s (cf. Sec. 4.4 in Ref. Reid et al. (2014)). We exclude from the full data set all objects in (or coincident to objects in) Reid+ ’14.
Stepanishchev & Bobylev ’11 Stepanishchev and
Bobylev (2011)
From Tab. 1, we take the parallaxes , heliocentric line-of-sight velocities and equatorial proper motions , corresponding to objects towards different equatorial coordinates (then converted to Galactic coordinates ) and add a virial motion of km/s in quadrature to following Sec. 3 in Ref. Reid et al. (2009). The parallaxes are converted to distances through , while the equatorial proper motions are converted to Galactic proper motions , . We further corrected for the mean peculiar motion of the masers, i.e. km/s (cf. Tab. 4, A5 in Ref. Reid et al. (2014)) and km/s (cf. Sec. 4.4 in Ref. Reid et al. (2014)). We exclude from the full data set all objects in (or coincident to objects in) Reid+ ’14; only one object remains: Ori GMRA.
Xu+ ’13 Xu et al. (2013)
From Tab. 4, we take the parallaxes , line-of-sight velocities and equatorial proper motions , corresponding to objects towards different Galactic coordinates and add a virial motion of km/s in quadrature to following Sec. 3 in Ref. Reid et al. (2009). The parallaxes are converted to distances through , while the equatorial proper motions are converted to Galactic proper motions , . The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We further corrected for the mean peculiar motion of the masers, i.e. km/s (cf. Tab. 4, A5 in Ref. Reid et al. (2014)) and km/s (cf. Sec. 4.4 in Ref. Reid et al. (2014)). We exclude from the full data set all objects in (or coincident to objects in) Reid+ ’14 and Honma+ ’12; only one maser remains: DoAr21/Ophiuchus.
Bobylev & Bajkova ’13 Bobylev and Bajkova (2013)
From Tab. 1, we take the parallaxes , line-of-sight velocities and equatorial proper motions , corresponding to objects towards different equatorial coordinates (then converted to Galactic coordinates ) and add a virial motion of km/s in quadrature to following Sec. 3 in Ref. Reid et al. (2009). The parallaxes are converted to distances through , while the equatorial proper motions are converted to Galactic proper motions , . The reported velocities correspond to a peculiar solar motion km/s, which is consistently subtracted off in our compilation. We further corrected for the mean peculiar motion of the masers, i.e. km/s (cf. Tab. 4, A5 in Ref. Reid et al. (2014)) and km/s (cf. Sec. 4.4 in Ref. Reid et al. (2014)). We exclude from the full data set all objects in (or coincident to objects in) Reid+ ’14, Honma+ ’12, Stepanishchev & Bobylev ’11 and Xu+ ’13.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Freeman (1970) K. C. Freeman, Ap J 160 , 811 (1970).
- 2Roberts and Rots (1973) M. S. Roberts and A. H. Rots, A&A 26 , 483 (1973).
- 3Bosma (1978) A. Bosma, Ph.D. thesis, Ph D Thesis, Groningen Univ., (1978) (1978).
- 4Rubin et al. (1978) V. C. Rubin, N. Thonnard, and W. K. Ford, Jr., Ap J 225 , L 107 (1978).
- 5Persic and Salucci (1995) M. Persic and P. Salucci, Astrophys. J. Suppl. 99 , 501 (1995), eprint astro-ph/9502091.
- 6Blitz (1979) L. Blitz, Ap J 231 , L 115 (1979).
- 7Burton and Gordon (1978) W. B. Burton and M. A. Gordon, A&A 63 , 7 (1978).
- 8Clemens (1985) D. P. Clemens, Ap J 295 , 422 (1985).
