Researcher profile

Garry W. Angus

Garry W. Angus contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 published item(s)

preprint2015arXiv

Dynamical measurement of the stellar surface density of face-on galaxies

The DiskMass survey recently provided measurements of the vertical velocity dispersions of disk stars in a sample of nearly face-on galaxies. By setting the disk scale-heights to be equal to those of edge-on galaxies with similar scale-lengths, it was found that these disks must be sub-maximal, with surprisingly low K-band mass-to-light ratios of the order of $M_\star/L_K \simeq 0.3 M_\odot/L_\odot$. This study made use of a simple relation between the disk surface density and the measured velocity dispersion and scale height of the disk, neglecting the shape of the rotation curve and the dark matter contribution to the vertical force, which can be especially important in the case of sub-maximal disks. Here, we point out that these simplifying assumptions led to an overestimation of the stellar mass-to-light ratios. Relaxing these assumptions, we compute even lower values than previously reported for the mass-to-light ratios, with a median $M_\star/L_K \simeq 0.18 M_\odot/L_\odot$, where 14 galaxies have $M_\star/L_K < 0.11$. Invoking prolate dark matter halos made only a small difference to the derived $M_\star/L_K$, although extreme prolate halos ($q>1.5$ for the axis ratios of the potential) might help. The cross-terms in the Jeans equation are also generally negligible. These deduced K-band stellar mass-to-light ratios are even more difficult to reconcile with stellar population synthesis models than the previously reported ones.

preprint2015arXiv

Mass models of disk galaxies from the DiskMass Survey in MOND

This article explores the agreement between the predictions of Modified Newtonian Dynamics (MOND) and the rotation curves and stellar velocity dispersion profiles measured by the DiskMass Survey. A bulge-disk decomposition was made for each of the thirty published galaxies, and a MOND Poisson solver was used to simultaneously compute, from the baryonic mass distributions, model rotation curves and vertical velocity dispersion profiles, which were compared to the measured values. The two main free parameters, the stellar disk&#39;s mass-to-light ratio ($M/L$) and its exponential scale-height ($h_z$), were estimated by Markov Chain Monte Carlo modelling. The average best-fit K-band stellar mass-to-light ratio was $M/L \simeq 0.55 \pm 0.15$. However, to match the DiskMass Survey data, the vertical scale-heights would have to be in the range $h_z=200$ to $400$ pc which is a factor of two lower than those derived from observations of edge-on galaxies with a similar scale-length. The reason is that modified gravity versions of MOND characteristically require a larger $M/L$ to fit the rotation curve in the absence of dark matter and therefore predict a stronger vertical gravitational field than Newtonian models. It was found that changing the MOND acceleration parameter, the shape of the velocity dispersion ellipsoid, the adopted vertical distribution of stars, as well as the galaxy inclination, within any realistic range, all had little impact on these results.

preprint2014arXiv

Modified Baryonic Dynamics: two-component cosmological simulations with light sterile neutrinos

In this article we continue to test cosmological models centred on Modified Newtonian Dynamics (MOND) with light sterile neutrinos, which could in principle be a way to solve the fine-tuning problems of the standard model on galaxy scales while preserving successful predictions on larger scales. Due to previous failures of the simple MOND cosmological model, here we test a speculative model where the modified gravitational field is produced only by the baryons and the sterile neutrinos produce a purely Newtonian field (hence Modified Baryonic Dynamics). We use two component cosmological simulations to separate the baryonic N-body particles from the sterile neutrino ones. The premise is to attenuate the over-production of massive galaxy cluster halos which were prevalent in the original MOND plus light sterile neutrinos scenario. Theoretical issues with such a formulation notwithstanding, the Modified Baryonic Dynamics model fails to produce the correct amplitude for the galaxy cluster mass function for any reasonable value of the primordial power spectrum normalisation.

preprint2014arXiv

N-body simulations of the Carina dSph in MOND

