Source author record

Carlos Frenk

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

41 published item(s)

preprint2026arXiv

Modelling the evolution and influence of dust in cosmological simulations that include the cold phase of the interstellar medium

While marginal in mass terms, dust grains play an outsized role in both the physics and observation of the interstellar medium (ISM). However, explicit modelling of this ISM constituent remains uncommon in large cosmological simulations. In this work, we present a model for the life-cycle of dust in the ISM that couples to the forthcoming COLIBRE galaxy formation model, which explicitly simulates the cold ISM. We follow 6 distinct grain types: 3 chemical species, including carbon and two silicate grains, with 2 size bins each. Our dust model accounts for seeding of grains from stellar ejecta, self-consistent element-by-element metal yields and growth by accretion, grain size transfer (shattering and coagulation) and destruction of dust by thermal sputtering in the ISM. We detail the calibration of this model, particularly the use of a clumping factor, to account for unresolved gas clouds in which dust readily evolves. We present a fiducial run in a 25$^3$~cMpc$^3$ cosmological volume that displays good agreement with observations of the cosmic evolution of dust density, as well as the $z=0$ galaxy dust mass function and dust scaling relations. We highlight known tensions between observational datasets of the dust-to-gas ratio as a function of metallicity depending on which metallicity calibrator is used; our model favours higher-normalisation metallicity calibrators, which agree with the observations within 0.1~dex for stellar masses $>10^9 \; {\rm M_\odot}$. We compare the grain size distribution to observations of local galaxies, and find that our simulation suggests a higher concentration of small grains, associated with more diffuse ISM and the warm-neutral medium (WNM), which both play a key role in boosting H$_2$ content. Putting these results and modelling approaches in context, we set the stage for upcoming insights into the dusty ISM of galaxies using the COLIBRE simulations.

preprint2022arXiv

Baryon-driven decontraction in Milky Way mass haloes

We select a sample of Milky Way (MW) mass haloes from a high-resolution version of the EAGLE simulation to study their inner dark matter (DM) content and how baryons alter it. As in previous studies, we find that all haloes are more massive at the centre compared to their DM-only (DMO) counterparts at the present day as a result of the dissipational collapse of baryons during the assembly of the galaxy. However, we identify two processes that can reduce the central halo mass during the evolution of the galaxy. Firstly, gas blowouts induced by AGN feedback can lead to a substantial decrease of the central DM mass. Secondly, the formation of a stellar bar and its interaction with the DM can induce a secular expansion of the halo; the rate at which DM is evacuated from the central region by this process is related to the average bar strength and the timescale on which it acts determines how much the halo has decontracted. Although the inner regions of the haloes we have investigated are still more massive than their DMO counterparts at $z = 0$, they are significantly less massive than in the past and less massive than expected from the classic adiabatic contraction model. Since the MW has both a central supermassive black hole and a bar, the extent to which its halo has contracted is uncertain. This may affect estimates of the mass of the MW halo and of the expected signals in direct and indirect DM detection experiments.

preprint2021arXiv

A high-resolution cosmological simulation of a strong gravitational lens

We present a cosmological hydrodynamical simulation of a 10^13 Msun galaxy group and its environment (out to 10 times the virial radius) carried out using the EAGLE model of galaxy formation. Exploiting a novel technique to increase the resolution of the dark matter calculation independently of that of the gas, the simulation resolves dark matter haloes and subhaloes of mass 5x10^6 Msun . It is therefore useful for studying the abundance and properties of the haloes and subhaloes targeted in strong lensing tests of the cold dark matter model. We estimate the halo and subhalo mass functions and discuss how they are affected both by the inclusion of baryons in the simulation and by the environment. We find that the halo and subhalo mass functions have lower amplitude in the hydrodynamical simulation than in its dark matter only counterpart. This reflects the reduced growth of haloes in the hydrodynamical simulation due to the early loss of gas by reionisation and galactic winds and, additionally, in the case of subhaloes, disruption by enhanced tidal effects within the host halo due to the presence of a massive central galaxy. The distribution of haloes is highly anisotropic reflecting the filamentary character of mass accretion onto the cluster. As a result, there is significant variation in the number of structures with viewing direction. The median number of structures near the centre of the halo, when viewed in projection, is reduced by a factor of two when baryons are included.

preprint2021arXiv

The effects of self-interacting dark matter on the stripping of galaxies that fall into clusters

We use the Cluster-EAGLE (C-EAGLE) hydrodynamical simulations to investigate the effects of self-interacting dark matter (SIDM) on galaxies as they fall into clusters. We find that SIDM galaxies follow similar orbits to their Cold Dark Matter (CDM) counterparts, but end up with ${\sim}$25 per cent less mass by the present day. One in three SIDM galaxies are entirely disrupted, compared to one in five CDM galaxies. However, the excess stripping will be harder to observe than suggested by previous DM-only simulations because the most stripped galaxies form cores and also lose stars: the most discriminating objects become unobservable. The best test will be to measure the stellar-to-halo mass relation (SHMR) for galaxies with stellar mass $10^{10-11}\,\mathrm{M}_{\odot}$. This is 8 times higher in a cluster than in the field for a CDM universe, but 13 times higher for an SIDM universe. Given intrinsic scatter in the SHMR, these models could be distinguished with noise-free galaxy-galaxy strong lensing of ${\sim}32$ cluster galaxies.

preprint2019arXiv

Cusp or core? Revisiting the globular cluster timing problem in Fornax

