Showing 1–50 of 207 results for author: Binney, J

Chemodynamical models of our Galaxy

Authors: James Binney, Eugene Vasiliev

Abstract: A chemodynamical model of our galaxy is fitted to data from DR17 of the APOGEE survey supplemented with data from the StarHorse catalogue and gaia DR3. Dynamically, the model is defined by action-based distribution functions for dark matter and six stellar components plus a gas disc. The gravitational potential jointly generated by the model's components is used to examine the galaxy's chemical co… ▽ More A chemodynamical model of our galaxy is fitted to data from DR17 of the APOGEE survey supplemented with data from the StarHorse catalogue and gaia DR3. Dynamically, the model is defined by action-based distribution functions for dark matter and six stellar components plus a gas disc. The gravitational potential jointly generated by the model's components is used to examine the galaxy's chemical composition within action space. The observational data probably cover all parts of action space that are populated by stars. The overwhelming majority of stars have angular momentum J_φふぁい>0 implying that they were born in the Galactic disc. High-alpha stars dominate in a region that is sharply bounded by J_φふぁい\la J_φふぁい(solar). Chemically the model is defined by giving each stellar component a Gaussian distribution in ([Fe/H],[Mg/Fe]) space about a mean that is a linear function of the actions. The model's 47 dynamical parameters are chosen to maximise the likelihood of the data given the model in 72 three-dimensional velocity spaces while its 70 chemical parameters are similarly chosen in five-dimensional chemo-dynamical space. The circular speed falls steadily from 237\kms at R=4\kpc to 218\kms at R=20\kpc. Dark matter contributes half the radial force on the Sun and has local density 0.011\msun\pc^{-3}, there being 24.5\msun\pc^{-2} in dark matter and 26.5\msun\pc^{-2} in stars within 1.1\kpc of the plane. △ Less

Submitted 27 October, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: 22 pp accepted by MNRAS

The Gaia-ESO Public Spectroscopic Survey: Motivation, implementation, GIRAFFE data processing, analysis, and final data products

Authors: G. Gilmore, S. Randich, C. C. Worley, A. Hourihane, A. Gonneau, G. G. Sacco, J. R. Lewis, L. Magrini, P. Francois, R. D. Jeffries, S. E. Koposov, A. Bragaglia, E. J. Alfaro, C. Allende Prieto, R. Blomme, A. J. Korn, A. C. Lanzafame, E. Pancino, A. Recio-Blanco, R. Smiljanic, S. Van Eck, T. Zwitter, T. Bensby, E. Flaccomio, M. J. Irwin , et al. (143 additional authors not shown)

Abstract: The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending a… ▽ More The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper (arXiv:2206.02901) introduces the survey results. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. The Gaia-ESO Survey obtained 202,000 spectra of 115,000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Comments: 38 pages. A&A in press

Journal ref: A&A 666, A120 (2022)

Self-consistent models of our Galaxy

Authors: James Binney, Eugene Vasiliev

Abstract: A new class of models of stellar discs is introduced and used to build a self-consistent model of our Galaxy. The model is defined by the parameters that specify the action-based distribution functions (DFs) f(J) of four stellar discs (three thin-disc age cohorts and a thick disc), spheroidal bulge and spheroidal stellar and dark haloes. From these DFs plus a specified distribution of gas, we solv… ▽ More A new class of models of stellar discs is introduced and used to build a self-consistent model of our Galaxy. The model is defined by the parameters that specify the action-based distribution functions (DFs) f(J) of four stellar discs (three thin-disc age cohorts and a thick disc), spheroidal bulge and spheroidal stellar and dark haloes. From these DFs plus a specified distribution of gas, we solve for the densities of stars and dark matter and the potential they generate. The principal observational constraints are the kinematics of stars with Gaia RVS data and the density of stars in the column above the Sun. The model predicts the density and kinematics of stars and dark matter throughout the Galaxy. We determine the structure of the dark halo prior to the infall of baryons. A simple extension of the DFs of stellar components to include chemistry allows the model to reproduce the way the Galaxy's chemistry is observed to vary in the (R,z) plane. Surprisingly, the data indicate that high-alpha stars are confined to orbits with J_z >~ 50 kpc km/s. The code used to create the model is available on Github. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 19 pages submitted to MNRAS

The Gaia-ESO Public Spectroscopic Survey: Implementation, data products, open cluster survey, science, and legacy

Authors: S. Randich, G. Gilmore, L. Magrini, G. G. Sacco, R. J. Jackson, R. D. Jeffries, C. C. Worley, A. Hourihane, A. Gonneau, C. Viscasillas Vàzquez, E. Franciosini, J. R. Lewis, E. J. Alfaro, C. Allende Prieto, T. Bensby R. Blomme, A. Bragaglia, E. Flaccomio, P. François, M. J. Irwin, S. E. Koposov, A. J. Korn, A. C. Lanzafame, E. Pancino, A. Recio-Blanco, R. Smiljanic , et al. (139 additional authors not shown)

Abstract: In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars… ▽ More In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars using FLAMES on the ESO VLT (both Giraffe and UVES spectrographs), covering all the Milky Way populations, with a special focus on open star clusters. This article provides an overview of the survey implementation (observations, data quality, analysis and its success, data products, and releases), of the open cluster survey, of the science results and potential, and of the survey legacy. A companion article (Gilmore et al.) reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. The GES has determined homogeneous good-quality radial velocities and stellar parameters for a large fraction of its more than 110,000 unique target stars. Elemental abundances were derived for up to 31 elements for targets observed with UVES. Lithium abundances are delivered for about 1/3 of the sample. The analysis and homogenisation strategies have proven to be successful; several science topics have been addressed by the Gaia-ESO consortium and the community, with many highlight results achieved. The final catalogue has been released through the ESO archive at the end of May 2022, including the complete set of advanced data products. In addition to these results, the Gaia-ESO Survey will leave a very important legacy, for several aspects and for many years to come. △ Less

Submitted 6 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Astronomy and Astrophysics. 30 pages, 30 figures, 4 tables

Our Galaxy's youngest disc

Authors: Chengdong Li, James Binney

Abstract: We investigate the structure of our Galaxy's young stellar disc by fitting the distribution functions (DFs) of a new family to five-dimensional Gaia data for a sample of $47\,000$ OB stars. Tests of the fitting procedure show that the young disc's DF would be strongly constrained by Gaia data if the distribution of Galactic dust were accurately known. The DF that best fits the real data accurately… ▽ More We investigate the structure of our Galaxy's young stellar disc by fitting the distribution functions (DFs) of a new family to five-dimensional Gaia data for a sample of $47\,000$ OB stars. Tests of the fitting procedure show that the young disc's DF would be strongly constrained by Gaia data if the distribution of Galactic dust were accurately known. The DF that best fits the real data accurately predicts the kinematics of stars at their observed locations, but it predicts the spatial distribution of stars poorly, almost certainly on account of errors in the best-available dust map. We argue that dust models could be greatly improved by modifying the dust model until the spatial distribution of stars predicted by a DF agreed with the data. The surface density of OB stars is predicted to peak at $R\simeq5.5\mbox{kpc}$, slightly outside the reported peak in the surface density of molecular gas; we suggest that the latter radius may have been under-estimated through the use of poor kinematic distances. The velocity distributions predicted by the best-fit DF for stars with measured line-of-sight velocities $v_\parallel$ reveal that the outer disc is disturbed at the level of 10 $\mbox{km}~\mbox{s}^{-1}$ in agreement with earlier studies, and that the measured values of $v_\parallel$ have significant contributions from the orbital velocities of binaries. Hence the outer disc is colder than it is sometimes reported to be. △ Less

Submitted 27 June, 2022; v1 submitted 3 May, 2022; originally announced May 2022.

Comments: 18 pages, 20 figures, accepted to publish on MNRAS

Modelling the stellar halo with RR-Lyrae stars

Authors: Chengdong Li, James Binney

Abstract: A seven-parameter distribution function (DF) is fitted to $20\,000$ RR-Lyrae stars for which only astrometric data are available. The observational data are predicted by the DF in conjunction with the gravitational potential of a self-consistent model Galaxy defined by DFs for the dark halo, the bulge and a four-component disc. Tests of the technique developed to deal with missing line-of-sight ve… ▽ More A seven-parameter distribution function (DF) is fitted to $20\,000$ RR-Lyrae stars for which only astrometric data are available. The observational data are predicted by the DF in conjunction with the gravitational potential of a self-consistent model Galaxy defined by DFs for the dark halo, the bulge and a four-component disc. Tests of the technique developed to deal with missing line-of-sight velocities show that adding such velocities tightens constraints on the DF only slightly. The recovered model of the RR-Lyrae population confirms that the population is flattened and has a strongly radially biased velocity distribution. At large radii its density profile tends to $ρろー\sim r^{-4.5}$ but no power law provides a good fit inside the solar sphere. The model is shown to provide an excellent fit to the data for stars brighter than $r=16.5$ but at certain longitudes it predicts too few faint stars at Galactocentric radii $\sim20\mbox{kpc}$, possibly signalling that the halo is not axisymmetric. The DF is used to predict the velocity distribution of BHB stars for which space velocities are available. The $z$ components are predicted successfully but too much anisotropy in the $v_Rv_φふぁい$ plane is expected. △ Less

Submitted 28 December, 2021; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 19 pages, 16 figures. Accepted by MNRAS

Probabilistic distribution functions

Authors: Jun Yan Lau, James Binney

Abstract: Observed clusters should be modelled by considering the distribution function to be a random variable that quantifies the degree of excitation of the system's normal modes. A system of canonical coordinates for the space of DFs is identified so DFs can be weighted in a consistent way. Observed clusters should be modelled by considering the distribution function to be a random variable that quantifies the degree of excitation of the system's normal modes. A system of canonical coordinates for the space of DFs is identified so DFs can be weighted in a consistent way. △ Less

