Source author record

N. W. Evans

N. W. Evans appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph astro-ph.SR astro-ph.HE math-ph math.MP physics.class-ph physics.hist-ph

Catalog footprint

What is connected

80works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A 6-d View of Stellar Shells

Stellar shells are low surface brightness features, created during nearly head-on galaxy mergers from the debris of the tidally disrupted satellite. Here, we investigate the formation and evolution mechanism of shells in six dimensions (3d positions and velocities). We propose a new description in action-angle coordinates which condenses the seemingly complex behaviour of an expanding shell system into a simple picture, and stresses the crucial role of the existence of different stripping episodes in the properties of shells. Based on our findings, we construct a method for constraining the potential of the host galaxy and the average epoch of stripping. The method is applicable even if the shells cannot be identified or isolated from the data, or if the data are heavily contaminated with additional foreground stars. These results open up a new possibility to study the ancient merger that built the Milky Way Galaxy's stellar halo.

preprint2020arXiv

Donald Lynden-Bell: A Biographical Memoir

Donald Lynden-Bell's many contributions to astrophysics encompass general relativity, galactic dynamics, telescope design and observational astronomy. In the 1960s, his papers on stellar dynamics led to fundamental insights into the equilibria of elliptical galaxies, the growth of spiral patterns in disc galaxies and the stability of differentially rotating, self-gravitating flows. Donald introduced the ideas of `violent relaxation' and `the gravothermal catastrophe' in pioneering work on the thermodynamics of galaxies and negative heat capacities. He shared the inaugural Kavli Prize in Astrophysics in 2008 for his contributions to our understanding of quasars. His prediction that dead quasars or supermassive black holes may reside in the nuclei of nearby galaxies has been confirmed by multiple pieces of independent evidence. His work on accretion discs led to new insights into their workings, as well as the realisation that the infrared excess in T Tauri stars was caused by protostellar discs around these young stars. He introduced the influential idea of monolithic collapse of a gas cloud as a formation mechanism for the Milky Way Galaxy. As this gave way to modern ideas of merging and accretion as drivers of galaxy formation, Donald was the first to realise the importance of tidal streams as measures of the past history and present day gravity field of the Galaxy. Though primarily a theorist, Donald participated in one of the first observational programs to measure the large-scale streaming of nearby galaxies. This led to the discovery of the `Great Attractor'. The depth and versatility of his contributions mark Donald out as one of the most influential and pre-eminent astronomers of his day.

preprint2020arXiv

The Early Merger that Made the Galaxy's Stellar Halo

The last two years have seen widespread acceptance of the idea that the Milky Way halo was largely created in an early (8-10 Gyr ago) and massive ($> 10^{10} M_\odot$) merger. The roots of this idea pre-date the Gaia mission, but the exquisite proper motions available from Gaia have made the hypothesis irresistible. We trace the history of this idea, reviewing the series of papers that led to our current understanding.

preprint2019arXiv

Age Demographics of the Milky Way Disk and Bulge

We use the extensive $Gaia$ Data Release 2 set of Long Period Variables to select a sample of Oxygen-rich Miras throughout the Milky Way disk and bulge for study. Exploiting the relation between Mira pulsation period and stellar age/chemistry, we slice the stellar density of the Galactic disk and bulge as a function of period. We find the morphology of both components evolves as a function of stellar age/chemistry with the stellar disk being stubby at old ages, becoming progressively thinner and more radially extended at younger stellar ages, consistent with the picture of inside-out and upside-down formation of the Milky Way's disk. We see evidence of a perturbed disk, with large-scale stellar over-densities visible both in and away from the stellar plane. We find the bulge is well modelled by a triaxial boxy distribution with an axis ratio of $\sim [1:0.4:0.3]$. The oldest of the Miras ($\sim$ 9-10 Gyr) show little bar-like morphology, whilst the younger stars appear inclined at a viewing angle of $\sim 21^{\circ}$ to the Sun-Galactic Centre line. This suggests that bar formation and buckling took place 8-9 Gyr ago, with the older Miras being hot enough to avoid being trapped by the growing bar. We find the youngest Miras to exhibit a strong peanut morphology, bearing the characteristic X-shape of an inclined bar structure.

preprint2019arXiv

Models of Bars with Exponential Density Profiles

We present a new model for galactic bars with exponentially falling major axis luminosity profiles and Gaussian cross-sections. This is based on the linear superposition of Gaussian potential-density pairs with an exponential weight function, using an extension of the method originally introduced by Long & Murali (1992). We compute the density, potential and forces, using Gaussian quadrature. These quantities are given as explicit functions of position. There are three independent scaled bar parameters that can be varied continuously to produce bespoke bars of a given mass and shape. We categorise the effective potential by splitting a reduced parameter space into six regions. Unusually, we find bars with three stable Lagrange points on the major axis are possible. Our model reveals a variety of unexpected orbital structure, including a bifurcating $x_1$ orbit coexisting with a stable $x_4$ orbit. Propeller orbits are found to play a dominant role in the orbital structure, and we find striking similarities between our bar configuration and the model of Kaufmann & Contopoulos (1996). We find a candidate orbital family, sired from the propeller orbits, that may be responsible for the observed high velocity peaks in the Milky Way's bar. As a cross-check, we inspect, for the first time, the proper motions of stars in the high velocity peaks, which also match our suggested orbital family well. This work adds to the increasing body of evidence that real galactic bars may be supported at least partly by propeller orbits rather than solely the $x_1$ family.

preprint2016arXiv

A Dipole on the Sky: Predictions for Hypervelocity Stars from the Large Magellanic Cloud

We predict the distribution of hypervelocity stars (HVSs) ejected from the Large Magellanic Cloud (LMC), under the assumption that the dwarf galaxy hosts a central massive black hole (MBH). For the majority of stars ejected from the LMC the orbital velocity of the LMC has contributed a significant fraction of their galactic rest frame velocity, leading to a dipole density distribution on the sky. We quantify the dipole using spherical harmonic analysis and contrast with the monopole expected for HVSs ejected from the Galactic Center. There is a tendril in the density distribution that leads the LMC which is coincident with the well-known and unexplained clustering of HVSs in the constellations of Leo and Sextans. Our model is falsifiable, since it predicts that Gaia will reveal a large density of HVSs in the southern hemisphere.

preprint2016arXiv

A Tail of Two Populations: Chemo-dynamics of the Sagittarius Stream and Implications for its Original Mass

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.

preprint2016arXiv

Chemodynamic subpopulations of the Carina dwarf galaxy

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.

preprint2016arXiv

Discy dwarf disruption and the shape of the Galactic halo

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.

preprint2016arXiv

Is the Dark Halo of the Milky Way Prolate?

We introduce the flattening equation, which relates the shape of the dark halo to the angular velocity dispersions and the density of a tracer population of stars. It assumes spherical alignment of the velocity dispersion tensor, as seen in the data on stellar halo stars in the Milky Way. The angular anisotropy and gradients in the angular velocity dispersions drive the solutions towards prolateness, whilst the gradient in the stellar density is a competing effect favouring oblateness. We provide an efficient numerical algorithm to integrate the flattening equation. Using tests on mock data, we show that the there is a strong degeneracy between circular speed and flattening, which can be circumvented with informative priors. Therefore, we advocate the use of the flattening equation to test for oblateness or prolateness, though the precise value of the flattening $q$ can only be measured with the addition of the radial Jeans equation. We apply the flattening equation to a sample extracted from the Sloan Digital Sky Survey of $\sim 15000$ halo stars with full phase space information and errors. We find that between Galactocentric radii of 5 and 10 kpc, the shape of the dark halo is prolate, whilst even mildly oblate models are disfavoured. Strongly oblate models are ruled out. Specifically, for a logarithmic halo model, if the asymptotic circular speed $v_0$ lies between $210$ and 250 kms$^{-1}$, then we find the axis ratio of the equipotentials $q$ satisfies $1.5 \lesssim q \lesssim 2$.

preprint2016arXiv

Is there substructure around M87?

We present a general method to identify infalling substructure in discrete datasets with position and line-of-sight velocity data. We exploit the fact that galaxies falling onto a brightest cluster galaxy (BCG) in a virialised cluster, or dwarf satellites falling onto a central galaxy like the Milky Way, follow nearly radial orbits. If the orbits are exactly radial, we show how to find the probability distribution for a satellite's energy, given a tracer density for the satellite population, by solving an Abel integral equation. This is an extension of Eddington (1916)'s classical formula for the isotropic distribution function. When applied to a system of galaxies, clustering in energy space can then be quantified using the Kullback-Leibler divergence, and groups of objects can be identified which, though separated in the sky, may be falling in on the same orbit. This method is tested using mock data and applied to the satellite galaxy population around M87, the BCG in Virgo, and a number of associations are found which may represent infalling galaxy groups.

preprint2016arXiv

Simple J-Factors and D-Factors for Indirect Dark Matter Detection

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 $γ$ and outer slope $β$. 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.

preprint2016arXiv

The Gaia-ESO survey: Metal-rich bananas in the bulge

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.

preprint2015arXiv

Dipping Our Toes in the Water: First Models of GD-1 as a Stream

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.

preprint2015arXiv

Haloes light and dark: dynamical models of the stellar halo and constraints on the mass of the Galaxy