We use N-body simulations to revisit the globular cluster (GC) ``timing problem'' in the Fornax dwarf spheroidal (dSph). In agreement with earlier work, we find that, due to dynamical friction, GCs sink to the center of dark matter halos with a cuspy inner density profile but ``stall'' at roughly 1/3 of the core radius ($r_{\rm core}$) in halos with constant-density cores. The timescales to sink or stall depend strongly on the mass of the GC and on the initial orbital radius, but are essentially the same for either cuspy (NFW) or cored halos normalized to have the same total mass within $r_{\rm core}$. Arguing against a cusp on the basis that GCs have not sunk to the center is thus no different from arguing against a core, unless all clusters are today at $\sim (1/3)\, r_{\rm core}$. This would imply a core radius exceeding $\sim 3$ kpc, much larger than seems plausible in any core-formation scenario. (The average projected distance of Fornax GCs is $\langle R_{\rm GC,Fnx}\rangle\sim 1$ kpc and its effective radius is $\sim 700$ pc.) A simpler explanation is that Fornax GCs have only been modestly affected by dynamical friction, as expected if clusters started orbiting at initial radii of order $\sim 1$-$2$ kpc, just outside Fornax's present-day half-light radius but well within the tidal radius imprinted by Galactic tides. This is not entirely unexpected. Fornax GCs are significantly older and more metal-poor than most Fornax stars, and such populations in dSphs tend to be more spatially extended than their younger and more metal-rich counterparts. Contrary to some earlier claims, our simulations further suggest that GCs do not truly ``stall'' at $\sim 0.3\, r_{\rm core}$, but rather continue decaying toward the center, albeit at reduced rates. We conclude that dismissing the presence of a cusp in Fornax based on the spatial distribution of its GC population is unwarranted.

preprint2019arXiv

Simulating cosmological substructure in the solar neighbourhood

We explore the predictive power of cosmological, hydrodynamical simulations for stellar phase space substructure and velocity correlations with the Auriga simulations and Aurigaia mock-Gaia catalogues. We show that at the solar circle the Auriga simulations commonly host phase space structures in the stellar component that have constant orbital energies and arise from accreted subhaloes. These structures can persist for a few Gyrs, even after coherent streams in position space have been erased. We also explore velocity two-point correlation functions and find this diagnostic is not deterministic for particular clustering patterns in phase space. Finally, we explore these structure diagnostics with the Aurigaia catalogues and show that the current catalogues have the ability to recover some structures in phase space but careful consideration is required to separate physical structures from numerical structures arising from catalogue generation methods.

preprint2019arXiv

Subhalo destruction in the Apostle and Auriga simulations

N-body simulations make unambiguous predictions for the abundance of substructures within dark matter halos. However, the inclusion of baryons in the simulations changes the picture because processes associated with the presence of a large galaxy in the halo can destroy subhalos and substantially alter the mass function and velocity distribution of subhalos. We compare the effect of galaxy formation on subhalo populations in two state-of-the-art sets of hydrodynamical CDM simulations of Milky Way mass halos, APOSTLE and AURIGA. We introduce a new method for tracking the orbits of subhalos between simulation snapshots that gives accurate results down to a few kiloparsecs from the centre of the halo. Relative to a dark matter-only simulation, the abundance of subhalos in APOSTLE is reduced by 50% near the centre and by 10% within r200. In AURIGA the corresponding numbers are 80% and 40%. The velocity distributions of subhalos are also affected by the presence of the galaxy, much more so in AURIGA than in APOSTLE . The differences on subhalo properties in the two simulations can be traced back to the mass of the central galaxies, which in AURIGA are typically twice as massive as those in APOSTLE . We show that some of the results from previous studies are inaccurate due to systematic errors in the modelling of subhalo orbits near the centre of halos.

preprint2015arXiv

Estimating the dark matter halo mass of our Milky Way using dynamical tracers

The mass of the dark matter halo of the Milky Way can be estimated by fitting analytical models to the phase-space distribution of dynamical tracers. We test this approach using realistic mock stellar halos constructed from the Aquarius N-body simulations of dark matter halos in the $Λ$CDM cosmology. We extend the standard treatment to include a Navarro-Frenk-White (NFW) potential and use a maximum likelihood method to recover the parameters describing the simulated halos from the positions and velocities of their mock halo stars. We find that the estimate of halo mass is highly correlated with the estimate of halo concentration. The best-fit halo masses within the virial radius, $R_{200}$, are biased, ranging from a 40\% underestimate to a 5\% overestimate in the best case (when the tangential velocities of the tracers are included). There are several sources of bias. Deviations from dynamical equilibrium can potentially cause significant bias; deviations from spherical symmetry are relatively less important. Fits to stars at different galactocentric radii can give different mass estimates. By contrast, the model gives good constraints on the mass within the half-mass radius of tracers even when restricted to tracers within 60kpc. The recovered velocity anisotropies of tracers, $β$, are biased systematically, but this does not affect other parameters if tangential velocity data are used as constraints.

preprint2015arXiv

Formation of In Situ Stellar Haloes in Milky Way-Mass Galaxies

We study the formation of stellar haloes in three Milky Way-mass galaxies using cosmological smoothed particle hydrodynamics simulations, focusing on the subset of halo stars that form in situ, as opposed to those accreted from satellites. In situ stars in our simulations dominate the stellar halo out to 20 kpc and account for 30-40 per cent of its total mass. We separate in situ halo stars into three straightforward, physically distinct categories according to their origin: stars scattered from the disc of the main galaxy ('heated disc'), stars formed from gas smoothly accreted on to the halo ('smooth' gas) and stars formed in streams of gas stripped from infalling satellites ('stripped' gas). We find that most belong to the stripped gas category. Those originating in smooth gas outside the disc tend to form at the same time and place as the stripped-gas population, suggesting that their formation is associated with the same gas-rich accretion events. The scattered disc star contribution is negligible overall but significant in the Solar neighbourhood, where ~90 per cent of stars on eccentric orbits once belonged to the disc. However, the distinction between halo and thick disc in this region is highly ambiguous. The chemical and kinematic properties of the different components are very similar at the present day, but the global properties of the in situ halo differ substantially between the three galaxies in our study. In our simulations, the hierarchical buildup of structure is the driving force behind not only the accreted stellar halo, but also those halo stars formed in situ.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): Panchromatic Data Release (far-UV --- far-IR) and the low-z energy budget

