BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Justin I. Read

Justin I. Read appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.SRastro-ph.HEastro-ph.IMhep-ph

Catalog footprint

What is connected

32works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

32 published item(s)

preprint2025arXiv

The radial acceleration relation at the EDGE of galaxy formation: testing its universality in low-mass dwarf galaxies

A tight correlation between the baryonic and observed acceleration of galaxies has been reported over a wide range of mass ($10^8 < M_{\rm bar}/{\rm M}_\odot < 10^{11}$) - the Radial Acceleration Relation (RAR). This has been interpreted as evidence that dark matter is actually a manifestation of some modified weak-field gravity theory. In this paper, we study the radially resolved RAR of 12 nearby dwarf galaxies, with baryonic masses in the range $10^4 < M_{\rm bar}/{\rm M}_\odot < 10^{7.5}$, using a combination of literature data and data from the MUSE-Faint survey. We use stellar line-of-sight velocities and the Jeans modelling code GravSphere to infer the mass distributions of these galaxies, allowing us to compute the RAR. We compare the results with the EDGE simulations of isolated dwarf galaxies with similar stellar masses in a $Λ$CDM cosmology. We find that most of the observed dwarf galaxies lie systematically above the low-mass extrapolation of the RAR. Each galaxy traces a locus in the RAR space that can have a multi-valued observed acceleration for a given baryonic acceleration, while there is significant scatter from galaxy to galaxy. Our results indicate that the RAR does not apply to low-mass dwarf galaxies and that the inferred baryonic acceleration of these dwarfs does not contain enough information, on its own, to derive the observed acceleration. The simulated EDGE dwarfs behave similarly to the real data, lying systematically above the extrapolated RAR. We show that, in the context of modified weak-field gravity theories, these results cannot be explained by differential tidal forces from the Milky Way, nor by the galaxies being far from dynamical equilibrium, since none of the galaxies in our sample seems to experience strong tides. As such, our results provide further evidence for the need for invisible dark matter in the smallest dwarf galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Constraining Ultra Light Dark Matter with the Galactic Nuclear Star Cluster

We use the Milky Way's nuclear star cluster (NSC) to test the existence of a dark matter 'soliton core', as predicted in ultra-light dark matter (ULDM) models. Since the soliton core size is proportional to mDM^{-1}, while the core density grows as mDM^{2}, the NSC (dominant stellar component within about 3 pc) is sensitive to a specific window in the dark matter particle mass, mDM. We apply a spherical isotropic Jeans model to fit the NSC line-of-sight velocity dispersion data, assuming priors on the precisely measured Milky Way's supermassive black hole (SMBH) mass and the well-measured NSC density profile. We find that the current observational data reject the existence of a soliton core for a single ULDM particle with mass in the range 10^{-20.4} < mDM < 10^{-18.5} eV, assuming that the soliton core structure is not affected by the Milky Way's SMBH. We test our methodology on mock data, confirming that we are sensitive to the same range in ULDM mass as for the real data. Dynamical modelling of a larger region of the Galactic centre, including the nuclear stellar disc, promises tighter constraints over a broader range of mDM. We will consider this in future work.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

EDGE: the puzzling ellipticity of Eridanus II's star cluster and its implications for dark matter at the heart of an ultra-faint dwarf

The Eridanus II (EriII) 'ultra-faint' dwarf has a large ($15\,\text{pc}$) and low mass ($4.3\times10^3\,\text{M}_\odot$) star cluster (SC) offset from its centre by $23\pm3\,\text{pc}$ in projection. Its size and offset are naturally explained if EriII has a central dark matter core, but such a core may be challenging to explain in a $Λ$CDM cosmology. In this paper, we revisit the survival and evolution of EriII's SC, focussing for the first time on its puzzlingly large ellipticity ($0.31^{+0.05}_{-0.06}$). We perform a suite of 960 direct $N$-body simulations of SCs, orbiting within a range of spherical background potentials fit to ultra-faint dwarf (UFD) galaxy simulations. We find only two scenarios that come close to explaining EriII's SC. In the first, EriII has a low density dark matter core (of size $\sim70\,\text{pc}$ and density $\lesssim2\times10^8\,\text{M}_{\odot}\,\text{kpc}^{-3}$). In this model, the high ellipticity of EriII's SC is set at birth, with the lack of tidal forces in the core allowing its ellipticity to remain frozen in for long times. In the second, EriII's SC orbits in a partial core, with its high ellipticity owing to its imminent tidal destruction. However, this latter model struggles to reproduce the large size of EriII's SC, and it predicts substantial tidal tails around EriII's SC that should have already been seen in the data. This leads us to favour the cored model. We discuss potential caveats to these findings, and the implications of the cored model for galaxy formation and the nature of dark matter.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

