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

Chris Power

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.HEastro-ph.IMphysics.data-an

Catalog footprint

What is connected

56works
5topics
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

56 published item(s)

preprint2025arXiv

The origin of the spectral versus dynamical age discrepancy in radio galaxies

We investigate the effect of turbulent magnetic fields on the observed spectral properties of synchrotron radio emission in large-scale radio galaxy lobes. We use three-dimensional relativistic magnetohydrodynamic simulations of fast, high-powered jets to study the structure of the lobe magnetic fields and how this structure affects the radio spectrum of the lobes. It has previously been argued that lobe ages inferred from radio spectra underestimate the true ages of radio galaxies due to re-acceleration of electrons in the lobe, mixing of electron populations, or the presence of turbulent magnetic fields in the lobes. We find that the spectral ages with and without accounting for the lobe magnetic field structure are consistent with each other, suggesting that mixing of radiating populations of different ages is the primary cause of the underestimation of radio lobe ages. By accounting for the structure of lobe magnetic fields, we find greater spectral steepening in the equatorial regions of the lobe. We demonstrate that the assumptions of the continuous injection, Jaffe-Perola, and Tribble models for radio lobe spectra do not hold in our simulations, and we show that young particles with high magnetic field strengths are the dominant contributors to the overall radio lobe spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

\textsc{The Three Hundred} project: The \textsc{Gizmo-Simba} run

We introduce \textsc{Gizmo-Simba}, a new suite of galaxy cluster simulations within \textsc{The Three Hundred} project. \textsc{The Three Hundred} consists of zoom re-simulations of 324 clusters with $M_{200}\gtrsim 10^{14.8}M_\odot$ drawn from the MultiDark-Planck $N$-body simulation, run using several hydrodynamic and semi-analytic codes. The \textsc{Gizmo-Simba} suite adds a state-of-the-art galaxy formation model based on the highly successful {\sc Simba} simulation, mildly re-calibrated to match $z=0$ cluster stellar properties. Comparing to \textsc{The Three Hundred} zooms run with \textsc{Gadget-X}, we find intrinsic differences in the evolution of the stellar and gas mass fractions, BCG ages, and galaxy colour-magnitude diagrams, with \textsc{Gizmo-Simba} generally providing a good match to available data at $z \approx 0$. \textsc{Gizmo-Simba}'s unique black hole growth and feedback model yields agreement with the observed BH scaling relations at the intermediate-mass range and predicts a slightly different slope at high masses where few observations currently lie. \textsc{Gizmo-Simba} provides a new and novel platform to elucidate the co-evolution of galaxies, gas, and black holes within the densest cosmic environments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A stochastic model to reproduce the star-formation history of individual galaxies in hydrodynamic simulations

The star formation history (SFH) of galaxies is critical for understanding galaxy evolution. Hydrodynamical simulations enable us to precisely reconstruct the SFH of galaxies and establish a link to the underlying physical processes. In this work, we present a model to describe individual galaxies' SFHs from three simulations: TheThreeHundred, Illustris-1 and TNG100-1. This model divides the galaxy SFH into two distinct components: the "main sequence" and the "variation". The "main sequence" part is generated by tracing the history of the $SFR-M_*$ main sequence of galaxies across time. The "variation" part consists of the scatter around the main sequence, which is reproduced by fractional Brownian motions. We find that: 1) The evolution of the main sequence varies between simulations; 2) fractional Brownian motions can reproduce many features of SFHs, however, discrepancies still exist; 3) The variations and mass-loss rate are crucial for reconstructing the SFHs of the simulations. This model provides a fair description of the SFHs in simulations. On the other hand, by correlating the fractional Brownian motion model to simulation data, we provide a 'standard' against which to compare simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An orbital perspective on the starvation, stripping, and quenching of satellite galaxies in the EAGLE simulations

Using the EAGLE suite of simulations, we demonstrate that both cold gas stripping {\it and} starvation of gas inflow play an important role in quenching satellite galaxies across a range of stellar and halo masses, $M_{\star}$ and $M_{200}$. By quantifying the balance between gas inflows, outflows, and star formation rates, we show that even at $z=2$, only $\approx30\%$ of satellite galaxies are able to maintain equilibrium or grow their reservoir of cool gas - compared to $\approx50\%$ of central galaxies at this redshift. We find that the number of orbits completed by a satellite is a very good predictor of its quenching, even more so than the time since infall. On average, we show that intermediate-mass satellites with $M_{\star}$ between $10^{9}{\rm M}_{\odot}-10^{10}{\rm M}_{\odot}$ will be quenched at first pericenter in massive group environments, $M_{200}>10^{13.5}{\rm M}_{\odot}$; and will be quenched at second pericenter in less massive group environments, $M_{200}<10^{13.5}{\rm M}_{\odot}$. On average, more massive satellites ($M_{\star}>10^{10}{\rm M}_{\odot}$) experience longer depletion time-scales, being quenched between first and second pericenters in massive groups; while in smaller group environments, just $\approx30\%$ will be quenched even after two orbits. Our results suggest that while starvation alone may be enough to slowly quench satellite galaxies, direct gas stripping, particularly at pericenters, is required to produce the short quenching time-scales exhibited in the simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spin transfer from dark matter to gas during halo formation

In the protogalactic density field, diffuse gas and collision-less cold dark matter (DM) are often assumed sufficiently mixed that both components experience identical tidal torques. However, haloes in cosmological simulations consistently end up with a higher specific angular momentum (sAM) in gas, even in simulations without radiative cooling and galaxy formation physics. We refine this result by analysing the spin distributions of gas and DM in $\sim$50,000 well-resolved haloes in a non-radiative cosmological simulation from the SURFS suite. The sAM of the halo gas on average ends up $\sim$40\% above that of the DM. This can be pinned down to an excess AM in the inner halo ($<$50\% virial radius), paralleled by a more coherent rotation pattern in the gas. We uncover the leading driver for this AM difference through a series of control simulations of a collapsing ellipsoidal top-hat, where gas and DM are initially well mixed. These runs reveal that the pressurised inner gas shells collapse more slowly, causing the DM ellipsoid to spin ahead of the gas ellipsoid. The arising torque generally transfers AM from the DM to the gas. The amount of AM transferred via this mode depends on the initial spin, the initial axes ratios and the collapse factor. These quantities can be combined in a single dimensionless parameter, which robustly predicts the AM transfer of the ellipsoidal collapse. This simplistic model can quantitatively explain the average AM excess of the gas found in the more complex non-radiative cosmological simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Three Hundred project: galaxy cluster mergers and their impact on the stellar component of brightest cluster galaxies

Using the data set of The Three Hundred project, i.e. a suite of 324 hydrodynamical resimulations of cluster-sized haloes, we study galaxy cluster mergers and their effect on colour and luminosity changes of their brightest cluster galaxies (BCG). We track the main progenitor of each halo at z=0 and search for merger situations based on its mass accretion history, defining mergers as very rapid increases in the halo mass. Based upon the evolution of the dynamical state of the cluster we define a pre- and post-merger phase. We create a list of all these events and statistically study their mass ratio and timescales, with the former verifying that all instances are in fact major mergers. By comparing to a control sample of clusters without mergers, we study the effect mergers have on the stellar component of the BCG. Analysing the mass, age and metallicity of the BCG stellar particles, we find that the stellar content of BCGs grows significantly during mergers and, even though the main growth mechanism is the accretion of older stars, there is even a burst in star formation induced by the merger. In our simulations, BCGs in mergers form in median around 70 per cent more stars than those normally growing, although this depends on the radius considered for defining the BCG. Regarding observable properties, we see an increase in SDSS-u luminosity of 20 per cent during mergers, accompanied by a slightly slower increase of the galaxy g-r colour as compared to the control sample.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Pre-processing, group accretion and the orbital trajectories of associated subhaloes