We present the GAMA Panchromatic Data Release (PDR) constituting over 230deg$^2$ of imaging with photometry in 21 bands extending from the far-UV to the far-IR. These data complement our spectroscopic campaign of over 300k galaxies, and are compiled from observations with a variety of facilities including: GALEX, SDSS, VISTA, WISE, and Herschel, with the GAMA regions currently being surveyed by VST and scheduled for observations by ASKAP. These data are processed to a common astrometric solution, from which photometry is derived for 221,373 galaxies with r<19.8 mag. Online tools are provided to access and download data cutouts, or the full mosaics of the GAMA regions in each band. We focus, in particular, on the reduction and analysis of the VISTA VIKING data, and compare to earlier datasets (i.e., 2MASS and UKIDSS) before combining the data and examining its integrity. Having derived the 21-band photometric catalogue we proceed to fit the data using the energy balance code MAGPHYS. These measurements are then used to obtain the first fully empirical measurement of the 0.1-500$μ$m energy output of the Universe. Exploring the Cosmic Spectral Energy Distribution (CSED) across three time-intervals (0.3-1.1Gyr, 1.1-1.8~Gyr and 1.8---2.4~Gyr), we find that the Universe is currently generating $(1.5 \pm 0.3) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$, down from $(2.5 \pm 0.2) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$ 2.3~Gyr ago. More importantly, we identify significant and smooth evolution in the integrated photon escape fraction at all wavelengths, with the UV escape fraction increasing from 27(18)% at z=0.18 in NUV(FUV) to 34(23)% at z=0.06. The GAMA PDR will allow for detailed studies of the energy production and outputs of individual systems, sub-populations, and representative galaxy samples at $z<0.5$. The GAMA PDR can be found at: http://gama-psi.icrar.org/

preprint2015arXiv

How well can cold-dark-matter substructures account for the observed radio flux-ratio anomalies?

Discrepancies between the observed and model-predicted radio flux ratios are seen in a number of quadruply-lensed quasars. The most favored interpretation of these anomalies is that CDM substructures present in lensing galaxies perturb the lens potentials and alter image magnifications and thus flux ratios. So far no consensus has emerged regarding whether or not the predicted CDM substructure abundance fully accounts for the lensing flux anomaly observations. Accurate modeling relies on a realistic lens sample in terms of both the lens environment and internal structures and substructures. In this paper we construct samples of generalised and specific lens potentials, to which we add (rescaled) subhalo populations from the galaxy-scale Aquarius and the cluster-scale Phoenix simulation suites. We further investigate the lensing effects from subhalos of masses several orders of magnitude below the simulation resolution limit. The resulting flux ratio distributions are compared to the currently best available sample of radio lenses. The observed anomalies in B0128+437, B0712+472 and B1555+375 are more likely to be caused by propagation effects or oversimplified lens modeling, signs of which are already seen in the data. Among the quadruple systems that have closely located image triplets/pairs, the anomalous flux ratios of MG0414+0534 can be reproduced by adding CDM subhalos to its macroscopic lens potential, with a probability of 5%-20%; for B0712+472, B1422+231, B1555+375 and B2045+265, these probabilities are only of a few percent. We hence find that CDM substructures are unlikely to be the whole reason for radio flux anomalies. We discuss other possible effects that might also be at work.

preprint2015arXiv

Particle tagging and its implications for stellar population dynamics

We establish a controlled comparison between the properties of galactic stellar halos obtained with hydrodynamical simulations and with `particle tagging'. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints' stars directly onto a suitably defined subset of dark matter particles in a collisionless, dark-matter-only simulation.Our study shows that there are conditions under which particle tagging generates good fits to the hydrodynamical stellar density profiles of a central Milky-Way-like galaxy and its most prominent substructure. Phase-space diffusion processes are crucial to reshaping the distribution of stars in infalling spheroidal systems and hence the final stellar halo. We conclude that the success of any particular tagging scheme hinges on this diffusion being taken into account, at a minimum by making use of `live' tagging schemes, in which particles are regularly tagged throughout the evolution of a galaxy.

preprint2015arXiv

The Milky Way system in LCDM cosmological simulations

We apply a semi-analytic galaxy formation model to two high resolution cosmological N-body simulations to investigate analogues of the Milky Way system. We select these according to observed properties of the Milky Way rather than by halo mass as in most previous work. For disk-dominated central galaxies with stellar mass (5--7) x 10d10Msun, the median host halo mass is 1.4 x 10d12Msun, with 1 sigma dispersion in the range [0.86, 3.1] x 10d12Msun, consistent with dynamical measurements of the Milky Way halo mass. For any given halo mass, the probability of hosting a Milky Way system is low, with a maximum of ~20% in haloes of mass ~10d12Msun. The model reproduces the V-band luminosity function and radial profile of the bright (MV < -9) Milky Way satellites. Galaxy formation in low mass haloes is found to be highly stochastic, resulting in an extremely large scatter in the relation between MV (or stellar mass) for satellites and the depth of the subhalo potential well in which they live, as measured by the maximum of the rotation curve, Vmax. We conclude that the "too big to fail" problem is an artifact of selecting satellites in N-body simulations according to subhalo properties: in 10% of cases we find that three or fewer of the brightest (or most massive) satellites have Vmax > 30 km/s. Our model predicts that around half of the dark matter subhaloes with Vmax > 20 km/s host satellites fainter than MV = -9 and so may be missing from existing surveys.

preprint2014arXiv

Creating mock catalogues of stellar haloes from cosmological simulations

