A dynamically young and perturbed Milky Way disk (original) (raw)
- Letter
- Published: 19 September 2018
- A. Helmi2,
- M. Romero-Gómez1,
- D. Katz3,
- C. Babusiaux3,4,
- R. Drimmel5,
- D. W. Evans6,
- F. Figueras1,
- E. Poggio5,7,
- C. Reylé8,
- A. C. Robin8,
- G. Seabroke9 &
- …
- C. Soubiran10
Nature volume 561, pages 360–362 (2018)Cite this article
- 8488 Accesses
- 348 Citations
- 664 Altmetric
- Metrics details
Subjects
Abstract
The evolution of the Milky Way disk, which contains most of the stars in the Galaxy, is affected by several phenomena. For example, the bar and the spiral arms of the Milky Way induce radial migration of stars1 and can trap or scatter stars close to orbital resonances2. External perturbations from satellite galaxies can also have a role, causing dynamical heating of the Galaxy3, ring-like structures in the disk4 and correlations between different components of the stellar velocity5. These perturbations can also cause ‘phase wrapping’ signatures in the disk6,7,8,9, such as arched velocity structures in the motions of stars in the Galactic plane. Some manifestations of these dynamical processes have already been detected, including kinematic substructure in samples of nearby stars10,11,12, density asymmetries and velocities across the Galactic disk that differ from the axisymmetric and equilibrium expectations13, especially in the vertical direction11,14,15,16, and signatures of incomplete phase mixing in the disk7,12,17,18. Here we report an analysis of the motions of six million stars in the Milky Way disk. We show that the phase-space distribution contains different substructures with various morphologies, such as snail shells and ridges, when spatial and velocity coordinates are combined. We infer that the disk must have been perturbed between 300 million and 900 million years ago, consistent with estimates of the previous pericentric passage of the Sagittarius dwarf galaxy. Our findings show that the Galactic disk is dynamically young and that modelling it as time-independent and axisymmetric is incorrect.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 51 print issues and online access
$199.00 per year
only $3.90 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Additional access options:
Similar content being viewed by others
Data availability
The datasets used and analysed for this study are derived from data available in the public Gaia archive (https://gea.esac.esa.int/archive). The Bayesian distances for the Gaia sources with radial velocity37are available at http://www.astro.lu.se/~paul/GaiaDR2_RV_star_distance.csv.gz. The rest of the relevant datasets and toy models are available from the corresponding author on reasonable request.
References
- Sellwood, J. A. & Binney, J. J. Radial mixing in galactic discs. Mon. Not. R. Astron. Soc. 336, 785–796 (2002).
Article ADS Google Scholar - Contopoulos, G. & Grosbol, P. Stellar dynamics of spiral galaxies: nonlinear effects at the 4/1 resonance. Astron. Astrophys. 155, 11–23 (1986).
ADS MATH Google Scholar - Quinn, P. J., Hernquist, L. & Fullagar, D. P. Heating of galactic disks by mergers. Astrophys. J. 403, 74–93 (1993).
Article ADS Google Scholar - Purcell, C. W., Bullock, J. S., Tollerud, E. J., Rocha, M. & Chakrabarti, S. The Sagittarius impact as an architect of spirality and outer rings in the Milky Way. Nature 477, 301–303 (2011).
Article ADS CAS Google Scholar - D’Onghia, E., Madau, P., Vera-Ciro, C., Quillen, A. & Hernquist, L. Excitation of coupled stellar motions in the Galactic disk by orbiting satellites. Astrophys. J. 823, 4 (2016).
Article ADS Google Scholar - Fux, R. Order and chaos in the local disc stellar kinematics induced by the Galactic bar. Astron. Astrophys. 373, 511–535 (2001).
Article ADS Google Scholar - Minchev, I. et al. Is the Milky Way ringing? The hunt for high-velocity streams. Mon. Not. R. Astron. Soc. 396, L56–L60 (2009).
Article ADS Google Scholar - Gómez, F. A., Minchev, I., Villalobos, Á., O’Shea, B. W. & Williams, M. E. K. Signatures of minor mergers in Milky Way like disc kinematics: ringing revisited. Mon. Not. R. Astron. Soc. 419, 2163–2172 (2012).
Article ADS Google Scholar - de la Vega, A., Quillen, A. C., Carlin, J. L., Chakrabarti, S. & D’Onghia, E. Phase wrapping of epicyclic perturbations in the Wobbly galaxy. Mon. Not. R. Astron. Soc. 454, 933–945 (2015).
Article ADS Google Scholar - Eggen, O. J. Star streams and Galactic structure. Astron. J. 112, 1595–1613 (1996).
Article ADS CAS Google Scholar - Dehnen, W. The distribution of nearby stars in velocity space inferred from HIPPARCOS data. Astron. J. 115, 2384–2396 (1998).
Article ADS Google Scholar - Gaia Collaboration. Gaia data release 2: mapping the Milky Way disc kinematics. Astron. Astrophys. 616, A11 (2018).
Article Google Scholar - Siebert, A. et al. Detection of a radial velocity gradient in the extended local disc with RAVE. Mon. Not. R. Astron. Soc. 412, 2026–2032 (2011).
Article ADS Google Scholar - Widrow, L. M., Gardner, S., Yanny, B., Dodelson, S. & Chen, H.-Y. Galactoseismology: discovery of vertical waves in the Galactic disk. Astrophys. J. 750, L41 (2012).
Article ADS Google Scholar - Schönrich, R. & Dehnen, W. Warp, waves, and wrinkles in the Milky Way. Mon. Not. R. Astron. Soc. 478, 3809–3824 (2018).
Article ADS Google Scholar - Quillen, A. C. et al. The GALAH survey: stellar streams and how stellar velocity distributions vary with Galactic longitude, hemisphere and metallicity. Mon. Not. R. Astron. Soc. 478, 228–254 (2018).
Article ADS CAS Google Scholar - Gómez, F. A. et al. Signatures of minor mergers in the Milky Way disc – I. The SEGUE stellar sample. Mon. Not. R. Astron. Soc. 423, 3727–3739 (2012).
Article ADS Google Scholar - Monari, G. et al. Coma Berenices: the first evidence for incomplete vertical phase-mixing in local velocity space with RAVE—confirmed with Gaia DR2. Res. Notes AAS 2, 32 (2018).
Article ADS Google Scholar - Gaia Collaboration. Gaia data release 2: summary of the contents and survey properties. Astron. Astrophys. 616, A1 (2018).
Article Google Scholar - Tremaine, S. The geometry of phase mixing. Mon. Not. R. Astron. Soc. 307, 877–883 (1999).
Article ADS Google Scholar - Afshordi, N., Mohayaee, R. & Bertschinger, E. Hierarchical phase space structure of dark matter haloes: tidal debris, caustics, and dark matter annihilation. Phys. Rev. D 79, 083526 (2009).
Article ADS Google Scholar - Candlish, G. N. et al. Phase mixing due to the Galactic potential: steps in the position and velocity distributions of popped star clusters. Mon. Not. R. Astron. Soc. 437, 3702–3717 (2014).
Article ADS Google Scholar - Manfredi, G. & Feix, R. M. Theory and simulation of classical and quantum echoes. Phys. Rev. E 53, 6460–6470 (1996).
Article ADS CAS Google Scholar - Law, D. R. & Majewski, S. R. The Sagittarius dwarf galaxy: a model for evolution in a triaxial Milky Way halo. Astrophys. J. 714, 229–254 (2010).
Article ADS Google Scholar - Laporte, C. F. P., Johnston, K. V., Gómez, F. A., Garavito-Camargo, N. & Besla, G. The influence of Sagittarius and the Large Magellanic Cloud on the Milky Way galaxy. Mon. Not. R. Astron. Soc. https://doi.org/10.1093/mnras/sty1574 (2018).
Article ADS CAS Google Scholar - Monari, G., Kawata, D., Hunt, J. A. S. & Famaey, B. Tracing the Hercules stream with Gaia and LAMOST: new evidence for a fast bar in the Milky Way. Mon. Not. R. Astron. Soc. 466, L113–L117 (2017).
Article ADS Google Scholar - Michtchenko, T. A., Lépine, J. R. D., Barros, D. A. & Vieira, R. S. S. Combined dynamical effects of the bar and spiral arms in a Galaxy model. Application to the solar neighbourhood. Astron. Astrophys. 615, A10 (2018).
Article ADS Google Scholar - Fouvry, J.-B., Binney, J. & Pichon, C. Self-gravity, resonances, and orbital diffusion in stellar disks. Astrophys. J. 806, 117 (2015).
Article ADS Google Scholar - Luri, X. et al. Gaia data release 2: using Gaia parallaxes. Astron. Astrophys. 616, A9 (2018).
Article Google Scholar - Chen, B. et al. Stellar population studies with the SDSS. I. The vertical distribution of stars in the Milky Way. Astrophys. J. 553, 184–197 (2001).
Article ADS Google Scholar - Reid, M. J. et al. Trigonometric parallaxes of high mass star forming regions: the structure and kinematics of the Milky Way. Astrophys. J. 783, 130 (2014).
Article ADS Google Scholar - Schönrich, R. Galactic rotation and solar motion from stellar kinematics. Mon. Not. R. Astron. Soc. 427, 274–287 (2012).
Article ADS Google Scholar - Reid, M. J. & Brunthaler, A. The proper motion of Sagittarius A*. II. The mass of Sagittarius A*. Astrophys. J. 616, 872–884 (2004).
Article ADS CAS Google Scholar - Lindegren, L. et al. Gaia data release 2: the astrometric solution. Astron. Astrophys. 616, A2 (2018).
Article Google Scholar - Taylor, M. B. TOPCAT & STIL: Starlink table/VOTable processing software. ASP Conf. Ser. 347, 29–33 (2005).
ADS Google Scholar - Astraatmadja, T. L. & Bailer-Jones, C. A. L. Estimating distances from parallaxes. II. Performance of Bayesian distance estimators on a Gaia-like catalogue. Astrophys. J. 832, 137 (2016).
Article ADS Google Scholar - McMillan, P. J. Simple distance estimates for Gaia DR2 stars with radial velocities. Res. Notes AAS 2, 51 (2018).
Article ADS Google Scholar - Binney, J. & Tremaine, S. Galactic Dynamics 2nd edn (Princeton Univ. Press, Princeton, 2008).
MATH Google Scholar - Miyamoto, M. & Nagai, R. Three-dimensional models for the distribution of mass in galaxies. Publ. Astron. Soc. Jpn 27, 533–543 (1975).
ADS Google Scholar - Allen, C. & Santillan, A. An improved model of the galactic mass distribution for orbit computations. Rev. Mex. Astron. Astrofis. 22, 255–263 (1991).
ADS Google Scholar - Irrgang, A., Wilcox, B., Tucker, E. & Schiefelbein, L. Milky Way mass models for orbit calculations. Astron. Astrophys. 549, A137 (2013).
Article ADS Google Scholar - McMillan, P. J. The mass distribution and gravitational potential of the Milky Way. Mon. Not. R. Astron. Soc. 465, 76–94 (2017).
Article ADS CAS Google Scholar - Romero-Gómez, M., Figueras, F., Antoja, T., Abedi, H. & Aguilar, L. The analysis of realistic stellar Gaia mock catalogues – I. Red clump stars as tracers of the central bar. Mon. Not. R. Astron. Soc. 447, 218–233 (2015).
Article ADS Google Scholar - Ferrers, N. On the potentials of ellipsoids, ellipsoidal shells, elliptic laminae and elliptic rings of variable densities. QJ Pure Appl. Math 14, 1–22 (1877).
MATH Google Scholar - Eggen, O. J. Stellar groups. II. The ζ Herculis, ϵ Indi and 61 Cygni groups of high-velocity stars. Mon. Not. R. Astron. Soc. 118, 154–160 (1958).
Article ADS Google Scholar - Blaauw, A. Remarks on Local Structure and Kinematics. Symp. IAU 38, 199–204 (1970).
ADS Google Scholar - Skuljan, J., Hearnshaw, J. B. & Cottrell, P. L. Velocity distribution of stars in the solar neighbourhood. Mon. Not. R. Astron. Soc. 308, 731–740 (1999).
Article ADS Google Scholar - Antoja, T., Figueras, F., Fernández, D. & Torra, J. Origin and evolution of moving groups. I. Characterization in the observational kinematic-age-metallicity space. Astron. Astrophys. 490, 135–150 (2008).
Article ADS CAS Google Scholar
Acknowledgements
This work made use of data from ESA mission Gaia (https://www.cosmos.esa.int/gaia), which was processed by the Gaia Data Processing and Analysis Consortium (DPAC; https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC is provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This project received funding from the European Union’s Horizon 2020 research and innovation programme under Marie Skłodowska-Curie grant agreement number 745617. This work was supported by the MDM-2014-0369 of ICCUB (Unidad de Excelencia ‘María de Maeztu’) and the European Community’s Seventh Framework Programme (FP7/2007-2013) under grant agreement GENIUS FP7-606740. A.H. acknowledges financial support from a VICI grant from the Netherlands Organisation for Scientific Research (NWO). We acknowledge the MINECO (Spanish Ministry of Economy) through grants ESP2016-80079-C2-1-R (MINECO/FEDER, UE) and ESP2014-55996-C2-1-R (MINECO/FEDER, UE). This work been funded in part by the Agenzia Spaziale Italiana (ASI) through contract 2014-025-R.1.2015 through the Italian Istituto Nazionale di Astrofisica (INAF). E.P. acknowledges the financial support of the 2014 PhD fellowship programme of INAF.
Author information
Authors and Affiliations
- Institut de Ciències del Cosmos, Universitat de Barcelona (IEEC-UB), Barcelona, Spain
T. Antoja, M. Romero-Gómez & F. Figueras - Kapteyn Astronomical Institute, University of Groningen, Groningen, The Netherlands
A. Helmi - GEPI, Observatoire de Paris, Université PSL, CNRS, Meudon, France
D. Katz & C. Babusiaux - Université Grenoble Alpes, CNRS, IPAG, Grenoble, France
C. Babusiaux - INAF—Osservatorio Astrofisico di Torino, Pino Torinese, Italy
R. Drimmel & E. Poggio - Institute of Astronomy, University of Cambridge, Cambridge, UK
D. W. Evans - Università di Torino, Dipartimento di Fisica, Torino, Italy
E. Poggio - Institut UTINAM, CNRS UMR6213, Université Bourgogne Franche-Comté, OSU THETA Franche-Comté Bourgogne, Observatoire de Besançon, Besançon, France
C. Reylé & A. C. Robin - Mullard Space Science Laboratory, University College London, Dorking, UK
G. Seabroke - Laboratoire d’astrophysique de Bordeaux, Université Bordeaux, CNRS, Pessac, France
C. Soubiran
Authors
- T. Antoja
You can also search for this author inPubMed Google Scholar - A. Helmi
You can also search for this author inPubMed Google Scholar - M. Romero-Gómez
You can also search for this author inPubMed Google Scholar - D. Katz
You can also search for this author inPubMed Google Scholar - C. Babusiaux
You can also search for this author inPubMed Google Scholar - R. Drimmel
You can also search for this author inPubMed Google Scholar - D. W. Evans
You can also search for this author inPubMed Google Scholar - F. Figueras
You can also search for this author inPubMed Google Scholar - E. Poggio
You can also search for this author inPubMed Google Scholar - C. Reylé
You can also search for this author inPubMed Google Scholar - A. C. Robin
You can also search for this author inPubMed Google Scholar - G. Seabroke
You can also search for this author inPubMed Google Scholar - C. Soubiran
You can also search for this author inPubMed Google Scholar
Contributions
T.A. contributed to the sample preparation, analysed and interpreted the data, performed most of the modelling and wrote the paper together with A.H. A.H. also provided interpretation of the findings. M.R.-G. performed the simulation with the barred potential and contributed to sample preparation. D.K., C.B., R.D., D.W.E., F.F., E.P., C.R., A.C.R., G.S. and C.S. contributed to sample preparation and validation of the Gaia data. All authors reviewed the manuscript.
Corresponding author
Correspondence toT. Antoja.
Ethics declarations
Competing interests
: The authors declare no competing interests.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Velocities of the stars at the solar Galactocentric radius.
Two-dimensional histograms of combinations of radial, azimuthal and vertical Galactic cylindrical velocities for the stars in our sample of Gaia data located at 8.24 kpc < R < 8.44 kpc, in bins of 1 km s−1. V R and V ϕ are positive towards the Galactic anticentre and the direction of Galactic rotation, respectively. The darkness is proportional to the number of counts. a, Although the bimodality seen here, separating the Hercules stream from the rest of the distribution, was known45,46, as well as some other elongated structures in this velocity projection11,47,48, the numerous and thin arches are a new phenomenon revealed by Gaia data12. The semi-circular dotted line marks an arbitrary line of constant kinetic energy in the plane \({E}_{k}=({V}_{R}^{2}+{V}_{\phi }^{2})/2\), as predicted for the substructure generated in horizontal phase mixing7,8. b, The data have a box-like appearance, where the extent in V Z of the arches varies with V ϕ (arrows), probably created by the correlation between the spiral shape and V ϕ seen in Fig. 1c. c, Although some velocity asymmetries where noticed before in the V ϕ_–_V Z projection11 and attributed to the Galaxy warp, the sharp shell-like features involving V Z, especially at V Z ≈ −30 km s−1 and V Z ≈ 25 km s−1, were not previously evident. These shells are different projections of the snail shell pattern of Fig. 1a.
Extended Data Fig. 2 Location of the stars in the sample.
a, b, Two-dimensional histograms with bins of 0.05 kpc in the X_–_Y (a) and X_–_Z (b) projections of our sample of Gaia data. The dotted lines mark the selection of stars in the solar Galactic ring between radii of 8.24 kpc and 8.44 kpc. The Sun is located at (X, Y, Z) = (−8.34, 0, 0.027) kpc and the Galactic centre (GC) is marked with a black dot.
Extended Data Fig. 3 Modelled vertical positions and velocities of stars with time.
The plots show the snail shells created in the phase space evolution under an anharmonic potential. a–c, Phase-space evolution at different times (t = 0, 10, 100, 200, 1,000 Myr) for an ensemble of particles at a fixed Galactocentric radius of R = 8.5 kpc with an initial Gaussian distribution in Z(t = 0) with mean of −0.1 kpc and dispersion of 0.04 kpc and in V Z(t = 0) with mean of −2 km s−1 and dispersion of 1 km s−1. d–f, Same as a–c, but for a skewed normal distribution of initial radius with skewness of 10, location parameter of 8.4 kpc and scale parameter of 0.2 kpc. In all cases, the evolution is under an anharmonic oscillator derived from the expansion of a Miyamoto–Nagai disk for small Z. In a and d the stars are colour-coded by vertical period.
Extended Data Fig. 4 Vertical frequency for orbits in a Galaxy model.
a, b, Frequencies as a function of Galactocentric radius R computed in the updated model from ref. 41, colour coded by the vertical amplitude (a) and by the vertical velocity amplitude (b) of the orbits.
Extended Data Fig. 5 Position of the spiral turns in the vertical positions and velocities.
The Z_–_V Z plane for stars at Galactocentric radii of 8.24 kpc to 8.44 kpc, coloured as a function of median guiding radius R_g in bins of Δ_Z = 0.04 kpc and Δ_V_ Z = 2 km s−1, with vertical and horizontal lines showing the approximate locations of the observed snail shell (turn-around and mid-plane points).
Extended Data Table 1 Time estimates from the turn-around points of the spiral
Extended Data Table 2 Time estimates from the mid-plane points of the spiral
Rights and permissions
About this article
Cite this article
Antoja, T., Helmi, A., Romero-Gómez, M. et al. A dynamically young and perturbed Milky Way disk.Nature 561, 360–362 (2018). https://doi.org/10.1038/s41586-018-0510-7
- Received: 26 April 2018
- Accepted: 26 July 2018
- Published: 19 September 2018
- Issue Date: 20 September 2018
- DOI: https://doi.org/10.1038/s41586-018-0510-7
Keywords
This article is cited by
The Radcliffe Wave is oscillating
- Ralf Konietzka
- Alyssa A. Goodman
- Núria Miret-Roig
Nature (2024)
A slightly oblate dark matter halo revealed by a retrograde precessing Galactic disk warp
- Yang Huang
- Qikang Feng
- Haibo Yuan
Nature Astronomy (2024)
Chronos - take the pulse of our galactic neighbourhood
- Eric Michel
- Misha Haywood
- Andy Moya
Experimental Astronomy (2021)
All-sky visible and near infrared space astrometry
- David Hobbs
- Anthony Brown
- Mattia Vaccari
Experimental Astronomy (2021)
The recurrent impact of the Sagittarius dwarf on the star formation history of the Milky Way
- Tomás Ruiz-Lara
- Carme Gallart
- Santi Cassisi
Nature Astronomy (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.