The classical dwarf spheroidals (dSphs) provide a critical test for Modified Newtonian Dynamics (MOND) because they are observable satellite galactic systems with low internal accelerations and low, but periodically varying, external acceleration. This varying external gravitational field is not commonly found acting on systems with low internal acceleration. Using Jeans modelling, Carina in particular has been demonstrated to require a V-band mass-to-light ratio greater than 5, which is the nominal upper limit for an ancient stellar population. We run MOND N-body simulations of a Carina-like dSph orbiting the Milky Way to test if dSphs in MOND are stable to tidal forces over the Hubble time and if those same tidal forces artificially inflate their velocity dispersions and therefore their apparent mass-to-light ratio. We run many simulations with various initial total masses for Carina, and Galactocentric orbits (consistent with proper motions), and compare the simulation line of sight velocity dispersions (losVDs) with the observed losVDs of Walker et al. (2007). We find that the dSphs are stable, but that the tidal forces are not conducive to artificially inflating the losVDs. Furthermore, the range of mass-to-light ratios that best reproduces the observed line of sight velocity dispersions of Carina is 5.3 to 5.7 and circular orbits are preferred to plunging orbits. Therefore, some tension still exists between the required mass-to-light ratio for the Carina dSph in MOND and those expected from stellar population synthesis models. It remains to be seen whether a careful treatment of the binary population or triaxiality might reduce this tension.

preprint2013arXiv

Cosmological simulations in MOND: the cluster scale halo mass function with light sterile neutrinos

We use our Modified Newtonian Dynamics (MOND) cosmological particle-mesh N-body code to investigate the feasibility of structure formation in a framework involving MOND and light sterile neutrinos in the mass range 11 - 300 eV: always assuming that Ω_{ν_s}=0.225 for H_o=72 \kms Mpc^{-1}. We run a suite of simulations with variants on the expansion history, cosmological variation of the MOND acceleration constant, different normalisations of the power spectrum of the initial perturbations and interpolating functions. Using various box sizes, but typically with ones of length 256 Mpc/h, we compare our simulated halo mass functions with observed cluster mass functions and show that (i) the sterile neutrino mass must be larger than 30 eV to account for the low mass (M_{200}<10^{14.6} solar masses) clusters of galaxies in MOND and (ii) regardless of sterile neutrino mass or any of the variations we mentioned above, it is not possible to form the correct number of high mass (M_{200}>10^{15.1} solar masses) clusters of galaxies: there is always a considerable over production. This means that the ansatz of considering the weak-field limit of MOND together with a component of light sterile neutrinos to form structure from z ~ 200 fails. If MOND is the correct description of weak-field gravitational dynamics, it could mean that subtle effects of the additional fields in covariant theories of MOND render the ansatz inaccurate, or that the gravity generated by light sterile neutrinos (or by similar hot dark matter particles) is different from that generated by the baryons.

preprint2013arXiv

The dynamics of the bulge dominated galaxy NGC 7814 in MOND

The bulge dominated galaxy NGC 7814 provides one of the strongest dynamical tests possible for Modified Newtonian Dynamics (MOND). Spitzer 3.6 micron photometry fixes the bulge parameterisation and strongly constrains the properties of the sub-dominant stellar disk. Furthermore, the distance is known to better than 5 percent, virtually eliminating it as a free parameter. The rotation curve is easily measured, since the H I (and stellar) disks are edge on, and both the receding and approaching sides agree very well. In this paper we explore the agreement between the model and observed rotation curves in MOND given that the only two free parameters available are the mass-to-light ratios of the bulge and disk. We use a grid based MOND Poisson solver that accurately solves for the MOND gravity and produces our model rotation curves from a given mass distribution. The input to the Poisson solver is a 3D distribution of N particles which is generated from modelling the observed distribution of stars and gas in the galaxy. By ensuring a superior fit to the radial surface brightness profile than previous works, by virtue of a double Sersic fit to the bulge, we were able to produce excellent fits to the rotation curve with typical values for both mass-to-light ratios. We conclude that the model rotation curve of a mass distribution in MOND is extremely sensitive to the bulge-disk decomposition and even slight deviation from the observed mass distribution can produce large differences in the model rotation curve.

preprint2012arXiv

A QUMOND galactic N-body code I: Poisson solver and rotation curve fitting