EDGE: The sensitivity of ultra-faint dwarfs to a metallicity-dependent initial mass function

Motivated by the observed bottom-light initial mass function (IMF) in faint dwarfs, we study how a metallicity-dependent IMF affects the feedback budget and observables of an ultra-faint dwarf galaxy. We model the evolution of a low-mass ($\approx 8 \, \times \, 10^{8} \, \rm M_{\odot}$) dark matter halo with cosmological, zoomed hydrodynamical simulations capable of resolving individual supernovae explosions, which we complement with an empirically motivated subgrid prescription for systematic IMF variations. In this framework, at the low gas metallicities typical of faint dwarfs, the IMF of newborn stellar populations becomes top-heavy, increasing the efficiency of supernova and photoionization feedback in regulating star formation. This results in a 100-fold reduction of the final stellar mass of the dwarf compared to a canonical IMF, at fixed dynamical mass. The increase in the feedback budget is nonetheless met by increased metal production from more numerous massive stars, leading to nearly constant iron content at $z=0$. A metallicity-dependent IMF therefore provides a mechanism to produce low-mass ($\rm M_{\star}\sim 10^3 \rm M_{\odot}$), yet enriched ($\rm [Fe/H]\approx -2$) field dwarf galaxies, thus opening a self-consistent avenue to populate the plateau in $\rm [Fe/H]$ at the faintest end of the mass--metallicity relation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

EDGE: What shapes the relationship between HI and stellar observables in faint dwarf galaxies?

We show how the interplay between feedback and mass-growth histories introduces scatter in the relationship between stellar and neutral gas properties of field faint dwarf galaxies ($M_{\star} \lessapprox 10^{6} M_{\odot}$). Across a suite of cosmological, high-resolution zoomed simulations, we find that dwarf galaxies of stellar masses $10^5 \leq M_{\star} \leq 10^{6} M_{\odot}$ are bimodal in their cold gas content, being either HI-rich or HI-deficient. This bimodality is generated through the coupling between (i) the modulation of HI contents by the background of ultraviolet radiation (UVB) at late times and (ii) the significant scatter in the stellar-mass-halo-mass relationship induced by reionization. Furthermore, our HI-rich dwarfs exhibit disturbed and time-variable neutral gas distributions primarily due to stellar feedback. Over the last four billion years, we observe order-of-magnitude changes around the median $M_{HI}$, factor-of-a-few variations in HI spatial extents, and spatial offsets between HI and stellar components regularly exceeding the galaxies' optical sizes. Time variability introduces further scatter in the $M_{\star}-M_{HI}$ relation and affects a galaxy's detectability in HI at any given time. These effects will need to be accounted for when interpreting observations of the population of faint, HI-bearing dwarfs by the combination of optical and radio wide, deep surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Observational constraints on stellar feedback in dwarf galaxies

Feedback to the interstellar medium (ISM) from ionising radiation, stellar winds and supernovae is central to regulating star formation in galaxies. Due to their low mass ($M_{*} < 10^{9}$\,M$_\odot$), dwarf galaxies are particularly susceptible to such processes, making them ideal sites to study the detailed physics of feedback. In this perspective, we summarise the latest observational evidences for feedback from star forming regions and how this drives the formation of 'superbubbles' and galaxy-wide winds. We discuss the important role of external ionising radiation -- 'reionisation' -- for the smallest galaxies. And, we discuss the observational evidences that this feedback directly impacts galaxy properties such as their star formation histories, metal content, colours, sizes, morphologies and even their inner dark matter densities. We conclude with a look to the future, summarising the key questions that remain unanswered and listing some of the outstanding challenges for galaxy formation theories.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Formation of the largest galactic cores through binary scouring and gravitational wave recoil