We present a new technique for creating mock catalogues of the individual stars that make up the accreted component of stellar haloes in cosmological simulations and show how the catalogues can be used to test and interpret observational data. The catalogues are constructed from a combination of methods. A semi-analytic galaxy formation model is used to calculate the star formation history in haloes in an N-body simulation and dark matter particles are tagged with this stellar mass. The tags are converted into individual stars using a stellar population synthesis model to obtain the number density and evolutionary stage of the stars, together with a phase-space sampling method that distributes the stars while ensuring that the phase-space structure of the original N-body simulation is maintained. A set of catalogues based on the $Λ$CDM Aquarius simulations of Milky Way mass haloes have been created and made publicly available on a website. Two example applications are discussed that demonstrate the power and flexibility of the mock catalogues. We show how the rich stellar substructure that survives in the stellar halo precludes a simple measurement of its density profile and demonstrate explicitly how pencil-beam surveys can return almost any value for the slope of the profile. We also show that localized variations in the abundance of particular types of stars, a signature of differences in the composition of stellar populations, allow streams to be easily identified.

preprint2014arXiv

Herschel-ATLAS/GAMA:How does the far-IR luminosity function depend on galaxy group properties?

We use the Herschel ATLAS (H-ATLAS) Phase I data to study the conditional luminosity function of far-IR (250 um) selected galaxies in optically-selected galaxy groups from the GAMA spectroscopic survey, as well as environmental effects on the far-IR-to-optical colour. We applied two methods, which gave consistent results for the far-IR conditional luminosity functions. The direct matching method matches H-ATLAS sources to GAMA/SDSS galaxies, then links the optical counterparts to GAMA groups. The stacking method counts the number of far-IR sources within the projected radii of GAMA groups, subtracting the local background. We investigated the dependence of the far-IR (250 um) luminosity function on group mass in the range 10d12 < M_h < 10d14Msun/h and on redshift in the range 0 < z < 0.4, using a sample of 3000 groups containing H-ATLAS sources with GAMA redshifts over an area of 126 sq. deg. We find that the characteristic 250 um luminosity, L^*(250), increases with group mass up to M_h ~ 10d13Msun/h, but is roughly constant above this, while it increases with redshift at high group masses, but less so at low masses. We also find that the group far-IR luminosity-to-mass ratio L(250)/M_h increases with redshift and is higher in low-mass groups. We estimate that around 70% of the 250 um luminosity density in the local universe is contributed by groups with M_h > 10d12Msun/h. We also find that the far-IR-to-optical colours of H-ATLAS galaxies are independent of group mass over the range 10d12 < M_h <10d14Msun/h in the local universe. We also compare our observational results with recent semi-analytical models, and find that none of these galaxy formation model can reproduce the conditional far-IR luminosity functions of galaxy groups.

preprint2014arXiv

How well can cold-dark-matter substructures account for the observed lensing flux-ratio anomalies?

Lensing flux-ratio anomalies are most likely caused by gravitational lensing by small-scale dark matter structures. These anomalies offer the prospect of testing a fundamental prediction of the cold dark matter (CDM) cosmological model: the existence of numerous substructures that are too small to host visible galaxies. In two previous studies we found that the number of subhalos in the six high-resolution simulations of CDM galactic halos of the Aquarius project is not sufficient to account for the observed frequency of flux ratio anomalies seen in selected quasars from the CLASS survey. These studies were limited by the small number of halos used, their narrow range of masses (1-2E12 solar masses) and the small range of lens ellipticities considered. We address these shortcomings by investigating the lensing properties of a large sample of halos with a wide range of masses in two sets of high resolution simulations of cosmological volumes and comparing them to a currently best available sample of radio quasars. We find that, as expected, substructures do not change the flux-ratio probability distribution of image pairs and triples with large separations, but they have a significant effect on the distribution at small separations. For such systems, CDM substructures can account for a substantial fraction of the observed flux-ratio anomalies. For large close-pair separation systems, the discrepancies existing between the observed flux ratios and predictions from smooth halo models are attributed to simplifications inherent in these models which do not take account of fine details in the lens mass distributions.

preprint2013arXiv

Constraining the WDM Particle Mass with Milky Way Satellites

Well-motivated particle physics theories predict the existence of particles (such as sterile neutrinos) which acquire non-negligible thermal velocities in the early universe. These particles could behave as warm dark matter (WDM) and generate a small-scale cutoff in the linear density power spectrum which scales approximately inversely with the particle mass. If this mass is of order a keV, the cutoff occurs on the scale of dwarf galaxies. Thus, in WDM models the abundance of small galaxies, such as the satellites that orbit in the halo of the Milky Way, depends on the mass of the warm particle. The abundance also scales with the mass of the host galactic halo. We use the \galform semi-analytic model of galaxy formation to calculate the properties of galaxies in universes in which the dark matter is warm. Using this method, we can compare the predicted satellite luminosity functions to the observed data for the Milky Way dwarf spheroidals, and determine a lower bound on the thermally produced WDM particle mass. This depends strongly on the value of the Milky Way halo mass and, to some extent, on the baryonic physics assumed; we examine both of these dependencies. For our fiducial model we find that for a particle mass of 3.3 keV (the 2$σ$ lower limit found by Viel et al. from a recent analysis of the Lyman-$α$ forest) the Milky Way halo mass is required to be $> 1.4 \times 10^{12}$ \msun. For this same fiducial model, we also find that all WDM particle masses are ruled out (at 95% confidence) if the halo of the Milky Way has a mass smaller than $1.1 \times 10^{12}$ \msun, while if the mass of the Galactic halo is greater than 1.8 $\times 10^{12}$ \msun, only WDM particle masses larger than 2 keV are allowed.

preprint2013arXiv

Spatial and luminosity distributions of galactic satellites

We investigate the luminosity functions (LFs) and projected number density profiles of galactic satellites around isolated primaries of different luminosities. We measure these quantities for model satellites placed into the Millennium and Millennium II dark matter simulations by the GALFORM semi-analytic galaxy formation model for different bins of primary galaxy magnitude and we investigate their dependence on satellite luminosity. We compare our model predictions to the data of Guo et al. from the Sloan Digital Sky Survey Data Release 8 (SDSS DR8). First, we use a mock light-cone catalogue to verify that the method we used to count satellites in the SDSS DR8 is unbiased. We find that the radial distributions of model satellites are similar to those around comparable primary galaxies in the SDSS DR8, with only slight differences at low luminosities and small projected radii. However, when splitting the satellites by colour, the model and SDSS satellite systems no longer resemble one another, with many red model satellites, in contrast to the dominant blue fraction at similar luminosity in SDSS. The few model blue satellites are also significantly less centrally concentrated in the halo of their stacked primary than their SDSS counterparts. The implications of this result for the GALFORM model are discussed.