We use a high-resolution cosmological dark matter-only simulation to study the orbital trajectories of haloes and subhaloes in the environs of isolated hosts. We carefully tally all apsis points and use them to distinguish haloes that are infalling for the first time from those that occupy more evolved orbits. We find that roughly 21 per cent of subhaloes within a host's virial radius are currently on first infall, and have not yet reached their first orbital pericentre; roughly 44 per cent are still approaching their first apocentre after infall. For the range of host masses studied, roughly half of all accreted systems were pre-processed prior to infall, and about 20 per cent were accreted in groups. We confirm that the entire population of accreted subhaloes -- often referred to as "associated" subhaloes -- extend far beyond the virial radii of their hosts, with roughly half currently residing at distances that exceed $\approx 1.2\times r_{200}$. Many of these backsplash haloes have gained orbital energy since infall, and occupy extreme orbits that carry them well past their initial turnaround radii. Such extreme orbits are created during the initial accretion and dissolution of loosely bound groups, but also through penetrating encounters between subhaloes on subsequent orbits. The same processes may also give rise to unexpectedly abrupt losses of orbital energy. These effects combine, giving rise to a large variation in the ratio of sequent apocentres for accreted systems. We find that, within 2 virial radii from host centres, the concentrations of first-infall halos are remarkably similar those of isolated field halos, whereas backsplash haloes, as well as systems that were pre-processed, are considerably more concentrated.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The distribution and properties of DLAs at z $\leq$ 2 in the EAGLE simulations

Determining the spatial distribution and intrinsic physical properties of neutral hydrogen on cosmological scales is one of the key goals of next-generation radio surveys. We use the EAGLE galaxy formation simulations to assess the properties of damped Lyman-alpha absorbers (DLAs) that are associated with galaxies and their underlying dark matter haloes between 0 $\leq$ z $\leq$ 2. We find that the covering fraction of DLAs increases at higher redshift; a significant fraction of neutral atomic hydrogen (HI) resides in the outskirts of galaxies with stellar mass greater than or equal to 10$^{10}$ M$_\odot$; and the covering fraction of DLAs in the circumgalactic medium (CGM) is enhanced relative to that of the interstellar medium (ISM) with increasing halo mass. Moreover, we find that the mean density of the HI in galaxies increases with increasing stellar mass, while the DLAs in high- and low-halo-mass systems have higher column densities than those in galaxies with intermediate halo masses (~ 10$^{12}$ M$_\odot$ at z = 0). These high-impact CGM DLAs in high-stellar-mass systems tend to be metal-poor, likely tracing smooth accretion. Overall, our results point to the CGM playing an important role in DLA studies at high redshift (z $\geq$ 1). However, their properties are impacted both by numerical resolution and the detailed feedback prescriptions employed in cosmological simulations, particularly that of AGN.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Hierarchical Structure of Galactic Haloes: Classification and characterisation with Halo-OPTICS

We build upon Ordering Points To Identify Clustering Structure (OPTICS), a hierarchical clustering algorithm well-known to be a robust data-miner, in order to produce Halo-OPTICS, an algorithm designed for the automatic detection and extraction of all meaningful clusters between any two arbitrary sizes. We then apply Halo-OPTICS to the 3D spatial positions of halo particles within four separate synthetic Milky Way type galaxies, classifying the stellar and dark matter structural hierarchies. Through visualisation of the Halo-OPTICS output, we compare its structure identification to the state-of-the-art galaxy/(sub)halo finder VELOCIraptor, finding excellent agreement even though Halo-OPTICS does not consider kinematic information in this current implementation. We conclude that Halo-OPTICS is a robust hierarchical halo finder, although its determination of lower spatial-density features such as the tails of streams could be improved with the inclusion of extra localised information such as particle kinematics and stellar metallicity into its distance metric.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Impact of Realistic Foreground and Instrument Models on 21cm Epoch of Reionization Experiments

Predictions for the ability of 21-cm interferometric experiments to discriminate Epoch of Reionization (EoR) signal models are typically limited by the simplicity of data models, whereby foreground signals and characteristics of the instrument are often simplified or neglected.To move towards more realistic scenarios, we explore the effects of applying more realistic foreground and instrument models to the 21cm signal, and the ability to estimate astrophysical parameters with these additional complexities. We use a highly-optimized version of \textsc{21cmFAST}, integrated into \textsc{21cmMC}, to generate lightcones of the brightness temperature fluctuation for Bayesian parameter estimation. We include a statistical point-source foreground model and an instrument model based on the Murchison Widefield Array (MWA) scaled in observation time to have an effective sensitivity similar to the future Square Kilometre Array (SKA). We also extend the current likelihood prescription to account for the presence of beam convolution and foregrounds, the 2-Dimensional Power Spectrum (PS), and the correlation of PS modes. We use frequency bands between 150 and 180 MHz to constrain the ionizing efficiency ($ζ$), the minimum virial temperature of halos ($T_{\mathrm{vir}}$), the soft X-ray emissivity per unit Star Formation Rate (SFR) ($L_X/SFR$ ), and the X-ray energy threshold ($E_0$). We find that the inclusion of realistic foregrounds and instrumental components biases the parameter constraints due to unaccounted for cross-power between the EoR signal, foregrounds and thermal noise. This causes estimates of $ζ$ to be biased by up to $5σ$ but the estimates of $T_{vir}$, L$_X$/SFR and E$_0$ remain unaffected and are all within $1σ$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The impact of stellar and AGN feedback on halo-scale baryonic and dark matter accretion in the EAGLE simulations

We use the EAGLE suite of hydrodynamical simulations to analyse accretion rates (and the breakdown of their constituent channels) onto haloes over cosmic time, comparing the behaviour of baryons and dark matter (DM). We also investigate the influence of sub-grid baryon physics on halo-scale inflow, specifically the consequences of modelling radiative cooling, as well as feedback from stars and active galactic nuclei (AGN). We find that variations in halo baryon fractions at fixed mass (particularly their circum-galactic medium gas content) are very well correlated with variations in the baryon fraction of accreting matter, which we show to be heavily suppressed by stellar feedback in low-mass haloes, $M_{\rm halo}\lesssim10^{11.5}M_{\odot}$. Breaking down accretion rates into first infall, recycled, transfer and merger components, we show that baryons are much more likely to be smoothly accreted than to have originated from mergers when compared to DM, finding (averaged across halo mass) a merger contribution of $\approx6\%$ for baryons, and $\approx15\%$ for DM at $z\approx0$. We also show that the breakdown of inflow into different channels is strongly dependent on sub-grid physics, particularly the contribution of recycled accretion (accreting matter that has been previously ejected from progenitor haloes). Our findings highlight the dual role that baryonic feedback plays in regulating the evolution of galaxies and haloes: by (i) directly removing gas from haloes, and (ii) suppressing gas inflow to haloes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The physical drivers of the atomic hydrogen-halo mass relation

We use SHARK, a semi-analytic galaxy formation model, to investigate the physical processes involved in dictating the shape, scatter and evolution of the HI-halo mass relation at $0\leq z \leq 2$. We compare SHARK with HI clustering and spectral stacking of the HI-halo mass relation derived from observations finding excellent agreement with the former and a deficiency of HI in SHARK at $M_{\rm vir}\approx 10^{12-13} M_{\odot}$ in the latter, but otherwise great agreement below and above that mass threshold. In SHARK, we find that the HI mass increases with the halo mass up to a critical mass of $\approx 10^{11.8} M_{\odot}$; between $\sim 10^{11.8}-10^{13}M_{\odot}$, the scatter in the relation increases by 0.7 dex and the HI mass decreases with the halo mass on average; at $M_{\rm vir} \geq 10^{13} M_{\odot}$, the HI content continues to increase with halo mass. We find that the critical halo mass of $\approx 10^{12} M_{\odot}$ is largely set by feedback from Active Galactic Nuclei (AGN), and the exact shape and scatter of the HI-halo mass relation around that mass is extremely sensitive to how AGN feedback is modelled, with other physical processes playing a less significant role. We determine the main secondary parameters responsible for the scatter of the HI-halo mass relation, namely the halo spin parameter at $M_{\rm vir}\leq 10^{11.8} M_{\odot}$, and the fractional contribution from substructure to the total halo mass for $M_{\rm vir}\geq 10^{13} M_{\odot}$. The scatter at $10^{11.8}<M_{\rm vir}<10^{13} M_{\odot}$ is best described by the black-hole-to-stellar mass ratio of the central galaxy, reflecting the AGN feedback relevance. We present a numerical model to populate dark matter-only simulations with HI at $0\leq z \leq 2$ based solely on halo parameters that are measurable in such simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Three Hundred project: shapes and radial alignment of satellite, infalling, and backsplash galaxies