Massive elliptical galaxies are typically observed to have central cores in their projected radial light profiles. Such cores have long been thought to form through `binary scouring' as supermassive black holes (SMBHs), brought in through mergers, form a hard binary and eject stars from the galactic centre. However, the most massive cores, like the ~3kpc core in A2261-BCG, remain challenging to explain in this way. In this paper, we run a suite of dry galaxy merger simulations to explore three different scenarios for central core formation in massive elliptical galaxies: `binary scouring', `tidal deposition' and `gravitational wave (GW) induced recoil'. Using the Griffin code, we self-consistently model the stars, dark matter and SMBHs in our merging galaxies, following the SMBH dynamics through to the formation of a hard binary. We find that we can only explain the large surface brightness core of A2261-BCG with a combination of a major merger that produces a small ~1kpc core through binary scouring, followed by the subsequent GW recoil of its SMBH that acts to grow the core size. We show that this same model can also explain the bright `knots' observed in the core region of A2261-BCG. Key predictions of this scenario are an offset SMBH surrounded by a compact cluster of bound stars and a non-divergent central density profile. We show that the bright `knots' observed in the core region of A2261-BCG are best explained as stalled perturbers resulting from minor mergers, though the brightest may also represent ejected SMBHs surrounded by a stellar cloak of bound stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dark matter constraints from dwarf galaxies with data-driven J-factors

We present an updated analysis of the gamma-ray flux from the directions of classical dwarf spheroidal galaxies, deriving new constraints on WIMP dark matter (DM) annihilation using a decade of Fermi-LAT data. Among the major novelties, we infer the dwarfs' J-factors by including new observations without imposing any a priori parametric profile for the DM distribution. While statistically compatible with results obtained from more conventional parameterisations, this procedure reduces the theoretical bias imposed on the data. Furthermore, we retain the full data-driven shape of the J-factors' empirical probability distributions when setting limits on DM, without imposing log-normality as is typically done. In conjunction with the data-driven J-factors, we use a new method for estimating the probability distribution function of the astrophysical background at the dwarf position, fully profiling over background uncertainties. We show that, for most classical dwarfs, the background systematic uncertainty dominates over the uncertainty on their J-factors. Raw distributions of J- and D-factors (the latter being the analogous of J-factors for decaying DM) are available upon request.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

EDGE: From quiescent to gas-rich to star-forming low-mass dwarf galaxies

We study how star formation is regulated in low-mass field dwarf galaxies ($10^5 \leq M_{\star} \leq 10^6 \, \text{M}_{\odot}$), using cosmological high-resolution ($3 \, \text{pc}$) hydrodynamical simulations. Cosmic reionization quenches star formation in all our simulated dwarfs, but three galaxies with final dynamical masses of $3 \times 10^{9} \,\text{M}_{\odot}$ are subsequently able to replenish their interstellar medium by slowly accreting gas. Two of these galaxies re-ignite and sustain star formation until the present day at an average rate of $10^{-5} \, \text{M}_{\odot} \, \text{yr}^{-1}$, highly reminiscent of observed low-mass star-forming dwarf irregulars such as Leo T. The resumption of star formation is delayed by several billion years due to residual feedback from stellar winds and Type Ia supernovae; even at $z=0$, the third galaxy remains in a temporary equilibrium with a large gas content but without any ongoing star formation. Using the "genetic modification'' approach, we create an alternative mass growth history for this gas-rich quiescent dwarf and show how a small $(0.2\,\mathrm{dex})$ increase in dynamical mass can overcome residual stellar feedback, re-igniting star formation. The interaction between feedback and mass build-up produces a diversity in the stellar ages and gas content of low-mass dwarfs, which will be probed by combining next-generation HI and imaging surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

New Constraints on the Mass of Fermionic Dark Matter from Dwarf Spheroidal Galaxies

Dwarf spheroidal galaxies are excellent systems to probe the nature of fermionic dark matter due to their high observed dark matter phase-space density. In this work, we review, revise and improve upon previous phase-space considerations to obtain lower bounds on the mass of fermionic dark matter particles. The refinement in the results compared to previous works is realised particularly due to a significantly improved Jeans analysis of the galaxies. We discuss two methods to obtain phase-space bounds on the dark matter mass, one model-independent bound based on Pauli's principle, and the other derived from an application of Liouville's theorem. As benchmark examples for the latter case, we derive constraints for thermally decoupled particles and (non-)resonantly produced sterile neutrinos. Using the Pauli principle, we report a model-independent lower bound of $m \geq 0.18\,\mathrm{keV}$ at 68% CL and $m \geq 0.13\,\mathrm{keV}$ at 95% CL. For relativistically decoupled thermal relics, this bound is strengthened to $m \geq 0.59\,\mathrm{keV}$ at 68% CL and $m \geq 0.41\,\mathrm{keV}$ at 95% CL, whilst for non-resonantly produced sterile neutrinos the constraint is $m \geq 2.80\,\mathrm{keV}$ at 68% CL and $m \geq 1.74\,\mathrm{keV}$ at 95% CL. Finally, the phase-space bounds on resonantly produced sterile neutrinos are compared with complementary limits from X-ray, Lyman-$α$ and Big Bang Nucleosynthesis observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Revealing the tidal scars of the Small Magellanic Cloud

