Galactic rotation from Cepheids with Gaia DR2 and effects of non-axisymmetry
01 Jan 2019-Monthly Notices of the Royal Astronomical Society (Oxford Academic)-Vol. 482, Iss: 1, pp 40-51
TL;DR: In this paper, a simple axisymmetric disc model was applied to 218 Galactic Cepheids whose accurate measurements of the distance and velocities are obtained by cross-matching an existing CEPheids catalogue with the Gaia DR2 data.
Abstract: We apply a simple axisymmetric disc model to 218 Galactic Cepheids whose accurate measurements of the distance and velocities are obtained by cross-matching an existing Cepheids catalogue with the Gaia DR2 data. Our model fit determines the "local centrifugal speed", $V_\mathrm{c}$ $-$ defined as the rotation speed required to balance the local radial gravitational force $-$ at the Sun's location to be $V_{c}(R_0)=236\pm 3$ km s$^{-1}$ and the Sun's azimuthal and radial peculiar motions to be $V_{\odot}=12.4\pm0.7$ km s$^{-1}$ and $U_{\odot}=7.7\pm0.9$ km s$^{-1}$, respectively. These results are obtained with strong priors on the solar radius, $R_0=8.2\pm0.1$ kpc, and Sun's angular rotation velocity, $\Omega_{\odot}=30.24\pm0.12$ km s$^{-1}$ kpc$^{-1}$. We also applied the axisymmetric model to mock data from an N-body/hydrodynamic simulation of a Milky Way-like galaxy with a bar and spiral arms. We find that our axisymmetric model fit to the young stars recovers the local centrifugal speed reasonably well, even in the face of significant non-axisymmetry. However, the local centrifugal speed determined from our Cepheid sample could suffer from systematic uncertainty as large as 6 km s$^{-1}$.
Citations
More filters
TL;DR: In this article, the origin of phase-space substructures revealed by the second Gaia data release in the disc of the Milky Way, such as the ridges in the $Vphi$-$r$-$V_r$ space and the streams in the$Vphi}$-V-r$ plane, was explored.
Abstract: We explore the origin of phase-space substructures revealed by the second Gaia data release in the disc of the Milky Way, such as the ridges in the $V_{\phi}$-$r$ plane, the undulations in the $V_{\phi}$-$r$-$V_r$ space and the streams in the $V_{\phi}$-$V_r$ plane. We use a collisionless N-body simulation with co-spatial thin and thick discs, along with orbit integration, to study the orbital structure close to the Outer Lindblad Resonance (OLR) of the bar. We find that a prominent, long-lived ridge is formed in the $V_{\phi}$-$r$ plane due to the OLR which translates to streams in the $V_{\phi}$-$V_r$ plane and examine which closed periodic and trapped librating orbits are responsible for these features. We find that orbits which carry out small librations around the $x_1(1)$ family are preferentially found at negative $V_r$, giving rise to a `horn'-like feature, while orbits with larger libration amplitudes, trapped around the $x_1(2)$ and $x_1(1)$ families, constitute the positive $V_r$ substructure, i.e. the Hercules-like feature. This changing libration amplitude of orbits will translate to a changing ratio of thin/thick disc stars, which could have implications on the metallicity distribution in this plane. We find that a scenario in which the Sun is placed close to the OLR gives rise to a strong asymmetry in $V_r$ in the $V_{\phi}$-$V_r$ plane (i.e. Hercules vs. `the horn') and subsequently to undulations in the $V_{\phi}$-$r$-$V_r$ space. We also explore a scenario in which the Sun is placed closer to the bar corotation and find that the bar perturbation $alone$ cannot give rise to the these features.
71 citations
TL;DR: In this paper, the authors used a sample of 773 Classical Cepheids with precise distances based on mid-infrared period-luminosity relations coupled with proper motions and radial velocities from Gaia to construct the accurate rotation curve of the Milky Way up to the distance of ~20 kpc from the Galactic center.
Abstract: Flat rotation curves of spiral galaxies are considered as an evidence for dark matter, but the rotation curve of the Milky Way is difficult to measure. Various objects were used to track the rotation curve in the outer parts of the Galaxy, but most studies rely on incomplete kinematical information and inaccurate distances. Here, we use a sample of 773 Classical Cepheids with precise distances based on mid-infrared period-luminosity relations coupled with proper motions and radial velocities from Gaia to construct the accurate rotation curve of the Milky Way up to the distance of ~20 kpc from the Galactic center. We use a simple model of Galactic rotation to measure the rotation speed of the Sun Theta_0 = 233.6 +/- 2.8 km/s, assuming a prior on the distance to the Galactic center R_0 = 8.122 +/- 0.031 kpc from the Gravity Collaboration. The rotation curve at Galactocentric distances 4 12 kpc constructed so far.
68 citations
Cites background from "Galactic rotation from Cepheids wit..."
...…and velocities of some tracers: H ii regions (Fich et al. 1989; Brand & Blitz 1993), Cepheids (Pont et al. 1994, 1997; Metzger et al. 1998; Kawata et al. 2018), open clusters (Hron 1987), or planetary nebulae (Durand et al. 1998), but the current uncertainties are considerable (see Figure…...
[...]
TL;DR: In this paper, the authors used a sample of 773 Classical Cepheids with precise distances based on mid-infrared period-luminosity relations coupled with proper motions and radial velocities from Gaia to construct the accurate rotation curve of the Milky Way up to the distance of ~20 kpc from the Galactic center.
Abstract: Flat rotation curves of spiral galaxies are considered as an evidence for dark matter, but the rotation curve of the Milky Way is difficult to measure. Various objects were used to track the rotation curve in the outer parts of the Galaxy, but most studies rely on incomplete kinematical information and inaccurate distances. Here, we use a sample of 773 Classical Cepheids with precise distances based on mid-infrared period-luminosity relations coupled with proper motions and radial velocities from Gaia to construct the accurate rotation curve of the Milky Way up to the distance of ~20 kpc from the Galactic center. We use a simple model of Galactic rotation to measure the rotation speed of the Sun Theta_0 = 233.6 +/- 2.8 km/s, assuming a prior on the distance to the Galactic center R_0 = 8.122 +/- 0.031 kpc from the Gravity Collaboration. The rotation curve at Galactocentric distances 4 12 kpc constructed so far.
56 citations
TL;DR: In this article, IB was partially supported by Science and Technology Facilities Council studentship 1506672, which was used to train a science and technology studentship student in the UK.
Abstract: Funding: IB was partially supported by Science and Technology Facilities Council studentship 1506672.
54 citations
TL;DR: In this article, two kinematical methods are adopted to measure the Galactic rotation curve in the Galactocentric distance range of $4\lesssim R_{\rm GC} \lesssim 19$ kpc.
Abstract: With the increasing numbers of large stellar survey projects, the quality and quantity of excellent tracers to study the Milky Way is rapidly growing, one of which is the classical Cepheids. Classical Cepheids are high precision standard candles with very low typical uncertainties ($<$ 3%) available via the mid-infrared period-luminosity relation. About 3500 classical Cepheids identified from OGLE, ASAS-SN, Gaia, WISE and ZTF survey data have been analyzed in this work, and their spatial distributions show a clear signature of Galactic warp. Two kinematical methods are adopted to measure the Galactic rotation curve in the Galactocentric distance range of $4\lesssim R_{\rm GC} \lesssim 19$ kpc. Gently declining rotation curves are derived by both the proper motion (PM) method and 3-dimensional velocity vector (3DV) method. The largest sample of classical Cepheids with most accurate 6D phase-space coordinates available to date are modeled in the 3DV method, and the resulting rotation curve is found to decline at the relatively smaller gradient of ($-1.33\pm0.1$) ${\rm km\,s^{-1}\,kpc^{-1}}$. Comparing to results from the PM method, a higher rotation velocity (($232.5\pm0.83$) ${\rm km\,s^{-1}}$) is derived at the position of Sun in the 3DV method. The virial mass and local dark matter density are estimated from the 3DV method which is the more reliable method, $M_{\rm vir} = (0.822\pm0.052)\times 10^{12}\,M_\odot$ and $\rho_{\rm DM,\odot} = 0.33\pm0.03$ GeV ${\rm cm^{-3}}$, respectively.
42 citations
Cites background or result from "Galactic rotation from Cepheids wit..."
...Thus, they are taken as excellent objects to study the structure, kine- matics and dynamics of the Galaxy, such as RR lyrae stars (Ablimit & Zhao 2017; Medina et al. 2018; Ablimit & Zhao 2018; Utkin et al. 2018; Wegg et al. 2018) and Cepheids (e.g., Kawata et al. 2018)....
[...]
...The effect of the asymmetric drift is not considered in the calculation of this work due to the very small systematic uncertainty it causes (e.g., estimated as ±0.28 km s−1 by Kawata et al. 2018)....
[...]
...The Sun’s rotation velocity obtained from the 3D velocity method is, within uncertainties, consistent with the relatively higher values reported by Metezger et al. (1998), Reid et al.(2014), Kawata et al.(2018) and Mróz et al. (2019)....
[...]
...It is well known that the motions of stars are affected by Galactic substructures (e.g. Grand, Kawata & Cropper 2014; Bovy et al.2015; Kawata et al. 2018; Martinez-Medina et al. 2019)....
[...]
References
More filters
TL;DR: The emcee algorithm as mentioned in this paper is a Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010).
Abstract: We introduce a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010). The code is open source and has already been used in several published projects in the astrophysics literature. The algorithm behind emcee has several advantages over traditional MCMC sampling methods and it has excellent performance as measured by the autocorrelation time (or function calls per independent sample). One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to ~N2 for a traditional algorithm in an N-dimensional parameter space. In this document, we describe the algorithm and the details of our implementation. Exploiting the parallelism of the ensemble method, emcee permits any user to take advantage of multiple CPU cores without extra effort. The code is available online at http://dan.iel.fm/emcee under the GNU General Public License v2.
8,805 citations
Anthony G. A. Brown1, Antonella Vallenari1, T. Prusti1, J. H. J. de Bruijne1 +449 more•Institutions (1)
TL;DR: The second Gaia data release, Gaia DR2 as mentioned in this paper, is a major advance with respect to Gaia DR1 in terms of completeness, performance, and richness of the data products.
Abstract: Context. We present the second Gaia data release, Gaia DR2, consisting of astrometry, photometry, radial velocities, and information on astrophysical parameters and variability, for sources brighter than magnitude 21. In addition epoch astrometry and photometry are provided for a modest sample of minor planets in the solar system. Aims: A summary of the contents of Gaia DR2 is presented, accompanied by a discussion on the differences with respect to Gaia DR1 and an overview of the main limitations which are still present in the survey. Recommendations are made on the responsible use of Gaia DR2 results. Methods: The raw data collected with the Gaia instruments during the first 22 months of the mission have been processed by the Gaia Data Processing and Analysis Consortium (DPAC) and turned into this second data release, which represents a major advance with respect to Gaia DR1 in terms of completeness, performance, and richness of the data products. Results: Gaia DR2 contains celestial positions and the apparent brightness in G for approximately 1.7 billion sources. For 1.3 billion of those sources, parallaxes and proper motions are in addition available. The sample of sources for which variability information is provided is expanded to 0.5 million stars. This data release contains four new elements: broad-band colour information in the form of the apparent brightness in the GBP (330-680 nm) and GRP (630-1050 nm) bands is available for 1.4 billion sources; median radial velocities for some 7 million sources are presented; for between 77 and 161 million sources estimates are provided of the stellar effective temperature, extinction, reddening, and radius and luminosity; and for a pre-selected list of 14 000 minor planets in the solar system epoch astrometry and photometry are presented. Finally, Gaia DR2 also represents a new materialisation of the celestial reference frame in the optical, the Gaia-CRF2, which is the first optical reference frame based solely on extragalactic sources. There are notable changes in the photometric system and the catalogue source list with respect to Gaia DR1, and we stress the need to consider the two data releases as independent. Conclusions: Gaia DR2 represents a major achievement for the Gaia mission, delivering on the long standing promise to provide parallaxes and proper motions for over 1 billion stars, and representing a first step in the availability of complementary radial velocity and source astrophysical information for a sample of stars in the Gaia survey which covers a very substantial fraction of the volume of our galaxy.
8,308 citations
TL;DR: This document introduces a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010).
Abstract: We introduce a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010). The code is open source and has already been used in several published projects in the astrophysics literature. The algorithm behind emcee has several advantages over traditional MCMC sampling methods and it has excellent performance as measured by the autocorrelation time (or function calls per independent sample). One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to $\sim N^2$ for a traditional algorithm in an N-dimensional parameter space. In this document, we describe the algorithm and the details of our implementation and API. Exploiting the parallelism of the ensemble method, emcee permits any user to take advantage of multiple CPU cores without extra effort. The code is available online at this http URL under the MIT License.
5,293 citations
[...]
T. Prusti1, J. H. J. de Bruijne1, Anthony G. A. Brown2, Antonella Vallenari3 +621 more•Institutions (93)
TL;DR: Gaia as discussed by the authors is a cornerstone mission in the science programme of the European Space Agency (ESA). The spacecraft construction was approved in 2006, following a study in which the original interferometric concept was changed to a direct-imaging approach.
Abstract: Gaia is a cornerstone mission in the science programme of the EuropeanSpace Agency (ESA). The spacecraft construction was approved in 2006, following a study in which the original interferometric concept was changed to a direct-imaging approach. Both the spacecraft and the payload were built by European industry. The involvement of the scientific community focusses on data processing for which the international Gaia Data Processing and Analysis Consortium (DPAC) was selected in 2007. Gaia was launched on 19 December 2013 and arrived at its operating point, the second Lagrange point of the Sun-Earth-Moon system, a few weeks later. The commissioning of the spacecraft and payload was completed on 19 July 2014. The nominal five-year mission started with four weeks of special, ecliptic-pole scanning and subsequently transferred into full-sky scanning mode. We recall the scientific goals of Gaia and give a description of the as-built spacecraft that is currently (mid-2016) being operated to achieve these goals. We pay special attention to the payload module, the performance of which is closely related to the scientific performance of the mission. We provide a summary of the commissioning activities and findings, followed by a description of the routine operational mode. We summarise scientific performance estimates on the basis of in-orbit operations. Several intermediate Gaia data releases are planned and the data can be retrieved from the Gaia Archive, which is available through the Gaia home page.
5,164 citations
Related Papers (5)
[...]