We develop a flexible set of action-based distribution functions (DFs) for stellar haloes. The DFs have five free parameters, controlling the inner and outer density slope, break radius, flattening and anisotropy respectively. The DFs generate flattened stellar haloes with a rapidly varying logarithmic slope in density, as well as a spherically aligned velocity ellipsoid with a long axis that points towards the Galactic centre - all attributes possessed by the stellar halo of the Milky Way. We use our action-based distribution function to model the blue horizontal branch stars extracted from the Sloan Digital Sky Survey as stellar halo tracers in a spherical Galactic potential. As the selection function is hard to model, we fix the density law from earlier studies and solve for the anisotropy and gravitational potential parameters. Our best fit model has a velocity anisotropy that becomes more radially anisotropic on moving outwards. It changes from $β\approx 0.4$ at Galactocentric radius of 15 kpc to $\approx 0.7$ at 60 kpc. This is a gentler increase than is typically found in simulations of stellar haloes built from the mutiple accretion of smaller systems. We find the potential corresponds to an almost flat rotation curve with amplitude of $\approx 200$ kms$^{-1}$ at these distances. This implies an enclosed mass of $\approx 4.5 \times 10^{11} M_\odot$ within a spherical shell of radius 50 kpc.

preprint2015arXiv

Made-to-Measure Dark Matter Haloes, Elliptical Galaxies and Dwarf Galaxies in Action Coordinates

We provide a family of action-based distribution functions (DFs) for the double-power law family of densities often used to model galaxies. The DF itself is a double-power law in combinations of the actions, and reduces to the known limits in the case of a pure power-law at small and large radii. Our method enables the velocity anisotropy of the model to be tuned, and so the anisotropy in the inner and outer parts can be specified for the application in hand. We provide self-consistent DFs for the Hernquist and Jaffe models - both with everywhere isotropic velocity dispersions, and with kinematics that gradually becomes more radially anisotropic on moving outwards. We also carry out this exercise for a cored dark-matter model. These are tailored to represent dark haloes and elliptical galaxies respectively with kinematic properties inferred from simulations or observational data. Finally, we relax a cored luminous component within a dark matter halo to provide a self-consistent model of a dwarf spheroidal embedded in dark matter. The DFs provide us with non-rotating spherical stellar systems, but one of the virtues of working with actions is the relative ease with which such models can be converted into axisymmetry and triaxiality.

preprint2015arXiv