Due to their close proximity, the Large and Small Magellanic Clouds (SMC/LMC) provide natural laboratories for understanding how galaxies form and evolve. With the goal of determining the structure and dynamical state of the SMC, we present new spectroscopic data for $\sim$ 3000 SMC red giant branch stars observed using the AAOmega spectrograph at the Anglo-Australian Telescope. We complement our data with further spectroscopic measurements from previous studies that used the same instrumental configuration and proper motions from the \textit{Gaia} Data Release 2 catalogue. Analysing the photometric and stellar kinematic data, we find that the SMC centre of mass presents a conspicuous offset from the velocity centre of its associated $\mbox{H\,{\sc i}}$ gas, suggesting that the SMC gas is likely to be far from dynamical equilibrium. Furthermore, we find evidence that the SMC is currently undergoing tidal disruption by the LMC within 2\,kpc of the centre of the SMC, and possibly all the way in to the very core. This is evidenced by a net outward motion of stars from the SMC centre along the direction towards the LMC and apparent tangential anisotropy at all radii. The latter is expected if the SMC is undergoing significiant tidal stripping, as we demonstrate using a suite of $N$-body simulations of the SMC/LMC system disrupting around the Milky Way. These results suggest that dynamical models for the SMC that assume a steady state will need to be revisited.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

To beta or not to beta: can higher-order Jeans analysis break the mass-anisotropy degeneracy in simulated dwarfs?

We test a non-parametric higher-order Jeans analysis method, GravSphere, on 32 simulated dwarf galaxies comparable to classical Local Group dwarfs like Fornax. The galaxies are selected from the APOSTLE suite of cosmological hydrodynamics simulations with Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, allowing us to investigate cusps and cores in density distributions. We find that, for CDM dwarfs, the recovered enclosed mass profiles have a bias of no more than 10 per cent, with a 50 per cent scatter in the inner regions and a 20 per cent scatter near the half-light radius, consistent with standard mass estimators. The density profiles are also recovered with a bias of no more than 10 per cent and a scatter of 30 per cent in the inner regions. For SIDM dwarfs, the mass and density profiles are recovered within our 95 per cent confidence intervals, but are biased towards cuspy dark matter distributions. This is mainly due to a lack of sufficient constraints from the data. We explore the sources of scatter in the accuracy of the recovered profiles and suggest a $χ^2$ statistic to separate successful models from biased ones. Finally, we show that the uncertainties on the mass profiles obtained with GravSphere are smaller than those for comparable Jeans methods, and that they can be further improved if stronger priors, motivated by cosmological simulations, are placed on the velocity anisotropy. We conclude that GravSphere is a promising Jeans-based approach for modelling dark matter distributions in dwarf galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the origin of the asymmetric dwarf galaxy distribution around Andromeda

The dwarf galaxy distribution surrounding M31 is significantly anisotropic in nature. Of the 30 dwarf galaxies in this distribution, 15 form a disc-like structure and 23 are contained within the hemisphere facing the Milky Way. Using a realistic local potential, we analyse the conditions required to produce and maintain these asymmetries. We find that some dwarf galaxies are required to have highly eccentric orbits in order to preserve the presence of the hemispherical asymmetry with an appropriately large radial dispersion. Under the assumption that the dwarf galaxies originate from a single association or accretion event, we find that the initial size and specific energy of that association must both be relatively large in order to produce the observed hemispherical asymmetry. However if the association was large in physical size, the very high-energy required would enable several dwarf galaxies to escape from the M31 and be captured by the Milky Way. Furthermore, we find that associations that result in this structure have total specific energies concentrated around $E = V_{esc}^{2} - V_{init}^{2} \sim 200^2$ -- $300^2\ \rm{km^2\ s^{-2}}$, implying that the initial velocity and initial position needed to produce the structure are strongly correlated. The overlap of initial conditions required to produce the radial dispersion, angular dispersion, and the planar structure is small and suggests that either they did not originate from a single accretion event, or that these asymmetric structures are short-lived.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A non-parametric method for measuring the local dark matter density

