-
Models of Bars I: Flattish Profiles for Early-Type Galaxies
Authors:
A. A. Williams,
N. W. Evans
Abstract:
We introduce a simple family of barred galaxy models with flat rotation curves. They are built by convolving the axisymmetric logarithmic model with a needle density. The density contours in the bar region are highly triaxial and elongated, but become spheroidal in the outer parts. The mass distribution differs markedly from the elliptical shape assumed in other analytical models, like Ferrers or…
▽ More
We introduce a simple family of barred galaxy models with flat rotation curves. They are built by convolving the axisymmetric logarithmic model with a needle density. The density contours in the bar region are highly triaxial and elongated, but become spheroidal in the outer parts. The mass distribution differs markedly from the elliptical shape assumed in other analytical models, like Ferrers or Freeman bars. The surface density profile along the bar major axis is flattish, as befits models for bars in early-type galaxies (SB0, SBa). The two-dimensional orbital structure of the models is analysed with surfaces of section and characteristic diagrams and it reveals qualitatively new features. For some pattern speeds, additional Lagrange points can exist along the major axis, and give rise to off-centered, trapped orbits. When the bar is weak, the orbital structure is very simple, comprising just prograde, aligned $x_1$ orbits and retrograde anti-aligned $x_4$ orbits. As the bar strength increases, the $x_1$ family becomes unstable and vanishes, with propeller orbits dominating the characteristic diagram.
△ Less
Submitted 12 May, 2017;
originally announced May 2017.
-
The Gaia-ESO Survey: Dynamical models of flattened, rotating globular clusters
Authors:
Sarah M. R. Jeffreson,
Jason L. Sanders,
N. W. Evans,
Angus A. Williams,
G. F. Gilmore,
Amelia Bayo,
Angela Bragaglia,
Ettore Flaccomio,
Richard Jackson,
Robert D. Jeffries,
Paul Jofre,
Carmela Lardo,
Lorenzo Morbidelli,
Elena Pancino,
Simone Zaggia
Abstract:
We present a family of self-consistent axisymmetric rotating globular cluster models which are fitted to spectroscopic data for NGC 362, NGC 1851, NGC 2808, NGC 4372, NGC 5927 and NGC 6752 to provide constraints on their physical and kinematic properties, including their rotation signals. They are constructed by flattening Modified Plummer profiles, which have the same asymptotic behaviour as clas…
▽ More
We present a family of self-consistent axisymmetric rotating globular cluster models which are fitted to spectroscopic data for NGC 362, NGC 1851, NGC 2808, NGC 4372, NGC 5927 and NGC 6752 to provide constraints on their physical and kinematic properties, including their rotation signals. They are constructed by flattening Modified Plummer profiles, which have the same asymptotic behaviour as classical Plummer models, but can provide better fits to young clusters due to a slower turnover in the density profile. The models are in dynamical equilibrium as they depend solely on the action variables. We employ a fully Bayesian scheme to investigate the uncertainty in our model parameters (including mass-to-light ratios and inclination angles) and evaluate the Bayesian evidence ratio for rotating to non-rotating models. We find convincing levels of rotation only in NGC 2808. In the other clusters, there is just a hint of rotation (in particular, NGC 4372 and NGC 5927), as the data quality does not allow us to draw strong conclusions. Where rotation is present, we find that it is confined to the central regions, within radii of $R \leq 2 r_h$. As part of this work, we have developed a novel q-Gaussian basis expansion of the line-of-sight velocity distributions, from which general models can be constructed via interpolation on the basis coefficients.
△ Less
Submitted 25 April, 2017;
originally announced April 2017.
-
On the Run: mapping the escape speed across the Galaxy with SDSS
Authors:
Angus A. Williams,
Vasily Belokurov,
Andrew R. Casey,
N. Wyn Evans
Abstract:
We measure the variation of the escape speed of the Galaxy across a range of $\sim$ 40 kpc in Galactocentric radius. The local escape speed is found to be $521^{+46}_{-30}\,\mathrm{km\,s^{-1}}$, in good agreement with other studies. We find that this has already fallen to $379^{+34}_{-28}\,\mathrm{km\,s^{-1}}$ at a radius of 50 kpc. Through measuring the escape speed and its variation, we obtain c…
▽ More
We measure the variation of the escape speed of the Galaxy across a range of $\sim$ 40 kpc in Galactocentric radius. The local escape speed is found to be $521^{+46}_{-30}\,\mathrm{km\,s^{-1}}$, in good agreement with other studies. We find that this has already fallen to $379^{+34}_{-28}\,\mathrm{km\,s^{-1}}$ at a radius of 50 kpc. Through measuring the escape speed and its variation, we obtain constraints on the Galactic potential as a whole. In particular, the gradient in the escape speed with radius suggests that the total mass contained within 50 kpc is $29^{+7}_{-5}\times10^{10}\mathrm{M}_\odot$, implying a relatively light dark halo for the Milky Way. Our method represents a novel way of estimating the mass of the Galaxy, and has very different systematics to more commonly used models of tracers, which are more sensitive to the central parts of the halo velocity distributions. Using our inference on the escape speed, we then investigate the orbits of high-speed Milky Way dwarf galaxies. For each dwarf we consider, we predict small pericenter radii and large orbital eccentricities. This naturally explains the large observed ellipticities of two of the dwarfs, which are likely to have been heavily disrupted as they passed through pericenter.
△ Less
Submitted 5 January, 2017;
originally announced January 2017.
-
The Gaia-ESO survey: Metal-rich bananas in the bulge
Authors:
Angus A. Williams,
N. W. Evans,
Matthew Molloy,
Georges Kordopatis,
M. C. Smith,
J. Shen,
G. Gilmore,
S. Randich,
T. Bensby,
P. Francois,
S. E Koposov,
A. Recio-Blanco,
A. Bayo,
G. Carraro,
A. Casey,
T. Costado,
E. Franciosini,
A. Hourihane,
P. de Laverny,
J. Lewis,
K. Lind,
L. Magrini,
L. Monaco,
L. Morbidelli,
G. G Sacco
, et al. (3 additional authors not shown)
Abstract:
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-ric…
▽ More
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.
△ Less
Submitted 31 May, 2016;
originally announced May 2016.
-
Is the Dark Halo of the Milky Way Prolate?
Authors:
A. Bowden,
N. W. Evans,
A. A. Williams
Abstract:
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 prolatene…
▽ More
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$.
△ Less
Submitted 23 April, 2016;
originally announced April 2016.
-
On the Alignment of the Stress Tensor in Galaxies
Authors:
N. W. Evans,
J. L. Sanders,
A. A. Williams,
J. An,
D. Lynden-Bell,
W. Dehnen
Abstract:
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 invar…
▽ More
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.
△ Less
Submitted 24 November, 2015; v1 submitted 29 September, 2015;
originally announced September 2015.
-
Haloes light and dark: dynamical models of the stellar halo and constraints on the mass of the Galaxy
Authors:
A. A. Williams,
N. W. Evans
Abstract:
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 poin…
▽ More
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.
△ Less
Submitted 21 August, 2015; v1 submitted 11 August, 2015;
originally announced August 2015.
-
The Doubloon Models of Dark Haloes
Authors:
N. W. Evans,
J. An,
A. Bowden,
A. A. Williams
Abstract:
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 (cusp…
▽ More
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.
△ Less
Submitted 12 March, 2015;
originally announced March 2015.
-
Made-to-Measure Dark Matter Haloes, Elliptical Galaxies and Dwarf Galaxies in Action Coordinates
Authors:
A. A. Williams,
N. W. Evans
Abstract:
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 inne…
▽ More
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.
△ Less
Submitted 14 January, 2015; v1 submitted 15 December, 2014;
originally announced December 2014.
-
A Very Simple Cusped Halo Model
Authors:
N. W. Evans,
A. A. Williams
Abstract:
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 readil…
▽ More
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.
△ Less
Submitted 14 June, 2014;
originally announced June 2014.
-
Hamiltonians of Spherically Symmetric, Scale-Free Galaxies in Action-Angle Coordinates
Authors:
A. A. Williams,
N. W. Evans,
A. Bowden
Abstract:
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 isotherm…
▽ More
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.
△ Less
Submitted 5 May, 2014;
originally announced May 2014.