On the Alignment of the Stress Tensor in Galaxies

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 ($β\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.

preprint2015arXiv

Resonant Clumping and Substructure in Galactic Discs

We describe a method to extract resonant orbits from N-body simulations exploiting the fact that they close in a frame rotating with a constant pattern speed. Our method is applied to the N-body simulation of the Milky Way by Shen et al. (2010). This simulation hosts a massive bar, which drives strong resonances and persistent angular momentum exchange. Resonant orbits are found throughout the disc, both close to the bar itself and out to the very edges of the disc. Using Fourier spectrograms, we demonstrate that the bar is driving kinematic substructure even in the very outer parts of the disc. We identify two major orbit families in the outskirts of the disc that make significant contributions to the kinematic landscape, namely the m:l = 3:-2 and 1:-1 families resonating with the pattern speed of the bar. A mechanism is described that produces bimodal distributions of Galactocentric radial velocities at selected azimuths in the outer disc. It occurs as a result of the temporal coherence of particles on the 3:-2 resonant orbits, which causes them to arrive simultaneously at pericentre or apocentre. This resonant clumping, due to the in-phase motion of the particles through their epicycle, leads to both inward and outward moving groups which belong to the same orbital family and consequently produce bimodal radial velocity distributions. This is a possible explanation of the bimodal velocity distributions observed towards the Galactic anti-Centre by Liu et al. (2012). Another consequence is that transient overdensities appear and dissipate (in a symmetric fashion) on timescales equal to the their epicyclic period resulting in a periodic pulsing of the disc's surface density.

preprint2015arXiv

The Doubloon Models of Dark Haloes

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.

preprint2015arXiv

What can the alignments of the velocity moments tell us about the nature of the potential?

We prove that, if the time-independent distribution function $F(v;x)$ of a steady-state stellar system is symmetric under velocity inversion such that $F(-v_1,v_2,v_3;x)=F(v_1,v_2,v_3;x)$ and the same for $v_2$ and $v_3$, where $(v_1,v_2,v_3)$ is the velocity component projected onto an orthogonal frame, then the potential within which the system is in equilibrium must be separable (i.e. the Staeckel potential). Furthermore, we find that the Jeans equations imply that, if all mixed second moments of the velocity vanish, that is, $\langle v_iv_j\rangle=0$ for any $i\ne j$, in some Staeckel coordinate system and the only non-vanishing fourth moments in the same coordinate are those in the form of $\langle v_i^4\rangle$ or $\langle v_i^2v_j^2\rangle$, then the potential must be separable in the same coordinates. Finally we also show that all second and fourth velocity moments of tracers with an odd power to the radial component $v_r$ being zero is a sufficient condition to guarantee the potential to be of the form $Φ=f(r)+r^{-2}g(θ,ϕ)$.

preprint2014arXiv

"Skinny Milky Way, Please", says Sagittarius

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".

preprint2014arXiv

A Very Simple Cusped Halo Model

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.

preprint2014arXiv

ATLAS lifts the Cup: Discovery of a New Milky Way satellite in Crater

We announce the discovery of a new Galactic companion found in data from the ESO VST ATLAS survey, and followed up with deep imaging on the 4m William Herschel Telescope. The satellite is located in the constellation of Crater (the Cup) at a distance of $\sim$ 170 kpc. Its half-light radius is $r_h=30$ pc and its luminosity is $M_V=-5.5$. The bulk of its stellar population is old and metal-poor. We would probably have classified the newly discovered satellite as an extended globular cluster were it not for the presence of a handful of Blue Loop stars and a sparsely populated Red Clump. The existence of the core helium burning population implies that star-formation occurred in Crater perhaps as recently as 400 Myr ago. No globular cluster has ever accomplished the feat of prolonging its star-formation by several Gyrs. Therefore, if our hypothesis that the blue bright stars in Crater are Blue Loop giants is correct, the new satellite should be classified as a dwarf galaxy with unusual properties. Note that only ten degrees to the North of Crater, two ultra-faint galaxies Leo IV and V orbit the Galaxy at approximately the same distance. This hints that all three satellites may once have been closely associated before falling together into the Milky Way halo.

preprint2014arXiv

Dynamical Models of Elliptical Galaxies -- I. Simple Methods

We study dynamical models for elliptical galaxies, deriving the projected kinematic profiles in a form that is valid for general surface-brightness profiles and (spherical) total mass profiles, without the need for any explicit deprojection. We then show that an almost flat rotation curve, combined with modest velocity anisotropy, is already sufficient to recover the kinematic profiles of nearby ellipticals. As an application, we provide two different sets of mass estimators for elliptical galaxies, based on either the velocity dispersion at a specific location near the effective radius, or the aperture-averaged velocity dispersion. In the large aperture (virial) limit, mass estimators are naturally independent of anisotropy. The spherical mass enclosed within the effective radius $R_{\rm e}$ can be estimated as $2.4 R_{\rm e} \langle σ^{2}_{\rm p} \rangle/ G$, where $\langle σ^2_{\rm p} \rangle$ is the average of the squared velocity dispersion over a finite aperture. This formula does not depend on assumptions such as mass-follows-light, and is a compromise between the cases of small and large apertures sizes. Its general agreement with results from other methods in the literature makes it a reliable means to infer masses in the absence of detailed kinematic information. If on the other hand the velocity dispersion profile is available, tight mass estimates can be found that are independent of the mass-model and anisotropy profile (within $\approx$ 10% accuracy). Explicit formulae are given for small anisotropy, large radii and/or power-law total densities. Motivated by recent observational claims, we also discuss the issue of weak homology of elliptical galaxies, emphasizing the interplay between morphology and orbital structure.

preprint2014arXiv

Dynamical Models of Elliptical Galaxies -- II. M87 and its Globular Clusters

We study the Globular Cluster (GC) system of the nearby elliptical galaxy M87 using the newly available dataset with accurate kinematics (Strader et al.2011). We find evidence for three distinct sub-populations of GCs in terms of colours, kinematics and radial profiles. We show that a decomposition into three populations (blue, intermediate and red GCs) is statistically preferred to one with two or four populations, and relate them to the stellar profile. We exploit the sub-populations to derive dynamical constraints on the mass and Dark Matter (DM) content of M87 out to $\sim100$ kpc. We use a class of global mass-estimators (from Paper I), obtaining mass measurements at different locations. M87's DM fraction changes from $\approx$0.2 at the starlight's effective radius (6 kpc) to $\approx$0.95 at the distance probed by the most extended, blue GCs (135 kpc). We supplement this with \textit{virial decompositions}, exploiting the dynamical model to produce a separation into multiple components. These yield the luminous mass as $5.5^{+1.5}_{-2.0}\times 10^{11}M_\odot$ and the DM within 135 kpc as $8.0^{+1.0}_{-4.0}\times 10^{12}M_\odot.$ The inner DM density behaves as $ρ\sim r^{-γ}$ with $γ\approx 1.6$. This is steeper than the cosmologically preferred cusp $ρ\sim r^{-1},$ providing evidence of DM contraction. Finally, we combine the GC separation into three sub-populations and the Jeans equations, obtaining information on the orbits of the GC system. The centrally concentrated red GCs exhibit tangential anisotropy, consistent with radial-orbit depletion by tidal shredding. The most extended blue GCs have an isotropic velocity dispersion tensor in the central parts, which becomes more tangential moving outwards, consistent with adiabatic contraction of the DM halo.

preprint2014arXiv

Extremely Flat Haloes and the Shape of the Galaxy

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.

preprint2014arXiv

Hamiltonians of Spherically Symmetric, Scale-Free Galaxies in Action-Angle Coordinates

We present a simple formula for the Hamiltonian in terms of the actions for spherically symmetric, scale-free potentials. The Hamiltonian is a power-law or logarithmic function of a linear combination of the actions. Our expression reduces to the well-known results for the familiar cases of the harmonic oscillator and the Kepler potential. For other power-laws, as well as for the singular isothermal sphere, it is exact for the radial and circular orbits, and very accurate for general orbits. Numerical tests show that the errors are always small, with mean errors across a grid of actions always less than 1 % and maximum errors less than 2.5 %. Simple first-order corrections can reduce mean errors to less than 0.6 % and maximum errors to less than 1 %. We use our new result to show that :[1] the misalignment angle between debris in a stream and a progenitor is always very nearly zero in spherical scale-free potentials, demonstrating that streams can be sometimes well approximated by orbits, [2] the effects of an adiabatic change in the stellar density profile in the inner regions of a galaxy weaken any existing 1/r density cusp, which is reduced to $1/r^{1/3}$. More generally, we derive the full range of adiabatic cusp transformations and show how to relate the starting cusp index to the final cusp index. It follows that adiabatic transformations can never erase a dark matter cusp.

preprint2014arXiv

Microlensing, Brown Dwarfs and GAIA

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.

preprint2014arXiv

On Asymmetric Distributions of Satellite Galaxies

We demonstrate that the asymmetric distribution of M31 satellites cannot be produced by tides from the Milky Way as such effects are too weak. However, loosely bound associations and groups of satellites can fall into larger haloes and give rise to asymmetries. We compute the survival times for such associations. We prove that the survival time is always shortest in Keplerian potentials, and can be ~ 3 times longer in logarithmic potentials. We provide an analytical formula for the dispersal time in terms of the size and velocity dispersion of the infalling structure. We show that, if an association of ~10 dwarfs fell into the M31 halo, its present aspect would be that of an asymmetric disk of satellites. We also discuss the case of cold substructure in the Andromeda II and Ursa Minor dwarfs.

preprint2014arXiv

The remnant of a merger between two dwarf galaxies in Andromeda II

Driven by gravity, massive structures like galaxies and clusters of galaxies are believed to grow continuously through hierarchical merging and accretion of smaller systems. Observational evidence of accretion events is provided by the coherent stellar streams crossing the outer haloes of massive galaxies, such as the Milky Way or Andromeda. At similar mass-scales, around $10^{11}$ solar masses in stars, further observational evidence of merging activity is also ample. Mergers of lower-mass galaxies are expected within the hierarchical process of galaxy formation, but have hitherto not been seen for galaxies with less than about $10^9$ solar masses in stars. Here, we report the kinematic detection of a stellar stream in one of the satellite galaxies of Andromeda, the dwarf spheroidal galaxy Andromeda II, which has a mass of only $10^7$ solar masses in stars. The properties of the stream show that we are observing the remnant of a merger between two dwarf galaxies. This had a dramatic influence on the dynamics of the remnant, which is now rotating around its projected major axis. The stellar stream in Andromeda II illustrates the scale-free character of the formation of galaxies, down to the lowest galactic mass scales.

preprint2013arXiv

The dark matter halo concentration and stellar initial mass function of a CASSOWARY group

We exploit the group environment of the CAmbridge Sloan Survey Of Wide ARcs in the skY (CASSOWARY) z=0.3 lens J2158+0257 to measure the group dynamical mass as a complement to the central dynamical and lensing mass constraints. Follow-up spectroscopy of candidate group members is performed using VLT/FORS2. From the resulting N=21 confirmed members we measure the group dynamical mass by calibrating an analytic tracer mass estimator with cosmological simulations. The luminosity weighted line-of-sight velocity dispersion and the Einstein radius of the lens are used as mass probes in the inner regions of the galaxy. Combining these three observational probes allows us to independently constrain the mass and concentration of the dark matter halo, in addition to the total stellar mass of the central galaxy. We find a dark matter halo in remarkably good agreement with simulations (log_10 M_200/M_sun = 14.2 +/- 0.2, c_200 = 4.4 (+1.6, -1.4)) and a stellar mass-to-light ratio which favors a Salpeter initial mass function ((M/L)* = 5.7 +/- 1.2). Our measurement of a normal halo concentration suggests that there is no discrepancy between simulations and observations at the group mass scale. This is in contrast to the cluster mass scale for which a number of studies have claimed over-concentrated halos. While the halo mass is robustly determined, and the halo concentration is not significantly affected by systematics, the resulting stellar mass-to-light ratio is sensitive to the choice of stellar parameters, such as density profile and velocity anisotropy.

preprint2013arXiv

Triaxial Cosmological Haloes and the Disc of Satellites

We construct simple triaxial generalisations of Navarro-Frenk-White haloes. The models have elementary gravitational potentials, together with a density that is cusped like 1/r at small radii and falls off like 1/r^3 at large radii. The ellipticity varies with radius in a manner that can be tailored to the user's specification. The closed periodic orbits in the planes perpendicular to the short and long axes of the model are well-described by epicyclic theory, and can be used as building blocks for long-lived discs. As an application, we carry out the simulations of thin discs of satellites in triaxial dark halo potentials. This is motivated by the recent claims of an extended, thin disc of satellites around the M31 galaxy with a vertical rms scatter of ~12 kpc and a radial extent of ~ 300 kpc (Ibata et al. 2013). We show that a thin satellite disc can persist over cosmological times if and only if it lies in the planes perpendicular to the long or short axis of a triaxial halo, or in the equatorial or polar planes of a spheroidal halo. In any other orientation, then the disc thickness doubles on ~5 Gyr timescales and so must have been born with an implausibly small vertical scaleheight.

preprint2012arXiv

A Troublesome Past: Chemodynamics of the Fornax dwarf spheroidal

We present compelling evidence for the complexity of the Fornax dwarf spheroidal. By disentangling three different stellar subpopulations in its red giant branch, we are able to study in detail the dependence between kinematics and metallicity. A well-defined ordering in velocity dispersion, spatial concentration, and metallicity is evident in the subpopulations. We also present evidence for a significant misalignment between the angular momentum vectors of the old and intermediate-age populations. According to the HST measurement of Fornax's proper motion, this corresponds to counter-rotation. These ingredients are used to construct a novel evolutionary history of the Fornax dwarf spheroidal, characterized as a late merger of a bound pair.

preprint2012arXiv

A Virial Core in the Sculptor Dwarf Spheroidal Galaxy

The projected virial theorem is applied to the case of multiple stellar populations in the nearby dwarf spheroidal galaxies. As each population must reside in the same gravitational potential, this provides strong constraints on the nature of the dark matter halo. We derive necessary conditions for two populations with Plummer or exponential surface brightnesses to reside in a cusped Navarro-Frenk-White (NFW) halo. We apply our methods to the Sculptor dwarf spheroidal, and show that there is no NFW halo compatible with the energetics of the two populations. The dark halo must possess a core radius of ~ 120 pc for the virial solutions for the two populations to be consistent. This conclusion remains true, even if the effects of flattening or self-gravity of the stellar populations are included.

preprint2012arXiv

Broken and unbroken: the Milky Way and M31 stellar haloes

We use the Bullock & Johnston suite of simulations to study the density profiles of L*-type galaxy stellar haloes. Observations of the Milky Way and M31 stellar haloes show contrasting results: the Milky Way has a `broken' profile, where the density falls off more rapidly beyond ~ 25 kpc, while M31 has a smooth profile out to 100 kpc with no obvious break. Simulated stellar haloes, built solely by the accretion of dwarf galaxies, also exhibit this behavior: some haloes have breaks, while others don't. The presence or absence of a break in the stellar halo profile can be related to the accretion history of the galaxy. We find that a break radius is strongly related to the build up of stars at apocentres. We relate these findings to observations, and find that the `break' in the Milky Way density profile is likely associated with a relatively early (~ 6-9 Gyr ago) and massive accretion event. In contrast, the absence of a break in the M31 stellar halo profile suggests that its accreted satellites have a wide range of apocentres. Hence, it is likely that M31 has had a much more prolonged accretion history than the Milky Way.

preprint2012arXiv

Broken degeneracies: the rotation curve and velocity anisotropy of the Milky Way halo

We use distant Blue Horizontal Branch stars with Galactocentric distances 16 < r/kpc < 48 as kinematic tracers of the Milky Way dark halo. We model the tracer density as an oblate, power-law embedded within a spherical power-law potential. Using a distribution function method, we estimate the overall power-law potential and the velocity anisotropy of the halo tracers. We measure the slope of the potential to be gamma ~ 0.4 and the overall mass within 50 kpc is ~ 4 x 10^11 M_sol. The tracer velocity anisotropy is radially biased with beta ~ 0.5, which is in good agreement with local solar neighbourhood studies. Our results provide an accurate outer circular velocity profile for the Milky Way and suggest a relatively high concentration dark matter halo (c_vir ~ 20).

preprint2012arXiv

Elliptical galaxy masses out to five effective radii: the realm of dark matter

We estimate the masses of elliptical galaxies out to five effective radii using planetary nebulae and globular clusters as tracers. A sample of 15 elliptical galaxies with a broad variation in mass is compiled from the literature. A distribution function-maximum likelihood analysis is used to estimate the overall potential slope, normalisation and velocity anisotropy of the tracers. We assume power-law profiles for the potential and tracer density and a constant velocity anisotropy. The derived potential power-law indices lie in between the isothermal and Keplerian regime and vary with mass: there is tentative evidence that the less massive galaxies have steeper potential profiles than the more massive galaxies. We use stellar mass-to-light ratios appropriate for either a Chabrier/KTG (Kroupa, Tout & Gilmore) or Salpeter initial mass function to disentangle the stellar and dark matter components. The fraction of dark matter within five effective radii increases with mass, in agreement with several other studies. We employ simple models to show that a combination of star formation efficiency and baryon extent are able to account for this trend. These models are in good agreement with both our measurements out to five effective radii and recent SLACS measurements within one effective radii when a universal Chabrier/KTG initial mass function is adopted.

preprint2012arXiv

Lensing and Dynamics in Two Simple Steps

We present a ready-to-use method to constrain the density distribution in early-type galaxy lenses. Assuming a power-law density profile, then joint use of the virial theorem and the lens equation yields simple formulae for the power-law index (or logarithmic density gradient). Any dependence on orbital anisotropy can be tightly constrained or even erased completely. Our results rely just on surface brightnesses and line-of-sight kinematics, making deprojection unnecessary. We revisit three systems that have already been examined in the literature (the Cosmic Horseshoe, the Jackpot and B1608+656) and provide our estimates. Finally, we show that the method yields a good approximation for the density profile even when the true profile is a broken power-law, albeit with a mild bias towards isothermality.

preprint2012arXiv

Line Profiles from Discrete Kinematic Data

We develop a method to extract the shape information of line profiles from discrete kinematic data. The Gauss-Hermite expansion, which is widely used to describe the line of sight velocity distributions extracted from absorption spectra of elliptical galaxies, is not readily applicable to samples of discrete stellar velocity measurements, accompanied by individual measurement errors and probabilities of membership. We introduce two parameter families of probability distributions describing symmetric and asymmetric distortions of the line profiles from Gaussianity. These are used as the basis of a maximum likelihood estimator to quantify the shape of the line profiles. Tests show that the method outperforms a Gauss-Hermite expansion for discrete data, with a lower limit for the relative gain of approx 2 for sample sizes N approx 800. To ensure that our methods can give reliable descriptions of the shape, we develop an efficient test to assess the statistical quality of the obtained fit. As an application, we turn our attention to the discrete velocity datasets of the dwarf spheroidals of the Milky Way. In Sculptor, Carina and Sextans the symmetric deviations are consistent with velocity distributions more peaked than Gaussian. In Fornax, instead, there is an evolution in the symmetric deviations of the line profile from a peakier to more flat-topped distribution on moving outwards. These results suggest a radially biased orbital structure for the outer parts of Sculptor, Carina and Sextans. On the other hand, tangential anisotropy is favoured in Fornax. This is all consistent with a picture in which Fornax may have had a different evolutionary history to Sculptor, Carina and Sextans.

preprint2012arXiv

Slicing and dicing the Milky Way disk in SDSS

We use the Stripe 82 proper motion catalogue of Bramich et al. (2008) to study the kinematics of Galactic disk stars in the solar neighborhood. We select samples of dwarf stars with reliable spectra and proper motions. They have cylindrical polar radius between 7 < R < 9 kpc, heights from the Galactic plane satisfying z < 2 kpc and span a range of metallicities -1.5 < [Fe/H] < 0. We develop a method for calculating and correcting for the halo contamination in our sample using the distribution of rotational velocities. Two Gaussians representing disk and halo populations are used to fit the radial and vertical velocity distributions via maximum likelihood methods. For the azimuthal velocities the same technique is used, except that a skewed non-Gaussian functional form now represents the disk velocity distribution. This enables us to compute the dispersions and cross-terms (the tilt and the vertex deviation) of the velocity ellipsoid as a function of height and metallicity. We also investigate the rotation lag of the disk, finding that the more metal-poor stars rotate significantly slower than the metal-rich stars. These samples provide important constraints on heating mechanisms in the Galactic disk and can be used for a variety of applications. We present one such application, employing the Jeans equations to provide a simple model of the potential close to the disk. Our model is in excellent agreement with others in the literature and provides an indication the disk, rather than the halo, dominates the circular speed at the solar neighborhood. We obtain a surface mass density within 1.1 kpc of around 66 M/pc^2 and estimate a local halo density of 0.015 M/pc^3 = 0.57 GeV/cm^3.

preprint2012arXiv

The cold veil of the Milky Way stellar halo

We build a sample of distant (D > 80 kpc) stellar halo stars with measured radial velocities. Faint ($20 < g <22$) candidate blue horizontal branch (BHB) stars were selected using the deep, but wide, multi-epoch Sloan Digital Sky Survey photometry. Follow-up spectroscopy for these A-type stars was performed using the VLT-FORS2 instrument. We classify stars according to their Balmer line profiles, and find 7 are bona fide BHB stars and 31 are blue stragglers (BS). Owing to the magnitude range of our sample, even the intrinsically fainter BS stars can reach out to D ~ 90 kpc. We complement this sample of A-type stars with intrinsically brighter, intermediate-age, asymptotic giant branch stars. A set of 4 distant cool carbon stars is compiled from the literature and we perform spectroscopic follow-up on a further 4 N-type carbon stars using the WHT-ISIS instrument. Altogether, this provides us with the largest sample to date of individual star tracers out to r ~ 150 kpc. We find that the radial velocity dispersion of these tracers falls rapidly at large distances and is surprisingly cold (sigma_r ~ 50-60 km/s) between 100-150 kpc. Relating the measured radial velocities to the mass of the Milky Way requires knowledge of the (unknown) tracer density profile and anisotropy at these distances. Nonetheless, by assuming the stellar halo stars between 50-150 kpc have a moderate density fall-off (with power-law slope alpha < 5) and are on radial orbits (sigma^2_t/sigma^2_r < 1), we infer that the mass within 150 kpc is less than 10^12 M_solar and suggest it probably lies in the range (5-10) x 10^11 M_solar. We discuss the implications of such a low mass for the Milky Way.

preprint2012arXiv

The core size of the Fornax dwarf Spheroidal

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.

preprint2012arXiv

The Sagittarius Streams in the Southern Galactic Hemisphere

The structure of the Sagittarius stream in the Southern Galactic hemisphere is analysed with the Sloan Digital Sky Survey Data Release 8. Parallel to the Sagittarius tidal track, but ~ 10deg away, there is another fainter and more metal-poor stream. We provide evidence that the two streams follow similar distance gradients but have distinct morphological properties and stellar populations. The brighter stream is broader, contains more metal-rich stars and has a richer colour-magnitude diagram with multiple turn-offs and a prominent red clump as compared to the fainter stream. Based on the structural properties and the stellar population mix, the stream configuration is similar to the Northern "bifurcation". In the region of the South Galactic Cap, there is overlapping tidal debris from the Cetus Stream, which crosses the Sagittarius stream. Using both photometric and spectroscopic data, we show that the blue straggler population belongs mainly to Sagittarius and the blue horizontal branch stars belong mainly to the Cetus stream in this confused location in the halo.

preprint2012arXiv

The stretching of Hercules

We present VLT/FORS2 spectroscopy of candidate blue horizontal branch (BHB) stars in the vicinity of the Hercules ultrafaint dwarf galaxy. We identify eight convincing Hercules BHB members, and a further five stars with similar systemic velocities to that of Hercules, but ~ 0.5 kpc from the centre of the galaxy along its major axis. It is likely that these stars once belonged to Hercules, but have been tidally stripped and are now unbound. We emphasise the usefulness of looking for any gradient in the systemic velocity of this stretched system, which would further support our interpretation of the origin of its elongated and distended morphology.

preprint2012arXiv

Virial Sequences for Thick Discs and Haloes: Flattening and Global Anisotropy

The virial theorem prescribes the ratio of the globally-averaged equatorial to vertical velocity dispersion of a tracer population in spherical and flattened dark haloes. This gives sequences of physical models in the plane of global anisotropy and flattening. The tracer may have any density, though there are particularly simple results for power-laws and exponentials. We prove the flattening theorem: for a spheroidally stratified tracer density with axis ratio q in a dark density potential with axis ratio g, the ratio of globally averaged equatorial to vertical velocity dispersion depends only on q/g. As the stellar halo density and velocity dispersion of the Milky Way are accessible to observations, this provides a new method for measuring the flattening of the dark matter. If the kinematics of the local halo subdwarfs are representative, then the Milky Way's dark halo is oblate with a flattening in the potential of g ~ 0.85, corresponding to a flattening in the dark matter density of ~ 0.7. The fractional pressure excess for power-law populations is roughly proportional to both the ellipticity and the fall-off exponent. Given the same pressure excess, if the density profile of one stellar population declines more quickly than that of another, then it must be rounder. This implies that the dual halo structure claimed by Carollo et al. (2007) for the Galaxy, a flatter inner halo and a rounder outer halo, is inconsistent with the virial theorem. For the thick disc, we provide formulae for the virial sequences of double-exponential discs in logarithmic and Navarro-Frenk-White (NFW) haloes. There are good matches to the observational data on the flattening and anisotropy of the thick disc if the thin disc is exponential with a short scalelength ~ 2.6 kpc and normalisation of 56 solar masses per square parsec, together with a logarithmic dark halo.

preprint2011arXiv

Accurate Masses for Navarro-Frenk-White Dark Matter Haloes

We consider the problem of estimating the virial mass of a dark halo from the positions and velocities of a tracer population. Although a number of general tools are available, more progress can be made if we are able to specify the functional form of the halo potential (although not its normalization). Here, we consider the particular case of the cosmologically motivated Navarro-Frenk-White (NFW) halo and develop two simple estimators. We demonstrate their effectiveness against numerical simulations and use them to provide new mass estimates of Carina, Fornax, Sculptor, and Sextans dSphs.

preprint2011arXiv

Chandrasekhar and Modern Stellar Dynamics

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.

preprint2011arXiv

Dark Matter Cores and Cusps: The Case of Multiple Stellar Populations in Dwarf Spheroidals

A number of dwarf spheroidal (dSph) galaxies are known to contain a more extended, metal-poor population with a flattish velocity dispersion profile, and a more concentrated, metal-rich population with a velocity dispersion declining with radius. The two populations can be modelled with Michie-King distribution functions (DFs) in the isothermal and in the sharply truncated limits, respectively. We argue that the truncation of the metal-rich population can be traced back to the spatial distribution of the star forming gas. Suppose D is the exponent of the first non-constant term in the Taylor expansion of the total potential at the center (D=1 for NFW halos, D=2 for cored halos). Then, we show that the ratio of the half-light radii of the populations raised to the D/2 power must be smaller than the ratio of the line-of-sight velocity dispersions. Specialising to the case of the Sculptor dSph, we develop a technique to fit simultaneously both populations with Michie-King DFs. This enables us to determine the mass profile of the Sculptor dSph with unprecedented accuracy in the radial range between 0.2 and 1.2 kpc. We show that cored halo models are preferred over cusped halo models, with a likelihood ratio test rejecting NFW models at any significance level higher than 0.05%. Even more worryingly, the best-fit NFW models require concentrations with c<20, which is not in the cosmologically preferred range for dwarf galaxies. We conclude that the kinematics of multiple populations in dSphs provides a substantial new challenge for theories of galaxy formation, with the weight of available evidence strongly against dark matter cusps at the centre.

preprint2011arXiv

Mass estimators in the Gaia era

Forthcoming astrometric missions such as the Gaia satellite will bring to the fore the problem of estimating the enclosed mass from a set of positions, radial velocities and proper motions of tracer stars. Here, we show how to construct the tracer mass estimator when the proper motion data are available in addition to the usual line-of-sight velocity data. Notably, the mass estimators do not require any assumption on the anisotropy, as it is naturally incorporated through the different components of proper motions. In addition, the separate treatment of the proper motions and the line-of-sight velocities is desirable because they are observationally independent and thus the propagation of the combined uncertainties is rather straightforward. The extension to projected data is also sketched, together with a possible application of measuring the masses of Galactic globular clusters.

preprint2011arXiv

Mismatch and Misalignment: Dark Haloes and Satellites of Disc Galaxies

(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.

preprint2011arXiv

Phase Space Models of the Dwarf Spheroidal Galaxies

This paper introduces new phase-space models of dwarf spheroidal galaxies (dSphs). The stellar component has an isotropic, lowered isothermal (or King) distribution function. A physical basis for the isotropization of stellar velocities is given by tidal stirring, whilst the isothermality of the distribution function guarantees the observed flatness of the velocity dispersion profile in the inner parts. Our models reproduce the data on the half-light radius and line of sight central velocity dispersion of the dSphs. We show that different dark halo profiles -- whether cored or cusped -- lead to very similar mass estimates within one particular radius, namely 1.7 half-light radii. Deviations between mass measures due to different density profiles are substantially smaller than the uncertainties propagated by the observational errors. We produce a mass measure for each of the Milky Way dSphs and find that the two most massive are the most luminous, namely Sagittarius (~ 2.8 x 10^8 solar masses) and Fornax (~ 1.3 x 10^8 solar masses). The least massive of the Milky Way satellites are Willman 1 (~ 4 x 10^5 solar masses) and Segue 1 (~ 6 x 10^5 solar masses).

preprint2011arXiv

The Milky Way stellar halo out to 40 kpc : Squashed, broken but smooth

We introduce a new maximum likelihood method to model the density profile of Blue Horizontal Branch and Blue Straggler stars and apply it to the Sloan Digital Sky Survey Data Release 8 (DR8) photometric catalogue. There are a large number (20,000) of these tracers available over an impressive 14,000 square degrees in both Northern and Southern Galactic hemispheres, and they provide a robust measurement of the shape of the Milky Way stellar halo. After masking out stars in the vicinity of the Virgo Overdensity and the Sagittarius stream, the data are consistent with a smooth, oblate stellar halo with a density that follows a broken power-law. The best fitting model has an inner slope = 2.3 and an outer slope = 4.6, together with a break radius occurring at 27 kpc and a constant halo flattening (that is, ratio of minor axis to major axis) of q = 0.6. Although a broken power-law describes the density fall-off most adequately, it is also well fit by an Einasto profile. There is no strong evidence for variations in flattening with radius, or for triaxiality of the stellar halo.

preprint2010arXiv

Big fish, small fish: Two New Ultra-Faint Satellites of the Milky Way

We report the discovery of two new Milky Way satellites in the neighboring constellations of Pisces and Pegasus identified in data from the Sloan Digital Sky Survey. Pisces II, an ultra-faint dwarf galaxy lies at the distance of ~180 kpc, some 15 degrees away from the recently detected Pisces I. Segue 3, an ultra-faint star cluster lies at the distance of 16 kpc. We use deep follow-up imaging obtained with the 4-m Mayall telescope at Kitt Peak National Observatory to derive their structural parameters. Pisces II has a half-light radius of ~60 pc, while Segue 3 is twenty times smaller at only 3pc.

preprint2010arXiv

Modified Virial Formulae and the Theory of Mass Estimators

We show how to estimate the enclosed mass from the observed motions of an ensemble of test particles. Traditionally, this problem has been attacked through virial or projected mass estimators. Here, we examine and extend these systematically, and show how to construct an optimal estimator for any given assumption as to the potential. The estimators do not explicitly depend on any properties of the density of the test objects, which is desirable as in practice such information is dominated by selection effects. As particular examples, we also develop estimators tailored for the problem of estimating the mass of the Hernquist or NFW dark matter haloes from the projected positions and velocities of stars.

preprint2010arXiv

Re-Assembling the Sagittarius Dwarf Galaxy

What is the mass of the progenitor of the Sagittarius (Sgr) dwarf galaxy? Here, we reassemble the stellar debris using SDSS and 2MASS data to find the total luminosity and likely mass. We find that the luminosity is in the range 9.6-13.2 x10^7 solar luminosities or M_V ~ -15.1 - 15.5, with 70% of the light residing in the debris streams. The progenitor is somewhat fainter than the present-day Small Magellanic Cloud, and comparable in brightness to the M31 dwarf spheroidals NGC 147 and NGC 185. Using cosmologically motivated models, we estimate that the mass of Sgr's dark matter halo prior to tidal disruption was ~10^10 solar masses.

preprint2010arXiv

Rotation of Halo Populations in the Milky Way and M31

We search for signs of rotation in the subsystems of the Milky Way and M31 that are defined by their satellite galaxies, their globular cluster populations, and their BHB stars. A set of simple distribution functions are introduced to describe anisotropic and rotating stellar populations embedded in dark haloes of approximate Navarro-Frenk-White form. The BHB stars in the Milky Way halo exhibit a dichotomy between a prograde rotating, comparatively metal-rich component ([Fe/H] > -2) and a retrograde rotating, comparatively metal-poor ([Fe/H] < -2) component. The prograde metal-rich population may be associated with the accretion of a massive satellite (~ 10^9 solar masses). The metal-poor population may characterise the primordial stellar halo and the net retrograde rotation could then reflect an underestimate in our adopted local standard of rest circular velocity Theta_0. If Theta_0 is ~ 240 km/s then the metal-poor component has no rotation and there is a net prograde rotation signal of ~45 km/s in the metal-rich component. There is reasonable evidence that the Milky Way globular cluster and satellite galaxy systems are rotating with <v_phi> ~50 km/s and <v_phi> ~ 40 km/s respectively. Furthermore, a stronger signal is found for the satellite galaxies when the angular momentum vector of the satellites is inclined with respect to the normal of the disc. The dwarf spheroidal satellites of M31 exhibit prograde rotation relative to the M31 disc with <v_phi> ~ 40 km/s. We postulate that this group of dwarf spheroidals may share a common origin. We also find strong evidence for systemic rotation in the globular clusters of M31 particularly for the most metal-rich.

preprint2010arXiv

The Masses of the Milky Way and Andromeda galaxies

We present a family of robust tracer mass estimators to compute the enclosed mass of galaxy haloes from samples of discrete positional and kinematical data of tracers, such as halo stars, globular clusters and dwarf satellites. The data may be projected positions, distances, line of sight velocities or proper motions. Forms of the estimator tailored for the Milky Way galaxy and for M31 are given. Monte Carlo simulations are used to quantify the uncertainty as a function of sample size. For the Milky Way, the satellite sample consists of 26 galaxies with line-of-sight velocities. We find that the mass of the Milky Way within 300 kpc is ~ 0.9 x 10^12 solar masses assuming velocity isotropy. However, the mass estimate is sensitive to the anisotropy and could plausibly lie between 0.7 - 3.4 x 10^12 solar masses. Incorporating the proper motions of 6 Milky Way satellites into the dataset, we find ~ 1.4 x 10^12 solar masses. The range here if plausible anisotropies are used is still broader, from 1.2 - 2.7 x 10^12 solar masses. For M31, there are 23 satellite galaxies with measured line-of-sight velocities, but only M33 and IC 10 have proper motions. We use the line of sight velocities and distances of the satellite galaxies to estimate the mass of M31 within 300 kpc as ~ 1.4 x 10^12 solar masses assuming isotropy. There is only a modest dependence on anisotropy, with the mass varying between 1.3 -1.6 x 10^12 solar masses. Given the uncertainties, we conclude that the satellite data by themselves yield no reliable insights into which of the two galaxies is actually the more massive.

preprint2009arXiv

A photometric and spectroscopic study of the new dwarf spheroidal galaxy in Hercules

Our aim is to provide as clean and as complete a sample as possible of red giant branch stars that are members of the Hercules dSph galaxy. With this sample we explore the velocity dispersion and the metallicity of the system. Stromgren photometry and multi-fibre spectroscopy are combined to provide information about the evolutionary state of the stars (via the Stromgren c_1 index) and their radial velocities. Based on this information we have selected a clean sample of red giant branch stars, and show that foreground contamination by Milky Way dwarf stars can greatly distort the results. Our final sample consists of 28 red giant branch stars in the Hercules dSph galaxy. Based on these stars we find a mean photometric metallicity of -2.35 dex which is consistent with previous studies. We find evidence for an abundance spread. Using those stars for which we have determined radial velocities we find a systemic velocity of 45.2 km/s with a dispersion of 3.72 km/s, this is lower than values found in the literature. Furthermore we identify the horizontal branch and estimate the mean magnitude of the horizontal branch of the Hercules dSph galaxy to be V_0=21.17, which corresponds to a distance of 147 kpc. We have shown that a proper cleaning of the sample results in a smaller value for the velocity dispersion of the system. This has implications for galaxy properties derived from such velocity dispersions.

preprint2009arXiv

Segue 2: A Prototype of the Population of Satellites of Satellites

We announce the discovery of a new Milky Way satellite Segue 2 found in the data of the Sloan Extension for Galactic Understanding and Exploration (SEGUE). We followed this up with deeper imaging and spectroscopy on the Multiple Mirror Telescope. From this, we derive a luminosity of M_v = -2.5, a half-light radius of 34 pc and a systemic velocity of -40$ km/s. Our MMT data also provides evidence for a stream around Segue 2 at a similar heliocentric velocity, and the SEGUE data show that it is also present in neighboring fields. We resolve the velocity dispersion of Segue 2 as 3.4 km/s and the possible stream as about 7 km/s. This object shows points of comparison with other recent discoveries, Segue 1, Boo II and Coma. We speculate that all four objects may be representatives of a population of satellites of satellites -- survivors of accretion events that destroyed their larger but less dense parents. They are likely to have formed at redshifts z > 10 and are good candidates for fossils of the reionization epoch.

preprint2009arXiv

The Origin of Segue 1

We apply the optimal filter technique to Sloan Digital Sky Survey photometry around Segue 1 and find that the outer parts of the cluster are distorted. There is strong evidence for ~ 1 degree elongations of extra-tidal stars, extending both eastwards and southwestwards of the cluster. The extensions have similar differential Hess diagrams to Segue 1. A Kolmogorov-Smirnov test suggests a high probability that both come from the same parent distribution. The location of Segue 1 is close to crossings of the tidal wraps of the Sagittarius stream. By extracting blue horizontal branch stars from Sloan's spectral database, two kinematic features are isolated and identified with different wraps of the Sagittarius stream. We show that Segue 1 is moving with a velocity that is close to one of the wraps. At this location, we estimate that there are enough Sagittarius stars, indistinguishable from Segue 1 stars, to inflate the velocity dispersion and hence the mass-to-light ratio. All the available evidence is consistent with the interpretation that Segue 1 is a star cluster, originally from the Sagittarius galaxy, and now dissolving in the Milky Way.

preprint2009arXiv

The Tilt of the Halo Velocity Ellipsoid and the Shape of the Milky Way Halo

A sample of roughly 1,800 halo subdwarf stars with radial velocities and proper motions is assembled, using the repeated multi-band Sloan Digital Sky Survey photometric measurements in Stripe 82. Our sample of halo subdwarfs is extracted via a reduced proper motion diagram and distances are obtained using photometric parallaxes, thus giving full phase space information. The tilt of the velocity ellipsoid with respect to the spherical polar coordinate system is computed and found to be consistent with zero for two of the three tilt angles, and very small for the third. We prove that if the inner halo is in a steady-state and the triaxial velocity ellipsoid is everywhere aligned in spherical polar coordinates, then the potential must be spherically symmetric. The detectable, but very mild, misalignment with spherical polars is consistent with the perturbative effects of the Galactic disk on a spherical dark halo. Banana orbits are generated at the 1:1 resonance (in horizontal and vertical frequency) by the disk. They populate Galactic potentials at the typical radii of our subdwarf sample, along with the much more dominant short-axis tubes. However, on geometric grounds alone, the tilt cannot vanish for the banana orbits and this leads to a slight, but detectable, misalignment. We argue that the tilt of the stellar halo velocity ellipsoid therefore provides a hitherto largely neglected but important line of argument that the Milky Way's dark halo, which dominates the potential, must be nearly spherical.

preprint2008arXiv

2MASS J01542930+0053266 : A New Eclipsing M-dwarf Binary System

We report on 2MASS J01542930+0053266, a faint eclipsing system composed of two M dwarfs. The variability of this system was originally discovered during a pilot study of the 2MASS Calibration Point Source Working Database. Additional photometry from the Sloan Digital Sky Survey yields an 8-passband lightcurve, from which we derive an orbital period of 2.6390157 +/- 0.0000016 days. Spectroscopic followup confirms our photometric classification of the system, which is likely composed of M0 and M1 dwarfs. Radial velocity measurements allow us to derive the masses (M_1 = 0.66 +/- 0.03 M_sun; M_2 = 0.62 +/- 0.03 M_sun) and radii (R_1 = 0.64 +/- 0.08 R_sun; R_2 = 0.61 +/- 0.09 R_sun) of the components, which are consistent with empirical mass-radius relationships for low-mass stars in binary systems. We perform Monte Carlo simulations of the lightcurves which allow us to uncover complicated degeneracies between the system parameters. Both stars show evidence of H-alpha emission, something not common in early-type M dwarfs. This suggests that binarity may influence the magnetic activity properties of low-mass stars; activity in the binary may persist long after the dynamos in their isolated counterparts have decayed, yielding a new potential foreground of flaring activity for next generation variability surveys.

preprint2007arXiv

Discovery of an Unusual Dwarf Galaxy in the Outskirts of the Milky Way

In this Letter, we announce the discovery of a new dwarf galaxy, Leo T, in the Local Group. It was found as a stellar overdensity in the Sloan Digital Sky Survey Data Release 5 (SDSS DR5). The color-magnitude diagram of Leo T shows two well-defined features, which we interpret as a red giant branch and a sequence of young, massive stars. As judged from fits to the color-magnitude diagram, it lies at a distance of about 420 kpc and has an intermediate-age stellar population with a metallicity of [Fe/H]= -1.6, together with a young population of blue stars of age of 200 Myr. There is a compact cloud of neutral hydrogen with mass roughly 10^5 solar masses and radial velocity 35 km/s coincident with the object visible in the HIPASS channel maps. Leo T is the smallest, lowest luminosity galaxy found to date with recent star-formation. It appears to be a transition object similar to, but much lower luminosity than, the Phoenix dwarf.

preprint2007arXiv

The discovery of two extremely low luminosity Milky Way globular clusters

We report the discovery of two extremely low luminosity globular clusters in the Milky Way Halo. These objects were detected in the Sloan Digital Sky Survey Data Release 5 and confirmed with deeper imaging at the Calar Alto Observatory. The clusters, Koposov 1 and Koposov 2, are located at $\sim 40-50$ kpc and appear to have old stellar populations and luminosities of only $M_V \sim -1$ mag. Their observed sizes of $\sim 3$ pc are well within the expected tidal limit of $\sim$10 pc at that distance. Together with Palomar 1, AM 4 and Whiting 1, these new clusters are the lowest luminosity globulars orbiting the Milky Way, with Koposov 2 the most extreme. Koposov 1 appears to lie close to distant branch of the Sagittarius stream. The half-mass relaxation times of Koposov 1 and 2 are only $\sim 70$ and $\sim 55$ Myr respectively (2 orders of magnitude shorter than the age of the stellar populations), so it would seem that they have undergone drastic mass segregation. Since they do not appear to be very concentrated, their evaporation timescales may be as low as $\sim 0.1 t_{\rm Hubble}$. These discoveries show that the structural parameter space of globular clusters in the Milky Way halo is not yet fully explored. They also add, through their short remaining survival times, significant direct evidence for a once much larger population of globular clusters.

preprint2006arXiv

A Curious New Milky Way Satellite in Ursa Major

In this Letter, we study a localized stellar overdensity in the constellation of Ursa Major, first identified in Sloan Digital Sky Survey (SDSS) data and subsequently followed up with Subaru imaging. Its color-magnitude diagram (CMD) shows a well-defined sub-giant branch, main sequence and turn-off, from which we estimate a distance of ~30 kpc and a projected size of ~250 x 125 pc. The CMD suggests a composite population with some range in metallicity and/or age. Based on its extent and stellar population, we argue that this is a previously unknown satellite galaxy of the Milky Way, hereby named Ursa Major II (UMa II) after its constellation. Using SDSS data, we find an absolute magnitude of M_V \~ -3.8, which would make it the faintest known satellite galaxy. UMa II's isophotes are irregular and distorted with evidence for multiple concentrations; this suggests that the satellite is in the process of disruption.

preprint2006arXiv

A Faint New Milky Way Satellite in Bootes

In this Letter, we announce the discovery of a new satellite of the Milky Way in the constellation of Bootes at a distance of 60 kpc. It was found in a systematic search for stellar overdensities in the North Galactic Cap using Sloan Digital Sky Survey Data Release 5 (SDSS DR5). The color-magnitude diagram shows a well-defined turn-off, red giant branch, and extended horizontal branch. Its absolute magnitude is -5.8, which makes it one of the faintest galaxies known. The half-light radius is 220 pc. The isodensity contours are elongated and have an irregular shape, suggesting that Boo may be a disrupted dwarf spheroidal galaxy.

preprint2006arXiv

A New Milky Way Dwarf Satellite in Canes Venatici

In this Letter, we announce the discovery of a new dwarf satellite of the Milky Way, located in the constellation Canes Venatici. It was found as a stellar overdensity in the North Galactic Cap using Sloan Digital Sky Survey Data Release 5 (SDSS DR5). The satellite's color-magnitude diagram shows a well-defined red giant branch, as well as a horizontal branch. As judged from the tip of the red giant branch, it lies at a distance of ~220 kpc. Based on the SDSS data, we estimate an absolute magnitude of Mv ~ -7.9, a central surface brightness of mu_0,V ~ 28 mag arcsecond^-2, and a half-light radius of \~ 8.5' (~ 550 pc at the measured distance). The outer regions of Canes Venatici appear extended and distorted. The discovery of such a faint galaxy in proximity to the Milky Way strongly suggests that more such objects remain to be found.

preprint2006arXiv

The Field of Streams: Sagittarius and its Siblings

We use Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5) u,g,r,i,z photometry to study Milky Way halo substructure in the area around the North Galactic Cap. A simple color cut (g-r < 0.4) reveals the tidal stream of the Sagittarius dwarf spheroidal, as well as a number of other stellar structures in the field. Two branches (A and B) of the Sagittarius stream are clearly visible in an RGB-composite image created from 3 magnitude slices, and there is also evidence for a still more distant wrap behind the A branch. A comparison of these data with numerical models suggests that the shape of the Galactic dark halo is close to spherical.

preprint2004arXiv

RIP: The Macho Era (1974-2004)

This article reviews the life and death of a scientific theory.

preprint2003arXiv

The Anomaly in the Candidate Microlensing Event PA-99-N2

The lightcurve of PA-99-N2, one of the recently announced microlensing candidates towards M31, shows small deviations from the standard Paczynski form. We explore a number of possible explanations, including correlations with the seeing, the parallax effect and a binary lens. We find that the observations are consistent with an unresolved RGB or AGB star in M31 being microlensed by a binary lens. We find that the best fit binary lens mass ratio is about one hundredth, which is one of most extreme values found for a binary lens so far. If both the source and lens lie in the M31 disk, then the standard M31 model predicts the probable mass range of the system to be 0.02-3.6 solar masses (95 % confidence limit). In this scenario, the mass of the secondary component is therefore likely to be below the hydrogen-burning limit. On the other hand, if a compact halo object in M31 is lensing a disk or spheroid source, then the total lens mass is likely to lie between 0.09-32 solar masses, which is consistent with the primary being a stellar remnant and the secondary a low mass star or brown dwarf. The optical depth (or alternatively the differential rate) along the line of sight toward the event indicates that a halo lens is more likely than a stellar lens provided that dark compact objects comprise no less than 15 per cent (or 5 per cent) of haloes.

preprint2003arXiv

Theory of pixel lensing towards M31 -- II. The velocity anisotropy and flattening of the MACHO distribution

The POINT-AGAPE collaboration is currently searching for massive compact halo objects (MACHOs) towards the Andromeda galaxy (M31). The survey aims to exploit the high inclination of the M31 disk, which causes an asymmetry in the spatial distribution of M31 MACHOs. Here, we investigate the effects of halo velocity anisotropy and flattening on the asymmetry signal using simple halo models. For a spherically symmetric and isotropic halo, we find that the underlying pixel-lensing rate in far-disk M31 MACHOs is more than 5 times the rate of near-disk events. We find that the asymmetry is increased further by about 30% if the MACHOs occupy radial orbits rather than tangential orbits, but is substantially reduced if the MACHOs lie in a flattened halo. However, even for haloes with a minor-to-major axis ratio q = 0.3, the numbers of M31 MACHOs in the far-side outnumber those in the near-side by a factor of ~2. We show that, if positional information is exploited in addition to number counts, then the number of candidate events required to confirm asymmetry for a range of flattened and anisotropic halo models is achievable, even with significant contamination by variable stars and foreground microlensing events. For pixel-lensing surveys which probe a representative portion of the M31 disk, a sample of around 50 candidates is likely to be sufficient to detect asymmetry within spherical haloes, even if half the sample is contaminated, or to detect asymmetry in haloes as flat as q = 0.3 provided less than a third of the sample comprises contaminants. We also argue that, provided its mass-to-light ratio is less than 100, the recently observed stellar stream around M31 is not problematic for the detection of asymmetry. (Abstract slightly abridged.)

preprint2003arXiv

Tracing the Star Stream Through M31 Using Planetary Nebula Kinematics

We present a possible orbit for the Southern Stream of stars in M31, which connects it to the Northern Spur. Support for this model comes from the dynamics of planetary nebulae (PNe) in the disk of M31: analysis of a new sample of 2611 PNe obtained using the Planetary Nebula Spectrograph reveals ~20 objects whose kinematics are inconsistent with the normal components of the galaxy, but which lie at the right positions and velocities to connect the two photometric features via this orbit. The satellite galaxy M32 is coincident with the stream both in position and velocity, adding weight to the hypothesis that the stream comprises its tidal debris.

preprint2002arXiv

New Mass Estimators for Tracer Populations

We introduce the tracer mass estimator. This is a new and simple way to estimate the enclosed mass from the projected positions and line of sight velocities of a tracer population (such as globular clusters, halo stars and planetary nebulae). Like the projected mass estimator, it works by averaging (projected distance) x (radial velocity squared)/ G over the sample. However, it applies to the commonest case of all, when the tracer population does not follow the overall dark matter density. The method is verified against simulated datasets drawn from Monte Carlo realisations of exact solutions of the collisionless Boltzmann equation and applied to the recent M31 globular cluster data set of Perrett et al. (2002), as well as to M31's satellite galaxies.

preprint2001arXiv

No Need for Dark Matter in Galaxies?

Unhappily, there has been a maelstrom of problems for dark matter theories over the last few years and many serious difficulties still have no resolution in sight. This article reviews the evidence for dark matter in galaxies. The haloes built up by hierarchical merging in dark matter cosmogonies are cusped and dominated by dark matter at the center. Evidence from the microlensing optical depth towards Baade's Window and from dynamical modelling of the Galactic bar already suggests that the Galactic halo is not cusped. Similarly, evidence from the stability of unbarred disk galaxies, as well as the survival of fast bars in barred galaxies, suggests that the this result holds good more generally. Judged on the data from galactic scales alone, the case for dark matter is weak and non-standard theories of gravity provide a better description. Of course, non-standard theories of gravity have their own problems, but not on galactic scales.

preprint1998arXiv

Stability of Power-Law Disks I -- The Fredholm Integral Equation

The power-law disks are a family of infinitesimally thin, axisymmetric stellar disks of infinite extent. The rotation curve can be rising, falling or flat. The self-consistent power-law disks are scale-free, so that all physical quantities vary as a power of radius. They possess simple equilibrium distribution functions depending on the two classical integrals, energy and angular momentum. While maintaining the scale-free equilibrium force law, the power-law disks can be transformed into cut-out disks by preventing stars close to the origin (and sometimes also at large radii) from participating in any disturbance. This paper derives the homogeneous Fredholm integral equation for the in-plane normal modes in the self-consistent and the cut-out power-law disks. This is done by linearising the collisionless Boltzmann equation to find the response density corresponding to any imposed density and potential. The normal modes -- that is, the self-consistent modes of oscillation -- are found by requiring the imposed density to equal the response density. In practice, this scheme is implemented in Fourier space, by decomposing both imposed and response densities in logarithmic spirals. The Fredholm integral equation then relates the transform of the imposed density to the transform of the response density. Numerical strategies to solve the integral equation and to isolate the growth rates and the pattern speeds of the normal modes are discussed.

preprint1996arXiv

Simple Three-Integral Scale-Free Galaxy Models

The Jeans equations give the second moments or stresses required to support a stellar population against the gravity field. A general solution of the Jeans equations for arbitrary axisymmetric scale-free densities in flattened scale-free potentials is given. A two-parameter subset of the solution for the second moments for the self-consistent density of the power-law models, which have exactly spheroidal equipotentials, is examined in detail. In the spherical limit, the potential of these models reduces to that of the singular power-law spheres. We build the physical three-integral distribution functions that correspond to the flattened stellar components. Next, we attack the problem of finding distribution functions associated with the Jeans solutions in flattened scale-free potentials. The third or partial integral introduced by de Zeeuw, Evans and Schwarzschild for Binney's model is generalised to thin and near-thin orbits moving in arbitrary axisymmetric scale-free potentials. The partial integral is a modification of the total angular momentum. For the self-consistent power-law models, we show how this enables the construction of simple three-integral distribution functions. The connexion between these approximate distribution functions and the Jeans solutions is discussed in some detail.

preprint1996arXiv

The Capture and Escape of Stars

The shape of galaxies depends on their orbital populations. These populations change through capture into and escape from resonance. Capture problems fall into distinct cases depending upon the shape of the potential well. To visualise the effective potential well for orbital capture, a diagrammatic approach to the resonant perturbation theory of Born is presented. These diagrams we call equiaction sections. To illustrate their use, we present examples drawn from both galactic and Solar System dynamics. The probability of capture for generic shapes of the potential well is calculated. A number of predictions are made. First, there are barred galaxies that possess two outer rings of gas and stars (type R1R2). We show how to relate changes in the pattern speed and amplitude of the bar to the strength of the two rings. Secondly, under certain conditions, small disturbances can lead to dramatic changes in orbital shape. This can be exploited as a mechanism to pump counter-rotating stars and gas into the nuclei of disk galaxies. Tidal resonant forcing of highly inclined orbits around a central mass causes a substantial increase in the likelihood of collision. Thirdly, the angular momentum of a potential well is changed by the passage of stars across or capture into the well. This can lead to the creation of holes, notches and high velocity tails in the stellar distribution function, whose form we explicitly calculate.

preprint1995arXiv

Jeans and Boltzmann Solutions for Oblate Galaxies with Flat Rotation Curves

A general solution of the Jeans equations for oblate scale-free logarithmic potentials is given. This provides all possible second velocity moments that can hold up a stellar population of flattened scale-free density against the gravity field. A two-parameter subset of second moments for the self-consistent density of Binney's model is examined in detail. These solutions have the desirable property that the observable dispersions in the radial and proper motions can be given explicitly. In the spherical limit, the potential of these models reduces to that of the singular isothermal sphere. The Jeans solutions for scale-free densities of arbitrary flattening that can correspond to physical three-integral distribution functions are identified. The problem of finding distribution functions associated with the Jeans solutions in flattened scale-free logarithmic potentials is then investigated for Binney's model. An approximate solution of the collisionless Boltzmann equation is found which provides a third (partial) integral of good accuracy for thin and near-thin tube orbits. It is a modification of the total angular momentum. This enables the construction of many simple three-integral distribution functions. The kinematic properties of these approximate DFs are shown to agree with a subset of the Jeans solutions --- which are thereby confirmed as good approximations to physical solutions.

preprint1995arXiv

The Nature of the Galactic Dark Matter

Disk-halo models of the Galaxy and LMC are constructed and used to analyse the microlensing data-set. Deflectors in the LMC bar, disk and halo provide an optical depth to microlensing of $\sim 2 \times 10^{-7}$. Deflectors in the Galactic disk and halo contribute $\sim 5 \times 10^{-7}$. The extent, flattening and velocity anisotropy of the dark objects in the Galactic halo are unknown. So, it is crucial to analyse the microlensing data-set with families of models that span the viable ranges of these structural parameters. Also uncertain is the contribution of the Galactic disk to the local circular speed, which affects the normalisation and size of the Galactic halo. Despite all the unknowns, a robust conclusion is that the Galactic and LMC haloes cannot be primarily built from objects in the mass range $10^{-7} \, \msun$ - $0.1 \, \msun$. By contrast, calculations of the baryon mass fraction of the Galactic and LMC haloes are very sensitive to details of the adopted models. This parameter is not constrained by the existing data-set. In particular, it is still possible for the halo to be entirely baryonic and composed of high mass compact objects, like $10^6\, \msun$ clusters or black holes.

preprint1994arXiv

Microlensing by the Milky Way Halo

The measurements of the possible gravitational microlensing events are analysed with a simple yet accurate disc--halo model of the Milky Way Galaxy. This comprises a luminous exponential disc embedded in a flattened dark matter halo with density varying like ${\rm distance}^{-1.8}$. Including a disc has the important effect of lowering the implied masses of the dark matter objects. For the possible detection reported by Alcock et al (1993), the inferred mass of the lens lies in the range $\sim 0.01 - 0.15 M_\odot$. The candidate events of Aubourg et al (1993) have slightly larger implied masses of $\sim 0.025 - 0.35\ M_\odot$ and $\sim 0.03 - 0.45\ M_\odot$ respectively. These are consistent with the deflecters being either brown dwarfs or low mass stars. If there is no disc dark matter and the halo is completely composed of baryonic dark objects of typical mass $ \sim 0.08\ M_\odot$, then the monitoring of $1.8 \times 10^6$ stars in the Large Magellanic Cloud will provide at least $8-9$ detections a year, in the limit of $100 \%$ efficiency.

N. W. Evans

What is connected

Connect this record

See the researcher in context

Building this map preview

80 published item(s)

A 6-d View of Stellar Shells

Donald Lynden-Bell: A Biographical Memoir

The Early Merger that Made the Galaxy's Stellar Halo

Age Demographics of the Milky Way Disk and Bulge

Models of Bars with Exponential Density Profiles

A Dipole on the Sky: Predictions for Hypervelocity Stars from the Large Magellanic Cloud

A Tail of Two Populations: Chemo-dynamics of the Sagittarius Stream and Implications for its Original Mass

Chemodynamic subpopulations of the Carina dwarf galaxy

Discy dwarf disruption and the shape of the Galactic halo

Is the Dark Halo of the Milky Way Prolate?

Is there substructure around M87?

Simple J-Factors and D-Factors for Indirect Dark Matter Detection

The Gaia-ESO survey: Metal-rich bananas in the bulge

Dipping Our Toes in the Water: First Models of GD-1 as a Stream

Haloes light and dark: dynamical models of the stellar halo and constraints on the mass of the Galaxy

Made-to-Measure Dark Matter Haloes, Elliptical Galaxies and Dwarf Galaxies in Action Coordinates

On the Alignment of the Stress Tensor in Galaxies

Resonant Clumping and Substructure in Galactic Discs

The Doubloon Models of Dark Haloes

What can the alignments of the velocity moments tell us about the nature of the potential?

"Skinny Milky Way, Please", says Sagittarius

A Very Simple Cusped Halo Model

ATLAS lifts the Cup: Discovery of a New Milky Way satellite in Crater

Dynamical Models of Elliptical Galaxies -- I. Simple Methods

Dynamical Models of Elliptical Galaxies -- II. M87 and its Globular Clusters

Extremely Flat Haloes and the Shape of the Galaxy

Hamiltonians of Spherically Symmetric, Scale-Free Galaxies in Action-Angle Coordinates

Microlensing, Brown Dwarfs and GAIA

On Asymmetric Distributions of Satellite Galaxies

The remnant of a merger between two dwarf galaxies in Andromeda II

The dark matter halo concentration and stellar initial mass function of a CASSOWARY group

Triaxial Cosmological Haloes and the Disc of Satellites

A Troublesome Past: Chemodynamics of the Fornax dwarf spheroidal

A Virial Core in the Sculptor Dwarf Spheroidal Galaxy

Broken and unbroken: the Milky Way and M31 stellar haloes

Broken degeneracies: the rotation curve and velocity anisotropy of the Milky Way halo

Elliptical galaxy masses out to five effective radii: the realm of dark matter

Lensing and Dynamics in Two Simple Steps

Line Profiles from Discrete Kinematic Data

Slicing and dicing the Milky Way disk in SDSS

The cold veil of the Milky Way stellar halo

The core size of the Fornax dwarf Spheroidal

The Sagittarius Streams in the Southern Galactic Hemisphere

The stretching of Hercules

Virial Sequences for Thick Discs and Haloes: Flattening and Global Anisotropy

Accurate Masses for Navarro-Frenk-White Dark Matter Haloes

Chandrasekhar and Modern Stellar Dynamics

Dark Matter Cores and Cusps: The Case of Multiple Stellar Populations in Dwarf Spheroidals

Mass estimators in the Gaia era

Mismatch and Misalignment: Dark Haloes and Satellites of Disc Galaxies

Phase Space Models of the Dwarf Spheroidal Galaxies

The Milky Way stellar halo out to 40 kpc : Squashed, broken but smooth

Big fish, small fish: Two New Ultra-Faint Satellites of the Milky Way

Modified Virial Formulae and the Theory of Mass Estimators

Re-Assembling the Sagittarius Dwarf Galaxy

Rotation of Halo Populations in the Milky Way and M31

The Masses of the Milky Way and Andromeda galaxies

A photometric and spectroscopic study of the new dwarf spheroidal galaxy in Hercules

Segue 2: A Prototype of the Population of Satellites of Satellites

The Origin of Segue 1

The Tilt of the Halo Velocity Ellipsoid and the Shape of the Milky Way Halo

2MASS J01542930+0053266 : A New Eclipsing M-dwarf Binary System

Discovery of an Unusual Dwarf Galaxy in the Outskirts of the Milky Way

The discovery of two extremely low luminosity Milky Way globular clusters

A Curious New Milky Way Satellite in Ursa Major

A Faint New Milky Way Satellite in Bootes

A New Milky Way Dwarf Satellite in Canes Venatici

The Field of Streams: Sagittarius and its Siblings

RIP: The Macho Era (1974-2004)

The Anomaly in the Candidate Microlensing Event PA-99-N2

Theory of pixel lensing towards M31 -- II. The velocity anisotropy and flattening of the MACHO distribution

Tracing the Star Stream Through M31 Using Planetary Nebula Kinematics

New Mass Estimators for Tracer Populations

No Need for Dark Matter in Galaxies?