We present a new method for determining the local dark matter density using kinematic data for a population of tracer stars. The Jeans equation in the $z$-direction is integrated to yield an equation that gives the velocity dispersion as a function of the total mass density, tracer density, and the tilt term that describes the coupling of vertical and radial motions. We then fit a dark matter mass profile to tracer density and velocity dispersion data to derive credible regions on the vertical dark matter density profile. Our method avoids numerical differentiation, leading to lower numerical noise, and is able to deal with the tilt term while remaining one dimensional. In this study we present the method and perform initial tests on idealised mock data. We also demonstrate the importance of dealing with the tilt term for tracers that sample $\gtrsim 1$ kpc above the disc plane. If ignored, this results in a systematic underestimation of the dark matter density.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A semi-analytic dynamical friction model for cored galaxies

We present a dynamical friction model based on Chandrasekhar's formula that reproduces the fast inspiral and stalling experienced by satellites orbiting galaxies with a large constant density core. We show that the fast inspiral phase does not owe to resonance. Rather, it owes to the background velocity distribution function for the constant density cores being dissimilar from the usually-assumed Maxwellian distribution. Using the correct background velocity distribution function and the semi-analytic model from Petts et al. (2015), we are able to correctly reproduce the infall rate in both cored and cusped potentials. However, in the case of large cores, our model is no longer able to correctly capture core-stalling. We show that this stalling owes to the tidal radius of the satellite approaching the size of the core. By switching off dynamical friction when rt(r) = r (where rt is the tidal radius at the satellite's position) we arrive at a model which reproduces the N-body results remarkably well. Since the tidal radius can be very large for constant density background distributions, our model recovers the result that stalling can occur for Ms/Menc << 1, where Ms and Menc are the mass of the satellite and the enclosed galaxy mass, respectively. Finally, we include the contribution to dynamical friction that comes from stars moving faster than the satellite. This next-to-leading order effect becomes the dominant driver of inspiral near the core region, prior to stalling.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An efficient positive potential-density pair expansion for modelling galaxies

We present a novel positive potential-density pair expansion for modelling galaxies, based on the Miyamoto-Nagai (MN) disc. By using three sets of such discs, each one of them aligned along each symmetry axis, we are able to reconstruct a broad range of potentials that correspond to density profiles from exponential discs to 3D power law models with varying triaxiality (henceforth simply "twisted" models). We increase the efficiency of our expansion by allowing the scale length parameter of each disc to be negative. We show that, for suitable priors on the scale length and height parameters, these "MNn discs" have just one negative density minimum. This allows us to ensure global positivity by demanding that the total density at the global minimum is positive. We find that at better than 10\% accuracy in our density reconstruction, we can represent a radial and vertical exponential disc over $0.1-10$ scale lengths/heights with 4 MNn discs, an NFW profile over $0.1-10$ scale lengths with 4 MNn discs, and a twisted triaxial NFW profile with 3 MNn discs per symmetry axis. Our expansion is efficient, fully analytic, and well-suited to reproducing the density distribution and gravitational potential of galaxies from discs to ellipsoids.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Collisionless loss-cone refilling: there is no final parsec problem

Coalescing massive black hole binaries, formed during galaxy mergers, are expected to be a primary source of low frequency gravitational waves. Yet in isolated gas-free spherical stellar systems, the hardening of the binary stalls at parsec-scale separations owing to the inefficiency of relaxation-driven loss-cone refilling. Repopulation via collisionless orbit diffusion in triaxial systems is more efficient, but published simulation results are contradictory. While sustained hardening has been reported in simulations of galaxy mergers with $N \sim 10^6$ stars and in early simulations of rotating models, in isolated non-rotating triaxial models the hardening rate continues to fall with increasing N, a signature of spurious two-body relaxation. We present a novel approach for studying loss cone repopulation in galactic nuclei. Since loss cone repopulation in triaxial systems owes to orbit diffusion, it is a purely collisionless phenomenon and can be studied with an approximated force calculation technique, provided the force errors are well behaved and sufficiently small. We achieve this using an accurate fast multipole method and define a proxy for the hardening rate that depends only on stellar angular momenta. We find that the loss cone is efficiently replenished even in very mildly triaxial models (with axis ratios 1 : 0.9 : 0.8). Such triaxiality is unavoidable following galactic mergers and can drive binaries into the gravitational wave regime. We conclude that there is no final parsec problem.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

nIFTy Galaxy Cluster simulations IV: Quantifying the Influence of Baryons on Halo Properties

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A low pre-infall mass for the Carina dwarf galaxy from disequilibrium modelling