Submitted 13 July, 2021; originally announced July 2021.

Comments: 5 pages. Submitted to MNRAS

Modes of a stellar system II: non-ergodic systems

Authors: Jun Yan Lau, James Binney

Abstract: An equation is derived for the energy of a small disturbance in a system that is generated by a distribution function (DF) of the form $f(\vJ)$ -- most galaxies and star clusters can be closely approximated by such a DF. The theory of van Kampen modes is extended to such general systems. An inner product on the space of DFs is defined such that the energy of a disturbance is its norm under this pr… ▽ More An equation is derived for the energy of a small disturbance in a system that is generated by a distribution function (DF) of the form $f(\vJ)$ -- most galaxies and star clusters can be closely approximated by such a DF. The theory of van Kampen modes is extended to such general systems. An inner product on the space of DFs is defined such that the energy of a disturbance is its norm under this product. It is shown that van Kampen modes that differ in frequency are then orthogonal, with the consequence that the energies of van Kampen modes are additive. Consequently, most of the insight into the dynamics of ergodic systems that was gained in a recent paper on the van Kampen modes of ergodic systems applies to real clusters and galaxies. △ Less

Submitted 14 July, 2021; v1 submitted 18 June, 2021; originally announced June 2021.

Comments: 6 pages, submitted to MNRAS

Modes of a stellar system I: ergodic systems

Authors: Jun Yan Lau, James Binney

Abstract: The excursions of star clusters and galaxies around statistical equilibria are studied. For a stable ergodic model Antonov's Hermitian operator on six-dimensional phase space has the normal modes as its eigenfunctions. The excitation energy of the system is just the sum of the (positive) energies associated with each normal mode. Formulae are given for the DFs of modes, which are of the type first… ▽ More The excursions of star clusters and galaxies around statistical equilibria are studied. For a stable ergodic model Antonov's Hermitian operator on six-dimensional phase space has the normal modes as its eigenfunctions. The excitation energy of the system is just the sum of the (positive) energies associated with each normal mode. Formulae are given for the DFs of modes, which are of the type first described by van Kampen rather than Landau, and Landau `modes' can be expressed as sums of van Kampen modes. Each van Kampen mode comprises the response of non-resonant stars to driving by the gravitational field of stars on a group of resonant tori, so its structure is sensitive to the degree of self gravity. The emergence of global distortions in N-body models when particles are started from an analytical equilibrium is explained in terms of the interplay of normal modes. The positivity of modal energies opens the way to modelling the thermal properties of clusters in close analogy with those of crystals. △ Less

Submitted 18 June, 2021; v1 submitted 14 April, 2021; originally announced April 2021.

Comments: 13 pages, Submitted to MNRAS

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- II: Stellar Atmospheric Parameters, Chemical Abundances and Distances

Authors: Matthias Steinmetz, Guillaume Guiglion, Paul J. McMillan, Gal Matijevic, Harry Enke, Georges Kordopatis, Tomaz Zwitter, Marica Valentini, Cristina Chiappini, Luca Casagrande, Jennifer Wojno, Borja Anguiano, Olivier Bienayme, Albert Bijaoui, James Binney, Donna Burton, Paul Cass, Patrick de Laverny, Kristin Fiegert, Kenneth Freeman, Jon P. Fulbright, Brad K. Gibson, Gerard Gilmore, Eva K. Grebel, Amina Helmi , et al. (36 additional authors not shown)

Abstract: We present part 2 of the 6th and final Data Release (DR6 or FDR) of the Radial Velocity Experiment (RAVE), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A) and span the complete time frame from the start of RAVE observations on 12 April 2003 to their… ▽ More We present part 2 of the 6th and final Data Release (DR6 or FDR) of the Radial Velocity Experiment (RAVE), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A) and span the complete time frame from the start of RAVE observations on 12 April 2003 to their completion on 4 April 2013. In the second of two publications, we present the data products derived from 518387 observations of 451783 unique stars using a suite of advanced reduction pipelines focussing on stellar atmospheric parameters, in particular purely spectroscopically derived stellar atmospheric parameters (Teff, log(g), and the overall metallicity), enhanced stellar atmospheric parameters inferred via a Bayesian pipeline using Gaia DR2 astrometric priors, and asteroseismically calibrated stellar atmospheric parameters for giant stars based on asteroseismic observations for 699 K2 stars. In addition, we provide abundances of the elements Fe, Al, and Ni, as well as an overall [alpha/Fe] ratio obtained using a new pipeline based on the GAUGUIN optimization method that is able to deal with variable signal-to-noise ratios. The RAVE DR6 catalogs are cross matched with relevant astrometric and photometric catalogs, and are complemented by orbital parameters and effective temperatures based on the infrared flux method. The data can be accessed via the RAVE Web site (http://rave-survey.org) or the Vizier database. △ Less

Submitted 9 June, 2020; v1 submitted 11 February, 2020; originally announced February 2020.

Comments: 65 pages, 33 figures, accepted for publication to AJ

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- I: Survey Description, Spectra and Radial Velocities

Authors: Matthias Steinmetz, Gal Matijevic, Harry Enke, Tomaz Zwitter, Guillaume Guiglion, Paul J. McMillan, Georges Kordopatis, Marica Valentini, Cristina Chiappini, Luca Casagrande, Jennifer Wojno, Borja Anguiano, Olivier Bienayme, Albert Bijaoui, James Binney, Donna Burton, Paul Cass, Patrick de Laverny, Kristin Fiegert, Kenneth Freeman, Jon P. Fulbright, Brad K. Gibson, Gerard Gilmore, Eva K. Grebel, Amina Helmi , et al. (37 additional authors not shown)

Abstract: The Radial Velocity Experiment (RAVE) is a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A). The 6th and final data release (DR6 or FDR) is based on 518387 observations of 451783 unique stars. RAVE observations were taken between 12 April 2003 and 4 Ap… ▽ More The Radial Velocity Experiment (RAVE) is a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A). The 6th and final data release (DR6 or FDR) is based on 518387 observations of 451783 unique stars. RAVE observations were taken between 12 April 2003 and 4 April 2013. Here we present the genesis, setup and data reduction of RAVE as well as wavelength-calibrated and flux-normalized spectra and error spectra for all observations in RAVE DR6. Furthermore, we present derived spectral classification and radial velocities for the RAVE targets, complemented by cross matches with Gaia DR2 and other relevant catalogs. A comparison between internal error estimates, variances derived from stars with more than one observing epoch and a comparison with radial velocities of Gaia DR2 reveals consistently that 68% of the objects have a velocity accuracy better than 1.4 km/s, while 95% of the objects have radial velocities better than 4.0 km/s. Stellar atmospheric parameters, abundances and distances are presented in subsequent publication. The data can be accessed via the RAVE Web (http://rave-survey.org) or the Vizier database. △ Less

Submitted 9 June, 2020; v1 submitted 11 February, 2020; originally announced February 2020.

Comments: 32 pages, 11 figures, accepted for publication to AJ

Angle-action variables for orbits trapped at a Lindblad resonance

Authors: James Binney

Abstract: The conventional approach to orbit trapping at Lindblad resonances via a pendulum equation fails when the parent of the trapped orbits is too circular. The problem is explained and resolved in the context of the Torus Mapper and a realistic Galaxy model. Tori are computed for orbits trapped at both the inner and outer Lindblad resonances of our Galaxy. At the outer Lindblad resonance, orbits are q… ▽ More The conventional approach to orbit trapping at Lindblad resonances via a pendulum equation fails when the parent of the trapped orbits is too circular. The problem is explained and resolved in the context of the Torus Mapper and a realistic Galaxy model. Tori are computed for orbits trapped at both the inner and outer Lindblad resonances of our Galaxy. At the outer Lindblad resonance, orbits are quasiperiodic and can be accurately fitted by torus mapping. At the inner Lindblad resonance, orbits are significantly chaotic although far from ergodic, and each orbit explores a small range of tori obtained by torus mapping. △ Less

Submitted 15 December, 2019; originally announced December 2019.

Comments: Submitted to MNRAS as companion paper to "Trapped orbits and solar-neighbourhood kinematics"

Trapped orbits and solar-neighbourhood kinematics

Authors: James Binney

Abstract: Torus mapping yields constants of motion for stars trapped at a resonance. Each such constant of motion yields a system of contours in velocity space at the Sun and neighbouring points. If Jeans' theorem applied to resonantly trapped orbits, the density of stars in velocity space would be equal at all intersections of any two contours. A quantitative measure of the violation of this principal is d… ▽ More Torus mapping yields constants of motion for stars trapped at a resonance. Each such constant of motion yields a system of contours in velocity space at the Sun and neighbouring points. If Jeans' theorem applied to resonantly trapped orbits, the density of stars in velocity space would be equal at all intersections of any two contours. A quantitative measure of the violation of this principal is defined and used to assess various pattern speeds for a model of the bar recently fitted to observations of interstellar gas. Trapping at corotation of a bar with pattern speed in the range 33-36 /Gyr is favoured and trapping at the outer Lindblad resonance is disfavoured. As one moves around the Sun the structure of velocity space varies quite rapidly, both as regards the observed star density and the zones of trapped orbits. The data seem consistent with trapping at corotation. △ Less

Submitted 15 December, 2019; originally announced December 2019.

Comments: Submitted to MNRAS

Modelling our Galaxy

Authors: James Binney

Abstract: Dynamical models will be key to exploitation of the incoming flood data for our Galaxy. Modelling techniques are reviewed with an emphasis on f(J) modelling. Dynamical models will be key to exploitation of the incoming flood data for our Galaxy. Modelling techniques are reviewed with an emphasis on f(J) modelling. △ Less

Submitted 5 September, 2019; originally announced September 2019.

Comments: 8 page text of invited review at IAU Symposium 353 `Galactic dynamics in the era of large surveys'

