• ### A two-parameter family of double-power-law biorthonormal potential-density expansions(1804.11190)

April 26, 2018 astro-ph.CO, astro-ph.GA
Biorthonormal basis function expansions are widely used in galactic dynamics, both to study problems in galactic stability and to provide numerical algorithms to evolve collisionless stellar systems. They also provide a compact and efficient description of the structure of numerical dark matter haloes in cosmological simulations. We present a two-parameter family of biorthonormal double-power-law potential-density expansions. Both the potential and density are given in closed analytic form and may be rapidly computed via recurrence relations. We show that this family encompasses all the known analytic biorthonormal expansions: the Zhao expansions (themselves generalizations of ones found earlier by Hernquist & Ostriker and by Clutton-Brock) and the recently discovered Lilley, Sanders, Evans & Erkal expansion. Our new two-parameter family includes expansions based around many familiar spherical density profiles as zeroth-order models, including the $\gamma$ models and the Jaffe model. It also contains a basis expansion that reproduces the famous Navarro-Frenk-White (NFW) profile at zeroth order. The new basis expansions have been found via a systematic methodology which has wide applications in finding further examples. In the process, we also uncovered a novel integral transform solution to Poisson's equation.
• ### Age Gradients Throughout the Galaxy with O-Miras(1804.09186)

April 24, 2018 astro-ph.GA
We assemble the largest sample of oxygen rich Miras to date and highlight their importance for age-dating the components of the Galaxy. Using data from the Catalina Rapid Transient Survey and the All Sky Automated Survey for Supernovae, we extract a clean sample of $\sim 2,400$ O-Miras, stretching from the Galactic Bulge to the distant halo. Given that the period of O-Miras correlates with age, this offers a new way of determining age gradients throughout the Galaxy. We use our sample to show (i) disk O-Miras have periods increasing on moving outwards from ~ 3 to 15 kpc, so the outer disk O-Miras are younger than the inner disk, (ii) the transition from younger disk to halo O-Miras occurs at r ~ 15 kpc and is marked by a plummeting in period, (iii) there exists a population of young O-Miras likely kicked from the disk to heights of order of |Z| ~ 10 kpc, (iv) great circle counts of old Miras show strong evidence for distant debris agglomeration associated with the Magellanic Clouds, (v) seven stars in our samples are located at distances between 200 and 500 kpc surpassing all previously established records, and, finally, (vi) O-Miras may be present in the Fornax, Sculptor, Sextans and Leo II Galactic dwarf spheroidals, as well as the distant globular cluster Pal 4. We spotlight the importance of O-Mira in the Era of Gaia as universal chronometers of the Galactic populations.
• ### Co-formation of the Galactic disc and the stellar halo(1802.03414)

Feb. 9, 2018 astro-ph.GA
Using a large sample of Main Sequence stars with 7-D measurements supplied by Gaia and SDSS, we study the kinematic properties of the local (within ~10 kpc from the Sun) stellar halo. We demonstrate that the halo's velocity ellipsoid evolves strongly with metallicity. At the low [Fe/H] end, the orbital anisotropy (the amount of motion in the radial direction compared to the tangential one) is mildly radial with 0.2<beta<0.4. However, for stars with [Fe/H]>-1.7 we measure extreme values of beta~0.9. Across the metallicity range considered, i.e. -3<[Fe/H]<-1, the stellar halo's spin is minimal, at the level of 20<v_theta (km/s) <30. Using a suite of cosmological zoom-in simulations of halo formation, we deduce that the observed acute anisotropy is inconsistent with the continuous accretion of dwarf satellites. Instead, we argue, the stellar debris in the inner halo were deposited in a major accretion event by a satellite with Mvir>10^10 Msun around the epoch of the Galactic disc formation, i.e. between 8 and 11 Gyr ago. The radical halo anisotropy is the result of the dramatic radialisation of the massive progenitor's orbit, amplified by the action of the growing disc.
• ### The Super-NFW model: an analytic dynamical model for cold dark matter haloes and elliptical galaxies(1802.03349)

