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

Andrew J. Benson

Andrew J. Benson 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-phastro-ph.IMhep-phastro-ph.HE

Catalog footprint

What is connected

47works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

47 published item(s)

preprint2022arXiv

A universal model for the evolution of tidally stripped systems

Accurate models of the structural evolution of dark matter subhaloes, as they orbit within larger systems, are fundamental to understanding the detailed distribution of dark matter at the present day. Numerical simulations of subhalo evolution support the idea that the mass loss associated with tidal stripping is most naturally understood in energy space, with the particles that are the least bound being removed first. Starting from this premise, we recently proposed a zero-parameter "energy-truncation model" for subhalo evolution. We tested this model with simulations of tidal stripping of satellites with initial NFW profiles, and showed that the energy-truncation model accurately predicts both the mass loss and density profiles. In this work, we apply the model to a variety of Hernquist, Einasto and King profiles. We show that it matches the simulation results quite closely in all cases, indicating that it may serve as a universal model to describe tidally stripped collisionless systems. A key prediction of the energy-truncation model is that the central density of dark matter subhaloes is conserved as they lose mass; this has important implications for dark matter annihilation calculations, and for other observational tests of dark matter.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Validating Synthetic Galaxy Catalogs for Dark Energy Science in the LSST Era

Large simulation efforts are required to provide synthetic galaxy catalogs for ongoing and upcoming cosmology surveys. These extragalactic catalogs are being used for many diverse purposes covering a wide range of scientific topics. In order to be useful, they must offer realistically complex information about the galaxies they contain. Hence, it is critical to implement a rigorous validation procedure that ensures that the simulated galaxy properties faithfully capture observations and delivers an assessment of the level of realism attained by the catalog. We present here a suite of validation tests that have been developed by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). We discuss how the inclusion of each test is driven by the scientific targets for static ground-based dark energy science and by the availability of suitable validation data. The validation criteria that are used to assess the performance of a catalog are flexible and depend on the science goals. We illustrate the utility of this suite by showing examples for the validation of cosmoDC2, the extragalactic catalog recently released for the LSST DESC second Data Challenge.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Random Walk Model for Dark Matter Halo Concentrations

For idealized (spherical, smooth) dark matter halos described by single-parameter density profiles (such as the NFW profile) there exists a one-to-one mapping between the energy of the halo and the scale radius of its density profile. The energy therefore uniquely determines the concentration parameter of such halos. We exploit this fact to predict the concentrations of dark matter halos via a random walk in halo energy space. Given a full merger tree for a halo, the total internal energy of each halo in that tree is determined by summing the internal and orbital energies of progenitor halos. We show that, when calibrated, this model can accurately reproduce the mean of the concentration--mass relation measured in N-body simulations, and reproduces more of the scatter in that relation than previous models. We further test this model by examining both the autocorrelation of scale radii across time, and the correlations between halo concentration and spin, and comparing to results measured from cosmological N-body simulations. In both cases we find that our model closely matches the N-body results. Our model is implemented within the open source Galacticus toolkit.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Probing Hot Gas Components of Circumgalactic Medium in Cosmological Simulations with the Thermal Sunyaev-Zel'dovich Effect

The thermal Sunyaev-Zel'dovich (tSZ) effect is a powerful tool with the potential for constraining directly the properties of the hot gas that dominates dark matter halos because it measures pressure and thus thermal energy density. Studying this hot component of the circumgalactic medium (CGM) is important because it is strongly impacted by star-formation and active galactic nuclei (AGN) activity in galaxies, participating in the feedback loop that regulates star and black hole mass growth in galaxies. We study the tSZ effect across a wide halo mass range using three cosmological hydrodynamical simulations: Illustris-TNG, EAGLE, and FIRE-2. Specifically, we present the scaling relation between tSZ signal and halo mass and radial profiles of gas density, temperature, and pressure for all three simulations. The analysis includes comparisons to Planck tSZ observations and to the thermal pressure profile inferred from the Atacama Cosmology Telescope (ACT) measurements. We compare these tSZ data to the simulations to interpret the measurements in terms of feedback and accretion processes in the CGM. We also identify as-yet unobserved potential signatures of these processes that may be visible in future measurements, which will have the capability of measuring tSZ signals to even lower masses. We also perform internal comparisons between runs with different physical assumptions. We conclude: (1) there is strong evidence for the impact of feedback at $R_{500}$ but that this impact decreases by $5R_{500}$, and (2) the thermodynamic profiles of the CGM are highly dependent on the implemented model, such as cosmic-ray or AGN feedback prescriptions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gravity and the Nonlinear Growth of Structure in the Carnegie-Spitzer-IMACS Redshift Survey

A key obstacle to developing a satisfying theory of galaxy evolution is the difficulty in extending analytic descriptions of early structure formation into full nonlinearity, the regime in which galaxy growth occurs. Extant techniques, though powerful, are based on approximate numerical methods whose Monte Carlo-like nature hinders intuition building. Here, we develop a new solution to this problem and its empirical validation. We first derive closed-form analytic expectations for the evolution of fixed percentiles in the real-space cosmic density distribution, {\it averaged over representative volumes observers can track cross-sectionally\}. Using the Lagrangian forms of the fluid equations, we show that percentiles in $δ$---the density relative to the median---should grow as $δ(t)\proptoδ_{0}^α\,t^β$, where $α\equiv2$ and $β\equiv2$ for Newtonian gravity at epochs after the overdensities transitioned to nonlinear growth. We then use 9.5 sq. deg. of Carnegie-Spitzer-IMACS Redshift Survey data to map {\it galaxy\} environmental densities over $0.2<z<1.5$ ($\sim$7 Gyr) and infer $α=1.98\pm0.04$ and $β=2.01\pm0.11$---consistent with our analytic prediction. These findings---enabled by swapping the Eulerian domain of most work on density growth for a Lagrangian approach to real-space volumetric averages---provide some of the strongest evidence that a lognormal distribution of early density fluctuations indeed decoupled from cosmic expansion to grow through gravitational accretion. They also comprise the first exact, analytic description of the nonlinear growth of structure extensible to (arbitrarily) low redshift. We hope these results open the door to new modeling of, and insight-building into, the galaxy growth and its diversity in cosmological contexts.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

LATIS: The Ly$α$ Tomography IMACS Survey

