A dynamically young and perturbed Milky Way disk - Nature
- ️Soubiran, C.
- ️Wed Sep 19 2018
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).
Contopoulos, G. & Grosbol, P. Stellar dynamics of spiral galaxies: nonlinear effects at the 4/1 resonance. Astron. Astrophys. 155, 11–23 (1986).
Quinn, P. J., Hernquist, L. & Fullagar, D. P. Heating of galactic disks by mergers. Astrophys. J. 403, 74–93 (1993).
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).
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).
Fux, R. Order and chaos in the local disc stellar kinematics induced by the Galactic bar. Astron. Astrophys. 373, 511–535 (2001).
Minchev, I. et al. Is the Milky Way ringing? The hunt for high-velocity streams. Mon. Not. R. Astron. Soc. 396, L56–L60 (2009).
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).
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).
Eggen, O. J. Star streams and Galactic structure. Astron. J. 112, 1595–1613 (1996).
Dehnen, W. The distribution of nearby stars in velocity space inferred from HIPPARCOS data. Astron. J. 115, 2384–2396 (1998).
Gaia Collaboration. Gaia data release 2: mapping the Milky Way disc kinematics. Astron. Astrophys. 616, A11 (2018).
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).
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).
Schönrich, R. & Dehnen, W. Warp, waves, and wrinkles in the Milky Way. Mon. Not. R. Astron. Soc. 478, 3809–3824 (2018).
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).
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).
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).
Gaia Collaboration. Gaia data release 2: summary of the contents and survey properties. Astron. Astrophys. 616, A1 (2018).
Tremaine, S. The geometry of phase mixing. Mon. Not. R. Astron. Soc. 307, 877–883 (1999).
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).
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).
Manfredi, G. & Feix, R. M. Theory and simulation of classical and quantum echoes. Phys. Rev. E 53, 6460–6470 (1996).
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).
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).
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).
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).
Fouvry, J.-B., Binney, J. & Pichon, C. Self-gravity, resonances, and orbital diffusion in stellar disks. Astrophys. J. 806, 117 (2015).
Luri, X. et al. Gaia data release 2: using Gaia parallaxes. Astron. Astrophys. 616, A9 (2018).
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).
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).
Schönrich, R. Galactic rotation and solar motion from stellar kinematics. Mon. Not. R. Astron. Soc. 427, 274–287 (2012).
Reid, M. J. & Brunthaler, A. The proper motion of Sagittarius A*. II. The mass of Sagittarius A*. Astrophys. J. 616, 872–884 (2004).
Lindegren, L. et al. Gaia data release 2: the astrometric solution. Astron. Astrophys. 616, A2 (2018).
Taylor, M. B. TOPCAT & STIL: Starlink table/VOTable processing software. ASP Conf. Ser. 347, 29–33 (2005).
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).
McMillan, P. J. Simple distance estimates for Gaia DR2 stars with radial velocities. Res. Notes AAS 2, 51 (2018).
Binney, J. & Tremaine, S. Galactic Dynamics 2nd edn (Princeton Univ. Press, Princeton, 2008).
Miyamoto, M. & Nagai, R. Three-dimensional models for the distribution of mass in galaxies. Publ. Astron. Soc. Jpn 27, 533–543 (1975).
Allen, C. & Santillan, A. An improved model of the galactic mass distribution for orbit computations. Rev. Mex. Astron. Astrofis. 22, 255–263 (1991).
Irrgang, A., Wilcox, B., Tucker, E. & Schiefelbein, L. Milky Way mass models for orbit calculations. Astron. Astrophys. 549, A137 (2013).
McMillan, P. J. The mass distribution and gravitational potential of the Milky Way. Mon. Not. R. Astron. Soc. 465, 76–94 (2017).
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).
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).
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).
Blaauw, A. Remarks on Local Structure and Kinematics. Symp. IAU 38, 199–204 (1970).
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).
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).
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 to T. 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. VR 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 VZ 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ϕ–VZ projection11 and attributed to the Galaxy warp, the sharp shell-like features involving VZ, especially at VZ ≈ −30 km s−1 and VZ ≈ 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 VZ(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–VZ plane for stars at Galactocentric radii of 8.24 kpc to 8.44 kpc, coloured as a function of median guiding radius Rg in bins of ΔZ = 0.04 kpc and ΔVZ = 2 km s−1, with vertical and horizontal lines showing the approximate locations of the observed snail shell (turn-around and mid-plane points).
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