Journal ref: Proc. IAU 14 (2019) 101-108

Flattened stellar systems based on distribution functions depending on actions

Authors: Raffaele Pascale, James Binney, Carlo Nipoti

Abstract: We address an issue that arises when self-consistently flattened dynamical stellar systems are constructed by adopting a distribution function (DF) that depends on the action integrals. The velocity distribution at points on the symmetry axis is controlled by the the dependence of the DF on just one action, while at points off the symmetry axis two actions are involved. Consequently, the physical… ▽ More We address an issue that arises when self-consistently flattened dynamical stellar systems are constructed by adopting a distribution function (DF) that depends on the action integrals. The velocity distribution at points on the symmetry axis is controlled by the the dependence of the DF on just one action, while at points off the symmetry axis two actions are involved. Consequently, the physical requirement that the velocity distribution evolves continuously in the neighbourhood of the symmetry axis restricts the functional forms of acceptable DFs. An algorithm for conforming to this restriction is presented and used to construct a variety of flattened models. △ Less

Submitted 22 July, 2019; originally announced July 2019.

Comments: 9 pages, 7 figures. Submitted to MNRAS. Comments welcome

The shearing sheet and swing amplification revisited

Authors: James Binney

Abstract: The principal results of the classic analysis of the shearing sheet and swing amplification by Julian & Toomre (1966) are re-derived in a more accessible way and then used to gain a better quantitative understanding of the dynamics of stellar discs. The axisymmetric limit of the shearing sheet is derived and used to re-derive Kalnajs' 1965 dispersion relation and Toomre's 1964 stability criterion… ▽ More The principal results of the classic analysis of the shearing sheet and swing amplification by Julian & Toomre (1966) are re-derived in a more accessible way and then used to gain a better quantitative understanding of the dynamics of stellar discs. The axisymmetric limit of the shearing sheet is derived and used to re-derive Kalnajs' 1965 dispersion relation and Toomre's 1964 stability criterion for axisymmetric disturbances. Using the shearing sheet to revisit Toomre's important 1969 paper on the group velocity implied by Lin-Shu-Kalnajs dispersion relation, we discover that two rather than one wavepackets emerges inside corotation: one each side of the inner Lindblad resonance. Although LSK dispersion relation provides useful interpretations of both wavepackets, the shearing sheet highlights the limitations of the LSK approach to disc dynamics. Disturbances by no means avoid an annulus around corotation, as the LSK dispersion relation implies. While disturbances of the shearing sheet have a limited life in real space, they live on much longer in velocity space, which Gaia allows us to probe extensively. C++ code is provided to facilitate applications of winding spiral waves. △ Less

Submitted 27 June, 2019; originally announced June 2019.

Comments: 15 pp submitted to MNRAS

Relaxation of spherical stellar systems

Authors: Jun Yan Lau, James Binney

Abstract: 10,000 simulations of 1000-particle realisations of the same cluster are computed by direct force summation. Over three crossing times the original Poisson noise is amplified more than tenfold by self-gravity. The cluster's fundamental dipole mode is strongly excited by Poisson noise, and this mode makes a major contribution to driving diffusion of stars in energy. The diffusive flow through actio… ▽ More 10,000 simulations of 1000-particle realisations of the same cluster are computed by direct force summation. Over three crossing times the original Poisson noise is amplified more than tenfold by self-gravity. The cluster's fundamental dipole mode is strongly excited by Poisson noise, and this mode makes a major contribution to driving diffusion of stars in energy. The diffusive flow through action space is computed for the simulations and compared with the predictions of both local-scattering theory and the Balescu-Lenard (BL) equation. The predictions of local-scattering theory are qualitatively wrong because the latter neglects self-gravity. These results imply that local-scattering theory is of little value. Future work on cluster evolution should employ either N-body simulation or the BL equation. However, significant code development will be required to make use of the BL equation practicable. △ Less

Submitted 27 June, 2019; originally announced June 2019.

Comments: Submitted to MNRAS

Action-based models for dwarf spheroidal galaxies and globular clusters

Authors: Raffaele Pascale, James Binney, Carlo Nipoti, Lorenzo Posti

Abstract: A new family of self-consistent DF-based models of stellar systems is explored. The stellar component of the models is described by a distribution function (DF) depending on the action integrals, previously used to model the Fornax dwarf spheroidal galaxy (dSph). The stellar component may cohabit with either a dark halo, also described by a DF, or with a massive central black hole. In all cases we… ▽ More A new family of self-consistent DF-based models of stellar systems is explored. The stellar component of the models is described by a distribution function (DF) depending on the action integrals, previously used to model the Fornax dwarf spheroidal galaxy (dSph). The stellar component may cohabit with either a dark halo, also described by a DF, or with a massive central black hole. In all cases we solve for the model's self-consistent potential. Focussing on spherically symmetric models, we show how the stellar observables vary with the anisotropy prescribed by the DF, with the dominance and nature of the dark halo, and with the mass of the black hole. We show that precise fits to the observed surface brightness profiles of four globular clusters can be obtained for a wide range of prescribed velocity anisotropies. We also obtain precise fits to the observed projected densities of four dSphs. Finally, we present a three-component model of the Scupltor dSph with distinct DFs for the red and blue horizontal branch stars and the dark matter halo. △ Less

Submitted 17 April, 2019; originally announced April 2019.

Comments: submitted to MNRAS

Models of rotating coronae

Authors: Mattia C. Sormani, Emanuele Sobacchi, Gabriele Pezzulli, James Binney, Ralf S. Klessen

Abstract: Fitting equilibrium dynamical models to observational data is an essential step in understanding the structure of the gaseous hot haloes that surround our own and other galaxies. However, the two main categories of models that are used in the literature are poorly suited for this task: (i) simple barotropic models are analytic and can therefore be adjusted to match the observations, but are clearl… ▽ More Fitting equilibrium dynamical models to observational data is an essential step in understanding the structure of the gaseous hot haloes that surround our own and other galaxies. However, the two main categories of models that are used in the literature are poorly suited for this task: (i) simple barotropic models are analytic and can therefore be adjusted to match the observations, but are clearly unrealistic because the rotational velocity $v_φふぁい(R,z)$ does not depend on the distance $z$ from the galactic plane, while (ii) models obtained as a result of cosmological galaxy formation simulations are more realistic, but are impractical to fit to observations due to high computational cost. Here we bridge this gap by presenting a general method to construct axisymmetric baroclinic equilibrium models of rotating galactic coronae in arbitrary external potentials. We consider in particular a family of models whose equipressure surfaces in the $(R,z)$ plane are ellipses of varying axis ratio. These models are defined by two one-dimensional functions, the axial ratio of pressure $q_{\rm axis}(z)$ and the value of the pressure $P_{\rm axis}(z)$ along the galaxy's symmetry axis. These models can have a rotation speed $v_φふぁい(R,z)$ that realistically decreases as one moves away from the galactic plane, and can reproduce the angular momentum distribution found in cosmological simulations. The models are computationally cheap to construct and can thus be used in fitting algorithms. We provide a python code that given $q_{\rm axis}(z)$, $P_{\rm axis}(z)$ and $Φふぁい(R,z)$ returns $ρろー(R,z)$, $T(R,z)$, $P(R,z)$, $v_φふぁい(R,z)$. We show a few examples of these models using the Milky Way as a case study. △ Less

Submitted 10 September, 2018; originally announced September 2018.

Comments: Accepted for publication in MNRAS

The GALAH survey and Gaia DR2: dissecting the stellar disc's phase space by age, action, chemistry and location

Authors: Joss Bland-Hawthorn, Sanjib Sharma, Thor Tepper-Garcia, James Binney, Ken Freeman, Michael Hayden, Janez Kos, Gayandhi De Silva, Simon Ellis, Shourya Khanna, Geraint Lewis, Martin Asplund, Sven Buder, Andrew Casey, Valentina D'Orazi, Ly Duong, Jane Lin, Karin Lind, Sarah Martell, Melissa Ness, Jeff Simpson, Daniel Zucker, Tomaz Zwitter, Prajwal Kafle, Alice Quillen , et al. (2 additional authors not shown)

Abstract: We use the second data releases of the ESA Gaia astrometric survey and the high-resolution GALAH spectroscopic survey to analyse the structure of our Galaxy's disc components. With GALAH, we separate the alpha-rich and alpha-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars… ▽ More We use the second data releases of the ESA Gaia astrometric survey and the high-resolution GALAH spectroscopic survey to analyse the structure of our Galaxy's disc components. With GALAH, we separate the alpha-rich and alpha-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars in the zV_z phase plane, for both V_phi and V_R, and recover the remarkable "phase spiral" discovered by Gaia. We identify the anticipated quadrupole signature in zV_z of a tilted velocity ellipsoid for stars above and below the Galactic plane. By connecting our work with earlier studies, we show that the phase spiral is likely to extend well beyond the narrow solar neighbourhood cylinder in which it was found. The phase spiral is a signature of corrugated waves that propagate through the disc, and the associated non-equilibrium phase mixing. The radially asymmetric distribution of stars involved in the phase spiral reveals that the corrugation, which is mostly confined to the alpha-poor disc, grows in z-amplitude with increasing radius. We present new simulations of tidal disturbance of the Galactic disc by the Sagittarius (Sgr) dwarf. The effect on the zV_z phase plane lasts >2 Gyr but a subsequent disc crossing wipes out the coherent structure. We find that the phase spiral was excited < 0.5 Gyr ago by an object like Sgr with total mass 3 x 10^10 Msun (stripped down from 5 x 10^10 Msun when it first entered the halo) passing through the plane. △ Less

Submitted 17 January, 2019; v1 submitted 7 September, 2018; originally announced September 2018.

Comments: 27 pages, 24 figures, MNRAS (accepted, Jan. 2019)

The origin of the Gaia phase-plane spiral