preprint2013arXiv

The Mass Profile and Accretion History of Cold Dark Matter Halos

We use the Millennium Simulation series to study the relation between the accretion history (MAH) and mass profile of cold dark matter halos. We find that the mean density within the scale radius, r_{-2} (where the halo density profile has isothermal slope), is directly proportional to the critical density of the Universe at the time when the main progenitor's virial mass equals the mass enclosed within r_{-2}. Scaled to these characteristic values of mass and density, the mean MAH, expressed in terms of the critical density of the Universe, M(ρ_{crit}(z)), resembles that of the enclosed density profile, M(<ρ>), at z=0. Both follow closely the NFW profile, suggesting that the similarity of halo mass profiles originates from the mass-independence of halo MAHs. Support for this interpretation is provided by outlier halos whose accretion histories deviate from the NFW shape; their mass profiles show correlated deviations from NFW and are better approximated by Einasto profiles. Fitting both M(<ρ>) and M(ρ_{crit}) with either NFW or Einasto profiles yield concentration and shape parameters that are correlated, confirming and extending earlier work linking the concentration of a halo with its accretion history. These correlations also confirm that halo structure is insensitive to initial conditions: only halos whose accretion histories differ greatly from the NFW shape show noticeable deviations from NFW in their mass profiles. As a result, the NFW profile provides acceptable fits to hot dark matter halos, which do not form hierarchically, and for fluctuation power spectra other than CDM. Our findings, however, predict a subtle but systematic dependence of mass profile shape on accretion history which, if confirmed, would provide strong support for the link between accretion history and halo structure we propose here.

preprint2013arXiv

The Mass-Concentration-Redshift Relation of Cold Dark Matter Halos

We use the Millennium Simulation series to investigate the mass and redshift dependence of the concentration of equilibrium cold dark matter (CDM) halos. We extend earlier work on the relation between halo mass profiles and assembly histories to show how the latter may be used to predict concentrations for halos of all masses and at any redshift. Our results clarify the link between concentration and the ``collapse redshift'' of a halo as well as why concentration depends on mass and redshift solely through the dimensionless ``peak height'' mass parameter, $ν(M,z)=δ_{\rm crit}(z)/σ(M,z)$. We combine these results with analytic mass accretion histories to extrapolate the $c(M,z)$ relations to mass regimes difficult to reach through direct simulation. Our model predicts that, at given $z$, $c(M)$ should deviate systematically from a simple power law at high masses, where concentrations approach a constant value, and at low masses, where concentrations are substantially lower than expected from extrapolating published empirical fits. This correction may reduce the expected self-annihilation boost factor from substructure by about one order of magnitude. The model also reproduces the $c(M,z)$ dependence on cosmological parameters reported in earlier work, and thus provides a simple and robust account of the relation between cosmology and the mass-concentration-redshift relation of CDM halos.

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.

preprint2012arXiv

On the Effects of Line-of-Sight Structures on Lensing Flux-ratio Anomalies in a LCDM Universe

The flux-ratio anomalies observed in multiply-lensed quasar images are most plausibly explained as the result of perturbing structures superposed on the underlying smooth matter distribution of the primary lens. The cold dark matter cosmological model predicts that a large number of substructures should survive inside larger halos but, surprisingly, this population alone has been shown to be insufficient to explain the observed distribution of the flux ratios of quasar's multiple images. Other halos (and their own subhalos) projected along the line of sight to the primary lens have been considered as additional source of perturbation. In this work, we use ray tracing through the Millennium II simulation to investigate the importance of projection effects due to halos and subhalos of mass m>1E8 Msun/h and extend our analysis to lower masses, m>1E6 Msun/h, using Monte-Carlo halo distributions. We find that the magnitude of the violation depends strongly on the density profile and concentration of the intervening halos, but clustering plays only a minor role. For a typical lensing geometry (lens at redshift 0.6 and source at redshift 2), background haloes (behind the main lens) are more likely to cause a violation than foreground halos. We conclude that line-of-sight structures can be as important as intrinsic substructures in causing flux-ratio anomalies. The combined effect of perturbing structures within the lens and along the line of sight in the LCDM universe results in a cusp-violation probability of 20-30%. This alleviates the discrepancy between models and current data, but a larger observational sample is required for a stronger test of the theory.

preprint2012arXiv

Satellite Galaxy Number Density Profiles in the Sloan Digital Sky Survey

We study the spatial distribution of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy catalogues. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For specific luminosity primaries we obtain robust statistical results by stacking as many as ~50, 000 galaxy systems. We find no evidence for any anistropy in the satellite galaxy distribution relative to the major axes of the primaries. We derive accurate projected number density profiles of satellites down to 4 magnitudes fainter than their primaries. We find the normalized satellite profiles generally have a universal form and can be well fitted by projected NFW profiles. The NFW concentration parameter increases with decreasing satellite luminosity while being independent of the luminosity of the primary except for very bright primaries. The profiles of the faintest satellites show deviations from the NFW form with an excess at small galactocentric projected distances. In addition, we quantify how the radial distribution of satellites depends on the colour of the satellites and on the colour and concentration of their primaries.

preprint2012arXiv

Space-quality data from balloon-borne telescopes: the High Altitude Lensing Observatory (HALO)