Dark matter only simulations of galaxy formation predict many more subhalos around a Milky Way like galaxy than the number of observed satellites. Proposed solutions require the satellites to inhabit dark matter halos with masses between one to ten billion solar masses at the time they fell into the Milky Way. Here we use a modelling approach, independent of cosmological simulations, to obtain a preinfall mass of 360 (+380,-230) million solar masses for one of the Milky Way's satellites: Carina. This determination of a low halo mass for Carina can be accommodated within the standard model only if galaxy formation becomes stochastic in halos below ten billion solar masses. Otherwise Carina, the eighth most luminous Milky Way dwarf, would be expected to inhabit a significantly more massive halo. The implication of this is that a population of "dark dwarfs" should orbit the Milky Way: halos devoid of stars and yet more massive than many of their visible counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Semi-Analytic dynamical friction model that reproduces core stalling

We present a new semi-analytic model for dynamical friction based on Chandrasekhar's formalism. The key novelty is the introduction of physically motivated, radially varying, maximum and minimum impact parameters. With these, our model gives an excellent match to full N-body simulations for isotropic background density distributions, both cuspy and shallow, without any fine-tuning of the model parameters. In particular, we are able to reproduce the dramatic core-stalling effect that occurs in shallow/constant density cores, for the first time. This gives us new physical insight into the core-stalling phenomenon. We show that core stalling occurs in the limit in which the product of the Coulomb logarithm and the local fraction of stars with velocity lower than the infalling body tends to zero. For cuspy backgrounds, this occurs when the infalling mass approaches the enclosed background mass. For cored backgrounds, it occurs at larger distances from the centre, due to a combination of a rapidly increasing minimum impact parameter and a lack of slow moving stars in the core. This demonstrates that the physics of core-stalling is likely the same for both massive infalling objects and low-mass objects moving in shallow density backgrounds. We implement our prescription for dynamical friction in the direct summation code NBODY6 as an analytic correction for stars that remain within the Roche volume of the infalling object. This approach is computationally efficient, since only stars in the inspiralling system need to be evolved with direct summation. Our method can be applied to study a variety of astrophysical systems, including young star clusters orbiting near the Galactic Centre; globular clusters moving within the Galaxy; and dwarf galaxies orbiting within dark matter halos.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiple populations in globular clusters: the distinct kinematic imprints of different formation scenarios

Several scenarios have been proposed to explain the presence of multiple stellar populations in globular clusters. Many of them invoke multiple generations of stars to explain the observed chemical abundance anomalies, but it has also been suggested that self-enrichment could occur via accretion of ejecta from massive stars onto the circumstellar disc of low-mass pre-main sequence stars. These scenarios imply different initial conditions for the kinematics of the various stellar populations. Given some net angular momentum initially, models for which a second generation forms from gas that collects in a cooling flow into the core of the cluster predict an initially larger rotational amplitude for the polluted stars compared to the pristine stars. This is opposite to what is expected from the accretion model, where the polluted stars are the ones crossing the core and are on preferentially radial (low-angular momentum) orbits, such that their rotational amplitude is lower. Here we present the results of a suite of $N$-body simulations with initial conditions chosen to capture the distinct kinematic properties of these pollution scenarios. We show that initial differences in the kinematics of polluted and pristine stars can survive to the present epoch in the outer parts of a large fraction of Galactic globular clusters. The differential rotation of pristine and polluted stars is identified as a unique kinematic signature that could allow us to distinguish between various scenarios, while other kinematic imprints are generally very similar from one scenario to the other.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

nIFTy galaxy cluster simulations I: dark matter & non-radiative models

We have simulated the formation of a galaxy cluster in a $Λ$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Magellanic Inter-Cloud Project (MAGIC) II: Slicing up the Bridge

The origin of the gas in between the Magellanic Clouds (MCs), known as the Magellanic Bridge (MB), has always been the subject of controversy. To shed light into this, we present the results from the MAGIC II project aimed at probing the stellar populations in ten large fields located perpendicular to the main ridge-line of HI in the Inter-Cloud region. We secured these observations of the stellar populations in between the MCs using the WFI camera on the 2.2 m telescope in La Silla. Using colour-magnitude diagrams (CMDs), we trace stellar populations across the Inter-Cloud region. In good agreement with MAGIC I, we find significant intermediate-age stars in the Inter-Cloud region as well as young stars of a similar age to the last pericentre passage in between the MCs (~200 Myr ago). We show here that the young, intermediate-age and old stars have distinct spatial distributions. The young stars correlate well with the HI gas suggesting that they were either recently stripped from the SMC or formed in-situ. The bulk of intermediate-age stars are located mainly in the bridge region where the HI column density is higher, but they are more spread out than the young stars. They have very similar properties to stars located ~2 Kpc from the SMC centre, suggesting that they were tidally stripped from this region. Finally, the old stars extend to some 8 Kpc from the SMC supporting the idea that all galaxies have a large extended metal poor stellar halo.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravitational Lens Recovery with GLASS: Measuring the mass profile and shape of a lens