Here we present a new particle-mesh galactic N-body code that uses the full multigrid algorithm for solving the modified Poisson equation of the Quasi Linear formulation of Modified Newtonian Dynamics (QUMOND). A novel approach for handling the boundary conditions using a refinement strategy is implemented and the accuracy of the code is compared with analytical solutions of Kuzmin disks. We then employ the code to compute the predicted rotation curves for a sample of five spiral galaxies from the THINGS sample. We generated static N-body realisations of the galaxies according to their stellar and gaseous surface densities and allowed their distances, mass-to-light ratios (M/L) and both the stellar and gas scale-heights to vary in order to estimate the best fit parameters. We found that NGC 3621, NGC 3521 and DDO 154 are well fit by MOND using expected values of the distance and M/L. NGC 2403 required a moderately larger $M/L$ than expected and NGC 2903 required a substantially larger value. The surprising result was that the scale-height of the dominant baryonic component was well constrained by the rotation curves: the gas scale-height for DDO 154 and the stellar scale-height for the others. In fact, if the suggested stellar scale-height (one-fifth the stellar scale-length) was used in the case of NGC 3621 and NGC 3521 it would not be possible to produce a good fit to the inner rotation curve. For each of the four stellar dominated galaxies, we calculated the vertical velocity dispersions which we found to be, on the whole, quite typical compared with observed stellar vertical velocity dispersions of face on spirals. We conclude that modelling the gas scale-heights of the gas rich dwarf spiral galaxies will be vital in order to make precise conclusions about MOND.

preprint2012arXiv

The Acceleration Scale, Modified Newtonian Dynamics, and Sterile Neutrinos

General Relativity is able to describe the dynamics of galaxies and larger cosmic structures only if most of the matter in the Universe is dark, namely it does not emit any electromagnetic radiation. Intriguingly, on the scale of galaxies, there is strong observational evidence that the presence of dark matter appears to be necessary only when the gravitational field inferred from the distribution of the luminous matter falls below an acceleration of the order of 10^(-10) m/s^2. In the standard model, which combines Newtonian gravity with dark matter, the origin of this acceleration scale is challenging and remains unsolved. On the contrary, the full set of observations can be neatly described, and were partly predicted, by a modification of Newtonian dynamics, dubbed MOND, that does not resort to the existence of dark matter. On the scale of galaxy clusters and beyond, however, MOND is not as successful as on the scale of galaxies, and the existence of some dark matter appears unavoidable. A model combining MOND with hot dark matter made of sterile neutrinos seems to be able to describe most of the astrophysical phenomenology, from the power spectrum of the cosmic microwave background anisotropies to the dynamics of dwarf galaxies. Whether there exists a yet unknown covariant theory that contains General Relativity and Newtonian gravity in the weak field limit, and MOND as the ultra-weak field limit is still an open question.

preprint2011arXiv

The abundance of galaxy clusters in MOND: Cosmological simulations with massive neutrinos

We present a new Particle-Mesh cosmological N-body code for accurately solving the modified Poisson equation of the Quasi Linear formulation of MOND. We generate initial conditions for the Angus (2009) cosmological model, which is identical to LCDM except that the cold dark matter is switched for a single species of thermal sterile neutrinos. We set the initial conditions at z=250 for a (512 Mpc/h)^3 box with 256^3 particles and we evolve them down to z=0. We clearly demonstrate the necessity of MOND for developing the large scale structure in a hot dark matter cosmology and contradict the naive expectation that MOND cannot form galaxy clusters. We find that the correct order of magnitude of X-ray clusters (with T_X > 4.5 keV) can be formed, but that we overpredict the number of very rich clusters and seriously underpredict the number of lower mass clusters. The latter is a shortcoming of the resolution of our simulations, whereas we suggest that the over production of very rich clusters might be prevented by incorporating a MOND acceleration constant that varies with redshift and an expansion history that cannot be described by the usual Friedmann models. We present evidence that suggests the density profiles of our simulated clusters are compatible with those of observed X-ray clusters in MOND. It remains to be seen if the low mass end of the cluster mass function can be reproduced and if the high densities of dark matter in the central 20 kpc of groups and clusters of galaxies, measured in the MOND framework, can be achieved. As a last test, we computed the relative velocity between pairs of halos within 10 Mpc and find that pairs with velocities larger than 3000 km/s like the bullet cluster, can form without difficulty.