Feb. 9, 2018 astro-ph.CO, astro-ph.GA
An analytic galaxy model with $\rho \sim r^{-1}$ at small radii and $\rho \sim r^{-3.5}$ at large radii is presented. The asymptotic density fall-off is slower than the Hernquist model, but faster than the Navarro-Frenk-White (NFW) profile for dark matter haloes, and so in accord with recent evidence from cosmological simulations. The model provides the zeroth-order term in a biorthornomal basis function expansion, meaning that axisymmetric, triaxial and lopsided distortions can easily be added (much like the Hernquist model itself which is the zeroth-order term of the Hernquist-Ostriker expansion). The properties of the spherical model, including analytic distribution functions which are either isotropic, radially anisotropic or tangentially anisotropic, are discussed in some detail. The analogue of the mass-concentration relation for cosmological haloes is provided.
• ### Unmixing the Galactic Halo with RR Lyrae tagging(1711.07485)

Nov. 20, 2017 astro-ph.GA
We show that tagging RR Lyrae stars according to their location in the period-amplitude diagram can be used to shed light on the genesis of the Galactic stellar halo. The mixture of RR Lyrae of ab type, separated into classes along the lines suggested by Oosterhoff, displays a strong and coherent evolution with Galactocentric radius. The change in the RR Lyrae composition appears to coincide with the break in the halo's radial density profile at ~25 kpc. Using simple models of the stellar halo, we establish that at least three different types of accretion events are necessary to explain the observed RRab behavior. Given that there exists a correlation between the RRab class fraction and the total stellar content of a dwarf satellite, we hypothesize that the field halo RRab composition is controlled by the mass of the progenitor contributing the bulk of the stellar debris at the given radius. This idea is tested against a suite of cosmological zoom-in simulations of Milky Way-like stellar halo formation. Finally, we study some of the most prominent stellar streams in the Milky Way halo and demonstrate that their RRab class fractions follow the trends established previously.
• Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.
• Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.
• ### A White Paper on keV Sterile Neutrino Dark Matter(1602.04816)

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.
• ### Gaia data release 1, the photometric data(1612.02952)

Dec. 9, 2016 astro-ph.IM
Context. This paper presents an overview of the photometric data that are part of the first Gaia data release. Aims. The principles of the processing and the main characteristics of the Gaia photometric data are presented. Methods. The calibration strategy is outlined briefly and the main properties of the resulting photometry are presented. Results. Relations with other broadband photometric systems are provided. The overall precision for the Gaia photometry is shown to be at the milli-magnitude level and has a clear potential to improve further in future releases.
• ### A Tail of Two Populations: Chemo-dynamics of the Sagittarius Stream and Implications for its Original Mass(1607.00803)

July 4, 2016 astro-ph.CO, astro-ph.GA
We use the SDSS/SEGUE spectroscopic sample of stars in the leading and trailing streams of the Sagittarius (Sgr) to demonstrate the existence of two sub-populations with distinct chemistry and kinematics. The metallicity distribution function (MDF) of the trailing stream is decomposed into two Gaussians describing a metal-rich sub-population with means and dispersions (-0.74, 0.18) dex and a metal-poor with (-1.33, 0.27) dex. The metal-rich sub-population has a velocity dispersion ~ 8 km/s, whilst the metal-poor is nearly twice as hot ~13 km/s. For the leading stream, the MDF is again well-described by a superposition of two Gaussians, though somewhat shifted as compared to the trailing stream. The metal-rich has mean and dispersion (-1.00, 0.34) dex, the metal-poor (-1.39, 0.22) dex. The velocity dispersions are inflated by projection effects to give 15 - 30 km/s for the metal-poor, and 6 - 20 km/s for the metal-rich, depending on longitude. We infer that, like many dwarf spheroidals, the Sgr progenitor possessed a more extended, metal-poor stellar component and less extended, metal-rich one. We study the implications of this result for the progenitor mass by simulating the disruption of the Sgr, represented as King light profiles in dark halos of masses between $10^{10}$ and $10^{11} M_\odot$, in a three-component Milky Way whose halo is a live Truncated Flat potential in the first phase of accretion and a triaxial Law \& Majewski model in the second phase. We show that that the dark halo of the Sgr must have been $\gtrsim 6 \times 10^{10} M_\odot$ to reproduce the run of velocity dispersion with longitude for the metal-rich and metal-poor sub-populations in the tails.
• ### The Gaia-ESO survey: Metal-rich bananas in the bulge(1605.09684)

