Professor James Binney FRS

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 fit both AVRs and heating histories with formulae ∝ tβ and determine the exponents β bar R and β bar z derived from in-plane and vertical AVRs and βR and βz from heating histories. Values of βz are in almost all simulations larger than values of βbarz , wheras values of βR are similar to or mildly larger than values of β bar R. Moreover, values of βz (β bar z) are generally larger than values of βR (β bar 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 in uence the measured AVR: they tend to decrease β values by smearing out ages and thus measured dispersions. We compare AVRs and velocity ellipsoid shapes σz/σ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.

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 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φ) as a function of iron abundance ([Fe/H]). We measure a positive gradient ∂Vφ/∂[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 data sets 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.

Abstract:

We fit an extended distribution function (EDF) to K giants in the Sloan Extension for Galactic Understanding and Exploration 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−1, a small gradient in the dispersion in [Fe/H] of ∼0.001 dex kpc−1, and a higher degree of radial anisotropy 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.

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 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 multicomponent Galactic potential. We introduce a new method for estimating actions that builds on the adiabatic approximation of Schönrich and Binney 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.