We introduce LATIS, the Ly$α$ Tomography IMACS Survey, a spectroscopic survey at Magellan designed to map the z=2.2-2.8 intergalactic medium (IGM) in three dimensions by observing the Ly$α$ forest in the spectra of galaxies and QSOs. Within an area of 1.7 deg${}^2$, we will observe approximately half of $\gtrsim L^*$ galaxies at z=2.2-3.2 for typically 12 hours, providing a dense network of sightlines piercing the IGM with an average transverse separation of 2.5 $h^{-1}$ comoving Mpc (1 physical Mpc). At these scales, the opacity of the IGM is expected to be closely related to the dark matter density, and LATIS will therefore map the density field in the $z \sim 2.5$ universe at $\sim$Mpc resolution over the largest volume to date. Ultimately LATIS will produce approximately 3800 spectra of z=2.2-3.2 galaxies that probe the IGM within a volume of $4 \times 10^6 h^{-3}$ Mpc${}^3$, large enough to contain a representative sample of structures from protoclusters to large voids. Observations are already complete over one-third of the survey area. In this paper, we describe the survey design and execution. We present the largest IGM tomographic maps at comparable resolution yet made. We show that the recovered matter overdensities are broadly consistent with cosmological expectations based on realistic mock surveys, that they correspond to galaxy overdensities, and that we can recover structures identified using other tracers. LATIS is conducted in Canada-France-Hawaii Telescope Legacy Survey fields, including COSMOS. Coupling the LATIS tomographic maps with the rich data sets collected in these fields will enable novel studies of environment-dependent galaxy evolution and the galaxy-IGM connection at cosmic noon.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mass loss in tidally stripped systems; the energy-based truncation method

The ability to accurately predict the evolution of tidally stripped haloes is important for understanding galaxy formation and testing the properties of dark matter. Most studies of substructure evolution make predictions based on empirical models of tidal mass loss that are calibrated using numerical simulations. This approach can be accurate in the cases considered, but lacks generality and does not provide a physical understanding of the processes involved. Recently, we demonstrated that truncating NFW distribution functions sharply in energy results in density profiles that resemble those of tidally stripped systems, offering a path to constructing physically motivated models of tidal mass loss. In this work, we review calculations of mass loss based on energy truncation alone, and then consider what secondary effects may modulate mass loss beyond this. We find that a combination of dependence on additional orbital parameters and variations in individual particle energies over an orbit results in a less abrupt truncation in energy space as a subhalo loses mass. Combining the energy truncation approach with a simple prediction for the mass-loss rate, we construct a full model of mass loss that can accurately predict the evolution of a subhalo in terms of a single parameter $η_{\rm eff}$. This parameter can be fully determined from the initial orbital and halo properties, and does not require calibration with numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Normalization and Slope of the Dark Matter (Sub-)Halo Mass Function on Sub-Galactic Scales

Simulations of cold dark matter make robust predictions about the slope and normalization of the dark matter halo and subhalo mass functions on small scales. Recent observational advances utilizing strong gravitational lensing have demonstrated the ability of this technique to place constraints on these quantities on subgalactic scales corresponding to dark matter halo masses of $10^6$--$10^9\mathrm{M}_\odot$. On these scales the physics of baryons, which make up around 17% of the matter content of the Universe but which are not included in pure dark matter N-body simulations, are expected to affect the growth of structure and the collapse of dark matter halos. In this work we develop a semi-analytic model to predict the amplitude and slope of the dark matter halo and subhalo mass functions on subgalactic scales in the presence of baryons. We find that the halo mass function is suppressed by up to 25%, and the slope is modified, ranging from $-1.916$ to $-1.868$ in this mass range. These results are consistent with current measurements, but differ sufficiently from the expectations for a dark matter only universe that it may be testable in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A unified multi-wavelength model of galaxy formation

We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy disks based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850 mu. This was not previously achieved by a single physical model within the LambdaCDM framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z=0 to z=6, and is found to be generally successful. These observations include the optical and near-IR luminosity functions, HI mass function, fraction of early type galaxies, Tully-Fisher, metallicity-luminosity and size-luminosity relations at z=0, as well as far-IR number counts, and far-UV luminosity functions at z ~ 3-6. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Achieving Convergence in Galaxy Formation Models by Augmenting N-body Merger Trees

Accurate modeling of galaxy formation in a hierarchical, cold dark matter universe requires the use of sufficiently high-resolution merger trees to obtain convergence in the predicted properties of galaxies. When semi-analytic galaxy formation models are applied to cosmological N-body simulation merger trees, it is often the case that those trees have insufficient resolution to give converged galaxy properties. We demonstrate a method to augment the resolution of N-body merger trees by grafting in branches of Monte Carlo merger trees with higher resolution, but which are consistent with the pre-existing branches in the N-body tree. We show that this approach leads to converged galaxy properties.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Origin and Evolution of the Galaxy Stellar Mass Function

Here we explore the evolution of galaxy ensembles at early times by writing the in situ stellar mass growth of galaxies purely as a stationary stochastic (e.g., quasi-steady state) process. By combining the mathematics of such processes with Newtonian gravity and a mean local star formation efficiency, we show that the stellar mass evolution of galaxy ensembles is directly related to the average acceleration of baryons onto dark matter halos at the onset of star formation, with explicit dependencies on initial local matter densities and halo mass. The density term specifically implies more rapid average rates of growth in higher density regions of the universe compared to low density regions, i.e., assembly bias. With this framework, using standard cosmological parameters, a mean star formation efficiency derived by other authors, and knowledge of the shape of the cosmological matter power spectrum at small scales, we analytically derive (1) the characteristic stellar masses of galaxies (M*), (2) the power-law low-mass slope (alpha) and normalization (phi*) of the stellar mass function, and (3) the evolution of the stellar mass function in time over 12.5 > z > 2. Correspondingly, the rise in the cosmic star formation rate density over these epochs, while the universe can sustain unabated fueling of star formation, also emerges naturally. All of our findings are consistent with the deepest available data, including the expectation of alpha~-7/5; i.e., a stellar mass function low-mass slope that is notably shallower than that of the halo mass function, and with no systematic deviations from a mean star formation efficiency with density or mass, nor any explicit, additional feedback mechanisms. These derivations yield a compelling richness and complexity but also show that very few astrophysical details are required to understand the evolution of cosmic ensemble of galaxies at early times.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The scaling relations of early-type dwarf galaxies across a range of environments