May 31, 2016 astro-ph.GA
We analyse the kinematics of $\sim 2000$ giant stars in the direction of the Galactic bulge, extracted from the Gaia-ESO survey in the region $-10^\circ \lesssim \ell \lesssim 10^\circ$ and $-11^\circ \lesssim b \lesssim -3^\circ$. We find distinct kinematic trends in the metal rich ($\mathrm{[M/H]}>0$) and metal poor ($\mathrm{[M/H]}<0$) stars in the data. The velocity dispersion of the metal-rich stars drops steeply with latitude, compared to a flat profile in the metal-poor stars, as has been seen previously. We argue that the metal-rich stars in this region are mostly on orbits that support the boxy-peanut shape of the bulge, which naturally explains the drop in their velocity dispersion profile with latitude. The metal rich stars also exhibit peaky features in their line-of-sight velocity histograms, particularly along the minor axis of the bulge. We propose that these features are due to stars on resonant orbits supporting the boxy-peanut bulge. This conjecture is strengthened through the comparison of the minor axis data with the velocity histograms of resonant orbits generated in simulations of buckled bars. The 'banana' or 2:1:2 orbits provide strongly bimodal histograms with narrow velocity peaks that resemble the Gaia-ESO metal-rich data.
• ### Simple J-Factors and D-Factors for Indirect Dark Matter Detection(1604.05599)

J-factors (or D-factors) describe the distribution of dark matter in an astrophysical system and determine the strength of the signal provided by annihilating (or decaying) dark matter respectively. We provide simple analytic formulae to calculate the J-factors for spherical cusps obeying the empirical relationship between enclosed mass, velocity dispersion and half-light radius. We extend the calculation to the spherical Navarro-Frenk-White (NFW) model, and demonstrate that our new formulae give accurate results in comparison to more elaborate Jeans models driven by Markov Chain Monte Carlo methods. Of the known ultrafaint dwarf spheroidals, we show that Ursa Major II, Reticulum II, Tucana II and Horologium I have the largest J-factors and so provide the most promising candidates for indirect dark matter detection experiments. Amongst the classical dwarfs, Draco, Sculptor and Ursa Minor have the highest J-factors. We show that the behaviour of the J-factor as a function of integration angle can be inferred for general dark halo models with inner slope $\gamma$ and outer slope $\beta$. The central and asymptotic behaviour of the J-factor curves are derived as a function of the dark halo properties. Finally, we show that models obeying the empirical relation on enclosed mass and velocity dispersion have J-factors that are most robust at the integration angle equal to the projected half-light radius of the dSph divided by heliocentric distance. For most of our results, we give the extension to the D-factor which is appropriate for the decaying dark matter picture.
• ### Discy dwarf disruption and the shape of the Galactic halo(1601.07720)

Feb. 17, 2016 astro-ph.GA
The shape of the Galactic dark halo can, in principle, be inferred through modelling of stellar tidal streams in the Milky Way halo. The brightest and the longest of these, the Sagittarius stream, reaches out to large Galactocentric distances and hence can deliver the tightest constraints on the Galaxy's potential. In this contribution, we revisit the idea that the Sagittarius Stream was formed from a rotating progenitor. We demonstrate that the angle between the disk's angular momentum and the progenitor's orbital angular momentum does not remain constant throughout the disruption. Instead, it undergoes a dramatic evolution caused, in part, by the changes in the progenitor's moment of inertia tensor. We show that, even in a spherical potential, the streams produced as a result of a disky dwarf disruption appear to be "precessing". Yet, this debris plane evolution is illusory as it is solely caused by the swaying and wobbling of the progenitor's disk. Stream plane precession is therefore not an unambiguous indicator of asphericity of the dark halo.
• ### Chemodynamic subpopulations of the Carina dwarf galaxy(1601.01589)