We present a method for attaining sub-arcsecond pointing stability during sub- orbital balloon flights, as designed for in the High Altitude Lensing Observatory (HALO) concept. The pointing method presented here has the potential to perform near-space quality optical astronomical imaging at 1-2% of the cost of space-based missions. We also discuss an architecture that can achieve sufficient thermomechanical stability to match the pointing stability. This concept is motivated by advances in the development and testing of Ultra Long Duration Balloon (ULDB) flights which promise to allow observation campaigns lasting more than three months. The design incorporates a multi-stage pointing architecture comprising: a gondola coarse azimuth control system, a multi-axis nested gimbal frame structure with arcsecond stability, a telescope de-rotator to eliminate field rotation, and a fine guidance stage consisting of both a telescope mounted angular rate sensor and guide CCDs in the focal plane to drive a fast-steering mirror. We discuss the results of pointing tests together with a preliminary thermo-mechanical analysis required for sub-arcsecond pointing at high altitude. Possible future applications in the areas of wide-field surveys and exoplanet searches are also discussed.

preprint2012arXiv

The haloes of bright satellite galaxies in a warm dark matter universe

High resolution N-body simulations of galactic cold dark matter haloes indicate that we should expect to find a few satellite galaxies around the Milky Way whose haloes have a maximum circular velocity in excess of 40 kms. Yet, with the exception of the Magellanic Clouds and the Sagittarius dwarf, which likely reside in subhaloes with significantly larger velocities than this, the bright satellites of the Milky Way all appear to reside in subhaloes with maximum circular velocities below 40 kms. As recently highlighted by Boylan-Kolchin et al., this discrepancy implies that the majority of the most massive subhaloes within a cold dark matter galactic halo are much too concentrated to be consistent with the kinematic data for the bright Milky Way satellites. Here we show that no such discrepancy exists if haloes are made of warm, rather than cold dark matter because these haloes are less concentrated on account of their typically later formation epochs. Warm dark matter is one of several possible explanations for the observed kinematics of the satellites.

preprint2011arXiv

A Halo Expansion (HEX) Technique for Approximating Simulated Dark Matter Haloes

We apply a basis function expansion method to create a time-evolving density/potential approximation of the late growth of simulated N-body dark matter haloes. We demonstrate how the potential of a halo from the Aquarius Project can be accurately represented by a small number of basis functions, and show that the halo expansion (HEX) method provides a way to replay simulations. We explore the level of accuracy of the technique as well as some of its limitations. We find that the number of terms included in the expansion must be large enough to resolve the large-scale distribution and shape of the halo but, beyond this, additional terms result in little further improvement. Particle and subhalo orbits can be integrated in this realistic, time-varying halo potential approximation, at much lower cost than the original simulation, with high fidelity for many individual orbits, and a good match to the distributions of orbital energy and angular momentum. Statistically, the evolution of structural subhalo properties, such as mass, half-mass radius and characteristic circular velocity, are very well reproduced in the halo expansion approximation over several gigayears. We demonstrate an application of the technique by following the evolution of an orbiting subhalo at much higher resolution than can be achieved in the original simulation. Our method represents a significant improvement over commonly used techniques based on static analytical descriptions of the halo potential.

preprint2011arXiv

The formation of shell galaxies similar to NGC 7600 in the cold dark matter cosmogony