Using 324 numerically modelled galaxy clusters we investigate the radial and galaxy-halo alignment of dark matter subhaloes and satellite galaxies orbiting within and around them. We find that radial alignment depends on distance to the centre of the galaxy cluster but appears independent of the dynamical state of the central host cluster. Furthermore, we cannot find a relation between radial alignment of the halo or galaxy shape with its own mass. We report that backsplash galaxies, i.e. objects that have already passed through the cluster radius but are now located in the outskirts, show a stronger radial alignment than infalling objects. We further find that there exists a population of well radially aligned objects passing very close to the central cluster's centre which were found to be on highly radial orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

xGASS: Robust quantification of asymmetries in global HI spectra and their relationship to environmental processes

We present an analysis of asymmetries in global HI spectra from the extended GALEX Arecibo SDSS Survey (xGASS), a stellar mass-selected and gas fraction-limited survey which is representative of the HI properties of galaxies in the local Universe. We demonstrate that the asymmetry in a HI spectrum is strongly linked to its signal-to-noise meaning that, contrary to what was done in previous works, asymmetry distributions for different samples cannot be compared at face value. We develop a method to account for noise-induced asymmetry and find that the typical galaxy detected by xGASS exhibits higher asymmetry than what can be attributed to noise alone, with 37% of the sample showing asymmetry greater than 10% at an 80% confidence level. We find that asymmetric galaxies contain, on average, 29% less HI mass compared to their symmetric counterparts matched in both stellar mass and signal-to-noise. We also present clear evidence that satellite galaxies, as a population, exhibit more asymmetric HI spectra than centrals and that group central galaxies show a slightly higher rate of HI asymmetries compared to isolated centrals. All these results support a scenario in which environmental processes, in particular those responsible for gas removal, are the dominant driver of asymmetry in xGASS.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

From Stellar Halos to Intracluster Light: the physics of the Intra-Halo Stellar Component in cosmological hydrodynamical simulations

We study the Intra-Halo Stellar Component (IHSC) of Milky Way-mass systems up to galaxy clusters in the Horizon-AGN cosmological hydrodynamical simulation. We identify the IHSC using an improved phase-space galaxy finder algorithm which provides an adaptive, physically motivated and shape-independent definition of this stellar component, that can be applied to halos of arbitrary masses. We explore the IHSC mass fraction-total halo's stellar mass, $f_{M*,IHSC}-M*$, relation and the physical drivers of its scatter. We find that on average the $f_{M*,IHSC}$ increases with $M_{*,tot}$, with the scatter decreasing strongly with mass from 2 dex at $M_{*,tot}\sim10^{11}M_\odot$ to 0.3 dex at group masses. At high masses, $M_{*,tot}>10^{11.5}M_\odot$, $f_{M*,IHSC}$ increases with the number of substructures, and with the mass ratio between the central galaxy and largest satellite, at fixed $M_{*,tot}$. From mid-size groups and systems below $M_{*,tot}<10^{12}M_\odot$, we find that the central galaxy's stellar rotation-to-dispersion velocity ratio, V/σ, displays the strongest (anti)-correlation with $f_{M*,IHSC}$ at fixed $M_{*,tot}$ of all the galaxy and halo properties explored, transitioning from $f_{M*,IHSC}$<0.1% for high V/σ, to $f_{M*,IHSC}\sim5$% for low V/σ galaxies. By studying the $f_{M*,IHSC}$ temporal evolution, we find that, in the former, mergers not always take place, but if they did, they happened early (z>1), while the high $f_{M*,IHSC}$ population displays a much more active merger history. In the case of massive groups and galaxy clusters, $M_{*,tot}>10^{12}M_\odot$, a fraction $f_{M*,IHSC}\sim$10-20% is reached at $z\sim1$ and then they evolve across lines of constant $f_{M*,IHSC}$ modulo some small perturbations. Because of the limited simulation's volume, the latter is only tentative and requires a larger sample of simulated galaxy clusters to confirm.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dark-ages Reionization & Galaxy Formation Simulation I: The dynamical lives of high redshift galaxies

We present the Dark-ages Reionization and Galaxy-formation Observables from Numerical Simulations (DRAGONS) program and Tiamat, the collisionless N-body simulation program upon which DRAGONS is built. The primary trait distinguishing Tiamat from other large simulation programs is its density of outputs at high redshift (100 from z=35 to z=5; roughly one every 10 Myr) enabling the construction of very accurate merger trees at an epoch when galaxy formation is rapid and mergers extremely frequent. We find that the friends-of-friends halo mass function agrees well with the prediction of Watson et al. at high masses, but deviates at low masses, perhaps due to our use of a different halo finder or perhaps indicating a break from "universal" behaviour. We then analyse the dynamical evolution of galaxies during the Epoch of Reionization finding that only a small fraction (~20%) of galactic halos are relaxed. We illustrate this using standard relaxation metrics to establish two dynamical recovery time-scales: i) halos need ~1.5 dynamical times following formation, and ii) ~2 dynamical times following a major (3:1) or minor (10:1) merger to be relaxed. This is remarkably consistent across a wide mass range. Lastly, we use a phase-space halo finder to illustrate that major mergers drive long-lived massive phase-space structures which take many dynamical times to dissipate. This can yield significant differences in the inferred mass build-up of galactic halos and we suggest that care must be taken to ensure a physically meaningful match between the galaxy-formation physics of semi-analytic models and the halo finders supplying their input.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Extended Stellar Components of Galaxies & the Nature of Dark Matter

Deep observations of galaxies reveal faint extended stellar components (hereafter ESCs) of streams, shells, and halos. These are a natural prediction of hierarchical galaxy formation, as accreted satellite galaxies are tidally disrupted by their host. We investigate whether or not global properties of the ESC could be used to test of dark matter, reasoning that they should be sensitive to the abundance of low-mass satellites, and therefore the underlying dark matter model. Using cosmological simulations of galaxy formation in the favoured Cold Dark Matter (CDM) and Warm Dark Matter (WDM) models ($m_{\rm WDM}$=0.5,1,2 keV/$c^2$), which suppress the abundance of low-mass satellites, we find that the kinematics and orbital structure of the ESC is consistent across models. However, we find striking differences in its spatial structure, as anticipated -- a factor of $\sim$10 drop in spherically averaged mass density between $\sim$10% and $\sim$75% of the virial radius in the more extreme WDM runs ($m_{\rm WDM}$=0.5, 1 keV/$c^2$) relative to the CDM run. These differences are consistent with the mass assembly histories of the different components, and are present across redshifts. However, even the least discrepant of the WDM models is incompatible with current observational limits on $m_{\rm WDM}$. Importantly, the differences we observe when varying the underlying dark matter are comparable to the galaxy-to-galaxy variation we expect within a fixed dark matter model. This suggests that it will be challenging to place limits on dark matter using only the unresolved spatial structure of the the ESC.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Matter in the beam: Weak lensing, substructures and the temperature of dark matter

Warm Dark Matter (WDM) models offer an attractive alternative to the current Cold Dark Matter (CDM) cosmological model. We present a novel method to differentiate between WDM and CDM cosmologies, namely using weak lensing; this provides a unique probe as it is sensitive to all the "matter in the beam", not just dark matter haloes and the galaxies that reside in them, but also the diffuse material between haloes. We compare the weak lensing maps of CDM clusters to those in a WDM model corresponding to a thermally produced $0.5$~keV dark matter particle. Our analysis clearly shows that the weak lensing magnification, convergence and shear distributions can be used to distinguish between CDM and WDM models. WDM models {\em increase} the probability of weak magnifications, with the differences being significant to $\gtrsim5σ$, while leaving no significant imprint on the shear distribution. WDM clusters analysed in this work are more homogeneous than CDM ones, and the fractional decrease in the amount of material in haloes is proportional to the average increase in the magnification. This difference arises from matter that would be bound in compact haloes in CDM being smoothly distributed over much larger volumes at lower densities in WDM. Moreover, the signature does not solely lie in the probability distribution function but in the full spatial distribution of the convergence field.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