We use a new non-parametric gravitational modelling tool -- \Glass{} -- to determine what quality of data (strong lensing, stellar kinematics, and/or stellar masses) are required to measure the circularly averaged mass profile of a lens and its shape. \Glass{} uses an under-constrained adaptive grid of mass pixels to model the lens, searching through thousands of models to marginalise over model uncertainties. Our key findings are as follows: (i) for pure lens data, multiple sources with wide redshift separation give the strongest constraints as this breaks the well-known mass-sheet or steepness degeneracy; (ii) a single quad with time delays also performs well, giving a good recovery of both the mass profile and its shape; (iii) stellar masses -- for lenses where the stars dominate the central potential -- can also break the steepness degeneracy, giving a recovery for doubles almost as good as having a quad with time delay data, or multiple source redshifts; (iv) stellar kinematics provide a robust measure of the mass at the half light radius of the stars $r_{1/2}$ that can also break the steepness degeneracy if the Einstein radius $r_E \neq r_{1/2}$; and (v) if $r_E \sim r_{1/2}$, then stellar kinematic data can be used to probe the stellar velocity anisotropy $β$ -- an interesting quantity in its own right. Where information on the mass distribution from lensing and/or other probes becomes redundant, this opens up the possibility of using strong lensing to constrain cosmological models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The hunt for the Milky Way's accreted disc

The Milky Way is expected to host an accreted disc of stars and dark matter. This forms as massive >1:10 mergers are preferentially dragged towards the disc plane by dynamical friction and then tidally shredded. The accreted disc likely contributes only a tiny fraction of the Milky Way's thin and thick stellar disc. However, it is interesting because: (i) its associated `dark disc' has important implications for experiments hoping to detect a dark matter particle in the laboratory; and (ii) the presence or absence of such a disc constrains the merger history of our Galaxy. In this work, we develop a chemo-dynamical template to hunt for the accreted disc. We apply our template to the high-resolution spectroscopic sample from Ruchti et al. (2011), finding at present no evidence for accreted stars. Our results are consistent with a quiescent Milky Way with no >1:10 mergers since the disc formed and a correspondingly light `dark disc'. However, we caution that while our method can robustly identify accreted stars, our incomplete stellar sample makes it more challenging to definitively rule them out. Larger unbiased stellar samples will be required for this.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The MAGellanic Inter-Cloud project (MAGIC) I: Evidence for intermediate-age stellar populations in between the Magellanic Clouds

The origin of the gas in between the Magellanic Clouds (MCs) -- known as the `Magellanic Bridge' (MB) -- is puzzling. Numerical simulations suggest that the MB formed from tidally stripped gas and stars in a recent interaction between the MCs. However, the apparent lack of stripped intermediate- or old-age stars associated with the MB is at odds with this picture. In this paper, we present the first results from the MAGellanic Inter-Cloud program (MAGIC) aimed at probing the stellar populations in the inter-cloud region. We present observations of the stellar populations in two large fields located in between the Large and Small Magellanic Clouds (LMC/SMC), secured using the WFI camera on the 2.2 m telescope in La Silla. Using a synthetic color-magnitude diagram (CMD) technique, we present the first quantitative evidence for the presence of intermediate-age and old stars in the inter-Cloud region. The intermediate-age stars -- that make up ~28% of all stars in the region -- are not present in fields at a similar distance from the SMC in a direction pointing away from the LMC. This provides potential evidence that these intermediate-age stars could have been tidally stripped from the SMC. However, spectroscopic studies will be needed to confirm or rule out the tidal origin for the inter-Cloud gas and stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A new determination of the local dark matter density from the kinematics of K dwarfs