preprint2011arXiv

Using dwarf satellite proper motions to determine their origin

The highly organised distribution of satellite galaxies surrounding the Milky Way is a serious challenge to the concordance cosmological model. Perhaps the only remaining solution, in this framework, is that the dwarf satellite galaxies fall into the Milky Way&#39;s potential along one or two filaments, which may or may not plausibly reproduce the observed distribution. Here we test this scenario by making use of the proper motions of the Fornax, Sculptor, Ursa Minor and Carina dwarf spheroidals, and trace their orbits back through several variations of the Milky Way&#39;s potential and account for dynamical friction. The key parameters are the proper motions and total masses of the dwarf galaxies. Using a simple model we find no tenable set of parameters that can allow Fornax to be consistent with filamentary infall, mainly because the 1 sigma error on its proper motion is relatively small. The other three must walk a tightrope between requiring a small pericentre (less than 20 kpc) to lose enough orbital energy to dynamical friction and avoiding being tidally disrupted. We then employed a more realistic model with host halo mass accretion, and found that the four dwarf galaxies must have fallen in at least 5 Gyrs ago. This time interval is longer than organised distribution is expected to last before being erased by the randomisation of the satellite orbits.

preprint2010arXiv

A lower limit on the dark particle mass from dSphs

We use dwarf spheroidal galaxies as a tool to attempt to put precise lower limits on the mass of the dark matter particle, assuming it is a sterile neutrino. We begin by making cored dark halo fits to the line of sight velocity dispersions as a function of projected radius (taken from Walker et al. 2007) for six of the Milky Way&#39;s dwarf spheroidal galaxies. We test Osipkov-Merritt velocity anisotropy profiles, but find that no benefit is gained over constant velocity anisotropy. In contrast to previous attempts, we do not assume any relation between the stellar velocity dispersions and the dark matter ones, but instead we solve directly for the sterile neutrino velocity dispersion at all radii by using the equation of state for a partially degenerate neutrino gas (which ensures hydrostatic equilibrium of the sterile neutrino halo). This yields a 1:1 relation between the sterile neutrino density and velocity dispersion, and therefore gives us an accurate estimate of the Tremaine-Gunn limit at all radii. By varying the sterile neutrino particle mass, we locate the minimum mass for all six dwarf spheroidals such that the Tremaine-Gunn limit is not exceeded at any radius (in particular at the centre). We find sizeable differences between the ranges of feasible sterile neutrino particle mass for each dwarf, but interestingly there exists a small range 270-280eV which is consistent with all dSphs at the 1-$σ$ level.

preprint2010arXiv

Implications for dwarf spheroidal mass content from interloper removal

Using the caustic method, we identify the member stars of five dwarf spheroidal (dSph) galaxies of the Milky Way, the smallest dark matter (DM) dominated systems in the Universe. After our interloper rejection, we compute line-of-sight velocity dispersion profiles that are substantially smoother than previous results. Moreover, two dSphs have line-of-sight velocity dispersions 20% smaller than previous calculations suggested. Our Jeans modelling confirms that the DM content interior to 300pc is roughly constant with satellite luminosity. Finally, if we assume that MOND provides the true law of gravity, our identification of interlopers implies that four dSphs have mass-to-light ratios in agreement with stellar population synthesis models, whereas Carina still has a mass-to-light ratio a factor of two too large and remains a problem for MOND.

preprint2009arXiv

Current data on the globular cluster Palomar 14 are not inconsistent with MOND

Certain types of globular clusters have the very important property that the predictions for their kinematics in the Newtonian and modified Newtonian dynamics (MOND) contexts are divergent. Here, we caution the recent claim that the stellar kinematics data (using 17 stars) of the globular cluster Palomar 14 are inconsistent with MOND. We compare the observations to the theoretical predictions using a Kolmogorov-Smirnov test, which is appropriate for small samples. We find that, with the currently available data, the MOND prediction for the velocity distribution can only be excluded with a very low confidence level, clearly insufficient to claim that MOND is falsified.