nIFTY galaxy cluster simulations III: The Similarity & Diversity of Galaxies & Subhaloes

We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) produced by hydrodynamical codes ranging from classic Smooth Particle Hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity only, non-radiative hydrodynamics and full feedback physics runs by looking at the overall subhalo/galaxy distribution and on an individual objects basis. We find the subhalo population is reproduced to within $\lesssim10\%$ for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim0.1$ dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by $\approx20\%$. The galaxy populations also show striking code-to-code variations. Although the Tully-Fisher relation is similar in almost all codes, the number of galaxies with $10^{9}M_\odot/h\lesssim M_*\lesssim 10^{12}M_\odot/h$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of $\sim0.5$ dex in stellar mass. Moreover, strong systematic differences exist, with some codes producing galaxies $70\%$ smaller than others. The diversity partially arises from the inclusion/absence of AGN feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies in all environments remains a challenge. We argue even basic galaxy properties, such as the stellar mass to halo mass, should be treated with errors bars of $\sim0.2-0.4$ dex.

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
preprint2016arXiv

nIFTy galaxy cluster simulations V: Investigation of the Cluster Infall Region

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $Λ$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the dynamical state of galaxy clusters: insights from cosmological simulations II

Using a suite of cosmology simulations of a sample of $> 120$ galaxy clusters with $\log(M_{DM, vir}) \le 14.5$. We compare clusters that form in purely dark matter run and their counterparts in hydro runs and investigate 4 independent parameters, that are normally used to classify dynamical state. We find that the virial ratio $η$ in hydro-dynamical runs is $\sim 10$ per cent lower than in the DM run, and there is no clear separation between the relaxed and unrelaxed clusters for any parameter. Further, using the velocity dispersion deviation parameter $ζ$, which is defined as the ratio between cluster velocity dispersion $σ$ and the theoretical prediction $σ_t = \sqrt{G M_{total}/R}$, we find that there is a linear correlation between the virial ratio $η$ and this $ζ$ parameter. We propose to use this $ζ$ parameter, which can be easily derived from observed galaxy clusters, as a substitute of the $η$ parameter to quantify the cluster dynamical state.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Stability of Satellite Planes I: Effects of Mass, Velocity, Halo Shape and Alignment

The recently discovered vast thin plane of dwarf satellites orbiting the Andromeda Galaxy (M31) adds to the mystery of the small scale distribution of the Local Group's galaxy population. Such well defined planar structures are apparently rare occurrences in cold dark matter cosmological simulations, and we lack a coherent explanation of their formation and existence. In this paper, we explore the long-term survivability of thin planes of dwarfs in galactic halos, focusing, in particular, on systems mimicking the observed Andromeda distribution. The key results show that, in general, planes of dwarf galaxies are fragile, sensitive to the shape of the dark matter halo and other perturbing effects. In fact, long lived planes of satellites only exist in polar orbits in spherical dark matter halos, presenting a challenge to the observed Andromeda plane which is significantly tilted with respect to the optical disk. Our conclusion is that, in standard cosmological models, planes of satellites are generally short lived, and hence we must be located at a relatively special time in the evolution of the Andromeda Plane, lucky enough to see its coherent pattern.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spurious Small-Scale Structure & Discreteness-Driven Relaxation in Cosmological Simulations

There is strong evidence that cosmological N-body simulations dominated by Warm Dark Matter (WDM) contain spurious or unphysical haloes, most readily apparent as regularly spaced low-mass haloes strung along filaments. We show that spurious haloes are a feature of traditional N-body simulations of cosmological structure formation models, including WDM and Cold Dark Matter (CDM) models, in which gravitational collapse proceeds in an initially anisotropic fashion, and arises naturally as a consequence of discreteness-driven relaxation. We demonstrate this using controlled N-body simulations of plane-symmetric collapse and show that spurious haloes are seeded at shell crossing by localised velocity perturbations induced by the discrete nature of the density field, and that their characteristic separation should be approximately the mean inter-particle separation of the N-body simulation, which is fixed by the mass resolution within the volume. Using cosmological N-body simulations in which particles are split into two collisionless components of fixed mass ratio, we find that the spatial distribution of the two components show signatures of discreteness-driven relaxation in their spatial distribution on both large and small scales. Adopting a spline kernel gravitational softening that is of order the comoving mean inter-particle separation helps to suppress the effect of discreteness-driven relaxation, but cannot eliminate it completely. These results provide further motivation for recent developments of new algorithms, which include, for example, revisions of the traditional N-body approach by means of spatially adaptive anistropric gravitational softenings or explicit solutions for the evolution of dark matter in phase space.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

HALOGEN: A tool for fast generation of mock halo catalogues

We present a simple method of generating approximate synthetic halo catalogues: HALOGEN. This method uses a combination of $2^{nd}$-order Lagrangian Perturbation Theory (2LPT) in order to generate the large-scale matter distribution, analytical mass functions to generate halo masses, and a single-parameter stochastic model for halo bias to position haloes. HALOGEN represents a simplification of similar recently published methods. Our method is constrained to recover the 2-point function at intermediate ($10Mpc/h<r<50Mpc/h$) scales, which we show is successful to within 2 per cent. Larger scales ($\sim100Mpc/h$) are reproduced to within 15 per cent. We compare several other statistics (e.g. power spectrum, point distribution function, redshift space distortions) with results from N-Body simulations to determine the validity of our method for different purposes. One of the benefits of HALOGEN is its flexibility, and we demonstrate this by showing how it can be adapted to varying cosmologies and simulation specifications. A driving motivation for the development of such approximate schemes is the need to compute covariance matrices and study the systematic errors for large galaxy surveys, which requires thousands of simulated realisations. We discuss the applicability of our method in this context, and conclude that it is well suited to mass production of appropriate halo catalogues. The code is publicly available at https://github.com/savila/halogen

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hidden from view: Coupled Dark Sector Physics and Small Scales

We study cluster mass dark matter haloes, their progenitors and surroundings in an coupled Dark Matter-Dark Energy model and compare it to quintessence and $Λ$CDM models with adiabatic zoom simulations. When comparing cosmologies with different expansions histories, growth functions & power spectra, care must be taken to identify unambiguous signatures of alternative cosmologies. Shared cosmological parameters, such as $σ_8$, need not be the same for optimal fits to observational data. We choose to set our parameters to $Λ$CDM $z=0$ values. We find that in coupled models, where DM decays into DE, haloes appear remarkably similar to $Λ$CDM haloes despite DM experiencing an additional frictional force. Density profiles are not systematically different and the subhalo populations have similar mass, spin, and spatial distributions, although (sub)haloes are less concentrated on average in coupled cosmologies. However, given the scatter in related observables ($V_{\rm max},R_{V_{\rm max}}$), this difference is unlikely to distinguish between coupled and uncoupled DM. Observations of satellites of MW and M31 indicate a significant subpopulation reside in a plane. Coupled models do produce planar arrangements of satellites of higher statistical significance than $Λ$CDM models, however, in all models these planes are dynamically unstable. In general, the nonlinear dynamics within and near large haloes masks the effects of a coupled dark sector. The sole environmental signature we find is that small haloes residing in the outskirts are more deficient in baryons than their $Λ$CDM counterparts. The lack of a pronounced signal for a coupled dark sector strongly suggests that such a phenomena would be effectively hidden from view.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

How does our choice of observable influence our estimation of the centre of a galaxy cluster? Insights from cosmological simulations

Galaxy clusters are an established and powerful test-bed for theories of both galaxy evolution and cosmology. Accurate interpretation of cluster observations often requires robust identification of the location of the centre. Using a statistical sample of clusters drawn from a suite of cosmological simulations in which we have explored a range of galaxy formation models, we investigate how the location of this centre is affected by the choice of observable - stars, hot gas, or the full mass distribution as can be probed by the gravitational potential. We explore several measures of cluster centre: the minimum of the gravitational potential, which would expect to define the centre if the cluster is in dynamical equilibrium; the peak of the density; the centre of BCG; and the peak and centroid of X-ray luminosity. We find that the centre of BCG correlates more strongly with the minimum of the gravitational potential than the X-ray defined centres, while AGN feedback acts to significantly enhance the offset between the peak X-ray luminosity and minimum gravitational potential. These results highlight the importance of centre identification when interpreting clusters observations, in particular when comparing theoretical predictions and observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