We present new deep observations of 'shell' structures in the halo of the nearby elliptical galaxy NGC 7600, alongside a movie of galaxy formation in a cold dark matter universe (available at http://www.virgo.dur.ac.uk/shell-galaxies). The movie, based on an ab initio cosmological simulation, shows how continuous accretion of clumps of dark matter and stars creates a swath of diffuse circumgalactic structures. The disruption of a massive clump on a near-radial orbit creates a complex system of transient concentric shells which bare a striking resemblance to those of NGC 7600. With the aid of the simulation we interpret NGC 7600 in the context of the CDM model.

preprint2011arXiv

The Link between Galactic Satellite Orbits and Subhalo Accretion

We calculate the orbital angular momentum of dark matter subhaloes in the Aquarius simulations of cold dark matter galactic haloes. We calculate the orientation of their angular momentum relative to that of the spin vector of their host halo and find a variety of different configurations. All six Aquarius haloes contain statistically significant populations of subhalo orbits that are aligned with the main halo spin. All haloes posses a population of subhaloes that rotates in the same direction as the main halo and three of them possess, in addition, a population that rotates in the opposite direction. These configurations arise from the filamentary accretion of subhaloes. Quasi-planar distributions of coherently rotating satellites, such as those inferred in the Milky Way and other galaxies, arise naturally in simulations of a $Λ$CDM universe.

preprint2011arXiv

The Satellite Luminosity Functions of Galaxies in SDSS

We study the luminosity function of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy samples. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For primaries of similar luminosity to the Milky Way and M31, we are able to stack as many as $\sim\negthinspace 20,000$ galaxy systems to obtain robust statistical results. We derive the satellite luminosity function extending almost 8 magnitudes fainter than the primary galaxy. We also determine how the satellite luminosity function varies with the luminosity, colour and concentration of the primary. We find that, in the mean, isolated primaries of comparable luminosity to the Milky Way and M31 contain about a factor of two fewer satellites brighter than $M_V=-14$ than the average of the Milky Way and M31.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2009arXiv

Physical constraints on the central mass and baryon content of satellite galaxies

Recent analysis of the Milky Way's satellite galaxies reveals that these objects share a common central mass density, even though their luminosities range over five orders of magnitude. This observation can be understood in the context of galaxy formation theory by quantifying the factors which restrict the central mass density to a small range. One limit is set by the maximum mass that can collapse into a given region by the hierarchical growth of structure in the standard cold dark matter cosmology. Another limit comes from the natural thresholds which exist for gas to be able to cool and form a galaxy. The wide range of luminosities in these satellites reflect the effects of supernova feedback on the fraction of cooled baryons which are retained.

preprint2009arXiv

Towards accurate modelling of the ISW effect, the non-linear contribution

In a universe with a cosmological constant, the large-scale gravitational potential varies in time and this is, in principle, observable. Using an N-body simulation of a $Λ$CDM universe, we show that linear theory is not sufficiently accurate to predict the power spectrum of the time derivative, $\dotΦ$, needed to compute the imprint of large-scale structure on the cosmic microwave background (CMB). The linear part of the $\dotΦ$ power spectrum (the integrated Sachs-Wolfe effect or ISW) drops quickly as the relative importance of $Ω_Λ$ diminishes at high redshift, while the non-linear part (the Rees-Sciama effect or RS) evolves more slowly with redshift. Therefore, the deviation of the total power spectrum from linear theory occurs at larger scales at higher redshifts. The deviation occurs at $k\sim 0.1 $ $h$ Mpc$^{-1}$ at $z=0$. The cross-correlation power spectrum of the density $δ$ with $\dotΦ$ behaves differently to the power spectrum of $\dotΦ$. Firstly, the deviation from linear theory occurs at smaller scales ($k\sim 1 $ $h$ Mpc$^{-1}$ at $z=0$). Secondly, the correlation becomes negative when the non-linear effect dominates. For the cross-correlation power spectrum of galaxy samples with the CMB, the non-linear effect becomes significant at $l\sim 500$ and rapidly makes the cross power spectrum negative. For high redshift samples, the cross-correlation is expected to be suppressed by $5-10%$ on arcminute scales. The RS effect makes a negligible contribution to the large-scale ISW cross-correlation measurement. However, on arc-minute scales it will contaminate the expected cross-correlation signal induced by the Sunyaev-Zel'dovich effect.

preprint2008arXiv

The Properties of Sub-mm Galaxies in Hierarchical Models

We use the combined GALFORM semi-analytical model of galaxy formation and GRASIL spectrophotometric code to investigate the properties of galaxies selected via their sub-mm emission. Our fiducial model has previously been shown to fit the properties of local ULIRGs, as well as the number counts of faint sub-mm galaxies. Here, we test the model in detail by comparing the SEDs and stellar, dynamical, gas and halo masses of sub-mm galaxies against observational data. We precisely mimic the sub-mm and radio selection function of the observations and show that the predicted far-infrared properties of model galaxies with S_850>5mJy and S_1.4>30uJy are in good agreement with observations. Although the dust emission model does not assume a single dust temperature, the far-infrared SEDs are well described by single component modified black-body spectrum with characteristic temperature 32+/-5K. We also find evidence that the observations may have uncovered evolution in the far-infrared--radio relation in ULIRGs out to z~2. We show that the predicted redshift distribution of sub-mm galaxies provides a reasonable fit to the observational data with a median redshift z=2.0, with the radio-selected subset predicted to make up approximately 75% of the population. However, the predicted K-band and mid-infrared (3--8um) flux densities of the sub-mm galaxies (and LBGs) are up to a factor 10x fainter than observed. This discrepancy may indicate that the stellar masses of the sub-mm galaxies in the model are too low: M~10^10Mo, while observations suggest more massive systems, M~10^11Mo. Finally, we discuss the potential modifications to the models which may improve the fit to the observational data. [Abridged]

preprint2006arXiv

A marked correlation function analysis of halo formation times in the Millennium Simulation

We study the environmental dependence of the formation epoch of dark matter haloes in the Millennium Simulation: a ten billion particle N-body simulation of standard LambdaCDM cosmology. A sensitive test of this dependence - the marked correlation function - reveals highly significant evidence that haloes of a given mass form earlier in denser regions. We define a marked cross-correlation function, which helps quantify how this effect depends upon the choice of the halo population used to define the environment. The mean halo formation redshift as a function of the local overdensity in dark matter is also well determined, and we see an especially clear dependence for galaxy-sized haloes. This contradicts one of the basic predictions of the excursion set model of structure formation, even though we see that this theory predicts other features of the distribution of halo formation epochs rather well. It also invalidates an assumption usually employed in the popular halo, or HOD, models of galaxy clustering, namely that the distribution of halo properties is a function of halo mass but not of halo environment.

preprint2004arXiv

Substructure Analysis of Selected Low Richness 2dFGRS Clusters of Galaxies

Complementary one-, two-, and three-dimensional tests for detecting the presence of substructure in clusters of galaxies are applied to recently obtained data from the 2dF Galaxy Redshift Survey. The sample of 25 clusters used in this study includes 16 clusters not previously investigated for substructure. Substructure is detected at or greater than the 99% CL level in at least one test for 21 of the 25 clusters studied here. From the results, it appears that low richness clusters commonly contain subclusters participating in mergers. About half of the clusters have two or more components within 0.5 h^{-1} Mpc of the cluster centroid, and at least three clusters (Abell 1139, Abell 1663, and Abell S333) exhibit velocity-position characteristics consistent with the presence of possible cluster rotation, shear, or infall dynamics. The geometry of certain features is consistent with influence by the host supercluster environments. In general, our results support the hypothesis that low richness clusters relax to structureless equilibrium states on very long dynamical time scales (if at all).

preprint2004arXiv

The 2dF Galaxy Redshift Survey: The blue galaxy fraction and implications for the Butcher-Oemler effect

We derive the fraction of blue galaxies in a sample of clusters at z < 0.11 and the general field at the same redshift. The value of the blue fraction is observed to depend on the luminosity limit adopted, cluster-centric radius and, more generally, local galaxy density, but it does not depend on cluster properties. Changes in the blue fraction are due to variations in the relative proportions of red and blue galaxies but the star formation rate for these two galaxy groups remains unchanged. Our results are most consistent with a model where the star formation rate declines rapidly and the blue galaxies tend to be dwarfs and do not favour mechanisms where the Butcher-Oemler effect is caused by processes specific to the cluster environment.

preprint2002arXiv

The 2dF Galaxy Redshift Survey: A targeted study of catalogued clusters of galaxies

We have carried out a study of known clusters within the 2dF Galaxy Redshift Survey (2dFGRS) observed areas and have identified 431 Abell, 173 APM and 343 EDCC clusters. Precise redshifts, velocity dispersions and new centroids have been measured for the majority of these objects, and this information has been used to study the completeness of these catalogues, the level of contamination from foreground and background structures along the cluster's line of sight, the space density of the clusters as a function of redshift, and their velocity dispersion distributions. We find that the Abell and EDCC catalogues are contaminated at the level of about 10%, whereas the APM catalogue suffers only 5% contamination. If we use the original catalog centroids, the level of contamination rises to approximately 15% for the Abell and EDCC catalogues, showing that the presence of foreground and background groups may alter the richness of clusters in these catalogues. There is a deficiency of clusters at $z \sim 0.05$ that may correspond to a large underdensity in the Southern hemisphere. From the cumulative distribution of velocity dispersions for these clusters, we derive an upper limit to the space density of $σ> 1000 \kms$ clusters of $3.6 \times 10^{-6} \hdens$. This result is used to constrain models for structure formation; our data favour low-density cosmologies, subject to the usual assumptions concerning the shape and normalization of the power spectrum.

preprint2002arXiv

The 2dF Galaxy Redshift Survey: The environmental dependence of galaxy star formation rates near clusters

We have measured the equivalent width of the H-alpha emission line for 11006 galaxies brighter than M_b=-19 (LCDM) at 0.05<z<0.1 in the 2dF Galaxy Redshift Survey (2dF), in the fields of seventeen known galaxy clusters. The limited redshift range ensures that our results are insensitive to aperture bias, and to residuals from night sky emission lines. We use these measurements to trace mustar, the star formation rate normalized to Lstar, as a function of distance from the cluster centre, and local projected galaxy density. We find that the distribution of mustar steadily skews toward larger values with increasing distance from the cluster centre, converging to the field distribution at distances greater than ~3 times the virial radius. A correlation between star formation rate and local projected density is also found, which is independent of cluster velocity dispersion and disappears at projected densities below ~1 galaxy (brighter than M_b=-19) per Mpc^2. This characteristic scale corresponds approximately to the mean density at the cluster virial radius. The same correlation holds for galaxies more than two virial radii from the cluster centre. We conclude that environmental influences on galaxy properties are not restricted to cluster cores, but are effective in all groups where the density exceeds this critical value. The present day abundance of such systems, and the strong evolution of this abundance, makes it likely that hierarchical growth of structure plays a significant role in decreasing the global average star formation rate. Finally, the low star formation rates well beyond the virialised cluster rule out severe physical processes, such as ram pressure stripping of disk gas, as being completely responsible for the variations in galaxy properties with environment.

preprint2001arXiv

The 2dF Galaxy Redshift Survey: Near Infrared Galaxy Luminosity Functions

We combine the 2MASS extended source catalogue and the 2dFGRS to produce an IR selected galaxy catalogue with 17,173 measured redshifts. We use this extensive dataset to estimate the J and K-band galaxy luminosity functions. The LFs are fairly well fit by Schechter functions with J: M*-5log h= -22.36+/-0.02, alpha= -0.93+/-0.04, Phi=0.0104+/-0.0016 h^3/Mpc^3 and K: M*-5log h= -23.44+/-0.03, alpha=-0.96+/-0.05, Phi=0.0108+/-0.0016 h^3/Mpc^3 (2MASS Kron magnitudes). These parameters assume a cosmological model with Omega=0.3 and Lambda=0.7. With datasets of this size, systematic rather than random errors are the dominant source of uncertainty in the determination of the LF. We carry out a careful investigation of possible systematic effects in our data. The surface brightness distribution of the sample shows no evidence that significant numbers of low surface brightness or compact galaxies are missed by the survey. We estimate the present-day distributions of B-K and J-K colours as a function of absolute magnitude and use models of the galaxy stellar populations, constrained by the observed optical and infrared colours, to infer the galaxy stellar mass function. Integrated over all galaxy masses, this yields a total mass fraction in stars (in units of the critical mass density) of Omega_*.h= (1.6+/-0.24)/10^3 for a Kennicutt IMF and Omega_*.h= (2.9+/-0.43)/10^3 for a Salpeter IMF. These values agree with those inferred from observational estimates of the star formation history of the universe provided that dust extinction corrections are modest.

preprint1994arXiv

The Clustering Properties of IRAS Galaxies

We investigate the clustering of galaxies in the QDOT redshift survey of IRAS galaxies. We find that the redshift space two-point correlation function is well approximated by a power-law of slope -1.11+\-0.09, with clustering length 3.87+\-0.32 Mpc/h out to pair separations of 25 Mpc\h. On scales larger than 40 Mpc/h, the correlation function is consistent with zero. The r.m.s. fluctuation in the count of QDOT galaxies above the Poisson level in spheres of radius 8 Mpc/h is σ_8^{IRAS}=0.58+\-0.14, showing that fluctuations in the distribution of IRAS galaxies on these scales are smaller than those of optical galaxies by a factor of about 0.65, with an uncertainty of 25%. We find no detectable difference between the correlation functions measured in redshift space and in real space, leading to a $2σ$ limit of b_{IRAS}/Ω^{0.6} >1.05, where b_{IRAS} is the bias factor for IRAS galaxies and Ωis the cosmological density parameter. The QDOT autocorrelation function calculated in concentric shells increases significantly with shell radius. This difference is more likely due to sampling fluctuations than to an increase of the clustering strength with galaxy luminosity, but the two effects are difficult to disentangle; our data allow at most an increase of \sim 20% in clustering strength for each decade in luminosity. For separations greater than \sim 3 Mpc/h, the cross-correlation function of Abell clusters (with Richness R\ge 1) and QDOT galaxies is well approximated by a power-law of slope -1.81+\-0.10, with clustering length 10.10+\-0.45 Mpc/h, and no significant signal beyond \gsim 50 Mpc\h. This cross-correlation depends only weakly on cluster richness.