Authors: James Binney, Ralph Schoenrich

Abstract: A simple model is presented of the formation of the spiral the (z,v_z) phase plane of solar-neighbourhood stars that was recently discovered in Gaia data. The key is that the frequency Omega_z at which stars oscillate vertically depends on angular momentum about the z axis in addition to the amplitude of the star's vertical oscillations. Spirals should form in both <v_phi> and <v_R> whenever a mas… ▽ More A simple model is presented of the formation of the spiral the (z,v_z) phase plane of solar-neighbourhood stars that was recently discovered in Gaia data. The key is that the frequency Omega_z at which stars oscillate vertically depends on angular momentum about the z axis in addition to the amplitude of the star's vertical oscillations. Spirals should form in both <v_phi> and <v_R> whenever a massive substructure, such as the Sgr dwarf galaxy, passes through the Galactic plane. The model yields similar spirals to those observed in both <v_phi> and <v_R>. The primary driver is the component of the tidal force that lies in the plane. We investigate the longevity of the spirals and the mass of the substructure, but the approximations inherent in the model make quantitative results unreliable. The work relies heavily on a self-consistent, multi-component model of our Galaxy produced by the AGAMA package for f(J) modelling. △ Less

Submitted 25 July, 2018; originally announced July 2018.

Comments: 6 pages submitted to MNRAS

Correlations between age, kinematics, and chemistry as seen by the RAVE survey

Authors: Jennifer Wojno, Georges Kordopatis, Matthias Steinmetz, Paul McMillan, James Binney, Benoit Famaey, Giacomo Monari, Ivan Minchev, Rosemary F. G. Wyse, Teresa Antoja, Arnaud Siebert, Ismael Carrillo, Joss Bland-Hawthorn, Eva K. Grebel, Tomaz Zwitter, Olivier Bienaymé, Brad Gibson, Andrea Kunder, Ulisse Munari, Julio Navarro, Quentin Parker, Warren Reid, George Seabroke

Abstract: We explore the connections between stellar age, chemistry, and kinematics across a Galactocentric distance of $7.5 < R\,(\mathrm{kpc}) < 9.0$, using a sample of $\sim 12\,000$ intermediate-mass (FGK) turnoff stars observed with the RAdial Velocity Experiment (RAVE) survey. The kinematics of this sample are determined using radial velocity measurements from RAVE, and parallax and proper motion meas… ▽ More We explore the connections between stellar age, chemistry, and kinematics across a Galactocentric distance of $7.5 < R\,(\mathrm{kpc}) < 9.0$, using a sample of $\sim 12\,000$ intermediate-mass (FGK) turnoff stars observed with the RAdial Velocity Experiment (RAVE) survey. The kinematics of this sample are determined using radial velocity measurements from RAVE, and parallax and proper motion measurements from the Tycho-Gaia Astrometric Solution (TGAS). In addition, ages for RAVE stars are determined using a Bayesian method, taking TGAS parallaxes as a prior. We divide our sample into young ($0 < τたう< 3$ Gyr) and old ($8 < τたう< 13$ Gyr) populations, and then consider different metallicity bins for each of these age groups. We find significant differences in kinematic trends of young and old, metal-poor and metal-rich, stellar populations. In particular, we find a strong metallicity dependence in the mean Galactocentric radial velocity as a function of radius ($\partial {V_{\rm R}}/\partial R$) for young stars, with metal-rich stars having a much steeper gradient than metal-poor stars. For $\partial {V_φふぁい}/\partial R$, young, metal-rich stars significantly lag the LSR with a slightly positive gradient, while metal-poor stars show a negative gradient above the LSR. We interpret these findings as correlations between metallicity and the relative contributions of the non-axisymmetries in the Galactic gravitational potential (the spiral arms and the bar) to perturb stellar orbits. △ Less

Submitted 17 April, 2018; originally announced April 2018.

Comments: 14 pages, 13 figures, accepted for publication in MNRAS

Revisiting relaxation in globular clusters

Authors: Chris Hamilton, Jean-Baptiste Fouvry, James Binney, Christophe Pichon

Abstract: The classical theory of cluster relaxation is unsatisfactory because it involves the Coulomb logarithm. The Balescu-Lenard (BL) equation provides a rigorous alternative that has no ill-defined parameter. Moreover, the BL equation, unlike classical theory, includes the cluster's self-gravity. A heuristic argument is given that indicates that relaxation does not occur predominantly through two-parti… ▽ More The classical theory of cluster relaxation is unsatisfactory because it involves the Coulomb logarithm. The Balescu-Lenard (BL) equation provides a rigorous alternative that has no ill-defined parameter. Moreover, the BL equation, unlike classical theory, includes the cluster's self-gravity. A heuristic argument is given that indicates that relaxation does not occur predominantly through two-particle scattering and is enhanced by self-gravity. The BL equation is adapted to a spherical system and used to estimate the flux through the action space of isochrone clusters with different velocity anisotropies. A range of fairly different secular behaviours is found depending on the fraction of radial orbits. Classical theory is also used to compute the corresponding classical fluxes. The BL and classical fluxes are very different because (a) the classical theory materially under-estimates the impact of large-scale collectively amplified fluctuations and (b) only the leading terms in an infinite sum for the BL flux are computed. A complete theory of cluster relaxation likely requires that the sum in the BL equation be decomposed into a sum over a finite number of small wavenumbers complemented by an integral over large wavenumbers analogous to classical theory. △ Less

Submitted 31 August, 2018; v1 submitted 29 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS. Replaced with updated manuscript including responses to referee's comments

Action-based dynamical models of dwarf spheroidal galaxies: application to Fornax

Authors: Raffaele Pascale, Lorenzo Posti, Carlo Nipoti, James Binney

Abstract: We present new dynamical models of dwarf spheroidal galaxies (dSphs) in which both the stellar component and the dark halo are described by analytic distribution functions that depend on the action integrals. In their most general form these distribution functions can represent axisymmetric and possibly rotating stellar systems. Here, as a first application, we model the Fornax dSph, limiting ours… ▽ More We present new dynamical models of dwarf spheroidal galaxies (dSphs) in which both the stellar component and the dark halo are described by analytic distribution functions that depend on the action integrals. In their most general form these distribution functions can represent axisymmetric and possibly rotating stellar systems. Here, as a first application, we model the Fornax dSph, limiting ourselves, for simplicity, to the non rotating, spherical case. The models are compared with state-of-the-art spectroscopic and photometric observations of Fornax, exploiting the knowledge of the line-of-sight velocity distribution of the models and accounting for the foreground contamination from the Milky Way. The model that best fits the structural and kinematic properties of Fornax has a cored dark halo, with core size $r_{\rm c}\simeq1.03$ kpc. The dark-to-luminous mass ratio is $(M_{\rm dm}/M_{\star})|_{R_{\rm eff}}\simeq9.6$ within the effective radius $R_{\rm eff} \simeq 0.62\,$kpc and $(M_{\rm dm}/M_{\star})|_{3 {\rm kpc}} \simeq 144$ within 3 kpc. The stellar velocity distribution is isotropic almost over the full radial range covered by the spectroscopic data and slightly radially anisotropic in the outskirts of the stellar distribution. The dark-matter annihilation $J$-factor and decay $D$-factor are, respectively, $\log_{10}(J$ $[$GeV$^2$ cm$^{-5}])\simeq18.34$ and $\log_{10}(D$ $[$GeV cm$^{-2}])\simeq18.55$, for integration angle $θしーた= 0.5^{\circ}$. This cored halo model of Fornax is preferred, with high statistical significance, to both models with a Navarro, Frenk and White dark halo and simple mass-follows-light models. △ Less

Submitted 9 July, 2018; v1 submitted 7 February, 2018; originally announced February 2018.

Comments: Accepted 2018 July 6; Received 2018 June 1; Submitted in original to MNRAS form 2018 February 5

Orbital tori for non-axisymmetric galaxies

Authors: James Binney

Abstract: Our Galaxy's bar makes the Galaxy's potential distinctly non-axisymmetric. All orbits are affected by non-axisymmetry, and significant numbers are qualitatively changed by being trapped at a resonance with the bar. Orbital tori are used to compute these effects. Thick-disc orbits are no less likely to be trapped by corotation or a Lindblad resonance than thin-disc orbits. Perturbation theory is us… ▽ More Our Galaxy's bar makes the Galaxy's potential distinctly non-axisymmetric. All orbits are affected by non-axisymmetry, and significant numbers are qualitatively changed by being trapped at a resonance with the bar. Orbital tori are used to compute these effects. Thick-disc orbits are no less likely to be trapped by corotation or a Lindblad resonance than thin-disc orbits. Perturbation theory is used to create non-axisymmetric orbital tori from standard axisymmetric tori, and both trapped and untrapped orbits are recovered to surprising accuracy. Code is added to the TorusModeller library that makes it as easy to manipulate non-axisymmetric tori as axisymmetric ones. The augmented TorusModeller is used to compute the velocity structure of the solar neighbourhood for bars of different pattern speeds and a simple action-based distribution function. The technique developed here can be applied to any non-axisymmetric potential that is stationary in a rotating from - hence also to classical spiral structure. △ Less

Submitted 31 October, 2017; originally announced October 2017.

Comments: 21 pp, 18 figs, accepted by MNRAS

Is the Milky Way still breathing? RAVE-Gaia streaming motions

Authors: I. Carrillo, I. Minchev, G. Kordopatis, M. Steinmetz, J. Binney, F. Anders, O. Bienaymé, J. Bland-Hawthorn, B. Famaey, K. C. Freeman, G. Gilmore, B. K. Gibson, E. K. Grebel, A. Helmi, A. Just, A. Kunder, P. McMillan, G. Monari, U. Munari, J. Navarro, Q. A. Parker, W. Reid, G. Seabroke, S. Sharma, A. Siebert , et al. (4 additional authors not shown)