Jan. 7, 2016 astro-ph.GA
We study the chemodynamical properties of the Carina dwarf spheroidal by combining an intermediate spectroscopic resolution dataset of more than 900 red giant and red clump stars, with high-precision photometry to derive the atmospheric parameters, metallicities and age estimates for our targets. Within the red giant branch population, we find evidence for the presence of three distinct stellar sub-populations with different metallicities, spatial distributions, kinematics and ages. As in the Fornax and Sculptor dwarf spheroidals, the subpopulation with the lowest average metallicity is more extended and kinematically hotter than all other populations. However, we identify an inversion in the parallel ordering of metallicity, kinematics and characteristic length scale in the two most metal rich subpopulations, which therefore do not contribute to a global negative chemical gradient. Contrary to common trends in the chemical properties with radius, the metal richest population is more extended and mildly kinematically hotter than the main component of intermediate metallicity. More investigations are required to ascertain the nature of this inversion, but we comment on the mechanisms that might have caused it.
• ### On the Alignment of the Stress Tensor in Galaxies(1509.08794)

Nov. 24, 2015 astro-ph.GA
We show that, provided the principal axes of the second velocity moment tensor of a stellar population are generally unequal and are oriented perpendicular to a set of orthogonal surfaces at each point, then those surfaces must be confocal quadric surfaces and the potential must be separable or Stackel. This is true under the mild assumption that the even part of the distribution function is invariant under time reversal $v_i \rightarrow -v_i$ of each velocity component. In particular, if the second velocity moment tensor is everywhere exactly aligned in spherical polar coordinates, then the potential must be of Stackel form (excepting degenerate cases where two or more of the semiaxes of ellipsoid are everywhere the same). The theorem also has consequences for alignment in cylindrical polar coordinates, which is used in the popular Jeans Anisotropic Models (JAM). We analyse data on the radial velocities and proper motions of a sample of $\sim 7400$ stars in the stellar halo of the Milky Way. We provide the distributions of the tilt angles or misalignments from the spherical polar coordinate systems. We show that in this sample the misalignment is always small (usually within $3^\circ$) for radii between 7 and 12 kpc. The velocity anisotropy is very radially biased ($\beta \approx 0.7$), and almost invariant across the volume in our study. Finally, we construct a triaxial stellar halo in a triaxial NFW dark matter halo using a made-to-measure method. Despite the triaxiality of the potential, the velocity ellipsoid of the stellar halo is nearly spherically aligned within $\sim6^\circ$ for large regions of space, particularly outside the scale radius of the stellar halo. We conclude that the second velocity moment ellipsoid can be close to spherically aligned for a much wider class of potentials than the strong constraints that arise from exact alignment might suggest.
• ### The Doubloon Models of Dark Haloes(1503.03689)

March 12, 2015 astro-ph.CO, astro-ph.GA
A family of spherical halo models with flat circular velocity curves is presented. This includes models in which the rotation curve has a finite central value but declines outwards (like the Jaffe model). It includes models in which the rotation curve is rising in the inner parts, but flattens asymptotically (like the Binney model). The family encompasses models with both finite and singular (cuspy) density profiles. The self-consistent distribution function depending on binding energy $E$ and angular momentum $L$ is derived and the kinematical properties of the models discussed. These really describe the properties of the total matter (both luminous and dark). For comparison with observations, it is better to consider tracer populations of stars. These can be used to represent elliptical galaxies or the spheroidal components of spiral galaxies. Accordingly, we study the properties of tracers with power-law or Einasto profiles moving in the doubloon potential. Under the assumption of spherical alignment, we provide a simple way to solve the Jeans equations for the velocity dispersions. This choice of alignment is supported by observations on the stellar halo of the Milky Way. Power-law tracers have prolate spheroidal velocity ellipsoids everywhere. However, this is not the case for Einasto tracers, for which the velocity ellipsoids change from prolate to oblate spheroidal near the pole. Asymptotic forms of the velocity distributions close to the escape speed are also derived, with an eye to application to the high velocity stars in the Milky Way. Power-law tracers have power-law or Maxwellian velocity distributions tails, whereas Einasto tracers have super-exponential cut-offs.
• ### Dipping Our Toes in the Water: First Models of GD-1 as a Stream(1502.00484)