We present the results of a Keck-ESI study of dwarf galaxies across a range of environment: the Perseus Cluster, the Virgo Cluster, the NGC 1407 group, and the NGC 1023 group. Eighteen dEs are targeted for spectroscopy, three for the first time. We confirm cluster membership for one Virgo dE, and group membership for one dE in the NGC 1023 group, and one dE in the NGC 1407 group for the first time. Regardless of environment, the dEs follow the same size-magnitude and $σ$-luminosity relation. Two of the Virgo dwarfs, VCC 1199 and VCC 1627, have among the highest central velocity dispersions ($σ_{0}$ = 58.4 km s$^{-1}$ and 49.2 km s$^{-1}$) measured for dwarfs of their luminosity ($M_{R}\approx -17$). Given their small sizes ($R_{e} < 300$ pc) and large central velocity dispersions, we classify these two dwarfs as compact ellipticals rather than dEs. Group dEs typically have higher mean dynamical-to-stellar mass ratios than the cluster dEs, with $M_{dyn}/M_{\star} = 5.1\pm0.6$ for the group dwarfs, vs. $M_{dyn}/M_{\star} = 2.2\pm0.5$ for the cluster sample, which includes two cEs. We also search for trends in $M_{dyn}/M_{\star}$ vs. distance from M87 for the Virgo Cluster population, and find no preference for dwarfs with high values of $M_{dyn}/M_{\star}$ to reside in the cluster outskirts vs. centre.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Galaxy Formation as a Cosmological Tool. I: The Galaxy Merger History as a Measure of Cosmological Parameters

As galaxy formation and evolution over long cosmic time-scales depends to a large degree on the structure of the universe, the assembly history of galaxies is potentially a powerful approach for learning about the universe itself. In this paper we examine the merger history of dark matter halos based on the Extended Press-Schechter formalism as a function of cosmological parameters, redshift and halo mass. We calculate how major halo mergers are influenced by changes in the cosmological values of $Ω_{\rm m}$, $Ω_Λ$, $σ_{8}$, the dark matter particle temperature (warm vs. cold dark matter), and the value of a constant and evolving equation of state parameter $w(z)$. We find that the merger fraction at a given halo mass varies by up to a factor of three for halos forming under the assumption of Cold Dark Matter, within different underling cosmological parameters. We find that the current measurements of the merger history, as measured through observed galaxy pairs as well as through structure, are in agreement with the concordance cosmology with the current best fit giving $1 - Ω_{\rm m} = Ω_{\rm Λ} = 0.84^{+0.16}_{-0.17}$. To obtain a more accurate constraint competitive with recently measured cosmological parameters from Planck and WMAP requires a measured merger accuracy of $δf_{\rm m} \sim 0.01$, implying surveys with an accurately measured merger history over 2 - 20 deg$^{2}$, which will be feasible with the next generation of imaging and spectroscopic surveys such as Euclid and LSST.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear evolution of dark matter subhalos and applications to warm dark matter

We describe the methodology to include nonlinear evolution, including tidal effects, in the computation of subhalo distribution properties in both cold (CDM) and warm (WDM) dark matter universes. Using semi-analytic modeling, we include effects from dynamical friction, tidal stripping, and tidal heating, allowing us to dynamically evolve the subhalo distribution. We calibrate our nonlinear evolution scheme to the CDM subhalo mass function in the Aquarius N-body simulation, producing a subhalo mass function within the range of simulations. We find tidal effects to be the dominant mechanism of nonlinear evolution in the subhalo population. Finally, we compute the subhalo mass function for $m_χ=1.5$ keV WDM including the effects of nonlinear evolution, and compare radial number densities and mass density profiles of subhalos in CDM and WDM models. We show that all three signatures differ between the two dark matter models, suggesting that probes of substructure may be able to differentiate between them.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Circular Velocity Function of Group Galaxies

A robust prediction of $Λ{\rm CDM}$ cosmology is the halo circular velocity function (CVF), a dynamical cousin of the halo mass function. The correspondence between theoretical and observed CVFs is uncertain, however: cluster galaxies are reported to exhibit a power-law CVF consistent with $N$-body simulations, but that of the field is distinctly Schechter-like, flattened relative to $Λ{\rm CDM}$ expectations at circular velocities $v_c \lesssim 200\, {\rm km\, s^{-1}}$. Groups offer a powerful probe of the role environment plays in this discrepancy as they bridge the field and clusters. Here, we construct the CVF for a large, mass- and multiplicity-complete sample of group galaxies from the Sloan Digital Sky Survey. Using independent photometric $v_c$ estimators, we find no transition from a field- to $Λ{\rm CDM}$-shaped CVF above $v_c = 50\, {\rm km\, s^{-1}}$ as a function of group halo mass. All groups with $12.4 \lesssim \log M_{\rm halo} / M_{\odot} \lesssim 15.1$ (Local Group analogs to rich clusters) display similar Schechter-like CVFs marginally suppressed at low-$v_c$ compared to that of the field. Conversely, some agreement with $N$-body results emerges for samples saturated with late-type galaxies, with isolated late-types displaying a CVF similar in shape to $Λ{\rm CDM}$ predictions. We conclude that the flattening of the low-$v_c$ slope in groups is due to their depressed late-type fractions -- environment affecting the CVF only to the extent that it correlates with this quantity -- and that previous cluster analyses may suffer from interloper contamination. These results serve as useful benchmarks for cosmological simulations of galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Research Note on the Implementation of Star Formation and Stellar Feedback in Semi-Analytic Models

We study the impact of star formation and stellar feedback prescriptions on galaxy properties predicted by means of "stripped-down" versions of independently developed semi-analytic models (SAMs). These include cooling, star formation, feedback from supernovae (SNe) and simplified prescriptions for galaxy merging, but no chemical evolution, disc instabilities or AGN feedback. We run these versions on identical samples of dark matter (DM) haloes extracted from high-resolution N-body simulations in order to perform both statistical analysis and object-by-object comparisons. We compare our results with previous work based on stripped-down versions of the same SAMs including only gas cooling, and show that all feedback models provide coherent modifications in the distribution of baryons between the various gas phases. In particular, we find that the predicted hot gas fractions are considerably increased by up to a factor of three, while the corresponding cold gas fractions are correspondingly decreased, and a significant amount of mass is ejected from the DM halo. Nonetheless, we also find relevant differences in the predicted properties of model galaxies among the three SAMs: these deviations are more relevant at mass scales comparable to that of our own Galaxy, and are reduced at larger masses, confirming the varying impact of stellar feedback at different mass scales. We also check the effect of enhanced star formation events (i.e. starbursts modes), defined in connection with galaxy mergers. We find that, in general, these episodes have a limited impact in the overall star formation histories of model galaxies, even in massive DM halos where merger-driven star formation has often been considered very important.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Excursion Set Theory for Correlated Random Walks