Abstract: 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 comp… ▽ More 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. △ Less

Submitted 15 February, 2018; v1 submitted 10 October, 2017; originally announced October 2017.

Comments: 19 p., 17 fig., accepted MNRAS

Distribution functions for resonantly trapped orbits in the Galactic disc

Authors: Giacomo Monari, Benoit Famaey, Jean-Baptiste Fouvry, James Binney

Abstract: The present-day response of a Galactic disc stellar population to a non-axisymmetric perturbation of the potential has previously been computed through perturbation theory within the phase-space coordinates of the unperturbed axisymmetric system. Such an Eulerian linearized treatment however leads to singularities at resonances, which prevent quantitative comparisons with data. Here, we manage to… ▽ More The present-day response of a Galactic disc stellar population to a non-axisymmetric perturbation of the potential has previously been computed through perturbation theory within the phase-space coordinates of the unperturbed axisymmetric system. Such an Eulerian linearized treatment however leads to singularities at resonances, which prevent quantitative comparisons with data. Here, we manage to capture the behaviour of the distribution function (DF) at a resonance in a Lagrangian approach, by averaging the Hamiltonian over fast angle variables and re-expressing the DF in terms of a new set of canonical actions and angles variables valid in the resonant region. We then follow the prescription of Binney (2016), assigning to the resonant DF the time average along the orbits of the axisymmetric DF expressed in the new set of actions and angles. This boils down to phase-mixing the DF in terms of the new angles, such that the DF for trapped orbits only depends on the new set of actions. This opens the way to quantitatively fitting the effects of the bar and spirals to Gaia data in terms of distribution functions in action space. △ Less

Submitted 17 July, 2017; originally announced July 2017.

Comments: 8 pages, 7 figures, accepted for publication in MNRAS

Improved distances and ages for stars common to TGAS and RAVE

Authors: Paul J. McMillan, Georges Kordopatis, Andrea Kunder, James Binney, Jennifer Wojno, Tomaž Zwitter, Matthias Steinmetz, Joss Bland-Hawthorn, Brad K. Gibson, Gerard Gilmore, Eva K. Grebel, Amina Helmi, Ulisse Munari, Julio F. Navarro, Quentin A. Parker, George Seabroke, Fred Watson, Rosemary F. G. Wyse

Abstract: We combine parallaxes from the first Gaia data release with the spectrophotometric distance estimation framework for stars in the fifth RAVE survey data release. The combined distance estimates are more accurate than either determination in isolation - uncertainties are on average two times smaller than for RAVE-only distances (three times smaller for dwarfs), and 1.4 times smaller than TGAS paral… ▽ More We combine parallaxes from the first Gaia data release with the spectrophotometric distance estimation framework for stars in the fifth RAVE survey data release. The combined distance estimates are more accurate than either determination in isolation - uncertainties are on average two times smaller than for RAVE-only distances (three times smaller for dwarfs), and 1.4 times smaller than TGAS parallax uncertainties (two times smaller for giants). We are also able to compare the estimates from spectrophotometry to those from Gaia, and use this to assess the reliability of both catalogues and improve our distance estimates. We find that the distances to the lowest log g stars are, on average, overestimated and caution that they may not be reliable. We also find that it is likely that the Gaia random uncertainties are smaller than the reported values. As a byproduct we derive ages for the RAVE stars, many with relative uncertainties less than 20 percent. These results for 219566 RAVE sources have been made publicly available, and we encourage their use for studies that combine the radial velocities provided by RAVE with the proper motions provided by Gaia. A sample that we believe to be reliable can be found by taking only the stars with the flag notification flag_any=0 △ Less

Submitted 18 April, 2018; v1 submitted 14 July, 2017; originally announced July 2017.

Comments: MNRAS accepted. Updated to final version, reflecting referee reports. Title changed to better reflect scientific content. Data available through RAVE website: https://www.rave-survey.org/project/

A theoretical explanation for the Central Molecular Zone asymmetry

Authors: Mattia C. Sormani, Robin G. Tress, Matthew Ridley, Simon C. O. Glover, Ralf S. Klessen, James Binney, John Magorrian, Rowan Smith

Abstract: It has been known for more than thirty years that the distribution of molecular gas in the innermost 300 parsecs of the Milky Way, the Central Molecular Zone, is strongly asymmetric. Indeed, approximately three quarters of molecular emission comes from positive longitudes, and only one quarter from negative longitudes. However, despite much theoretical effort, the origin of this asymmetry has rema… ▽ More It has been known for more than thirty years that the distribution of molecular gas in the innermost 300 parsecs of the Milky Way, the Central Molecular Zone, is strongly asymmetric. Indeed, approximately three quarters of molecular emission comes from positive longitudes, and only one quarter from negative longitudes. However, despite much theoretical effort, the origin of this asymmetry has remained a mystery. Here we show that the asymmetry can be neatly explained by unsteady flow of gas in a barred potential. We use high-resolution 3D hydrodynamical simulations coupled to a state-of-the-art chemical network. Despite the initial conditions and the bar potential being point-symmetric with respect to the Galactic Centre, asymmetries develop spontaneously due to the combination of a hydrodynamical instability known as the "wiggle instability" and the thermal instability. The observed asymmetry must be transient: observations made tens of megayears in the past or in the future would often show an asymmetry in the opposite sense. Fluctuations of amplitude comparable to the observed asymmetry occur for a large fraction of the time in our simulations, and suggest that the present is not an exceptional moment in the life of our Galaxy. △ Less

Submitted 15 December, 2017; v1 submitted 12 July, 2017; originally announced July 2017.

Comments: Accepted for publication in MNRAS. Videos of the simulations are available at http://www.ita.uni-heidelberg.de/~mattia/download.html

Migration and kinematics in growing disc galaxies with thin and thick discs

Authors: Michael Aumer, James Binney, Ralph Schönrich

Abstract: We analyse disc heating and radial migration in N-body models of growing disc galaxies with thick and thin discs. Similar to thin-disc-only models, galaxies with appropriate non-axisymmetric structures reproduce observational constraints on radial disc heating in and migration to the Solar Neighbourhood (Snhd). The presence of thick discs can suppress non-axisymmetries and thus higher baryonic-to-… ▽ More We analyse disc heating and radial migration in N-body models of growing disc galaxies with thick and thin discs. Similar to thin-disc-only models, galaxies with appropriate non-axisymmetric structures reproduce observational constraints on radial disc heating in and migration to the Solar Neighbourhood (Snhd). The presence of thick discs can suppress non-axisymmetries and thus higher baryonic-to-dark matter fractions are required than in models that only have a thin disc. Models that are baryon-dominated to roughly the Solar radius R_0 are favoured, in agreement with data for the Milky Way. For inside-out growing discs, today's thick-disc stars at R_0 are dominated by outwards migrators. Whether outwards migrators are vertically hotter than non-migrators depends on the radial gradient of the thick disc vertical velocity dispersion. There is an effective upper boundary in angular momentum that thick disc stars born in the centre of a galaxy can reach by migration, which explains the fading of the high-alpha sequence outside R_0. Our models compare well to Snhd kinematics from RAVE-TGAS. For such comparisons it is important to take into account the azimuthal variation of kinematics at R ~ R_0 and biases from survey selection functions. The vertical heating of thin disc stars by giant molecular clouds is only mildly affected by the presence of thick discs. Our models predict higher vertical velocity dispersions for the oldest stars than found in the Snhd age-velocity dispersion relation, possibly because of measurement uncertainties or an underestimation of the number of old cold stars in our models. △ Less

Submitted 13 June, 2017; originally announced June 2017.

Comments: accepted for publication in MNRAS, 22 pages, 11 figures, 1 Table, Appendix

PLATO as it is: a legacy mission for Galactic archaeology

Authors: A. Miglio, C. Chiappini, B. Mosser, G. R. Davies, K. Freeman, L. Girardi, P. Jofre, D. Kawata, B. M. Rendle, M. Valentini, L. Casagrande, W. J. Chaplin, G. Gilmore, K. Hawkins, B. Holl, T. Appourchaux, K. Belkacem, D. Bossini, K. Brogaard, M. -J. Goupil, J. Montalban, A. Noels, F. Anders, T. Rodrigues, G. Piotto , et al. (80 additional authors not shown)

Abstract: 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 s… ▽ More 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. △ Less

Submitted 7 July, 2017; v1 submitted 12 June, 2017; originally announced June 2017.

Comments: 17 pages, 9 figures, accepted for publication in Astronomical Notes

Self-consistent Modelling of the Milky Way using Gaia data

Authors: David R Cole, James Binney

Abstract: Angle/action based distribution function (DF) models can be optimised based on how well they reproduce observations thus revealing the current matter distribution in the Milky Way. Gaia data combined with data from other surveys, e.g. the RAVE/TGAS sample, and its full selection function will greatly improve their accuracy. Angle/action based distribution function (DF) models can be optimised based on how well they reproduce observations thus revealing the current matter distribution in the Milky Way. Gaia data combined with data from other surveys, e.g. the RAVE/TGAS sample, and its full selection function will greatly improve their accuracy. △ Less

Submitted 6 June, 2017; originally announced June 2017.

Comments: 4 pages, 1 figure, to appear in proceedings of IAU Symposium No. 330, 2017 'Astrometry and Astrophysics in the Gaia sky' eds. Recio-Blanco, de Laverny & Brown

Self-consistent modelling of our Galaxy with Gaia data

Authors: James Binney

Abstract: Galaxy models are fundamental to exploiting surveys of our Galaxy. There is now a significant body of work on axisymmetric models. A model can be defined by giving the DF of each major class of stars and of dark matter. Then the self-consistent gravitational potential is determined. Other modelling techniques are briefly considered before an overview of some early work on non-axisymmetric models. Galaxy models are fundamental to exploiting surveys of our Galaxy. There is now a significant body of work on axisymmetric models. A model can be defined by giving the DF of each major class of stars and of dark matter. Then the self-consistent gravitational potential is determined. Other modelling techniques are briefly considered before an overview of some early work on non-axisymmetric models. △ Less