Feb. 2, 2015 astro-ph.CO, astro-ph.GA
We present a model for producing tidal streams from disrupting progenitors in arbitrary potentials, utilizing the idea that the majority of stars escape from the progenitor's two Lagrange points. The method involves releasing test particles at the Lagrange points as the satellite orbits the host and dynamically evolving them in the potential of both host and progenitor. The method is sufficiently fast to allow large-dimensional parameter exploration using Monte Carlo methods. We provide the first direct modelling of 6-D stream observations -- assuming a stream rather than an orbit -- by applying our methods to GD-1. This is a kinematically cold stream spanning $60^{\circ}$ of the sky and residing in the outer Galaxy $\approx 15$ kpc distant from the centre. We assume the stream moves in a flattened logarithmic potential characterised by an asymptotic circular velocity $v_0$ and a flattening $q$. We recover values of normalisation $v_0$ = $227.2^{+15.6}_{-18.2}$ kms$^{-1}$ and flattening $q$ = $0.91^{+0.04}_{-0.1}$, if the stream is assumed to leading, and $v_0$ = $226.5^{+17.9}_{-17.0}$ kms$^{-1}$, $q$ = $0.90^{+0.05}_{-0.09}$, if it is assumed to be trailing. This can be compared to the values $v_0 = 224 \pm 13$ kms$^{-1}$ and $q= 0.87^{+0.07}_{-0.04}$ obtained by Koposov et al (2010) using the simpler technique of orbit fitting. Although there are differences between stream and orbit fitting, we conclude that orbit fitting can provide accurate results given the current quality of the data, at least for this kinematically cold stream in this logarithmic model of the Galaxy.
• ### "Skinny Milky Way, Please", says Sagittarius(1406.2243)

Sept. 24, 2014 astro-ph.GA
Motivated by recent observations of the Sagittarius stream, we devise a rapid algorithm to generate faithful representations of the centroids of stellar tidal streams formed in a disruption of a progenitor of an arbitrary mass in an arbitrary potential. Our method works by releasing swarms of test particles at the Lagrange points around the satellite and subsequently evolving them in a combined potential of the host and the progenitor. We stress that the action of the progenitor's gravity is crucial to making streams that look almost indistinguishable from the N-body realizations, as indeed ours do. The method is tested on mock stream data in three different Milky Way potentials with increasing complexity, and is shown to deliver unbiased inference on the Galactic mass distribution out to large radii. When applied to the observations of the Sagittarius stream, our model gives a natural explanation of the stream's apocentric distances and the differential orbital precession. We, therefore, provide a new independent measurement of the Galactic mass distribution beyond 50 kpc. The Sagittarius stream model favours a light Milky Way with the mass 4.1 +/- 0.4 x 10^11 M_sun at 100 kpc, which can be extrapolated to 5.6 +/- 1.2 x 10^11 M_sun at 200 kpc. Such a low mass for the Milky Way Galaxy is in good agreement with estimates from the kinematics of halo stars and from the satellite galaxies (once Leo I is removed from the sample). It entirely removes the "Too Big To Fail Problem".
• ### Microlensing, Brown Dwarfs and GAIA(1409.1788)

Sept. 5, 2014 astro-ph.GA, astro-ph.SR
The GAIA satellite can precisely measure the masses of nearby brown dwarfs and lower main sequence stars by the microlensing effect. The scientific yield is maximised if the microlensing event is also followed with ground-based telescopes to provide densely sampled photometry. There are two possible strategies. First, ongoing events can be triggered by photometric or astrometric alerts by GAIA. Second, events can be predicted using known high proper motion stars as lenses. This is much easier, as the location and time of an event can be forecast. Using the GAIA source density, we estimate that the sample size of high proper motion ($>300$ mas yr$^{-1}$) brown dwarfs needed to provide predictable events during the 5 year mission lifetime is surprisingly small, only of the order of a hundred. This is comparable to the number of high proper motion brown dwarfs already known from the work of the UKIDSS Large Area Survey and the all-sky WISE satellite. Provided the relative parallax of the lens and the angular Einstein radius can be recovered from astrometric data, then the mass of the lens can be found. Microlensing provides the only way of measuring the masses of individual objects irrespective oftheir luminosity. So, microlensing with GAIA is the best way to carry out an inventory of masses in the brown dwarf regime.
• ### Extremely Flat Haloes and the Shape of the Galaxy(1406.3710)