We apply a new method to determine the local disc matter and dark halo matter density to kinematic and position data for \sim2000 K dwarf stars taken from the literature. Our method assumes only that the disc is locally in dynamical equilibrium, and that the 'tilt' term in the Jeans equations is small up to \sim1 kpc above the plane. We present a new calculation of the photometric distances to the K dwarf stars, and use a Monte Carlo Markov Chain to marginalise over uncertainties in both the baryonic mass distribution, and the velocity and distance errors for each individual star. We perform a series of tests to demonstrate that our results are insensitive to plausible systematic errors in our distance calibration, and we show that our method recovers the correct answer from a dynamically evolved N-body simulation of the Milky Way. We find a local dark matter density of ρdm = 0.025+0.014-0.013 M\odotpc^{-3} (0.95+0.53-0.49 GeV cm^{-3}) at 90% confidence assuming no correction for the non-flatness of the local rotation curve, and ρdm = 0.022+0.015-0.013 M\odotpc^-3 (0.85+0.57-0.50 GeV cm^{-3}) if the correction is included. Our 90% lower bound on ρdm is larger than the canonical value typically assumed in the literature, and is at mild tension with extrapolations from the rotation curve that assume a spherical halo. Our result can be explained by a larger normalisation for the local Milky Way rotation curve, an oblate dark matter halo, a local disc of dark matter, or some combination of these.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The mass distribution of the Fornax dSph: constraints from its globular cluster distribution

Uniquely among the dwarf spheroidal (dSph) satellite galaxies of the Milky Way, Fornax hosts globular clusters. It remains a puzzle as to why dynamical friction has not yet dragged any of Fornax's five globular clusters to the centre, and also why there is no evidence that any similar star cluster has been in the past (for Fornax or any other dSph). We set up a suite of 2800 N-body simulations that sample the full range of globular-cluster orbits and mass models consistent with all existing observational constraints for Fornax. In agreement with previous work, we find that if Fornax has a large dark-matter core then its globular clusters remain close to their currently observed locations for long times. Furthermore, we find previously unreported behaviour for clusters that start inside the core region. These are pushed out of the core and gain orbital energy, a process we call 'dynamical buoyancy'. Thus a cored mass distribution in Fornax will naturally lead to a shell-like globular cluster distribution near the core radius, independent of the initial conditions. By contrast, CDM-type cusped mass distributions lead to the rapid infall of at least one cluster within Δt = 1-2Gyr, except when picking unlikely initial conditions for the cluster orbits (\sim 2% probability), and almost all clusters within Δt = 10Gyr. Alternatively, if Fornax has only a weakly cusped mass distribution, dynamical friction is much reduced. While over Δt = 10Gyr this still leads to the infall of 1-4 clusters from their present orbits, the infall of any cluster within Δt = 1-2Gyr is much less likely (with probability 0-70%, depending on Δt and the strength of the cusp). Such a solution to the timing problem requires that in the past the globular clusters were somewhat further from Fornax than today; they most likely did not form within Fornax, but were accreted.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Haloes gone MAD: The Halo-Finder Comparison Project

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Limits on the local dark matter density

We revisit systematics in determining the local dark matter density (rho_dm) from the vertical motion of stars in the Solar Neighbourhood. Using a simulation of a Milky Way-like galaxy, we determine the data-quality required to detect the dark matter density at its expected local value. We introduce a new method for recovering rho_dm that uses moments of the Jeans equations, combined with a Monte Carlo Markov Chain technique to marginalise over the unknown parameters. Given sufficiently good data, we show that our method can recover the correct local dark matter density even in the face of disc inhomogeneities, non-isothermal tracers and a non-separable distribution function. We illustrate the power of our technique by applying it to Hipparcos data [Holmberg & Flynn 2000,2004]. We first make the assumption that the A and F star tracer populations are isothermal. This recovers rho_dm=0.003^{+0.009}_{-0.007}Msun/pc^3 (with 90 per cent confidence), consistent with previous determinations. However, the vertical dispersion profile of these tracers is poorly known. If we assume instead a non-isothermal profile similar to the blue disc stars from SDSS DR-7 [Abazajian et al. 2009] measured by Bond et al. (2009), we obtain a fit with a very similar chi^2 value, but with rho_dm=0.033^{+0.008}_{-0.009}Msun/pc^3 (with 90 per cent confidence). This highlights that it is vital to measure the vertical dispersion profile of the tracers to recover an unbiased estimate of the local dark matter density.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Measuring the local dark matter density

We examine systematic problems in determining the local matter density from the vertical motion of stars, i.e. the 'Oort limit'. Using collisionless simulations and a Monte Carlo Markov Chain technique, we determine the data quality required to detect local dark matter at its expected density. We find that systematic errors are more important than observational errors and apply our technique to Hipparcos data to reassign realistic error bars to the local dark matter density.

0 citations0 reviewsNo code linkedOpen access