We present a new method to compute the first crossing distribution in excursion set theory for the case of correlated random walks. We use a combination of the path integral formalism of Maggiore & Riotto, and the integral equation solution of Zhang & Hui, and Benson et al. to find a numerically robust and convenient algorithm to derive the first crossing distribution in terms of a perturbative expansion around the limit of an uncorrelated random walk. We apply this methodology to the specific case of a Gaussian random density field filtered with a Gaussian smoothing function. By comparing our solutions to results from Monte Carlo calculations of the first crossing distribution we demonstrate that our method accurate for power spectra $P(k)\propto k^n$ for $n=1$, becoming less accurate for smaller values of $n$. It is therefore complementary to the method of Musso & Sheth, which will therefore be more useful for standard $Λ$CDM power spectra. Our approach is quite general, and can be adapted to other smoothing functions, and also to non-Gaussian density fields.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lightcone mock catalogues from semi-analytic models of galaxy formation - I. Construction and application to the BzK colour selection

We introduce a method for constructing end-to-end mock galaxy catalogues using a semi-analytical model of galaxy formation, applied to the halo merger trees extracted from a cosmological N-body simulation. The mocks that we construct are lightcone catalogues, in which a galaxy is placed according to the epoch at which it first enters the past lightcone of the observer, and incorporate the evolution of galaxy properties with cosmic time. We determine the position between the snapshot outputs at which a galaxy enters the observer's lightcone by interpolation. As an application, we consider the effectiveness of the BzK colour selection technique, which was designed to isolate galaxies in the redshift interval 1.4<z<2.5. The mock catalogue is in reasonable agreement with the observed number counts of all BzK galaxies, as well as with the observed counts of the subsample of BzKs that are star-forming galaxies. We predict that over 75 per cent of the model galaxies with K_{AB}<=23, and 1.4<z<2.5, are selected by the BzK technique. Interloper galaxies, outside the intended redshift range, are predicted to dominate bright samples of BzK galaxies (i.e. with K_{AB}<=21). Fainter K-band cuts are necessary to reduce the predicted interloper fraction. We also show that shallow B-band photometry can lead to confusion in classifying BzK galaxies as being star-forming or passively evolving. Overall, we conclude that the BzK colour selection technique is capable of providing a sample of galaxies that is representative of the 1.4<z<2.5 galaxy population.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dark Matter Halo Merger Histories Beyond Cold Dark Matter: I - Methods and Application to Warm Dark Matter

We describe a methodology to accurately compute halo mass functions, progenitor mass functions, merger rates and merger trees in non-cold dark matter universes using a self-consistent treatment of the generalized extended Press-Schechter formalism. Our approach permits rapid exploration of the subhalo population of galactic halos in dark matter models with a variety of different particle properties or universes with rolling, truncated, or more complicated power spectra. We make detailed comparisons of analytically derived mass functions and merger histories with recent warm dark matter cosmological N-body simulations, and find excellent agreement. We show that, once the accretion of smoothly distributed matter is accounted for, coarse-grained statistics such as the mass accretion history of halos can be almost indistinguishable between cold and warm dark matter cases. However, the halo mass function and progenitor mass functions differ significantly, with the warm dark matter cases being strongly suppressed below the free-streaming scale of the dark matter. We demonstrate the importance of using the correct solution for the excursion set barrier first-crossing distribution in warm dark matter - if the solution for a flat barrier is used instead the truncation of the halo mass function is much slower, leading to an overestimate of the number of low mass halos.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Threshing in Action - The tidal disruption of a dwarf galaxy by the Hydra I Cluster

We report on the discovery of strong tidal features around a dwarf spheroidal galaxy in the Hydra I galaxy cluster, indicating its ongoing tidal disruption. This very low surface brightness object, HCC-087, was originally classified as an early-type dwarf in the Hydra Cluster Catalogue (HCC), but our re-analysis of the ESO-VLT/FORS images of the HCC unearthed a clear indication of an S-shaped morphology and a large spatial extent. Its shape, luminosity (M_V=-11.6 mag), and physical size (at a half-light radius of 3.1 kpc and a full length of ~5.9 kpc) are comparable to the recently discovered NGC 4449B and the Sagittarius dwarf spheroidal, all of which are undergoing clear tidal disruption. Aided by N-body simulations we argue that HCC-087 is currently at its first apocenter, at 150 kpc, around the cluster center and that it is being tidally disrupted by the galaxy cluster's potential itself. An interaction with the near-by (50 kpc) S0 cluster galaxy HCC-005, at M* ~ 3 x 10^10 M_sun is rather unlikely, as this constellation requires a significant amount of dynamical friction and thus low relative velocities. The S-shaped morphology and large spatial extent of the satellite would, however, also appear if HCC-087 would orbit the cluster center. These features appear to be characteristic properties of satellites that are seen in the process of being tidally disrupted, independent of the environment of the destruction. An important finding of our simulations is an orientation of the tidal tails perpendicular to the orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Convergence of Galaxy Properties with Merger Tree Temporal Resolution

Dark matter halo merger trees are now routinely extracted from cosmological simulations of structure formation. These trees are frequently used as inputs to semi-analytic models of galaxy formation to provide the backbone within which galaxy formation takes place. By necessity, these merger trees are constructed from a finite set of discrete "snapshots" of the N-body simulation and so have a limited temporal resolution. To date, there has been little consideration of how this temporal resolution affects the properties of galaxies formed within these trees. In particular, the question of how many snapshots are needed to achieve convergence in galaxy properties has not be answered. Therefore, we study the convergence in the stellar and total baryonic masses of galaxies, distribution of merger times, stellar mass functions and star formation rates in the Galacticus model of galaxy formation as a function of the number of "snapshot" times used to represent dark matter halo merger trees. When utilizing snapshots between z=20 and z=0, we find that at least 128 snapshots are required to achieve convergence to within 5% for galaxy masses. This convergence is obtained for mean quantities averaged over large samples of galaxies - significant variance for individual galaxies remains even when using very large numbers of snapshots. We find only weak dependence of the rate of convergence on the distribution of snapshots in time - snapshots spaced uniformly in the expansion factor, uniformly in the logarithm of expansion factor or uniformly in the logarithm of critical overdensity for collapse work equally well in almost all cases. We provide input parameters to Galacticus which allow this type of convergence study to be tuned to other simulations and to be carried out for other galaxy properties.

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