June 14, 2014 astro-ph.GA
We present a set of highly flattened galaxy models with asymptotically constant rotation curves. The mass density in the equatorial plane falls like (distance)$^{-1}$ at large radii. Although the inner equidensity contours may be spherical, oblate or prolate, the outer parts are always severely flattened. The elongated shape is supported by rotation or tangential velocity anisotropy. The models are thickened Mestel discs, and form a previously undiscovered part of the Miyamoto & Nagai sequence of flattened galaxies. The properties of the models -- axis ratios, velocity dispersions, streaming velocities and distribution functions -- are all discussed in some detail. We pose the question: are extremely flattened or disk-like haloes possible for the Milky Way galaxy? This has never been examined before, as very flattened halo models were not available. We fit the rotation curve and the vertical kinematics of disc stars in the solar neighbourhood to constrain the overall shape of the Galaxy. Denoting the ratio of polar axis to major axis by $q$, we show that models with $q\lesssim 0.57$ cannot simultaneously reproduce the in-plane and out-of-plane constraints. The kinematics of the Sagittarius galaxy also strongly disfavour models with high flattening, as the orbital plane precession is too great and the height reached above the Galactic plane is too small. At least for our Galaxy, the dark halo cannot be flatter than E4 (or axis ratio $q \sim 0.57$) at the Solar circle. Models in which the dark matter is accounted for by a massive baryonic disc or by decaying neutrinos are therefore ruled out by constraints from the rotation curve and the vertical kinematics.
• ### A Very Simple Cusped Halo Model(1406.3730)

June 14, 2014 astro-ph.GA
We introduce a very simple model of a dark halo. It is a close relative of Hernquist's model, being generated by the same transformation but this time applied to the logarithmic potential rather than the point mass. The density is proportional to (distance)$^{-1}$ at small radii, whilst the rotation curve is flat at large radii. Isotropic and radially anisotropic distributions functions are readily found, and the intrinsic and line of sight kinematical quantities are available as simple formulae. We also provide an analytical approximation to the Hamiltonian as a function of the actions. As an application, we study the kinematic properties of stellar haloes and tracers in elliptical galaxies. We show that the radial velocity dispersion of a power-law population in a galaxy with a flat rotation curve always tends to the constant value. This holds true irrespective of the anisotropy or the lengthscales of the dark or luminous matter. An analogous result holds for the line of sight or projected velocity dispersion of a power-law surface brightness profile. The radial velocity dispersion of Population II stars in the Milky Way is a strongly declining function of Galactocentric radius. So, if the rotation curve is flat, we conclude that the stellar halo density cannot follow a power-law at large radii, but must decrease more sharply (like an Einasto profile) or be abruptly truncated at large radii. Both the starcount and kinematic data of the Milky Way stellar halo are well-represented by an Einasto profile with index $m\approx 2$ and effective radius $\approx 20$ kpc.
• ### Precession of the Sagittarius stream(1301.7069)