preprint2009arXiv

Equilibrium configurations of 11eV sterile neutrinos in MONDian galaxy clusters

(Abridged) Modified Newtonian Dynamics (MOND) can fit a broad range of galaxy kinematic data, but struggles with clusters of galaxies. MONDian clusters need dark matter, and here we test the 11 eV sterile neutrino - used to fit the first three acoustic peaks of the cosmic microwave background - by investigating their equilibrium distributions in 30 groups and clusters over a wide range of temperatures. Knowing both the sterile neutrino density and velocity dispersion tells us the Tremaine-Gunn phase-space limit at all radii. We find that all 30 systems serendipitously reach the Tremaine-Gunn limit by the centre, which means a portion of the dynamical mass must always be covered by the brightest cluster galaxy. Interestingly, the typical fitted K-band mass-to-light ratio is unity and at most 1.2, which is very consistent - although leaving no margin for error - with stellar population synthesis models. Amidst the sample there are several special cases including the Coma cluster (for which dark matter was first proposed), NGC 720 (where geometrical evidence for dark matter was found) and the bullet cluster (where dark matter - of some kind - in clusters was directly proven to exist). We demonstrate that 11 eV sterile neutrinos are unlikely to influence spiral galaxy rotation curves, as they don&#39;t influence even some very massive early-types (NGC 4125 and NGC 6482). Finally, we conclude that it is intriguing that the minimum mass of sterile neutrino particle that can match the cosmic microwave background is identical to the minimum mass found here to be consistent with equilibrium configurations of MONDian clusters of galaxies.

preprint2009arXiv

Resolving the timing problem of the globular clusters orbiting the Fornax dwarf galaxy

We re-investigate the old problem of the survival of the five globular clusters orbiting the Fornax dwarf galaxy in both standard and modified Newtonian dynamics. For the first time in the history of the topic, we use accurate mass models for the Fornax dwarf, obtained through Jeans modelling of the recently published line of sight velocity dispersion data, and we are also not resigned to circular orbits for the globular clusters. Previously conceived problems stem from fixing the starting distances of the globulars to be less than half the tidal radius. We relax this constraint since there is absolutely no evidence for it and show that the dark matter paradigm, with either cusped or cored dark matter profiles, has no trouble sustaining the orbits of the two least massive globular clusters for a Hubble time almost regardless of their initial distance from Fornax. The three most massive globulars can remain in orbit as long as their starting distances are marginally outside the tidal radius. The outlook for modified Newtonian dynamics is also not nearly as bleak as previously reported. Although dynamical friction inside the tidal radius is far stronger in MOND, outside dynamical friction is negligible due to the absence of stars. This allows highly radial orbits to survive, but more importantly circular orbits at distances more than 85% of Fornax&#39;s tidal radius to survive indefinitely. The probability of the globular clusters being on circular orbits at this distance compared with their current projected distances is discussed and shown to be plausible. Finally, if we ignore the presence of the most massive globular (giving it a large line of sight distance) we demonstrate that the remaining four globulars can survive within the tidal radius for the Hubble time with perfectly sensible orbits.

preprint2006arXiv

Tidal Destruction of The First Dark Microhalos

We point out that the usual self-similarity in cold dark matter models is broken by encounters with individual normal galactic stars on sub-pc scale. Tidal heating and stripping must have redefined the density and velocity structures of the population of the Earth-mass dark matter halos, which are likely to have been the first bound structures to form in the Universe. The disruption rate depends strongly on {\it galaxy types} and the orbital distribution of the microhalos; in the Milky Way, stochastic radial orbits are destroyed first by stars in the triaxial bulge, microhalos on non-planar retrograde orbits with large pericenters and/or apocenters survive the longest. The final microhalo distribution in the {\it solar neighborhood} is better described as a superposition of filamentry microstreams rather than as a set of discrete spherical clumps in an otherwise homogeneous medium. We discuss its important consequences to our detections of microhalos by direct recoil signal and indirect annihilation signal.