On the impact of empirical and theoretical star formation laws on galaxy formation

We investigate the consequences of applying different star formation laws in the galaxy formation model GALFORM. Three broad star formation laws are implemented: the empirical relations of Kennicutt and Schmidt and Blitz & Rosolowsky and the theoretical model of Krumholz, McKee & Tumlinson. These laws have no free parameters once calibrated against observations of the star formation rate (SFR) and gas surface density in nearby galaxies. We start from published models, and investigate which observables are sensitive to a change in the star formation law, without altering any other model parameters. We show that changing the star formation law (i) does not significantly affect either the star formation history of the universe or the galaxy luminosity functions in the optical and near-IR, due to an effective balance between the quiescent and burst star formation modes; (ii) greatly affects the cold gas contents of galaxies; (iii) changes the location of galaxies in the SFR versus stellar mass plane, so that a second sequence of "passive" galaxies arises, in addition to the known "active" sequence. We show that this plane can be used to discriminate between the star formation laws.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Fossil Record of Two-phase Galaxy Assembly: Kinematics and Metallicities in the Nearest S0 Galaxy

We present a global analysis of kinematics and metallicity in the nearest S0 galaxy, NGC 3115, along with implications for its assembly history. The data include high-quality wide-field imaging from Suprime-Cam on the Subaru telescope, and multi-slit spectra of the field stars and globular clusters (GCs) obtained using Keck-DEIMOS/LRIS and Magellan-IMACS. Within two effective radii, the bulge (as traced by the stars and metal-rich GCs) is flattened and rotates rapidly (v/sigma > 1.5). At larger radii, the rotation declines dramatically to v/sigma ~ 0.7, but remains well-aligned with the inner regions. The radial decrease in characteristic metallicity of both the metal-rich and metal-poor GC subpopulations produces strong gradients with power law slopes of -0.17 +/- 0.04 and -0.38 +/- 0.06 dex per dex, respectively. We argue that this pattern is not naturally explained by a binary major merger, but instead by a two-phase assembly process where the inner regions have formed in an early violent, dissipative phase, followed by the protracted growth of the outer parts via minor mergers with typical mass ratios of ~ 15-20:1.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The population of Milky Way satellites in the LambdaCDM cosmology

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

X-rays and hard UV radiation From the First Galaxies: Ionization Bubbles and 21 cm Observations

The first stars and quasars are known sources of hard ionizing radiation in the first billion years of the Universe. We examine the joint effects of X-rays and hard UV radiation from such first-light sources on the hydrogen and helium reionization of the intergalactic medium (IGM) at early times, and the associated heating. We study the growth and evolution of individual HII, HeII and HeIII regions around early galaxies with first stars and/or QSO populations. We find that in the presence of helium-ionizing radiation, X-rays may not dominate the ionization and thermal history of the IGM at redshifts, z, of 10-20, contributing relatively modest increases to IGM ionization, and heating up to about 10^3--10^5 K in IGM temperatures. We also calculate the 21 cm signal expected from a number of scenarios with metal-free starbursts and quasars at these redshifts. The peak values for the spin temperature reach about 10^4 to 10^5 K in such cases. The maximum values for the 21 cm brightness temperature are around 30-40 mK in emission, while the net values of the 21 cm absorption signal range from about a few to 60 mK on scales of 0.01-1 Mpc. We find that the 21 cm signature of X-ray versus UV ionization could be distinct, with the emission signal expected from X-rays alone occurring at smaller scales than that from UV radiation, resulting from the inherently different spatial scales at which X-ray and UV ionization/heating manifest. This difference is time-dependent, and becomes harder to distinguish with an increasing X-ray contribution to the total ionizing photon production. Such differing scale-dependent contributions from X-ray and UV photons may therefore "blur" the 21 cm signature of the percolation of ionized bubbles around early halos (depending on whether a cosmic X-ray or UV background built up first), and affect the interpretation of 21 cm data constraints on reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A semi-analytic model comparison - gas cooling and galaxy mergers

We use stripped-down versions of three semi-analytic galaxy formation models to study the influence of different assumptions about gas cooling and galaxy mergers. By running the three models on identical sets of merger trees extracted from high-resolution cosmological N-body simulations, we are able to perform both statistical analyses and halo-by-halo comparisons. Our study demonstrates that there is a good statistical agreement between the three models used here, when operating on the same merger trees, reflecting a general agreement in the underlying framework for semi-analytic models. We also show, however, that various assumptions that are commonly adopted to treat gas cooling and galaxy mergers can lead to significantly different results, at least in some regimes. In particular, we find that the different models adopted for gas cooling lead to similar results for mass scales comparable to that of our own Galaxy. Significant differences, however, arise at larger mass scales. These are largely (but not entirely) due to different treatments of the `rapid cooling' regime, and different assumptions about the hot gas distribution. At this mass regime, the predicted cooling rates can differ up to about one order of magnitude, with important implications on the relative weight that these models give to AGN feedback in order to counter-act excessive gas condensation in relatively massive haloes at low redshift. Different assumptions in the modelling of galaxy mergers can also result in significant differences in the timings of mergers, with important consequences for the formation and evolution of massive galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Angular Momentum Evolution in Dark Matter Halos

We have analyzed high resolution N-body simulations of dark matter halos, focusing specifically on the evolution of angular momentum. We find that not only is individual particle angular momentum not conserved, but the angular momentum of radial shells also varies over the age of the Universe by up to factors of a few. We find that torques from external structure are the most likely cause for this distribution shift. Since the model of adiabatic contraction that is often applied to model the effects of galaxy evolution on the dark-matter density profile in a halo assumes angular momentum conservation, this variation implies that there is a fundamental limit on the possible accuracy of the adiabatic contraction model in modeling the response of DM halos to the growth of galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cold Mode Accretion in Galaxy Formation

A generic expectation for gas accreted by high mass haloes is that it is shock heated to the virial temperature of the halo. In low mass haloes, or at high redshift, however, the gas cooling rate is sufficiently rapid that an accretion shock is unlikely to form. Instead, gas can accrete directly into the centre of the halo in a `cold mode' of accretion. Although semi-analytic models have always made a clear distinction between hydrostatic and rapid cooling they have not made a distinction between whether or not an accretion shock forms. Starting from the well-established Galform code, we investigate the effect of explicitly accounting for cold mode accretion using the shock stability model of Birnboim & Dekel. When we modify the code so that there is no effective feedback from galaxy formation, we find that cold mode accretion is the dominant channel for feeding gas into the galaxies at high redshifts. However, this does not translate into a significant difference in the star formation history of the universe compared to the previous code. When effective feedback is included in the model, we find that the the cold mode is much less apparent because of the presence of gas ejected from the galaxy. Thus the inclusion of the additional cold mode physics makes little difference to basic results from earlier semi-analytic models which used a simpler treatment of gas accretion. For more sophisticated predictions of its consequences, we require a better understanding of how the cold mode delivers angular momentum to galaxies and how it interacts with outflows.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraints on Decaying Dark Matter