Dec. 10, 2013 astro-ph.GA
Using a variety of stellar tracers -- blue horizontal branch stars, main-sequence turn-off stars and red giants -- we follow the path of the Sagittarius (Sgr) stream across the sky in Sloan Digital Sky Survey data. Our study presents new Sgr debris detections, accurate distances and line-of-sight velocities that together help to shed new light on the puzzle of the Sgr tails. For both the leading and the trailing tail, we trace the points of their maximal extent, or apo-centric distances, and find that they lie at $R^L$ = 47.8 $\pm$ 0.5 kpc and $R^T$ = 102.5 $\pm$ 2.5 kpc respectively. The angular difference between the apo-centres is 93.2 $\pm$ 3.5 deg, which is smaller than predicted for logarithmic haloes. Such differential orbital precession can be made consistent with models of the Milky Way in which the dark matter density falls more quickly with radius. However, currently, no existing Sgr disruption simulation can explain the entirety of the observational data. Based on its position and radial velocity, we show that the unusually large globular cluster NGC 2419 can be associated with the Sgr trailing stream. We measure the precession of the orbital plane of the Sgr debris in the Milky Way potential and show that, surprisingly, Sgr debris in the primary (brighter) tails evolves differently to the secondary (fainter) tails, both in the North and the South.
• ### The core size of the Fornax dwarf Spheroidal(1210.3157)

Oct. 11, 2012 astro-ph.CO
We exploit the detection of three distinct stellar subpopulations in the red giant branch of the Fornax dwarf Spheroidal to probe its density distribution. This allows us to resolve directly the evolution with radius of the dark matter mass profile. We find that a cored dark matter halo provides a perfect fit to the data, being consistent with all three stellar populations well within 1-sigma, and for the first time we are able to put constraints on the core size of such a halo. With respect to previous work, we do not strengthen the statistical exclusion of a dark matter cusp in Fornax, but we find that Navarro-Frenk-White haloes would be required to have unrealistically large scale radii in order to be compatible with the data, hence low values of the concentration parameter. We are then forced to conclude that the Fornax dwarf Spheroidal sits within a dark matter halo having a constant density core, with a core size of between 0.6 and 1.8 kpc.
• ### Mismatch and Misalignment: Dark Haloes and Satellites of Disc Galaxies(1101.0816)

June 13, 2011 astro-ph.CO, astro-ph.GA
(Abridged) We study the phase-space distribution of satellite galaxies associated with late-type galaxies in the GIMIC suite of simulations. GIMIC consists of re-simulations of 5 cosmologically representative regions from the Millennium simulation, which have higher resolution and incorporate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions (r ~ 0.1r_200) of the dark matter halo, both in shape and angular momentum, there can be substantial misalignments at larger radii (r ~ r_200). Misalignments of > 45 deg are seen in ~ 30% of our sample. We find that the satellite population aligns with the shape (and angular momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak owing to the mismatch between the disc and dark matter halo. Roughly 20% of the satellite systems with ten bright galaxies within r_200 exhibit a polar spatial alignment with respect to the galaxy --- an orientation reminiscent of the classical satellites of the Milky Way. We find that a small fraction (~10%) of satellite systems show evidence for rotational support which we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). This preference towards co-rotation is stronger in the inner regions of the halo where the most massive satellites accreted at relatively early times are located. We attribute the anisotropic spatial distribution and angular momentum bias of the satellites at z=0 to their directional accretion along the major axes of the dark matter halo. The satellite galaxies have been accreted relatively recently compared to the dark matter mass and have experienced less phase-mixing and relaxation --- the memory of their accretion history can remain intact to z=0.
• ### Chandrasekhar and Modern Stellar Dynamics(1102.4514)

Feb. 22, 2011 astro-ph.GA
Stellar dynamics occupied Chandrasekhar's interest for a brief interlude between his more prolonged studies of stellar structure and radiative transfer. This paper traces the history of one of his ideas -- namely, that the shape of the galactic potential controls the orientation of the stellar velocity dispersion tensor. It has its roots in papers by Eddington (1915) and Chandrasekhar (1939), and provoked a fascinating dispute between these two great scientists -- less well-known than their famous controversy over the white dwarf stars. In modern language, Eddington claimed that the integral curves of the eigenvectors of the velocity dispersion tensor provide a one-dimensional foliation into mutually orthogonal surfaces. Chandrasekhar challenged this, and explicitly constructed a counter-example. In fact, the work of neither of these great scientists was without flaws, though further developments in stellar dynamics were to ultimately draw more on Eddington's insight than Chandrasekhar's. We conclude with a description of modern attempts to measure the orientation of the velocity dispersion tensor for populations in the Milky Way Galaxy, a subject that is coming into its own with the dawning of the Age of Precision Astrometry.