Submitted 5 June, 2017; originally announced June 2017.

Comments: 8 pages, 1 figure to appear in proceedings of IAU Symposium 330, eds. Recio-Blanco, de Laverny & Brown

The structural evolution of galaxies with both thin and thick discs

Authors: Michael Aumer, James Binney

Abstract: We perform controlled N-body simulations of disc galaxies growing within live dark matter (DM) haloes to present-day galaxies that contain both thin and thick discs. We consider two types of models: a) thick-disc initial conditions to which stars on near-circular orbits are continuously added over ~10 Gyr and b) models in which the birth velocity dispersion of stars decreases continuously over the… ▽ More We perform controlled N-body simulations of disc galaxies growing within live dark matter (DM) haloes to present-day galaxies that contain both thin and thick discs. We consider two types of models: a) thick-disc initial conditions to which stars on near-circular orbits are continuously added over ~10 Gyr and b) models in which the birth velocity dispersion of stars decreases continuously over the same time-scale. We show that both schemes produce double-exponential vertical profiles similar to that of the Milky Way (MW). We indicate how the spatial age structure of galaxies can be used to discriminate between scenarios. We show that the presence of a thick disc significantly alters and delays bar formation and thus makes possible models with a realistic bar and a high baryon-to-DM mass ratio in the central regions, as required by microlensing constraints. We examine how the radial mass distribution in stars and DM is affected by disc growth and non-axisymmetries. We discuss how bar buckling shapes the vertical age distribution of thin- and thick-disc stars in the bar region. The extent to which the combination of observationally motivated inside-out growth histories and cosmologically motivated dark halo properties leads to the spontaneous formation of non-axisymmetries that steer the models towards present-day MW-like galaxies is noteworthy. △ Less

Submitted 13 June, 2017; v1 submitted 25 May, 2017; originally announced May 2017.

Comments: Accepted for publication in MNRAS, 22 pages, 13 figures, 2 tables

Modelling the Milky Way's globular cluster system

Authors: James Binney, Leong Khim Wong

Abstract: We construct a model for the Galactic globular cluster system based on a realistic gravitational potential and a distribution function (DF) analytic in the action integrals. The DF comprises disc and halo components whose functional forms resemble those recently used to describe the stellar discs and stellar halo. We determine the posterior distribution of our model parameters using a Bayesian app… ▽ More We construct a model for the Galactic globular cluster system based on a realistic gravitational potential and a distribution function (DF) analytic in the action integrals. The DF comprises disc and halo components whose functional forms resemble those recently used to describe the stellar discs and stellar halo. We determine the posterior distribution of our model parameters using a Bayesian approach. This gives us an understanding of how well the globular cluster data constrain our model. The favoured parameter values of the disc and halo DFs are similar to values previously obtained from fits to the stellar disc and halo, although the cluster halo system shows clearer rotation than does the stellar halo. Our model reproduces the generic features of the globular cluster system, namely the density profile, the mean rotation velocity. The fraction of disc clusters coincides with the observed fraction of metal-rich clusters. However, the data indicate either incompatibility between catalogued cluster distances and current estimates of distance to the Galactic Centre, or failure to identify clusters behind the bulge. As the data for our Galaxy's components increase in volume and precision over the next few years, it will be rewarding to revisit the present analysis. △ Less

Submitted 24 January, 2017; originally announced January 2017.

Comments: 13 pages accepted by MNRAS

The angular momentum of cosmological coronae and the inside-out growth of spiral galaxies

Authors: Gabriele Pezzulli, Filippo Fraternali, James Binney

Abstract: Massive and diffuse haloes of hot gas (coronae) are important intermediaries between cosmology and galaxy evolution, storing mass and angular momentum acquired from the cosmic web until eventual accretion on to star-forming discs. We introduce a method to reconstruct the rotation of a galactic corona, based on its angular momentum distribution (AMD). This allows us to investigate in what condition… ▽ More Massive and diffuse haloes of hot gas (coronae) are important intermediaries between cosmology and galaxy evolution, storing mass and angular momentum acquired from the cosmic web until eventual accretion on to star-forming discs. We introduce a method to reconstruct the rotation of a galactic corona, based on its angular momentum distribution (AMD). This allows us to investigate in what conditions the angular momentum acquired from tidal torques can be transferred to star forming discs and explain observed galaxy-scale processes, such as inside-out growth and the build-up of abundance gradients. We find that a simple model of an isothermal corona with a temperature slightly smaller than virial and a cosmologically motivated AMD is in good agreement with galaxy evolution requirements, supporting hot-mode accretion as a viable driver for the evolution of spiral galaxies in a cosmological context. We predict moderately sub-centrifugal rotation close to the disc and slow rotation close to the virial radius. Motivated by the observation that the Milky Way has a relatively hot corona (T ~ 2 x 10^6 K), we also explore models with a temperature larger than virial. To be able to drive inside-out growth, these models must be significantly affected by feedback, either mechanical (ejection of low angular momentum material) or thermal (heating of the central regions). However, the agreement with galaxy evolution constraints becomes, in these cases, only marginal, suggesting that our first and simpler model may apply to a larger fraction of galaxy evolution history. △ Less

Submitted 19 January, 2017; v1 submitted 5 January, 2017; originally announced January 2017.

Comments: MNRAS, in press. 21 pages, 10 figures, 2 tables

The selection function of the RAVE survey

Authors: Jennifer Wojno, Georges Kordopatis, Tilmann Piffl, James Binney, Matthias Steinmetz, Gal Matijevič, Joss Bland-Hawthorn, Sanjib Sharma, Paul McMillan, Fred Watson, Warren Reid, Andrea Kunder, Harry Enke, Eva K. Grebel, George Seabroke, Rosemary F. G. Wyse, Tomaž Zwitter, Olivier Bienaymé, Kenneth C. Freeman, Brad K. Gibson, Gerry Gilmore, Amina Helmi, Ulisse Munari, Julio F. Navarro, Quentin A. Parker , et al. (1 additional authors not shown)

Abstract: We characterize the selection function of RAVE using 2MASS as our underlying population, which we assume represents all stars which could have potentially been observed. We evaluate the completeness fraction as a function of position, magnitude, and color in two ways: first, on a field-by-field basis, and second, in equal-size areas on the sky. Then, we consider the effect of the RAVE stellar para… ▽ More We characterize the selection function of RAVE using 2MASS as our underlying population, which we assume represents all stars which could have potentially been observed. We evaluate the completeness fraction as a function of position, magnitude, and color in two ways: first, on a field-by-field basis, and second, in equal-size areas on the sky. Then, we consider the effect of the RAVE stellar parameter pipeline on the final resulting catalogue, which in principle limits the parameter space over which our selection function is valid. Our final selection function is the product of the completeness fraction and the selection function of the pipeline. We then test if the application of the selection function introduces biases in the derived parameters. To do this, we compare a parent mock catalogue generated using Galaxia with a mock-RAVE catalogue where the selection function of RAVE has been applied. We conclude that for stars brighter than I = 12, between $4000 \rm K < T_{\rm eff} < 8000 \rm K$ and $0.5 < \rm{log}\,g < 5.0$, RAVE is kinematically and chemically unbiased with respect to expectations from Galaxia. △ Less

Submitted 9 March, 2017; v1 submitted 2 November, 2016; originally announced November 2016.

Comments: 14 pages, 12 figures, accepted for publication in MNRAS

A centrally heated dark halo for our Galaxy

Authors: David R. Cole, James Binney

Abstract: We construct a new family of models of our Galaxy in which dark matter and disc stars are both represented by distribution functions that are analytic functions of the action integrals of motion. The potential that is self-consistently generated by the dark matter, stars and gas is determined, and parameters in the distribution functions are adjusted until the model is compatible with observationa… ▽ More We construct a new family of models of our Galaxy in which dark matter and disc stars are both represented by distribution functions that are analytic functions of the action integrals of motion. The potential that is self-consistently generated by the dark matter, stars and gas is determined, and parameters in the distribution functions are adjusted until the model is compatible with observational constraints on the circular-speed curve, the vertical density profile of the stellar disc near the Sun, the kinematics of nearly 200 000 giant stars within 2 kpc of the Sun, and estimates of the optical depth to microlensing of bulge stars. We find that the data require a dark halo in which the phase-space density is approximately constant for actions |J| \lesssim 140 kpc km ^-1. In real space these haloes have core radii ~ 2 kpc. △ Less

Submitted 25 October, 2016; originally announced October 2016.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

The Radial Velocity Experiment (RAVE): Fifth Data Release

Authors: Andrea Kunder, Georges Kordopatis, Matthias Steinmetz, Tomaz Zwitter, Paul McMillan, Luca Casagrande, Harry Enke, Jennifer Wojno, Marica Valentini, Cristina Chiappini, Gal Matijevic, Alessandro Siviero, Patrick de Laverny, Alejandra Recio-Blanco, Albert Bijaoui, Rosemary F. G. Wyse, James Binney, E. K. Grebel, Amina Helmi, Paula Jofre, Teresa Antoja, Gerard Gilmore, Arnaud Siebert, Benoit Famaey, Olivier Bienayme , et al. (29 additional authors not shown)