We explore a dark-matter model in which there are two dark-matter species nearly degenerate in mass, with epsilon = Delta M/M << 1. The heavier particle undergoes two-body decay with a half-life tau, to the lighter dark-matter particle and a noninteracting massless particle. Unlike previous work on decaying dark matter, we explore the regime tau > 100 Myr and non-relativistic kick speeds vk / c = epsilon. Using a set of N-body simulations of isolated dark-matter halos, we show how halos change as a function of tau and vk. We find that tau < 40 Gyr is ruled out for vk > 20 km s^{-1} (epsilon > 10^{-4}) when we compare the simulations to observations of dwarf-galaxy- to cluster-mass dark matter halos. We highlight which set of observations should provide better future constraints for decays and other types of dark-matter physics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dark-matter decays and Milky Way satellite galaxies

We consider constraints on a phenomenological dark-matter model consisting of two nearly degenerate particle species using observed properties of the Milky Way satellite galaxy population. The two parameters of this model, assuming the particle masses are >~ GeV, are v_k, the recoil speed of the daughter particle, and tau, the lifetime of the parent particle. The satellite constraint that spans the widest range of v_k is the number of satellites that have a mass within 300 pc M300 > 5 x 10^6 solar masses, although constraints based on M300 in the classical dwarfs and the overall velocity function are competitive for v_k >~ 50 km/s. In general, we find that tau <~ 30 Gyr is ruled out for 20 km/s <~ v_k <~ 200 km/s, although we find that the limits on tau for fixed v_k can change constraints by a factor of ~3 depending on the star-formation histories of the satellites. We advocate using the distribution of M300 in Milky Way satellites determined by next-generation all-sky surveys and follow-up spectroscopy as a probe of dark-matter properties.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Galacticus: A Semi-Analytic Model of Galaxy Formation

We describe a new, free and open source semi-analytic model of galaxy formation, Galacticus. The Galacticus model was designed to be highly modular to facilitate expansion and the exploration of alternative descriptions of key physical ingredients. We detail the Galacticus engine for evolving galaxies through a merging hierarchy of dark matter halos and give details of the specific implementations of physics currently available in Galacticus. Finally, we show results from an example model that is in reasonably good agreement with several observational datasets. We use this model to explore numerical convergence and to demonstrate the types of information which can be extracted from Galacticus.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Galaxy Formation Spanning Cosmic History

Over the past several decades, galaxy formation theory has met with significant successes. In order to test current theories thoroughly we require predictions for as yet unprobed regimes. To this end, we describe a new implementation of the Galform semi-analytic model of galaxy formation. Our motivation is the success of the model described by Bower et al. in explaining many aspects of galaxy formation. Despite this success, the Bower et al. model fails to match some observational constraints and certain aspects of its physical implementation are not as realistic as we would like. The model described in this work includes substantially updated physics, taking into account developments in our understanding over the past decade, and removes certain limiting assumptions made by this (and most other) semi-analytic models. This allows it to be exploited reliably in high-redshift and low mass regimes. Furthermore, we have performed an exhaustive search of model parameter space to find a particular set of model parameters which produce results in good agreement with a wide range of observational data (luminosity functions, galaxy sizes and dynamics, clustering, colours, metal content) over a wide range of redshifts. This model represents a solid basis on which to perform calculations of galaxy formation in as yet unprobed regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Galaxy Formation Theory

We review the current theory of how galaxies form within the cosmological framework provided by the cold dark matter paradigm for structure formation. Beginning with the pre-galactic evolution of baryonic material we describe the analytical and numerical understanding of how baryons condense into galaxies, what determines the structure of those galaxies and how internal and external processes (including star formation, merging, active galactic nuclei etc.) determine their gross properties and evolution. Throughout, we highlight successes and failings of current galaxy formation theory. We include a review of computational implementations of galaxy formation theory and assess their ability to provide reliable modelling of this complex phenomenon. We finish with a discussion of several "hot topics" in contemporary galaxy formation theory and assess future directions for this field.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The impact of dark matter cusps and cores on the satellite galaxy population around spiral galaxies

(Abridged) We use N-body simulations to study the effects that a divergent (i.e. "cuspy") dark matter (DM) profile introduces on the tidal evolution of dwarf spheroidal galaxies (dSphs). Our models assume cosmologically-motivated initial conditions where dSphs are DM-dominated systems on eccentric orbits about a host galaxy composed of a dark halo and a baryonic disc. We find that the resilience of dSphs to tidal stripping is extremely sensitive to the halo cuspiness; whereas dwarfs with a cored profile can be easily destroyed by the host disc, those with cusps always retain a bound remnant. For a given halo profile the evolution of the structural parameters as driven by tides is controlled solely by the total amount of mass lost. This information is used to construct a semi-analytic code that simulates the hierarchical build-up of spiral galaxies assuming different halo profiles and disc masses. We find that tidal encounters with discs tend to decrease the average mass of satellites at all galactocentric radii. Interestingly, satellites accreted before re-ionization (z>6), which may be singled out by anomalous metallicity patterns, survive only if haloes are cuspy. We show that the size-mass relation established from Milky Way (MW) dwarfs strongly supports the presence of cusps in the majority of these systems, as cored models systematically underestimate the masses of the known Ultra-Faint dSphs. Our models also indicate that a massive M31 disc may explain why many of its dSphs fall below the size-mass relationship derived from MW dSphs. We use our models to constrain the mass threshold below which star formation is suppressed in DM haloes, finding that luminous satellites must be accreted with masses above 10^8--10^9 M_sol in order to explain the size-mass relation observed in MW dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The mass assembly of galaxy groups and the evolution of the magnitude gap

