Jason L. Sanders

Stream-orbit misalignment II: A new algorithm to constrain the Galactic potential

In the first of these two papers we demonstrated that assuming streams delineate orbits can lead to order one errors in potential parameters for realistic Galactic potentials. Motivated by the need for an improvement on orbit-fitting, we now present an algorithm for constraining the Galactic potential using tidal streams without assuming that streams delineate orbits. This approach is independent of the progenitor mass so is valid for all observed tidal streams. The method makes heavy use of angle-action variables and seeks the potential which recovers the expected correlations in angle space. We demonstrate that the method can correctly recover the parameters of a simple two-parameter logarithmic potential by analysing an N-body simulation of a stream. We investigate the magnitude of the errors in observational data for which the method can still recover the correct potential and compare this to current and future errors in data. The errors in the observables of individual stars for current and near future data are shown to be too large for the direct use of this method, but when the data are averaged in bins on the sky, the resulting averaged data are accurate enough to constrain correctly the potential parameters for achievable observational errors. From pseudo-data with errors comparable to those that will be furnished in the era of Gaia (20 per cent distance errors, 1.2 mas/yr proper motion errors, and 10 km/s line-of-sight velocity errors) we recover the circular velocity, V_c=220 km/s, and the flattening of the potential, q=0.9, to be V_c=223+/-10km/s and q=0.91+/-0.09.

The number and size of subhalo-induced gaps in stellar streams

European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013) (Grant ID: 308024), Natural Sciences and Engineering Research Council of Canada, Science and Technology Facilities Council

Extended distribution functions for our Galaxy

We extend models of our Galaxy based on distribution functions (DFs) that are analytic functions of the action integrals to extended distribution functions (EDFs), which have an analytic dependence on metallicity as well. We use a simple, but physically-motivated, functional forms for the metallicity of the interstellar medium as a function of radius and time and for the star-formation rate, and a model for the diffusion of stars through phase space to suggest the required functional form of an EDF. We introduce a simple prescription for radial migration that preserves the overall profile of the disc while allowing individual stars to migrate throughout the disc. Our models explicitly consider the thin and thick discs as two distinct components separated in age. We show how an EDF can be used to incorporate realistic selection functions in models, and to construct mock catalogues of observed samples. We show that the selection function of the Geneva-Copenhagen Survey (GCS) biases in favour of young stars, which have atypically small random velocities. With the selection function taken into account our models produce good fits of the GCS data in chemo-dynamical space and the Gilmore and Reid (1983) density data. From our EDF, we predict the structure of the SEGUE G-dwarf sample. The kinematics are successfully predicted. The predicted metallicity distribution has too few stars with [Fe/H]

Linear perturbation theory for tidal streams and the small-scale CDM power spectrum

\simeq-0.5

Probabilistic model for constraining the Galactic potential using tidal streams

dex and too many metal-rich stars. A significant problem may be the lack of any chemical-kinematic correlations in our thick disc. We argue that EDFs will prove essential tools for the analysis of both observational data and sophisticated models of Galaxy formation and evolution.

The RAVE-on Catalog of Stellar Atmospheric Parameters and Chemical Abundances for Chemo-dynamic Studies in the Gaia Era

Tidal streams in the Milky Way are sensitive probes of the population of dark-matter subhalos predicted in cold-dark-matter (CDM) simulations. We present a new calculus for computing the effect of subhalo fly-bys on cold tidal streams based on the action-angle representation of streams. The heart of this calculus is a line-of-parallel-angle approach that calculates the perturbed distribution function of a given stream segment by undoing the effect of all impacts. This approach allows one to compute the perturbed stream density and track in any coordinate system in minutes for realizations of the subhalo distribution down to 10^5 Msun, accounting for the streams internal dispersion and overlapping impacts. We study the properties of density and track fluctuations with suites of simulations. The one-dimensional density and track power spectra along the stream trace the subhalo mass function, with higher-mass subhalos producing power only on large scales, while lower mass subhalos cause structure on smaller scales. The time-dependence of impacts and of the evolution of the stream after an impact gives rise to bispectra. We further find that tidal streams are essentially corrugated sheets in the presence of subhalo perturbations: different projections of the track all reflect the same pattern of perturbations, facilitating their observational measurement. We apply this formalism to density data for the Pal 5 stream and make a first rigorous determination of 10^{+11}_{-6} dark-matter subhalos with masses between 3x10^6 and 10^9 Msun within 20 kpc from the Galactic center (corresponding to 1.4^{+1.6}_{-0.9} times the number predicted by CDM-only simulations or to f_{sub}(r<20 kpc) ~ 0.2%). Improved data will allow measurements of the subhalo mass function down to 10^5 Msun, thus definitively testing whether dark matter clumps on the smallest scales relevant for galaxy formation.

Dynamics of stream–subhalo interactions

We present a generative probabilistic model for a tidal stream and demonstrate how this model is used to constrain the Galactic potential. The model takes advantage of the simple structure of a stream in angle and frequency space for the correct potential. We investigate how the method performs on full 6D mock stream data, and mock data with outliers included. As currently formulated the technique is computationally costly when applied to data with large observational errors, but we describe several modifications that promise to make the technique computationally tractable.

Actions, angles and frequencies for numerically integrated orbits

The orbits, atmospheric parameters, chemical abundances, and ages of individual stars in the Milky Way provide the most comprehensive illustration of galaxy formation available. The Tycho-Gaia Astrometric Solution (TGAS) will deliver astrometric parameters for the largest ever sample of Milky Way stars, though its full potential cannot be realized without the addition of complementary spectroscopy. Among existing spectroscopic surveys, the RAdial Velocity Experiment (RAVE) has the largest overlap with TGAS (

Angle-Action Estimation in a General Axisymmetric Potential

\gtrsim

A review of action estimation methods for galactic dynamics

200,000 stars). We present a data-driven re-analysis of 520,781 RAVE spectra using The Cannon. For red giants, we build our model using high-fidelity APOGEE stellar parameters and abundances for stars that overlap with RAVE. For main-sequence and sub-giant stars, our model uses stellar parameters from the K2/EPIC. We derive and validate effective temperature