Abstract: Data Release 5 (DR5) of the Radial Velocity Experiment (RAVE) is the fifth data release from a magnitude-limited (9< I < 12) survey of stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra ($R\sim7500$) covering the Ca-triplet region (8410-8795Å) span the complete time frame from the start of RAVE observations in 2003 to their completion in 2013. Radial velocities… ▽ More Data Release 5 (DR5) of the Radial Velocity Experiment (RAVE) is the fifth data release from a magnitude-limited (9< I < 12) survey of stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra ($R\sim7500$) covering the Ca-triplet region (8410-8795Å) span the complete time frame from the start of RAVE observations in 2003 to their completion in 2013. Radial velocities from 520,781 spectra of 457,588 unique stars are presented, of which 255,922 stellar observations have parallaxes and proper motions from the Tycho-Gaia astrometric solution (TGAS) in Gaia DR1. For our main DR5 catalog, stellar parameters (effective temperature, surface gravity, and overall metallicity) are computed using the RAVE DR4 stellar pipeline, but calibrated using recent K2 Campaign 1 seismic gravities and Gaia benchmark stars, as well as results obtained from high-resolution studies. Also included are temperatures from the Infrared Flux Method, and we provide a catalogue of red giant stars in the dereddened color $(J-Ks)_0$ interval (0.50,0.85) for which the gravities were calibrated based only on seismology. Further data products for sub-samples of the RAVE stars include individual abundances for Mg, Al, Si, Ca, Ti, Fe, and Ni, and distances found using isochrones. Each RAVE spectrum is complemented by an error spectrum, which has been used to determine uncertainties on the parameters. The data can be accessed via the RAVE Web site or the Vizier database. △ Less

Submitted 16 November, 2016; v1 submitted 11 September, 2016; originally announced September 2016.

Comments: accepted Astronomical Journal version, data tables on RAVE Web site: https://www.rave-survey.org/project/

Characterizing stellar halo populations II: The age gradient in blue horizontal-branch stars

Authors: Payel Das, Angus Williams, James Binney

Abstract: The distribution of Milky Way halo blue horizontal-branch (BHB) stars is examined using action-based extended distribution functions (EDFs) that describe the locations of stars in phase space, metallicity, and age. The parameters of the EDFs are fitted using stars observed in the Sloan Extension for Galactic Understanding and Exploration-II (SEGUE-II) survey that trace the phase-space kinematics a… ▽ More The distribution of Milky Way halo blue horizontal-branch (BHB) stars is examined using action-based extended distribution functions (EDFs) that describe the locations of stars in phase space, metallicity, and age. The parameters of the EDFs are fitted using stars observed in the Sloan Extension for Galactic Understanding and Exploration-II (SEGUE-II) survey that trace the phase-space kinematics and chemistry out to ~70 kpc. A maximum a posteriori probability (MAP) estimate method and a Markov Chain Monte Carlo method are applied, taking into account the selection function in positions, distance, and metallicity for the survey. The best-fit EDF declines with actions less steeply at actions characteristic of the inner halo than at the larger actions characteristic of the outer halo, and older ages are found at smaller actions than at larger actions. In real space, the radial density profile steepens smoothly from -2 at ~2 kpc to -4 in the outer halo, with an axis ratio ~0.7 throughout. There is no indication for rotation in the BHBs, although this is highly uncertain. A moderate level of radial anisotropy is detected, with $βべーた_s$ varying from isotropic to between ~0.1 and ~0.3 in the outer halo depending on latitude. The BHB data are consistent with an age gradient of -0.03 Gyr kpc$^{-1}$, with some uncertainty in the distribution of the larger ages. These results are consistent with a scenario in which older, larger systems contribute to the inner halo, whilst the outer halo is primarily comprised of younger, smaller systems. △ Less

Submitted 25 August, 2016; originally announced August 2016.

Comments: 18 pages, 15 figures, 3 tables, accepted for publication in MNRAS

Managing resonant trapped orbits in our Galaxy

Authors: James Binney