We investigate the assembly of groups and clusters of galaxies using the Millennium dark matter simulation and the associated gas simulations and semi-analytic catalogues of galaxies. In particular, in order to find an observable quantity that could be used to identify early-formed groups, we study the development of the difference in magnitude between their brightest galaxies to assess the use of magnitude gaps as possible indicators. We select galaxy groups and clusters at redshift z=1 with dark matter halo mass M(R200) > 1E13/h Msun, and trace their properties until the present time (z=0). We consider only the systems with X-ray luminosity L_X> 0.25E42/h^2 erg/s at z=0. While it is true that a large magnitude gap between the two brightest galaxies of a particular group often indicates that a large fraction of its mass was assembled at an early epoch, it is not a necessary condition. More than 90% of fossil groups defined on the basis of their magnitude gaps (at any epoch between 0<z<1) cease to be fossils within 4 Gyr, mostly because other massive galaxies are assembled within their cores, even though most of the mass in their haloes might have been assembled at early times. We show that, compared to the conventional definition of fossil galaxy groups based on the magnitude gap Delta m(12)> 2 (in the R-band, within 0.5R200 of the centre of the group), an alternative criterion Delta m(14)>2.5 (within the same radius) finds 50% more early-formed systems, and those that on average retain their fossil phase longer. However, the conventional criterion performs marginally better at finding early-formed groups at the high-mass end of groups. Nevertheless, both criteria fail to identify a majority of the early-formed systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Role of Ram Pressure Stripping in the Quenching of Cluster Star Formation

Recent observations of galaxy clusters have shown that environmental effects apparently associated with the cluster begin to lower the star formation rates of galaxies at distances as great as three times the cluster virial radius. These observations may indicate preprocessing of cluster galaxies in groups or in the cluster core for galaxies on highly elliptical orbits, but may also imply that the environmental effects due to the cluster are directly affecting galaxies on their first infall. To explore these issues, we investigate different models of ram pressure stripping as it acts on satellite galaxies in clusters, and compare to observations of the radial star formation gradient in clusters. We calculate the location of the accretion shock around model clusters, and use this as the radius of onset of ram pressure stripping in the GALFORM semi-analytic model of galaxy formation. Comparison of the results of our model, and previously considered, simpler ram pressure models, with recent observations indicates that current data is unable to strongly discriminate between models of ram pressure stripping due to the complex interplay of preprocessing effects at work. However, future observations of a larger sample of clusters will likely be able to place stronger constraints on the process of ram pressure stripping and its role in shaping radial trends in and around clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Maximum Spin of Black Holes Driving Jets

Unbounded outflows in the form of highly collimated jets and broad winds appear to be a ubiquitous feature of accreting black hole systems. The most powerful jets are thought to derive a significant fraction, if not the majority, of their power from the rotational energy of the black hole. Whatever the precise mechanism that causes them, these jets must therefore exert a braking torque on the black hole. We calculate the spin-up function for an accreting black hole, accounting for this braking torque. We find that the predicted black hole spin-up function depends only on the black hole spin and dimensionless parameters describing the accretion flow. Using recent relativistic magnetohydrodynamical numerical simulation results to calibrate the efficiency of angular momentum transfer in the flow, we find that an ADAF flow will spin a black hole up (or down) to an equilibrium value of about 96% of the maximal spin value in the absence of jets. Combining our ADAF system with a simple model for jet power, we demonstrate that an equilibrium is reached at approximately 93% of the maximal spin value, as found in the numerical simulation studies of the spin-up of accreting black holes, at which point the spin-up of the hole by accreted material is balanced by the braking torque arising from jet production. Our model also yields a relationship between jet efficiency and black hole spin that is in surprisingly good agreement with that seen in the simulation studies, indicating that our simple model is a useful and convenient description of ADAF inflow - jet outflow about a spinning black hole for incorporation in models of the formation and evolution of galaxies, groups and clusters of galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Origin of the Hubble Sequence in Lambda-CDM Cosmology

The Galform semi-analytic model of galaxy formation is used to explore the mechanisms primarily responsible for the three types of galaxies seen in the local universe: bulge, bulge+disk and disk, identified with the visual morphological types E, S0/a-Sbc, and Sc-Scd, respectively. With a suitable choice of parameters the Galform model can accurately reproduce the observed local K_s-band luminosity function (LF) for galaxies split by visual morphological type. The successful set of model parameters is used to populate the Millennium Simulation with 9.4 million galaxies and their dark matter halos. The resulting catalogue is then used to explore the evolution of galaxies through cosmic history. The model predictions concur with recent observational results including the galaxy merger rate, the star formation rate and the seemingly anti-hierarchical evolution of ellipticals. However, the model also predicts significant evolution of the elliptical galaxy LF that is not observed. The discrepancy raises the possibility that samples of z~1 galaxies which have been selected using colour and morphological criteria may be contaminated with galaxies that are not actually ellipticals.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Red Galaxy Growth and the Halo Occupation Distribution

We have traced the past 7 Gyr of red galaxy stellar mass growth within dark matter halos. We have determined the halo occupation distribution, which describes how galaxies reside within dark matter halos, using the observed luminosity function and clustering of 40,696 0.2<z<1.0 red galaxies in Bootes. Half of 10^{11.9} Msun/h halos host a red central galaxy, and this fraction increases with increasing halo mass. We do not observe any evolution of the relationship between red galaxy stellar mass and host halo mass, although we expect both galaxy stellar masses and halo masses to evolve over cosmic time. We find that the stellar mass contained within the red population has doubled since z=1, with the stellar mass within red satellite galaxies tripling over this redshift range. In cluster mass halos most of the stellar mass resides within satellite galaxies and the intra-cluster light, with a minority of the stellar mass residing within central galaxies. The stellar masses of the most luminous red central galaxies are proportional to halo mass to the power of a third. We thus conclude that halo mergers do not always lead to rapid growth of central galaxies. While very massive halos often double in mass over the past 7 Gyr, the stellar masses of their central galaxies typically grow by only 30%.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Colours of Satellite Galaxies in Groups and Clusters