nIFTy Cosmology: Comparison of Galaxy Formation Models

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics

We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.

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

nIFTy galaxy cluster simulations II: radiative models

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $Λ$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Novel Adaptive softening for collisionless N-body simulations: Eliminating spurious halos

We describe a NOVel form of Adaptive softening (NovA) for collisionless $N$-body simulations, implemented in the Ramses adaptive mesh refinement code. We introduce a refinement criterion that the particle distribution within each cell be sufficiently isotropic, as measured by its moment of inertia tensor. In this way, collapse is only refined if it occurs along all three axes, ensuring that the softening $ε$ is always of order twice the largest inter-particle spacing in a cell. This more conservative force softening criterion is designed to minimise spurious two-body effects, while maintaining high force resolution in collapsed regions of the flow. We test NovA using an antisymmetric perturbed plane wave collapse (`Valinia' test) before applying it to warm dark matter (WDM) simulations. For the Valinia test, we show that -- unlike the standard $N$-body method -- NovA produces no numerical fragmentation while still being able to correctly capture fine caustics and shells around the collapsing regions. For the WDM simulations, we find that NovA converges significantly more rapidly than standard $N$-body, producing little or no spurious halos on small scales. We show, however, that determining whether or not halos exist below the free streaming mass $M_{\rm fs}$ is complicated by the fact that our halo finder (AHF) likely incorrectly labels some caustics and criss-crossing filaments as halos, while one or two particularly massive filaments appear to fragment in any version of NovA where refinement is allowed. Such massive filaments may be physically unstable to collapse, as is the case for infinite, static, self-gravitating cylinders. We will use NovA in forthcoming papers to study the issue of halo formation below $M_{\rm fs}$; filament stability; and to obtain new constraints on the temperature of dark matter.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulating feedback from nuclear clusters: the impact of multiple sources

Nuclear star clusters (NCs) are found to exist in the centres of many galaxies and appear to follow scaling relations similar to those of super-massive black holes. Previous analytical work has suggested that such relations are a consequence of feedback regulated growth. We explore this idea using high resolution hydrodynamical simulations, focusing on the validity of the simplifying assumptions made in analytical models. In particular, we investigate feedback emanating from multiple stellar sources rather than from a single source, as is usually assumed, and show that collisions betweens shells of gas swept up by feedback leads to momentum cancellation and the formation of high density clumps and filaments. This high density material is resistant both to expulsion from the galaxy potential and to disruption by feedback; if it falls back onto the NC, we expect the gas to be available for further star formation or for feeding a central black hole. We also note our results may have implications for the evolution of globular clusters and stellar clusters in high redshift dark matter halos.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Gigaparsec WiggleZ Simulations: Characterising scale dependant bias and associated systematics in growth of structure measurements

We use the Gigaparsec WiggleZ (GiggleZ) simulations to characterise galaxy bias and its scale dependence for a range of redshifts and halo masses in a standard $Λ$LCDM cosmology. Assuming bias converges to a scale independent form at large scales, we develop a phenomenological model which fully expresses the mass and redshift dependence of bias and its scale dependence in real or redshift space. We then use this to illustrate how scale-dependent bias can systematically skew measurements of the growth-rate of cosmic structure obtained from redshift-space distortion measurements. When data is fit only to scales $k_{\rm max}{\le}0.1$ $[h^{-1} \rm{Mpc}]^{-1}$, we find that these effects are significant only for large biases ($b{\gtrsim}3$) at large redshifts ($z{\gtrsim}1$). However, when smaller scales are incorporated ($k_{\rm max}{\le}0.2$ $[h^{-1} \rm{Mpc}]^{-1}$) to increase measurement precision, the combination of reduced statistical uncertainty and increased scale dependent bias can result in highly significant systematics for most large halos across all redshifts. We identify several new interesting aspects of bias, including a significant large-scale bias boost for small halos at low-redshifts due to substructure effects ($\sim$20\% for Milky Way-like systems) and a nearly redshift-independent halo mass (corresponding to a redshift-space bias of ${\sim}1.5$) for which halo bias has little-or-no scale dependence on scales greater than $3$ $[h^{-1} {\rm Mpc}]$. This suggests an optimal strategy of targeting bias ${\sim}{1.5}$ systems for clustering studies which are dominated by systematic uncertainties in how observed halo (or galaxy) distributions map to their underlying mass distribution, such as cosmological measurements of neutrino masses. Code for generating our fitting formula is publicly available at http://gbpoole.github.io/Poole_2014a_code/ (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Wide Area VISTA Extra-galactic Survey (WAVES)

The "Wide Area VISTA Extra-galactic Survey" (WAVES) is a 4MOST Consortium Design Reference Survey which will use the VISTA/4MOST facility to spectroscopically survey ~2million galaxies to $r_{\rm AB} < 22$ mag. WAVES consists of two interlocking galaxy surveys ("WAVES-Deep" and "WAVES-Wide"), providing the next two steps beyond the highly successful 1M galaxy Sloan Digital Sky Survey and the 250k Galaxy And Mass Assembly survey. WAVES will enable an unprecedented study of the distribution and evolution of mass, energy, and structures extending from 1-kpc dwarf galaxies in the local void to the morphologies of 200-Mpc filaments at $z\sim1$. A key aim of both surveys will be to compare comprehensive empirical observations of the spatial properties of galaxies, groups, and filaments, against state-of-the-art numerical simulations to distinguish between various Dark Matter models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravitational lensing in WDM cosmologies: The cross section for giant arcs

The nature of the dark sector of the Universe remains one of the outstanding problems in modern cosmology, with the search for new observational probes guiding the development of the next generation of observational facilities. Clues come from tension between the predictions from ΛCDM and observations of gravitationally lensed galaxies. Previous studies showed that galaxy clusters in the ΛCDM are not strong enough to reproduce the observed number of lensed arcs. This work aims to constrain the warm dark matter cosmologies by means of the lensing efficiency of galaxy clusters drawn from these alternative models. The lensing characteristics of two samples of simulated clusters in the warm dark matter (ΛWDM) and cold dark matter (ΛCDM) cosmologies have been studied. The results show that even though the CDM clusters are more centrally concentrated and contain more substructures, the WDM clusters have slightly higher lensing efficiency than their CDM counterparts. The key difference is that WDM clusters have more extended and more massive subhaloes than CDM analogues. These massive substructures significantly stretch the critical lines and caustics and hence they boost the lensing efficiency of the host halo. Despite the increase in the lensing efficiency due to the contribution of massive substructures in the WDM clusters, this is not enough to resolve the arc statistics problem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Heating and Ionization of the Primordial Intergalactic Medium by High Mass X-ray Binaries

We investigate the influence of High Mass X-ray Binaries on their high redshift environments. Using a one-dimensional radiative transfer code, we predict the ionization and temperature profiles surrounding a coeval stellar population, composed of main sequence stars and HMXBs, at various times after its formation. We consider both uniform density surroundings, and a cluster embedded in a 10^8 solar mass NFW halo. HMXBs in a constant density environment produce negligible enhanced ionization because of their high-energy SEDs and short lifetimes. In this case, HMXBs only marginally contribute to the local heating rate. For NFW profiles, radiation from main sequence stars cannot prevent the initially ionized volume from recombining since it is unable to penetrate the high density galactic core. However, HMXB photons stall recombinations behind the front, keeping it partially ionized for longer. The increased electron density in these partially ionized regions promotes further cooling, resulting in lower IGM temperatures. In the context of this starburst model, we have shown that HMXBs do not make a major contribution to reionization or IGM heating. However, X-ray escape fractions are high in both density profile cases. Continuous star formation may result in the build up of X-rays over time, reducing the ionization timescale and potentially leading to low level ionization of the distant IGM.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stochastic Star Formation & Feedback: Mapping Low-Mass Galaxies to Dark Matter Haloes

Comparison of observed satellite galaxies of the Milky Way (hereafter MW) with dark matter subhaloes in cosmological $N$-body simulations of MW-mass haloes suggest that such subhaloes, if they exist, are occupied by satellites in a stochastic fashion. We examine how inefficient massive star formation and associated supernova feedback in high-redshift progenitors of present-day low-mass subhaloes might contribute to this stochasticity. Using a Monte Carlo approach to follow the assembly histories of present-day low-mass haloes with $10^7 \lesssim M \leq 10^{10}$ ${\rm M}_{\odot}$, we identify when cooling and star formation is likely to proceed, and observe that haloes with present-day masses $\lesssim 10^9 {\rm M}_{\odot}$ never grow sufficiently massive to support atomic hydrogen line cooling. Noting that the star formation timescale decreases sharply with stellar mass as $t_{\rm PMS} \propto m_{\ast}^{-2.5}$, we argue that, should the conditions for high mass star formation arise in low-mass haloes, the ensuing supernovae are likely to disrupt ongoing lower-mass star formation and unbind gas within the halo. This potentially star-forming gas is unlikely to be replenished in lower mass haloes because of, e.g. cosmological reionization, and so we expect galaxy formation to be stymied in a manner that depends on host halo assembly history and the efficiency and timing of star formation in proto-galaxies, which we illustrate using a Monte Carlo model. Based on these simple physical arguments, we assert that stochasticity of star formation and feedback is an essential but overlooked ingredient in modelling galaxy formation on the smallest scales.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Warm Dark Haloes Accretion Histories and their Gravitational Signatures

We study clusters in Warm Dark Matter (WDM) models of a thermally produced dark matter particle $0.5$ keV in mass. We show that, despite clusters in WDM cosmologies having similar density profiles as their Cold Dark Matter (CDM) counterparts, the internal properties, such as the amount of substructure, shows marked differences. This result is surprising as clusters are at mass scales that are {\em a thousand times greater} than that at which structure formation is suppressed. WDM clusters gain significantly more mass via smooth accretion and contain fewer substructures than their CDM brethren. The higher smooth mass accretion results in subhaloes which are physically more extended and less dense. These fine-scale differences can be probed by strong gravitational lensing. We find, unexpectedly, that WDM clusters have {\em higher} lensing efficiencies than those in CDM cosmologies, contrary to the naive expectation that WDM clusters should be less efficient due to the fewer substructures they contain. Despite being less dense, the larger WDM subhaloes are more likely to have larger lensing cross-sections than CDM ones. Additionally, WDM subhaloes typically reside at larger distances, which radially stretches the critical lines associated with strong gravitational lensing, resulting in excess in the number of clusters with large radial cross-sections at the $\sim2σ$ level. Though lensing profile for an individual cluster vary significantly with the line-of-sight, the radial arc distribution based on a sample of $\gtrsim100$ clusters may prove to be the crucial test for the presence of WDM.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Confronting Cold Dark Matter Predictions with Observed Galaxy Rotations

The rich statistics of galaxy rotations as captured by the velocity function (VF) provides invaluable constraints on galactic baryon physics and the nature of dark matter (DM). However, the comparison of observed galaxy rotations against cosmological models is prone to subtle caveats that can easily lead to misinterpretations. Our analysis reveals full statistical consistency between ~5000 galaxy rotations, observed in line-of-sight projection, and predictions based on the standard cosmological model (LCDM) at the mass-resolution of the Millennium simulation (HI line-based circular velocities above ~50 km/s). Explicitly, the HI linewidths in the HI Parkes All Sky Survey (HIPASS) are found consistent with those in S3-SAX, a post-processed semi-analytic model for the Millennium simulation. Previously found anomalies in the VF can be plausibly attributed to (1) the mass-limit of the Millennium simulation, (2) confused sources in HIPASS, (3) inaccurate inclination measurements for optically faint sources, and (4) the non-detectability of gas-poor early-type galaxies. These issues can be bypassed by comparing observations and models using linewidth source counts rather than VFs. We investigate if and how well such source counts can constrain the temperature of DM.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

HMFcalc: An Online Tool for Calculating Dark Matter Halo Mass Functions

The dark matter halo mass function (HMF) is a characteristic property of cosmological structure formation models, quantifying the number density of dark matter haloes per unit mass in the Universe. A key goal of current and planned large galaxy surveys is to measure the HMF and to use it to test theories of dark matter and dark energy. We present a new web application for calculating the HMF -- the frontend HMFcalc and the engine hmf. HMFcalc has been designed to be flexible, efficient and easy to use, providing observational and theoretical astronomers alike with the means to explore standard functional forms of the HMF or to tailor their own. We outline the theoretical background needed to compute the HMF, we show how it has been implemented in hmf, and finally we provide worked examples that illustrate HMFcalc's versatility as an analysis tool.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How well do we know the Halo Mass Function?

The parameters governing the standard ΛCold Dark Matter cosmological model have been constrained with unprecedented accuracy by precise measurements of the cosmic microwave background by the Wilkinson Microwave Anisotropy Probe (WMAP) and Planck satellites. Each new data release has refined further our knowledge of quantities -- such as the matter density parameter Ω_M -- that are imprinted on the dark matter halo mass function (HMF), a powerful probe of dark matter and dark energy models. In this letter we trace how changes in the cosmological parameters over the last decade have influenced uncertainty in our knowledge of the HMF. We show that this uncertainty has reduced significantly since the 3rd WMAP data release, but the rate of this reduction is slowing. This is limited by uncertainty in the normalisation σ_8, whose influence is most pronounced at the high mass end of the mass function. Interestingly, we find that the accuracy with which we can constrain the HMF in terms of the cosmological parameters has now reached the point at which it is comparable to the scatter in HMF fitting functions. This suggests that the power of the HMF as a precision probe of dark matter and dark energy hinges on more accurate determination of the theoretical HMF. Finally, we assess prospects of using the HMF to differentiate between Cold and Warm Dark Matter models based on ongoing improvements in measurements of Ω_M, and we comment briefly on optimal survey strategies for constraining dark matter and dark energy models using the HMF.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modelling the Growth of Supermassive Black Holes in Cosmological Simulations

There is strong evidence that supermassive black holes reside in all galaxies that contain a stellar spheroid and their mass is tightly correlated with properties such as stellar bulge mass and velocity dispersion. There are also strong theoretical arguments that feedback from supermassive black holes plays an important role in shaping the high mass end of the galaxy mass function, hence to accurately model galaxies we also need to model the black holes. We present a comparison of two black hole growth models implemented within a large-scale, cosmological SPH simulation including star formation and feedback. One model is a modified Bondi-Hoyle prescription that grows black holes based on the smooth density of local gas, while the other is the recently proposed Accretion Disc Particle (ADP) method. This model swallows baryonic particles that pass within an accretion radius of the black hole and adds them to a subgrid accretion disc. Black holes are then grown by material from this disc. We find that both models can reproduce local scaling relations, although the ADP model is offset from the observed relations at high black hole masses. The total black hole mass density agrees between models to within a factor of three, but both struggle to reproduce the black hole mass function. The simulated mass functions are too steep and underestimate the number of intermediate and high mass black holes. In addition, the ADP model swallows excessive amounts of material at the resolution of large-scale, cosmological simulations producing unrealistically large accretion discs. Future work needs to be performed to improve the black hole mass function within simulations. This should be done through the mass growth and feedback as they are strongly coupled and should not be treated as separate entities.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quenching star formation in cluster galaxies

In order to understand the processes that quench star formation within rich clusters, we construct a library of subhalo orbits drawn from $Λ$CDM cosmological N-body simulations of four rich clusters. The orbits are combined with models of star formation followed by quenching in the cluster environment. These are compared with observed bulge and disc colours and stellar absorption linestrength indices of satellite galaxies. Models in which the bulge stellar populations depend only on the galaxy subhalo mass while the disc quenching depends on the cluster environment are acceptable fits to the data. An exponential disc quenching timescale of 3 - 3.5 Gyr is preferred. Models with short ($\lesssim 1$ Gyr) quenching timescales yield cluster-centric gradients in disc colours and Balmer line indices that are too steep compared to observations. We also examine models in which there is quenching in lower mass groups prior to cluster infall ("pre-processing"), finding that such models are a better fit to the data than models without pre-processing and require similar quenching times. The data slightly prefer models where quenching occurs only for galaxies falling within about 0.5 $r_{200}$. Finally, we have examined models with short quenching timescales of 1 Gyr, but a long delay time of 3 Gyr prior to quenching. All models with short quenching timescales, even such "delayed-then-rapid" quenching models, produce excessively red galaxies near the cluster core and are strongly disfavoured by the data. These results imply that the environments of rich clusters must impact star formation rates of infalling galaxies on relatively long timescales -- several times longer than a typical halo spends within the virial radius of a cluster. This scenario favours gentler quenching mechanisms such as slow "strangulation" over more rapid ram-pressure stripping.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Seeking Observable Imprints of Small-Scale Structure on the Properties of Dark Matter Haloes

The characteristic prediction of the Cold Dark Matter (CDM) model of cosmological structure formation is that the Universe should contain a wealth of small-scale structure -- low-mass dark matter haloes and subhaloes. However, galaxy formation is inefficient in their shallow potential wells and so we expect these low-mass haloes and subhaloes to be dark. Can we tell the difference between a Universe in which low-mass haloes are present but dark and one in which they never formed, thereby providing a robust test of the CDM model? We address this question using cosmological simulations to examine how properties of low-mass haloes that are potentially accessible to observation, such as their spatial clustering, rate of accretions and mergers onto massive galaxies and the angular momentum content of massive galaxies, differ between a LCDM model and dark matter models in which low-mass halo formation is suppressed. Adopting an effective cut-off mass scale M_cut below which small-scale power is suppressed in the initial conditions, we study dark matter models in which M_cut varies between 5e9 Msol/h and 1e11 Msol/h, equivalent to the host haloes of dwarf and low mass galaxies. Our results show that both the clustering strength of low-mass haloes around galaxy-mass primaries and the rate at which they merge with these primaries is sensitive to the assumed value of M_cut; in contrast, suppressing low-mass halo formation has little influence on the angular momentum content of galaxy-mass haloes -- it is the quiescence or violence of a halo's assembly history that has a more marked effect. However, we expect that measuring the effect on spatial clustering or the merger rate is likely to be observationally difficult for realistic values of M_cut, and so isolating the effect of this small-scale structure would appear to be remarkably difficult to detect, at least in the present day Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Robust Measure of Cosmic Structure beyond the Power-Spectrum: Cosmic Filaments and the Temperature of Dark Matter

We discover that the mass of dark matter particles mDM is imprinted in phase-correlations of the cosmic density field more significantly than in the 2-point correlation. In particular, phase-correlations trace mDM out to scales about five times larger than the 2-point correlation. This result relies on a new estimator l(r) of pure phase-information in Fourier space, which can be interpreted as a parameter-free and scale-invariant tracer of filament-like structure. Based on simulated density fields we show how mDM can, in principle, be measured using l(r), given a suitably reconstructed density field.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Feedback from High-Mass X-Ray Binaries on the High Redshift Intergalactic Medium : Model Spectra

Massive stars at redshifts z > 6 are predicted to have played a pivotal role in cosmological reionization as luminous sources of ultra-violet (UV) photons. However, the remnants of these massive stars could be equally important as X-ray luminous (L_X 1e38 erg/s) high-mass X-ray binaries (HMXBs). Because the absorption cross section of neutral hydrogen decreases sharply with photon energy (proportional to the inverse cube), X-rays can escape more freely than UV photons from the star-forming regions in which they are produced, allowing HMXBs to make a potentially significant contribution to the ionizing X-ray background during reionization. In this paper, we explore the ionizing power of HMXBs at redshifts z > 6 using a Monte Carlo model for a coeval stellar population of main sequence stars and HMXBs. Using the archetypal Galactic HMXB Cygnus X-1 as our template, we propose a composite HMXB spectral energy distribution consisting of black-body and power-law components, whose contributions depend on the accretion state of the system. We determine the time-dependent ionizing power of a combined population of UV-luminous stars and X-ray luminous HMXBs, and deduce fitting formulae for the boost in the population's ionizing power arising from HMXBs; these fits allow for simple implementation of HMXB feedback in numerical simulations. Based on this analysis, we estimate the contribution of high redshift HMXBs to the present-day soft X-ray background, and we show that it is a factor of ~100-1000 smaller than the observed limit. Finally, we discuss the implications of our results for the role of HMXBs in reionization and in high redshift galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modelling supermassive black hole growth: towards an improved sub-grid prescription

Accretion onto supermassive black holes (SMBHs) in galaxy formation simulations is frequently modelled by the Bondi-Hoyle formalism. Here we examine the validity of this approach analytically and numerically. We argue that the character of the flow where one evaluates the gas properties is unlikely to satisfy the simple Bondi-Hoyle model. Only in the specific case of hot virialised gas with zero angular momentum and negligible radiative cooling is the Bondi-Hoyle solution relevant. In the opposite extreme, where the gas is in a state of free-fall at the evaluation radius due to efficient cooling and the dominant gravity of the surrounding halo, the Bondi-Hoyle formalism can be erroneous by orders of magnitude in either direction. This may impose artificial trends with halo mass in cosmological simulations by being wrong by different factors for different halo masses. We propose an expression for the sub-grid accretion rate which interpolates between the free-fall regime and the Bondi-Hoyle regime, therefore taking account of the contribution of the halo to the gas dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The observed $M - σ$ relations imply that SMBHs grow by cold chaotic accretion

We argue that current observations of $M - σ$ relations for galaxies can be used to constrain theories of super-massive black holes (SMBH) feeding. In particular, assuming that SMBH mass is limited only by the feedback on the gas that feeds it, we show that SMBHs fed via a planar galaxy scale gas flow, such as a disc or a bar, should be much more massive than their counterparts fed by quasi-spherical inflows. This follows from the relative inefficiency of AGN feedback on a flattened inflow. We find that even under the most optimistic conditions for SMBH feedback on flattened inflows, the mass at which the SMBH expels the gas disc and terminates its own growth is a factor of several higher than the one established for quasi-spherical inflows. Any beaming of feedback away from the disc and any disc self-shadowing strengthens this result further. Contrary to this theoretical expectation, recent observations have shown that SMBH in pseudobulge galaxies (which are associated with barred galaxies) are typically under- rather than over-massive when compared with their classical bulge counterparts at a fixed value of $σ$. We conclude from this that SMBHs are not fed by large (100 pc to many kpc) scale gas discs or bars, most likely because such planar flows are turned into stars too efficiently to allow any SMBH growth. Based on this and other related observational evidence, we argue that most SMBHs grow by chaotic accretion of gas clouds with a small and nearly randomly distributed direction of angular momentum.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermal instabilities in cooling galactic coronae: fuelling star formation in galactic discs

We investigate the means by which cold gas can accrete onto Milky Way mass galaxies from a hot corona of gas, using a new smoothed particle hydrodynamics code, 'SPHS'. We find that the 'cold clumps' seen in many classic SPH simulations in the literature are not present in our SPHS simulations. Instead, cold gas condenses from the halo along filaments that form at the intersection of supernovae-driven bubbles from previous phases of star formation. This positive feedback feeds cold gas to the galactic disc directly, fuelling further star formation. The resulting galaxies in the SPH and SPHS simulations differ greatly in their morphology, gas phase diagrams, and stellar content. We show that the classic SPH cold clumps owe to a numerical thermal instability caused by an inability for cold gas to mix in the hot halo. The improved treatment of mixing in SPHS suppresses this instability leading to a dramatically different physical outcome. In our highest resolution SPHS simulation, we find that the cold filaments break up into bound clumps that form stars. The filaments are overdense by a factor of 10-100 compared to the surrounding gas, suggesting that the fragmentation results from a physical non-linear instability driven by the overdensity. This 'fragmenting filament' mode of disc growth has important implications for galaxy formation, in particular the role of star formation in bringing cold gas into disc galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cosmic evolution of the atomic and molecular gas content of galaxies

We study the evolution of the cold gas content of galaxies by splitting the interstellar medium into its atomic and molecular hydrogen components, using the galaxy formation model GALFORM in the LCDM framework. We calculate the molecular-to-atomic hydrogen mass ratio, H2/HI, in each galaxy using two different approaches; the pressure-based empirical relation of Blitz & Rosolowsky and the theoretical model of Krumholz, McKeee & Tumlinson, and apply them to consistently calculate the star formation rates of galaxies. We find that the model based on the Blitz & Rosolowsky law predicts an HI mass function, CO(1-0) luminosity function, correlations between the H2/HI ratio and stellar and cold gas mass, and infrared-CO luminosity relation in good agreement with local and high redshift observations. The HI mass function evolves weakly with redshift, with the number density of high mass galaxies decreasing with increasing redshift. In the case of the H2 mass function, the number density of massive galaxies increases strongly from z=0 to z=2, followed by weak evolution up to z=4. We also find that the H2/HI ratio of galaxies is strongly dependent on stellar and cold gas mass, and also on redshift. The slopes of the correlations between H2/HI and stellar and cold gas mass hardly evolve, but the normalisation increases by up to two orders of magnitude from z=0-8. The strong evolution in the H2 mass function and the H2/HI ratio is primarily due to the evolution in the sizes of galaxies and secondarily, in the gas fractions. The predicted cosmic density evolution of HI agrees with the observed evolution inferred from DLAs, and is dominated by low/intermediate mass halos. We find that previous theoretical studies have largely overestimated the redshift evolution of the global H2/HI ratio due to limited resolution. We predict a maximum of rho_H2/rho_HI~1.2 at z~3.5.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Feeding SMBHs through supersonic turbulence and ballistic accretion

It has long been recognised that the main obstacle to accretion of gas onto supermassive black holes (SMBHs) is large specific angular momentum. It is feared that the gas settles in a large scale disc, and that accretion would then proceed too inefficiently to explain the masses of the observed SMBHs. Here we point out that, while the mean angular momentum in the bulge is very likely to be large, the deviations from the mean can also be significant. Indeed, cosmological simulations show that velocity and angular momentum fields of gas flows onto galaxies are very complex. Furthermore, inside bulges the gas velocity distribution can be further randomised by the velocity kicks due to feedback from star formation. We perform hydrodynamical simulations of gaseous rotating shells infalling onto an SMBH, attempting to quantify the importance of velocity dispersion in the gas at relatively large distances from the black hole. We implement this dispersion by means of a supersonic turbulent velocity spectrum. We find that, while in the purely rotating case the circularisation process leads to efficient mixing of gas with different angular momentum, resulting in a low accretion rate, the inclusion of turbulence increases this accretion rate by up to several orders of magnitude. We show that this can be understood based on the notion of "ballistic" accretion, whereby dense filaments, created by convergent turbulent flows, travel through the ambient gas largely unaffected by hydrodynamical drag. We derive a simple analytical formula that captures the numerical results to within a factor of a few. Rescaling our results to astrophysical bulges, we argue that this "ballistic" mode of accretion could provide the SMBHs with a sufficient supply of fuel without the need to channel the gas via large-scale discs or bars, and therefore that star formation in bulges can be a strong catalyst for SMBH accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Dynamical State of Dark Matter Haloes in Cosmological Simulations I: Correlations with Mass Assembly History

Using a statistical sample of dark matter haloes drawn from a suite of cosmological N-body simulations of the Cold Dark Matter (CDM) model, we quantify the impact of a simulated halo's mass accretion and merging history on two commonly used measures of its dynamical state, the virial ratio eta and the centre of mass offset Delta r. Quantifying this relationship is important because the degree to which a halo is dynamically equilibrated will influence the reliability with which we can measure characteristic equilibrium properties of the structure and kinematics of a population of haloes. We begin by verifying that a halo's formation redshift zform correlates with its virial mass Mvir and we show that the fraction of its recently accreted mass and the likelihood of it having experienced a recent major merger increases with increasing Mvir and decreasing zform. We then show that both eta and Delta r increase with increasing Mvir and decreasing zform, which implies that massive recently formed haloes are more likely to be dynamically unrelaxed than their less massive and older counterparts. Our analysis shows that both eta and Delta r are good indicators of a halo's dynamical state, showing strong positive correlations with recent mass accretion and merging activity, but we argue that Delta r provides a more robust and better defined measure of dynamical state for use in cosmological N-body simulations at z~0. We find that Delta r < 0.04 is sufficient to pick out dynamically relaxed haloes at z=0. Finally, we assess our results in the context of previous studies, and consider their observational implications.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Accretion Disc Particle Method for Simulations of Black Hole Feeding and Feedback

Black holes grow by accreting matter from their surroundings. However, angular momentum provides an efficient natural barrier to accretion and so only the lowest angular momentum material will be available to feed the black holes. The standard sub-grid model for black hole accretion in galaxy formation simulations - based on the Bondi-Hoyle method - does not account for the angular momentum of accreting material, and so it is unclear how representative the black hole accretion rate estimated in this way is likely to be. In this paper we introduce a new sub-grid model for black hole accretion that naturally accounts for the angular momentum of accreting material. Both the black hole and its accretion disc are modelled as a composite accretion disc particle. Gas particles are captured by the accretion disc particle if and only if their orbits bring them within its accretion radius R_acc, at which point their mass is added to the accretion disc and feeds the black hole on a viscous timescale t_visc. The resulting black hole accretion rate (dM/dt)_BH powers the accretion luminosity L_acc ~ (dM/dt)_BH, which drives black hole feedback. Using a series of controlled numerical experiments, we demonstrate that our new accretion disc particle method is more physically self-consistent than the Bondi-Hoyle method. We also discuss the physical implications of the accretion disc particle method for systems with a high degree of rotational support, and we argue that the M_BH-sigma relation in these systems should be offset from the relation for classical bulges and ellipticals, as appears to be observed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Accuracy of Subhalo Detection

With the ever increasing resolution of N-body simulations, accurate subhalo detection is becoming essential in the study of the formation of structure, the production of merger trees and the seeding of semi-analytic models. To investigate the state of halo finders, we compare two different approaches to detecting subhaloes; the first based on overdensities in a halo and the second being adaptive mesh refinement. A set of stable mock NFW dark matter haloes were produced and a subhalo was placed at different radii within a larger halo. SUBFIND (a Friends-of-Friends based finder) and AHF (an adaptive mesh based finder) were employed to recover the subhalo. As expected, we found that the mass of the subhalo recovered by SUBFIND has a strong dependence on the radial position and that neither halo finder can accurately recover the subhalo when it is very near the centre of the halo. This radial dependence is shown to be related to the subhalo being truncated by the background density of the halo and originates due to the subhalo being defined as an overdensity. If the subhalo size is instead determined using the peak of the circular velocity profile, a much more stable value is recovered. The downside to this is that the maximum circular velocity is a poor measure of stripping and is affected by resolution. For future halo finders to recover all the particles in a subhalo, a search of phase space will need to be introduced.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Simulations of momentum feedback by black hole winds

The observed super-massive black hole (SMBH) mass -- galaxy velocity dispersion ($M_{\rm bh} - σ$) correlation may be established when winds/outflows from the SMBH drive gas out of the potential wells of classical bulges. Here we present numerical simulations of this process in a static isothermal potential. Simple spherically symmetric models of SMBH feedback at the Eddington luminosity can successfully explain the $M_{\rm bh} - σ$ and nuclear cluster mass $M_{\rm NC}-σ$ correlations, as well as why larger bulges host SMBHs while smaller ones host nuclear star clusters. However these models do not specify how SMBHs feed on infalling gas whilst simultaneously producing feedback that drives gas out of the galaxy. More complex models with rotation and/or anisotropic feedback allow SMBHs to feed via a disc or regions not exposed to SMBH winds, but in these more realistic cases it is not clear why a robust $M_{\rm bh} - σ$ relation should be established. In fact, some of the model predictions contradict observations. For example, an isotropic SMBH wind impacting on a disc (rather than a shell) of aspect ratio $H/R \ll 1$ requires the SMBH mass to be larger by a factor $\sim R/H$, which is opposite to what is observed. We conclude that understanding how a SMBH feeds is as important a piece of the puzzle as understanding how its feedback affects its host galaxy. Finally, we note that in aspherical cases the SMBH outflows induce differential motions in the bulge. This may pump turbulence that is known to hinder star formation in star forming regions. SMBH feedback thus may not only drive gas out of the bulge but also reduce the fraction of gas turned into stars.

0 citations0 reviewsNo code linkedOpen access