Abstract: Galaxy modelling is greatly simplified by assuming the existence of a global system of angle-action coordinates. Unfortunately, global angle-action coordinates do not exist because some orbits become trapped by resonances, especially where the radial and vertical frequencies coincide. We show that in a realistic Galactic potential such trapping occurs only on thick-disc and halo orbits (speed rela… ▽ More Galaxy modelling is greatly simplified by assuming the existence of a global system of angle-action coordinates. Unfortunately, global angle-action coordinates do not exist because some orbits become trapped by resonances, especially where the radial and vertical frequencies coincide. We show that in a realistic Galactic potential such trapping occurs only on thick-disc and halo orbits (speed relative to the guiding centre >~80 km/s). We explain how the Torus Mapper code (TM) behaves in regions of phase space in which orbits are resonantly trapped, and we extend TM so trapped orbits can be manipulated as easily as untrapped ones. The impact that the resonance has on the structure of velocity space depends on the weights assigned to trapped orbits. The impact is everywhere small if each trapped orbit is assigned the phase space density equal to the time average along the orbit of the DF for untrapped orbits. The impact could be significant with a different assignment of weights to trapped orbits. △ Less

Submitted 4 August, 2016; originally announced August 2016.

Comments: 13 pages, 17 figures, accepted by MNRAS

Age velocity dispersion relations and heating histories in disc galaxies

Authors: Michael Aumer, James Binney, Ralph Schönrich

Abstract: We analyse the heating of stellar discs by non axisymmetric structures and giant molecular clouds (GMCs) in N-body simulations of growing disc galaxies. The analysis resolves long-standing discrepancies between models and data by demonstrating the importance of distinguishing between measured age-velocity dispersion relations (AVRs) and the heating histories of the stars that make up the AVR. We f… ▽ More We analyse the heating of stellar discs by non axisymmetric structures and giant molecular clouds (GMCs) in N-body simulations of growing disc galaxies. The analysis resolves long-standing discrepancies between models and data by demonstrating the importance of distinguishing between measured age-velocity dispersion relations (AVRs) and the heating histories of the stars that make up the AVR. We fit both AVRs and heating histories with formulae proportional to t^beta and determine the exponents beta_R and beta_z derived from in-plane and vertical AVRs and ~beta_R and ~beta_z from heating histories. Values of beta_z are in almost all simulations larger than values of ~beta_z, whereas values of beta_R are similar to or mildly larger than values of ~beta_R. Moreover, values of beta_z (~beta_z) are generally larger than values of beta_R (~beta_R). The dominant cause of these relations is the decline over the life of the disc in importance of GMCs as heating agents relative to spiral structure and the bar. We examine how age errors and biases in solar neighbourhood surveys influence the measured AVR: they tend to decrease beta values by smearing out ages and thus measured dispersions. We compare AVRs and velocity ellipsoid shapes sigma_z/sigma_R from simulations to Solar neighbourhood data. We conclude that for the expected disc mass and dark halo structure, combined GMC and spiral/bar heating can explain the AVR of the Galactic thin disc. Strong departures of the disc mass or the dark halo structure from expectation spoil fits to the data. △ Less

Submitted 7 July, 2016; originally announced July 2016.

Comments: Accepted for publication in MNRAS, 19 pages, 8 figures, 1 table

The quiescent phase of galactic disc growth

Authors: Michael Aumer, James Binney, Ralph Schönrich

Abstract: We perform a series of controlled N-body simulations of growing disc galaxies within non-growing, live dark matter haloes of varying mass and concentration. Our initial conditions include either a low-mass disc or a compact bulge. New stellar particles are continuously added on near-circular orbits to the existing disc, so spiral structure is continuously excited. To study the effect of combined s… ▽ More We perform a series of controlled N-body simulations of growing disc galaxies within non-growing, live dark matter haloes of varying mass and concentration. Our initial conditions include either a low-mass disc or a compact bulge. New stellar particles are continuously added on near-circular orbits to the existing disc, so spiral structure is continuously excited. To study the effect of combined spiral and giant molecular cloud (GMC) heating on the discs we introduce massive, short-lived particles that sample a GMC mass function. An isothermal gas component is introduced for a subset of the models. We perform a resolution study and vary parameters governing the GMC population, the histories of star formation and radial scale growth. Models with GMCs and standard values for the disc mass and halo density provide the right level of self-gravity to explain the age velocity dispersion relation of the Solar neighbourhood (Snhd). GMC heating generates remarkably exponential vertical profiles with scaleheights that are radially constant and agree with observations of galactic thin discs. GMCs are also capable of significantly delaying bar formation. The amount of spiral induced radial migration agrees with what is required for the metallicity distribution of the Snhd. However, in our standard models the outward migrating populations are not hot enough vertically to create thick discs. Thick discs can form in models with high baryon fractions, but the corresponding bars are too long, the young stellar populations too hot and the discs flare considerably. △ Less

Submitted 1 April, 2016; originally announced April 2016.

Comments: Accepted for publication in MNRAS, 26 pages, 16 figures, 4 tables

Chemical separation of disc components using RAVE

Authors: Jennifer Wojno, Georges Kordopatis, Matthias Steinmetz, Paul J. McMillan, Gal Matijevič, James Binney, Rosemary F. G. Wyse, Corrado Boeche, Andreas Just, Eva K. Grebel, Arnaud Siebert, Olivier Bienaymé, Brad K. Gibson, Tomaž Zwitter, Joss Bland-Hawthorn, Julio F. Navarro, Quentin A. Parker, Warren Reid, George Seabroke, Fred Watson

Abstract: We present evidence from the RAdial Velocity Experiment (RAVE) survey of chemically separated, kinematically distinct disc components in the solar neighbourhood. We apply probabilistic chemical selection criteria to separate our sample into $αあるふぁ$-low (`thin disc') and $αあるふぁ$-high (`thick disc') sequences. Using newly derived distances, which will be utilized in the upcoming RAVE DR5, we explore the kin… ▽ More We present evidence from the RAdial Velocity Experiment (RAVE) survey of chemically separated, kinematically distinct disc components in the solar neighbourhood. We apply probabilistic chemical selection criteria to separate our sample into $αあるふぁ$-low (`thin disc') and $αあるふぁ$-high (`thick disc') sequences. Using newly derived distances, which will be utilized in the upcoming RAVE DR5, we explore the kinematic trends as a function of metallicity for each of the disc components. For our $αあるふぁ$-low disc, we find a negative trend in the mean rotational velocity ($V_{\mathrmφふぁい}$) as a function of iron abundance ([Fe/H]). We measure a positive gradient $\partial V_{\mathrmφふぁい}$/$\partial$[Fe/H] for the $αあるふぁ$-high disc, consistent with results from high-resolution surveys. We also find differences between the $αあるふぁ$-low and $αあるふぁ$-high discs in all three components of velocity dispersion. We discuss the implications of an $αあるふぁ$-low, metal-rich population originating from the inner Galaxy, where the orbits of these stars have been significantly altered by radial mixing mechanisms in order to bring them into the solar neighbourhood. The probabilistic separation we propose can be extended to other data sets for which the accuracy in [$αあるふぁ$/Fe] is not sufficient to disentangle the chemical disc components a priori. For such datasets which will also have significant overlap with Gaia DR1, we can therefore make full use of the improved parallax and proper motion data as it becomes available to investigate kinematic trends in these chemical disc components. △ Less

Submitted 7 July, 2016; v1 submitted 30 March, 2016; originally announced March 2016.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS

Characterising stellar halo populations I: An extended distribution function for halo K giants

Authors: Payel Das, James Binney

Abstract: We fit an Extended Distribution Function (EDF) to K giants in the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey. These stars are detected to radii ~80 kpc and span a wide range in [Fe/H]. Our EDF, which depends on [Fe/H] in addition to actions, encodes the entanglement of metallicity with dynamics within the Galaxy's stellar halo. Our maximum-likelihood fit of the EDF t… ▽ More We fit an Extended Distribution Function (EDF) to K giants in the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey. These stars are detected to radii ~80 kpc and span a wide range in [Fe/H]. Our EDF, which depends on [Fe/H] in addition to actions, encodes the entanglement of metallicity with dynamics within the Galaxy's stellar halo. Our maximum-likelihood fit of the EDF to the data allows us to model the survey's selection function. The density profile of the K giants steepens with radius from a slope ~-2 to ~-4 at large radii. The halo's axis ratio increases with radius from 0.7 to almost unity. The metal-rich stars are more tightly confined in action space than the metal-poor stars and form a more flattened structure. A weak metallicity gradient ~-0.001 dex/kpc, a small gradient in the dispersion in [Fe/H] of ~0.001 dex/kpc, and a higher degree of radial anistropy in metal-richer stars result. Lognormal components with peaks at ~-1.5 and ~-2.3 are required to capture the overall metallicity distribution, suggestive of the existence of two populations of K giants. The spherical anisotropy parameter varies between 0.3 in the inner halo to isotropic in the outer halo. If the Sagittarius stream is included, a very similar model is found but with a stronger degree of radial anisotropy throughout. △ Less

Submitted 30 March, 2016; originally announced March 2016.

Comments: 15 pages, 15 figures, accepted for publication in MNRAS

Chemodynamics of the Milky Way and disc formation history: insight from the RAVE and Gaia-ESO surveys

Authors: G. Kordopatis, R. F. G. Wyse, J. Binney

Abstract: Multi-object spectrographs have opened a new window on the analyses of the chemo-dynamical properties of old Milky Way stars. These analyses allow us to trace back the internal mechanisms and the external factors that have influenced the evolution of our Galaxy, and therefore understand fundamental aspects of galaxy evolution in general. Here, we present recent results from the RAdial Velocity Exp… ▽ More Multi-object spectrographs have opened a new window on the analyses of the chemo-dynamical properties of old Milky Way stars. These analyses allow us to trace back the internal mechanisms and the external factors that have influenced the evolution of our Galaxy, and therefore understand fundamental aspects of galaxy evolution in general. Here, we present recent results from the RAdial Velocity Experiment (RAVE) and the Gaia-ESO survey. These surveys explore the Milky Way properties in different ways, in terms of sample size and selection, magnitude range, and spectral resolution. We focus here on (i) the first direct detection of evidence for radial migration within the thin disc, providing insight into the history of spiral structure of the Milky Way, and (ii) the chemo-dynamical characterisation of the metal-weak thick and thin discs, for which chemo-dynamical models still have difficulties in reproducing. △ Less

Submitted 11 January, 2016; originally announced January 2016.

Comments: 5 pages, to appear in Astronomische Nachrichten, special issue "Reconstruction the Milky Way's History: Spectroscopic surveys, Asteroseismology and Chemo-dynamical models", Guest Editors C. Chiappini, J. Montalban, and M. Steffen, AN 2016 (in press)"

Chemodynamical modelling of the Milky Way

Authors: James Binney, Jason L Sanders

Abstract: Chemodynamical models of our Galaxy that have analytic Extended Distribution Functions (EDFs) are likely to play a key role in extracting science from surveys in the era of Gaia. Chemodynamical models of our Galaxy that have analytic Extended Distribution Functions (EDFs) are likely to play a key role in extracting science from surveys in the era of Gaia. △ Less

Submitted 26 November, 2015; originally announced November 2015.

Comments: 6 pages to appear in "Reconstructing the Milky Way history: spectroscopic surveys, asteroseismology and chemodynamical models", eds C Chiappini, J Montalban & M Steffen

A review of action estimation methods for galactic dynamics

Authors: Jason L. Sanders, James Binney

Abstract: We review the available methods for estimating actions, angles and frequencies of orbits in both axisymmetric and triaxial potentials. The methods are separated into two classes. Unless an orbit has been trapped by a resonance, convergent, or iterative, methods are able to recover the actions to arbitrarily high accuracy given sufficient computing time. Faster non-convergent methods rely on the po… ▽ More We review the available methods for estimating actions, angles and frequencies of orbits in both axisymmetric and triaxial potentials. The methods are separated into two classes. Unless an orbit has been trapped by a resonance, convergent, or iterative, methods are able to recover the actions to arbitrarily high accuracy given sufficient computing time. Faster non-convergent methods rely on the potential being sufficiently close to a separable potential and the accuracy of the action estimate cannot be improved through further computation. We critically compare the accuracy of the methods and the required computation time for a range of orbits in an axisymmetric multi-component Galactic potential. We introduce a new method for estimating actions that builds on the adiabatic approximation of Schönich & Binney (2012) and discuss the accuracy required for the actions, angles and frequencies using suitable distribution functions for the thin and thick discs, the stellar halo and a star stream. We conclude that for studies of the disc and smooth halo component of the Milky Way the most suitable compromise between speed and accuracy is the Stäckel Fudge, whilst when studying streams the non-convergent methods do not offer sufficient accuracy and the most suitable method is computing the actions from an orbit integration via a generating function. All the software used in this study can be downloaded from https://github.com/jls713/tact. △ Less

Submitted 19 January, 2016; v1 submitted 25 November, 2015; originally announced November 2015.

Comments: 15 pages, 7 figures, accepted for publication in MNRAS (new subsection on resonant trapping), associated software downloadable at https://github.com/jls713/tact

Torus mapper: a code for dynamical models of galaxies

Authors: James Binney, Paul J McMillan

Abstract: We present a freely downloadable software package for modelling the dynamics of galaxies, which we call the Torus Mapper (TM). The package is based around `torus mapping', which is a non-perturbative technique for creating orbital tori for specified values of the action integrals. Given an orbital torus and a star's position at a reference time, one can compute its position at any other time, no m… ▽ More We present a freely downloadable software package for modelling the dynamics of galaxies, which we call the Torus Mapper (TM). The package is based around `torus mapping', which is a non-perturbative technique for creating orbital tori for specified values of the action integrals. Given an orbital torus and a star's position at a reference time, one can compute its position at any other time, no matter how remote. One can also compute the velocities with which the star will pass through any given point and the contribution it will make to the time-averaged density there. A system of angle-action coordinates for the given potential can be created by foliating phase space with orbital tori. Such a foliation is facilitated by the ability of TM to create tori by interpolating on a grid of tori. We summarise the advantages of using TM rather than a standard time-stepper to create orbits, and give segments of code that illustrate applications of TM in several contexts, including setting up initial conditions for an N-body simulation. We examine the precision of the orbital tori created by TM and the behaviour of the code when orbits become trapped by a resonance. △ Less

Submitted 24 November, 2015; originally announced November 2015.

Comments: 19 pages accepted for publication in MNRAS

The distribution function of the Galaxy's dark halo

Authors: James Binney, Tilmann Piffl

Abstract: Starting from the hypothesis that the Galaxy's dark halo responded adiabatically to the infall of baryons, we have constructed a self-consistent dynamical model of the Galaxy that satisfies a large number of observations, including measurements of gas terminal velocities and masers, the kinematics of a 180,000 giant stars from the RAVE survey, and star count data from the SDSS. The stellar disc an… ▽ More Starting from the hypothesis that the Galaxy's dark halo responded adiabatically to the infall of baryons, we have constructed a self-consistent dynamical model of the Galaxy that satisfies a large number of observations, including measurements of gas terminal velocities and masers, the kinematics of a 180,000 giant stars from the RAVE survey, and star count data from the SDSS. The stellar disc and the dark halo are both specified by distribution functions (DFs) of the action integrals. The model is obtained by extending the work of Piffl Penoyre & Binney (2015} from the construction of a single model to a systematic search of model space. Whereas the model of Piffl et al violated constraints on the terminal-velocity curve, our model respects these constraints by adopting a long scale length R_d=3.66 kpc for the thin and thick discs. The model is, however, inconsistent with the measured optical depth for microlensing of bulge stars because it attributes too large a fraction of the density at R <~ 3 kpc to dark matter rather than stars. Moreover, it now seems likely that the thick disc's scale-length is significantly shorter than the model implies. Shortening this scale-length would cause the constraints from the rotation curve to be violated anew. We conclude that we can now rule out adiabatic compression of our Galaxy's dark halo. △ Less

Submitted 23 September, 2015; originally announced September 2015.

Comments: 11 pages MNRAS accepted