Current models of galaxy formation predict satellite galaxies in groups and clusters that are redder than observed. We investigate the effect on the colours of satellite galaxies produced by the ram pressure stripping of their hot gaseous atmospheres as the satellites orbit within their parent halo. We incorporate a model of the stripping process based on detailed hydrodynamic simulations within the Durham semi-analytic model of galaxy formation. The simulations show that the environment in groups and clusters is less aggressive than previously assumed. The main uncertainty in the model is the treatment of gas expelled by supernovae. With reasonable assumptions for the stripping of this material, we find that satellite galaxies are able to retain a significant fraction of their hot gas for several Gigayears, thereby replenishing their reservoirs of cold, star forming gas and remaining blue for a relatively long period of time. A bimodal distribution of galaxy colours, similar to that observed in SDSS data, is established and the colours of the satellite galaxies are in good agreement with the data. In addition, our model naturally accounts for the observed dependence of satellite colours on environment, from small groups to high mass clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Evolving Luminosity Function of Red Galaxies

We trace the assembly history of red galaxies since z=1, by measuring their evolving space density with the B-band luminosity function. Our sample of 39599 red galaxies, selected from 6.96 square degrees of imaging from the NOAO Deep Wide-Field and Spitzer IRAC Shallow surveys, is an order of magnitude larger, in size and volume, than comparable samples in the literature. We measure a higher space density of z=0.9 red galaxies than some of the recent literature, in part because we account for the faint yet significant galaxy flux which falls outside of our photometric aperture. The B-band luminosity density of red galaxies, which effectively measures the evolution of ~L* galaxies, increases by only 36 percent from z=0 to z=1. If red galaxy stellar populations have faded by 1.24 B-band magnitudes since z=1, the stellar mass contained within the red galaxy population has roughly doubled over the past 8 Gyr. This is consistent with star-forming galaxies being transformed into ~L* red galaxies after a decline in their star formation rates. In contrast, the evolution of 4L* red galaxies differs only slightly from a model with negligible star formation and no galaxy mergers since z=1. If this model approximates the luminosity evolution of red galaxy stellar populations, then 80 percent of the stellar mass contained within today's 4L* red galaxies was already in place at z=0.7. While red galaxy mergers have been observed, such mergers do not produce rapid growth of 4L* red galaxy stellar masses between z=1 and the present day.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Orbital Parameters of Infalling Dark Matter Substructures

We present distributions of the orbital parameters of dark matter substructures at the time of merging into their host halo. Accurate knowledge of the orbits of dark matter substructures is a crucial input to studies which aim to assess the effects of the cluster environment on galaxies, the heating of galaxy disks and many other topics. Orbits are measured for satellites in a large number of N-body simulations. We focus on the distribution of radial and tangential velocities, but consider also distributions of orbital eccentricity and semi-major axis. We show that the distribution of radial and tangential velocities has a simple form and provide a fitting formula for this distribution. We also search for possible correlations between the infall directions of pairs of satellites, finding evidence for positive correlation at small angular separations as expected if some infall occurs along filaments. We also find (weak) evidence for correlations between the direction of the infall and infall velocity and the spin of the host halo.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The Halo Occupation Distribution and the Physics of Galaxy Formation

The halo occupation distribution (HOD) describes the bias between galaxies and dark matter by specifying (a) the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class and (b) the relative distributions of galaxies and dark matter within halos. We calculate predicted HODs for a Lambda-CDM cosmological model using an SPH hydrodynamic simulation and a semi-analytic (SA) galaxy formation model. Although the two methods predict different galaxy mass functions, their HOD predictions agree remarkably well. For mass-selected samples, the mean occupation <N(M)> exhibits a sharp cutoff at low halo masses, a slowly rising plateau for <N>~1-2, and a more steeply rising high occupancy regime. At low <N>, the mean pair and triple counts are well below Poisson expectations, with important consequences for small scale behavior of 2- and 3-point correlation functions. The HOD depends strongly on galaxy age, with high mass halos populated mainly by old galaxies and low mass halos by young galaxies. The SPH simulation supports several simplifying assumptions about HOD bias: the most massive galaxy in a halo usually lies close to the center and moves near the halo's mean velocity; satellite galaxies have the same radial profile and velocity dispersion as the dark matter; and the mean occupation at fixed halo mass is independent of the halo's larger scale environment. By applying the SPH and SA HODs to a large volume N-body simulation, we show that both methods predict slight, observable departures from a power-law galaxy correlation function. The predicted HODs are closely tied to the underlying galaxy formation physics, they offer useful guidance to theoretical models of galaxy clustering, and they will be tested empirically by ongoing analyses of galaxy redshift surveys. (Shortened)

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Polarization of the Cosmic Microwave Background from Non-Uniform Reionization

We study the signal in the Cosmic Microwave Background (CMB) polarization anisotropy resulting from patchy reionization. It is well known that the primordial polarization of the CMB is very sensitive to the details of reionization. Combining a semi-analytic model of galaxy formation, in which the optical depth to the reionization epoch is in the range 0.014 to 0.048, with a high resolution N-body simulation we find that reionization generates a peak with amplitude 0.05~0.15 μK at large angular scales. The position of this peak reveals the size of the horizon at reionization, whilst its amplitude is a measure of the optical depth to reionization. On small scales, ionized patches prior to full reionization create a second order polarization signal due to the coupling of the free electron density fluctuation with the quadrupole moment of the temperature anisotropy. Careful study reveals that the coupling generates the same power spectra for electric and magnetic modes, whose amplitudes of polarization anisotropies from this process are predicted to be ~10nK. The amplitude depends strongly on the total baryon density $Ω_b$ and on the spatial correlations of the free electron density in the ionized regions, and weakly on the fraction of ionizing photons able to escape their source galaxy, $f_{esc}$. The first- and second-order signals are therefore sensitive to the details of how the reionization occurred. Detection of these signals will place important constraints on the reionization history of the Universe.

0 citations0 reviewsNo code linkedOpen access