Search

A new generation of interferometric instruments is emerging which aim to use intensity mapping of redshifted $21\,$cm radiation to measure the large-scale structure of the Universe at $z\simeq 1-2$ over wide areas of sky. While these instruments typically have limited angular resolution, they cover huge volumes and thus can be used to provide large samples of rare objects. In this paper we study how well such instruments could find spatially extended large-scale structures, such as cosmic voids, using a matched filter formalism. Such a formalism allows us to work in Fourier space, the natural space for interferometers, and to study the impact of finite $u-v$ coverage, noise and foregrounds on our ability to recover voids. We find that in the absence of foregrounds such instruments would provide enormous catalogs of voids, with high completeness, but that control of foregrounds is key to realizing this goal., Comment: 14 pages, 8 figures, minor revisions to match version accepted by MNRAS

We present the CAlibrating LYMan-$\alpha$ with H$\alpha$ (CALYMHA) pilot survey and new results on Lyman-$\alpha$ (Lya) selected galaxies at z~2. We use a custom-built Lya narrow-band filter at the Isaac Newton Telescope, designed to provide a matched volume coverage to the z=2.23 Ha HiZELS survey. Here we present the first results for the COSMOS and UDS fields. Our survey currently reaches a 3$\sigma$ line flux limit of ~4x10$^{-17}$ erg/s/cm$^{2}$, and a Lya luminosity limit of ~10$^{42.3}$ erg/s. We find 188 Lya emitters over 7.3x10$^5$ Mpc$^{3}$, but also find significant numbers of other line emitting sources corresponding to HeII, CIII] and CIV emission lines. These sources are important contaminants, and we carefully remove them, unlike most previous studies. We find that the Lya luminosity function at z=2.23 is very well described by a Schechter function up to L~10$^{43}$ erg/s with L$^*=10^{42.59+-0.05}$ erg/s, $\phi^*=10^{-3.09+-0.08}$ Mpc$^{-3}$ and $\alpha$=-1.75+-0.15. Above L~10$^{43}$ erg/s the Lya luminosity function becomes power-law like, driven by X-ray AGN. We find that Lya-selected emitters have a high escape fraction of 37+-7%, anti-correlated with Lya luminosity and correlated with Lya equivalent width. Lya emitters have ubiquitous large (~40 kpc) Lya haloes, 2x larger than their Ha extents. By directly comparing our Lya and Ha luminosity functions we find that the global/overall escape fraction of Lya photons (within a 13 kpc radius) from the full population of star-forming galaxies is 5.1+-0.2% at the peak of the star formation history. An extra 3.3+-0.3% of Lya photons likely still escape, but at larger radii., Comment: 19 pages, 11 figures. Accepted for publication in MNRAS. Catalogue of line emitters available with MNRAS refereed paper

Mid-infrared (IR) colour selection techniques have proved to be very efficient in finding active galactic nuclei (AGN). This is because the AGN heats the surrounding dust producing warm mid-IR colours. Using the Wide-field Infrared Survey Explorer 3.6, 4.5 and 12 μm colours, the largest sample of IR selected AGN has already been produced containing 1.4 million AGN over the whole sky. Here, we explore the X-ray properties of this AGN sample by cross-correlating it with the subsample of the 3XMM X-ray catalogue that has available X-ray spectra and at the same time optical spectroscopy from Sloan Digital Sky Survey (SDSS). Our goal is to find rare luminous obscured AGN. Our final sample contains 65 QSOs with log νLν ≥ 46.2 erg s−1. This IR luminosity cut corresponds to log LX ≈ 45 erg s−1, at the median redshift of our sample (z = 2.3), that lies at the bright end of the X-ray luminosity function at z > 2. The X-ray spectroscopic analysis reveals seven obscured AGN having a column density NH > 1022 cm−2. Six of them show evidence for broad [C IV] absorption lines and five are classified as broad absorption line QSOs. We fit the optical spectra of our X-ray absorbed sources to estimate the optical reddening. We find that none of these show any obscuration according to the optical continuum. These sources add to the growing evidence for populations of luminous QSOs with evidence for substantial absorption by outflowing ionized material, similar to those expected to be emerging from their absorbing cocoons in the framework of AGN/galaxy co-evolution., The research leading to these results has received funding from the European Union’s Horizon 2020 Programme under the AHEAD project (grant agreement no. 654215). GM acknowledges financial support from the AHEAD project that is funded by the European Union as Research and Innovation Action under Grant no. 654215. FJC and IO-P acknowledge financial support through grant AYA2015-64346-C2-1-P (MINECO/FEDER). SM acknowledges financial support by the Spanish Ministry of Economy and Competitiveness through Grant no. AYA2016-76730-P. Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the US Department of Energy Office of Science.

A prediction of the standard λ cold dark matter cosmology is that dark matter (DM) haloes are teeming with numerous self-bound substructure or subhaloes. The precise properties of these subhaloes represent important probes of the underlying cosmological model. We use data from Via Lactea II and Exploring the Local Volume in Simulations N-body simulations to learn about the structure of subhaloes with masses 10-10 hM. Thanks to a superb subhalo statistics, we study subhalo properties as a function of distance to host halo centre and subhalo mass, and provide a set of fits that accurately describe the subhalo structure. We also investigate the role of subhaloes on the search for DM annihilation. Previous work has shown that subhaloes are expected to boost the DM signal of their host haloes significantly. Yet, these works traditionally assumed that subhaloes exhibit similar structural properties than those of field haloes, while it is known that subhaloes are more concentrated. Building upon our N-body data analysis, we refine the substructure boost model of SA & Acute;nchez-Conde & Prada (2014), and find boosts that are a factor 2-3 higher. We further refine the model to include unavoidable tidal stripping effects on the subhalo population. For field haloes, this introduces a moderate (20-30 per cent) suppression. Yet, for subhaloes like those hosting dwarf galaxy satellites, tidal stripping plays a critical role, the boost being at the level of a few tens of percent at most. We provide a parametrization of the boost for field haloes that can be safely applied over a wide halo mass range.© 2017 The Authors., AM has been supported by the Fundacao para a Ciencia e a Tecnologia (FCT) of Portugal and thanks SLAC National Accelerator Laboratory and Instituto de Fisica Corpuscular in Valencia for hospitality. AM has been also partially supported by the FCT through the project PTDC/FIS- NUC/0548/2012. MASC is a Wenner-Gren Fellow and acknowledges the support of the WennerGren Foundations to develop his research. He also acknowledges the support from the MULTIDARK project of Spanish MCINN Consolider-Ingenio: CSD2009-00064. MASC is very grateful to SLAC and the Instituto de Fisica Teorica in Madrid, where part of this work was done. SPR is supported by a Ramon y Cajal contract, by the Spanish MINECO under Grants No. FPA2014-54459-P and SEV-2014-0398 and by the Generalitat Valenciana under Grant No. PROMETEOII/2014/049. AM and SPR are also partially supported by the Portuguese FCT through the CFTP-FCT Unit 777 (PEst-OE/FIS/UI0777/2013).

It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional semi-analytic galaxy formation models. We use the Eagle simulations to assess assumptions built into the Galform semi-analytic model, focussing on those relating to baryon cycling, angular momentum and feedback. We show that the assumption in Galform that newly formed stars have the same specific angular momentum as the total disc leads to a significant overestimate of the total stellar specific angular momentum of disc galaxies. In Eagle, stars form preferentially out of low specific angular momentum gas in the interstellar medium (ISM) due to the assumed gas density threshold for stars to form, leading to more realistic galaxy sizes. We find that stellar mass assembly is similar between Galform and Eagle but that the evolution of gas properties is different, with various indications that the rate of baryon cycling in Eagle is slower than is assumed in Galform. Finally, by matching individual galaxies between Eagle and Galform, we find that an artificial dependence of AGN feedback and gas infall rates on halo mass doubling events in Galform drives most of the scatter in stellar mass between individual objects. Put together our results suggest that the Galform semi-analytic model can be significantly improved in light of recent advances., Comment: some references corrected, 31 pages, 19 figures, submitted to MNRAS after responding to a first referee report

Large-scale surveys of galaxies require accurate theoretical predictions of the dark matter clustering for thousands of mock galaxy catalogues. We demonstrate that this goal can be achieved with the new Parallel Particle-Mesh (PM) N-body code GLAM at a very low computational cost. We run ~22 000 simulations with ~2 billion particles that provide ~1 per cent accuracy of the dark matter power spectra P(k) for wavenumbers up to k ~ 1 h Mpc. Using this large data set we study the power spectrum covariance matrix. In contrast to many previous analytical and numerical results, we find that the covariance matrix normalized to the power spectrum C(k, k')/P(k)P(k') has a complex structure of non-diagonal components: an upturn at small k, followed by a minimum at k ≈ 0.1-0.2 h Mpc, and a maximum at k ≈ 0.5-0.6 h Mpc. The normalized covariance matrix strongly evolves with redshift: C(k, k') ∝ δα(t)P(k)P(k'), where δ is the linear growth factor and α ≈ 1-1.25, which indicates that the covariance matrix depends on cosmological parameters. We also show that waves longer than 1 h Gpc have very little impact on the power spectrum and covariance matrix. This significantly reduces the computational costs and complexity of theoretical predictions: relatively small volume ~(1 h Gpc) simulations capture the necessary properties of dark matter clustering statistics. As our results also indicate, achieving ~1 per cent errors in the covariance matrix for k < 0.50 h Mpc requires a resolution better than ε ~ 0.5 h Mpc.© 2018 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society., AK acknowledges support of the Fulbright Foundation and support of the Instituto de Fisica Teorica (IFT), Madrid, Spain. FP acknowledges support from the Spanish Ministerio de Economia, Industria y Competitividad grant AYA2010-21231-C02-01.

We use the SDSS/DR8 galaxy sample to study the radial distribution of satellite galaxies around isolated primaries, comparing to semi-analytic models of galaxy formation based on the Millennium and Millennium-II simulations. SDSS satellites behave differently around high- and low-mass primaries: those orbiting objects with $M_*>10^{11}M_\odot$ are mostly red and are less concentrated towards their host than the inferred dark matter halo, an effect that is very pronounced for the few blue satellites. On the other hand, less massive primaries have steeper satellite profiles that agree quite well with the expected dark matter distribution and are dominated by blue satellites, even in the inner regions where strong environmental effects are expected. In fact, such effects appear to be strong only for primaries with $M_* > 10^{11}M_\odot$. This behaviour is not reproduced by current semi-analytic simulations, where satellite profiles always parallel those of the dark matter and satellite populations are predominantly red for primaries of all masses. The disagreement with SDSS suggests that environmental effects are too efficient in the models. Modifying the treatment of environmental and star formation processes can substantially increase the fraction of blue satellites, but their radial distribution remains significantly shallower than observed. It seems that most satellites of low-mass primaries can continue to form stars even after orbiting within their joint halo for 5 Gyr or more., 19 pages, 11 figures, accepted by MNRAS

We present the first results from the Juropa Hubble Volume project, based a large N-body, dark matter-only cosmological simulation with a volume of V = (6 h-1 Gpc)3, containing 60003 particles, performed within the concordance cold dark matter cosmological model. The simulation volume is sufficient to probe extremely large length-scales in the universe, whilst at the same time the particle count is high enough so that dark matter haloes down to 1.5 × 1012 h-1M⊙ can be resolved. At z = 0, we identify over 400 million haloes. The cluster mass function is derived using three different halofinders and compared to fitting functions in the literature. The distribution of clusters of maximal mass across redshifts agrees well with predicted masses of extreme objects, and we explicitly confirm that the Poisson distribution is very good at describing the distribution of rare clusters. The Poisson distribution also matches well the level to which cosmic variance can be expected to affect number counts of high-mass clusters. We find that objects like the Bullet cluster exist in the far-tail of the distribution of mergers in terms of relative collisional speed. We also derive the number counts of voids in the simulation box for z = 0, 0.5 and 1. © 2013 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society., ITI was supported by the SEPNet and the Science and Technology Facilities Council grants ST/F002858/1 and ST/I000976/1. AK was supported by the Spanish Ministerio de Ciencia e Innovación (MICINN) in Spain through the Ramón y Cajal programme as well as the grants AYA 2009-13875-C03-02, AYA2009-12792-C03-03, CSD2009-00064, CAM S2009/ESP-1496 and the Ministerio de Economía y Competitividad (MINECO) through grantAYA2012-31101.GYacknowledges support fromMINECO(Spain) under research grantsAYA2009-13875-C03-02, FPA2009-08958, AYA2012-31101 and Consolider Ingenio SyeC CSD2007-0050. He also thanks Comunidad de Madrid for support under PRICIT project ASTROMADRID (S2009/ESP-146). JMD and EMG acknowledge support of the consolider project CSD2010-00064 funded by the Ministerio de Economía y Competitividad and AYA2012-39475-C02-01.

Quasars represent the brightest Active Galactic Nuclei (AGN) in the Universe and are thought to indicate the location of prodigiously growing Black Holes (BHs), with luminosities as high as 10^48 erg/sec. It is often expected though that such an extremely energetic process will take place in the most massive bound structures in the dark matter (DM) distribution. We show that in contrast to this expectation, in a galaxy formation model which includes AGN feedback, quasars are predicted to live in average DM halo environments with typical masses of a few times 10^12 Msun. This fundamental prediction arises from the fact that quasar activity (i.e., BH accretion with luminosity greater than 10^46 erg/sec) is inhibited in DM haloes where AGN feedback operates. The galaxy hosts of quasars in our simulations are identified with over massive (in gas and stars) spheroidal galaxies, in which BH accretion is triggered via a galaxy merger or secular processes. We further show that the z=0 descendants of high redshift (z~6) QSOs span a wide range of morphologies, galaxy and halo masses. The z~6 BHs typically grow only by a modest factor by the present day. Remarkably, high redshift QSOs never inhabit the largest DM haloes at that time and their descendants are very seldom found in the most massive haloes at z=0. We also show that observationally it is very likely to find an enhancement in the abundance of galaxies around quasars at z~5. However, these enhancements are considerably weaker compared to the overdensities expected at the extreme peaks of the DM distribution. Thus, it is very unlikely that a quasar detected in the $z\gtrsim5$ Universe pinpoints the location of the progenitors of superclusters in the local Universe., 11 pages, 6 figures, submitted to MNRAS (comments are welcome)

We present Gemini griz' photometry for 521 globular cluster (GC) candidates in a 5.5 × 5.5 arcmin2 field centred 3.8 arcmin to the south and 0.9arcmin to the west of the centre of the giant elliptical galaxy NGC 4486. All these objects have previously published (C - T1) photometry. We also present new (C - T1) photometry for 338 globulars, within 1.7 arcmin in galactocentric radius, which have (g - z) colours in the photometric system adopted by the Virgo Cluster Survey of the Advanced Camera for Surveys of the Hubble Space Telescope (HST). These photometric data are used to define a self-consistent multicolour grid (avoiding polynomial fits) and preliminarily calibrated in terms of two chemical abundance scales. The resulting multicolour colour-chemical abundance relations are used to test GC chemical abundance distributions. This is accomplished by modelling the 10 GC colour histograms that can be defined in terms of the Cgriz' bands. Our results suggest that the best fit to the GC observed colour histograms is consistent with a genuinely bimodal chemical abundance distribution NGC(Z). On the other side, each ('blue' and 'red') GC subpopulation follows a distinct colour-colour relation., Facultad de Ciencias Astronómicas y Geofísicas, Instituto de Astrofísica de La Plata

CH, MV, MC, HHo, CS and AC acknowledge support from the European Research Council under FP7 grant number 279396 (MV, MC, CS, HHo), grant number 240185 (AC and CH) and grant number G47112 (CH). MV acknowledges support from the Netherlands Organisation for Scientific Research (NWO) through grants 614.001.103. HHi is supported by an Emmy Noether grant (No. Hi 1495/2-1) of the Deutsche Forschungsgemeinschaft. RN acknowledges support from the German Federal Ministry for Economic Affairs and Energy (BMWi) provided via DLR under project no. 50QE1103. TM-R and PN acknowledge support from an European Research Council Starting Grant (DEGAS-259586). This work is supported by the Deutsche Forschungsgemeinschaft in the framework of the TR33 ‘The Dark Universe’. GVK acknowledges financial support from the Netherlands Research School for Astronomy (NOVA) and Target. Target is supported by Samenwerkingsverband Noord Nederland, European fund for regional development, Dutch Ministry of economic affairs, Pieken in de Delta, Provinces of Groningen and Drenthe. EvU acknowledges support from an STFC Ernest Rutherford Research Grant, grant reference ST/L00285X/1. Galaxies and their dark matter haloes are part of a complex network of mass structures, collectively called the cosmic web. Using the tidal tensor prescription these structures can be classified into four cosmic environments: voids, sheets, filaments and knots. As the cosmic web may influence the formation and evolution of dark matter haloes and the galaxies they host, we aim to study the effect of these cosmic environments on the average mass of galactic haloes. To this end we measure the galaxy-galaxy lensing profile of 91 195 galaxies, within 0.039 < z

We study the contribution of the first galaxies to the far-infrared/sub-millimeter (FIR/sub-mm) extragalactic background light (EBL) by implementing an analytical model for dust emission. We explore different dust models, assuming different grain size distributions and chemical compositions. According to our findings, observed re-radiated emission from dust in dwarf-size galaxies at $z \sim 10$ would peak at a wavelength of $\sim 500 \mu {\rm m}$ with observed fluxes of $\sim 10^{-3} - 10^{-2}$ nJy, which is below the capabilities of current observatories. In order to be detectable, model sources at these high redshifts should exhibit luminosities of $\gtrsim 10^{12} L_{\odot}$, comparable to that of local ultra-luminous systems. The FIR/sub-mm EBL generated by primeval galaxies peaks at $\sim 500 \mu {\rm m}$, with an intensity ranging from $\sim 10^{-4}$ to $10^{-3} {\rm nW \ m^{-2} \ sr^{-1}}$, depending on dust properties. These values are $\sim 3 - 4$ orders of magnitude below the absolute measured cosmic background level, suggesting that the first galaxies would not contribute significantly to the observed FIR/sub-mm EBL. Our model EBL exhibits a strong correlation with the dust-to-metal ratio, where we assume a fiducial value of $D = 0.005$, increasing almost proportionally to it. Thus, measurements of the FIR/sub-mm EBL could provide constraints on the amount of dust in the early Universe. Even if the absolute signal from primeval dust emission may be undetectable, it might still be possible to obtain information about it by exploring angular fluctuations at $\sim 500 \mu {\rm m}$, close to the peak of dust emission from the first galaxies., Comment: 14 pages, 7 figures. Accepted for publication in MNRAS

© 2016 The Authors. Current and future large redshift surveys, as the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), will use emission-line galaxies (ELGs) to probe cosmological models by mapping the large-scale structure of the Universe in the redshift range 0.6 < z < 1.7. With current data, we explore the halo-galaxy connection by measuring three clustering properties of g-selected ELGs as matter tracers in the redshift range 0.6 < z < 1: (i) the redshift-space two-point correlation function using spectroscopic redshifts from the BOSS ELG sample and VIPERS; (ii) the angular two-point correlation function on the footprint of the CFHT-LS; (iii) the galaxy-galaxy lensing signal around the ELGs using the CFHTLenS. We interpret these observations by mapping them on to the latest high-resolution MultiDark Planck N-body simulation, using a novel (Sub)Halo-Abundance Matching technique that accounts for the ELG incompleteness. ELGs at z ~ 0.8 live in haloes of (1 ± 0.5) × 1012 h-1M⊙ and 22.5 ± 2.5 per cent of them are satellites belonging to a larger halo. The halo occupation distribution of ELGs indicates that we are sampling the galaxies in which stars form in the most efficient way, according to their stellar-to-halo mass ratio.

We present a direct measurement of the mean halo occupation distribution (HOD) of galaxies taken from the eleventh data release (DR11) of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey (BOSS). The HOD of BOSS low-redshift (LOWZ: $0.2 < z < 0.4$) and Constant-Mass (CMASS: $0.43, Comment: Accepted for publication in MNRAS; 16 pages, 9 figures, 6 tables (1 electronic)

© 2016 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society. We use subhalo abundance matching (SHAM) to model the stellar mass function (SMF) and clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) 'CMASS' sample at z ~ 0.5. We introduce a novel method which accounts for the stellar mass incompleteness of CMASS as a function of redshift, and produce CMASS mock catalogues which include selection effects, reproduce the overall SMF, the projected two-point correlation function wp, the CMASS dn/dz, and are made publicly available. We study the effects of assembly bias above collapse mass in the context of 'age matching' and show that these effects are markedly different compared to the ones explored by Hearin et al. at lower stellar masses. We construct two models, one in which galaxy colour is stochastic ('AbM' model) as well as a model which contains assembly bias effects ('AgM' model). By confronting the redshift dependent clustering of CMASS with the predictions from our model, we argue that that galaxy colours are not a stochastic process in high-mass haloes. Our results suggest that the colours of galaxies in high-mass haloes are determined by other halo properties besides halo peak velocity and that assembly bias effects play an important role in determining the clustering properties of this sample.

We present a study of the clustering and halo occupation distribution of Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a halo abundance matching (HAM) clustering model that assigns galaxies to dark matter haloes selected from the large BigMultiDark N-body simulation of a flat Λ cold dark matter Planck cosmology. We compare the observational data with the simulated ones on a light cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fibre collisions are included in the model, which reproduces within 1σ errors the observed monopole of the two-point correlation function at all relevant scales: from the smallest scales, 0.5 h-1 Mpc, up to scales beyond the baryon acoustic oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to k ~ 1 h Mpc-1), and the three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar-to-halo mass relation for the CMASS sample measured using weak lensing in the Canada-France-Hawaii Telescope Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1σ. The BigMD-BOSS light cone including properties of BOSS galaxies and halo properties is made publicly available.

This paper presents Gemini-gri' high-quality photometry for cluster candidates in the field of NGC 1316 (Fornax A) as part of a study that also includes GMOS spectroscopy. A preliminary discussion of the photometric data indicates the presence of four stellar cluster populations with distinctive features in terms of age, chemical abundance and spatial distribution. Two of them seem to be the usually old (metal poor and metal rich) populations typically found in elliptical galaxies. In turn, an intermediate-age (5 Gyr) globular cluster population is the dominant component of the sample (as reported by previous papers). We also find a younger cluster population with a tentative age of ≈ 1 Gyr., Facultad de Ciencias Astronómicas y Geofísicas, Instituto de Astrofísica de La Plata

Using realistic cosmological simulations of Milky Way sized haloes, we study their dynamical state and the accuracy of inferring their mass profiles with steady-state models of dynamical tracers. We use a new method that describes the phase-space distribution of a steady-state tracer population in a spherical potential without any assumption regarding the distribution of their orbits. Applying the method to five haloes from the Aquarius $\Lambda$CDM N-body simulation, we find that dark matter particles are an accurate tracer that enables the halo mass and concentration parameters to be recovered with an accuracy of $5\%$. Assuming a potential profile of the NFW form does not significantly affect the fits in most cases, except for halo A whose density profile differs significantly from the NFW form, leading to a $30\%$ bias in the dynamically fitted parameters. The existence of substructures in the dark matter tracers only affects the fits by $\sim 1\%$. Applying the method to mock stellar haloes generated by a particle-tagging technique, we find the stars are farther from equilibrium than dark matter particles, yielding a systematic bias of $\sim 20\%$ in the inferred mass and concentration parameter. The level of systematic biases obtained from a conventional distribution function fit to stars is comparable to ours, while similar fits to DM tracers are significantly biased in contrast to our fits. In line with previous studies, the mass bias is much reduced near the tracer half-mass radius., Comment: MNRAS accepted. minor changes to match published version

We develop two general methods to infer the gravitational potential of a system using steady-state tracers, i.e., tracers with a time-independent phase-space distribution. Combined with the phase-space continuity equation, the time independence implies a universal Orbital Probability Density Function (oPDF) $\mathrm{d} P(\lambda|{\rm orbit})\propto \mathrm{d} t$, where $\lambda$ is the coordinate of the particle along the orbit. The oPDF is equivalent to Jeans theorem, and is the key physical ingredient behind most dynamical modelling of steady-state tracers. In the case of a spherical potential, we develop a likelihood estimator that fits analytical potentials to the system, and a non-parametric method ("phase-mark") that reconstructs the potential profile, both assuming only the oPDF. The methods involve no extra assumptions about the tracer distribution function and can be applied to tracers with any arbitrary distribution of orbits, with possible extension to non-spherical potentials. The methods are tested on Monte Carlo samples of steady-state tracers in dark matter haloes to show that they are unbiased as well as efficient. A fully documented \textsc{C/Python} code implementing our method is freely available at a GitHub repository linked from \url{http://icc.dur.ac.uk/data/#oPDF}., Comment: Added more references and discussions. MNRAS in press

We study the physical properties of a homogeneous sample of 157 optically-thick absorption line systems at redshifts ~1.8-4.4, selected from a high-dispersion spectroscopic survey of Lyman limit systems (LLSs). By means of multiple ionisation models and Bayesian techniques, we derive the posterior probability distribution functions for the density, metallicity, temperature, and dust content of the absorbing gas. We find that z>2 LLSs are highly ionised with ionisation parameters between -32 are characterised by a broad unimodal distribution over >4 orders of magnitude, with a peak at log Z/Zsun~-2. LLSs are metal poor, significantly less enriched than DLAs, with ~70% of the metallicity PDF below log Z/Zsun19 rapidly evolves with redshift, with a ten-fold increase between z~2.1-3.6 (~1.5 Gyr). Based on this sample, we find that LLSs at z=2.5-3.5 account for ~15% of all the metals produced by UV-selected galaxies. The implications for theories of cold gas accretion and metal ejection from galaxies are also discussed., Comment: 23 pages, 21 figures. MNRAS in press

We study the stellar-to-halo mass relation of central galaxies in the range 9.7 < log 10(M*/h- 2 M☉) < 11.7 and z < 0.4, obtained from a combined analysis of the Kilo Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey. We use ∼100 deg2 of KiDS data to study the lensing signal around galaxies for which spectroscopic redshifts and stellar masses were determined by GAMA. We show that lensing alone results in poor constraints on the stellar-to-halo mass relation due to a degeneracy between the satellite fraction and the halo mass, which is lifted when we simultaneously fit the stellar mass function. At M* > 5 × 1010 h- 2 M☉, the stellar mass increases with halo mass as {∼ }M_h^{0.25}. The ratio of dark matter to stellar mass has a minimum at a halo mass of 8 × 1011 h-1 M☉ with a value of M_h/M_{*}=56_{-10}^{+16} [h]. We also use the GAMA group catalogue to select centrals and satellites in groups with five or more members, which trace regions in space where the local matter density is higher than average, and determine for the first time the stellar-to-halo mass relation in these denser environments. We find no significant differences compared to the relation from the full sample, which suggests that the stellar-to-halo mass relation does not vary strongly with local density. Furthermore, we find that the stellar-to-halo mass relation of central galaxies can also be obtained by modelling the lensing signal and stellar mass function of satellite galaxies only, which shows that the assumptions to model the satellite contribution in the halo model do not significantly bias the stellar-to-halo mass relation. Finally, we show that the combination of weak lensing with the stellar mass function can be used to test the purity of group catalogues.

We explore the massive bluer star-forming population of the Sloan Digital Sky Survey (SDSS) III/BOSS CMASS DR11 galaxies at z > 0.55 to quantify their differences, in terms of redshiftspace distortions and large-scale bias, with respect to the luminous red galaxy sample. We perform a qualitative analysis to understand the significance of these differences and whether we can model and reproduce them in mock catalogues. Specifically, we measure galaxy clustering in CMASS on small and intermediate scales (0.1 ≲ r ≲ 50 h Mpc) by computing the two-point correlation function - both projected and redshift-space - of these galaxies, and a new statistic, Σ(π), able to separate the coherent and dispersed redshift-space distortion contributions and the large-scale bias. We interpret our clustering measurements by adopting a Halo Occupation Distribution (HOD) scheme that maps them on to high-resolution N-body cosmological simulations to produce suitable mock galaxy catalogues. The traditional HOD prescription can be applied to the red and the blue samples, independently, but this approach is unphysical since it allows the same mock galaxies to be either red or blue. To overcome this ambiguity, we modify the standard formulation and infer the red and the blue models by splitting the full mock catalogue into two complementary and non-overlapping submocks. This separation is performed by constraining the HOD with the observed CMASS red and blue galaxy fractions and produces reliable and accurate models. © 2016 The Authors., GF is supported by the Ministerio de Educacion y Ciencia of the Spanish Government through FPI grant AYA2010-2131-C02-01 and wishes to thank the Smithsonian Center for Astrophysics and the Harvard University, Astronomy Department, for the hospitality during the creation of this work. GF acknowledges support from the Spanish Government through EEBB-I-13-07167 and EEBB-I-12-05220 grants. GF, FP, and CC acknowledge support from the Spanish MICINNs Consolider-Ingenio 2010 Programme under grant MultiDark CSD2009-00064, MINECO Centro de Excelencia Severo Ochoa Programme under grant SEV-2012-0249, and MINECO grant AYA2014-60641-C2-1-P. Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation and the U.S. Department of Energy.

Current and future large redshift surveys, as the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), will use emission-line galaxies (ELGs) to probe cosmological models by mapping the large-scale structure of the Universe in the redshift range 0.6 < z < 1.7. With current data, we explore the halo-galaxy connection by measuring three clustering properties of g-selected ELGs as matter tracers in the redshift range 0.6 < z < 1: (i) the redshift-space two-point correlation function using spectroscopic redshifts from the BOSS ELG sample and VIPERS; (ii) the angular two-point correlation function on the footprint of the CFHT-LS; (iii) the galaxy-galaxy lensing signal around the ELGs using the CFHTLenS. We interpret these observations by mapping them on to the latest high-resolution MultiDark Planck N-body simulation, using a novel (Sub)Halo-Abundance Matching technique that accounts for the ELG incompleteness. ELGs at z ~ 0.8 live in haloes of (1 ± 0.5) × 10 hM⊙ and 22.5 ± 2.5 per cent of them are satellites belonging to a larger halo. The halo occupation distribution of ELGs indicates that we are sampling the galaxies in which stars form in the most efficient way, according to their stellar-to-halo mass ratio. © 2016 The Authors., GF is supported by the Ministerio de Educacion y Ciencia of the Spanish Government through FPI grant AYA2010-2131-C02-01. JC acknowledges financial support from MINECO (Spain) under project number AYA2012-31101. GY acknowledges financial support from MINECO (Spain) under project number AYA2012-31101 and grant FPA2012-34694. EJ acknowledges the support of CNRS, and the Labex OCEVU. SEN acknowledges support by the Deutsche Forschungsgemeinschaft under the grant NU 332/21. GF, FP, SART, AK, SN and CC acknowledge financial support from the Spanish MICINN Consolider-Ingenio 2010 Programme under grant MultiDark CSD2009-00064, MINECO Centro de Excelencia Severo Ochoa Programme under grant SEV-2012-0249, and MINECO grant AYA2014-60641-C2-1-P. GF, JC and FP wish to thank the Lawrence Berkeley National Laboratory for the hospitality during the creation of this work. FP acknowledges the spanish MEC 'Salvador de Madariaga' program, Ref. PRX14/00444.

Using data from four deep fields (COSMOS, AEGIS, ECDFS, and CDFN), we study the correlation between the position of galaxies in the star formation rate (SFR) versus stellar mass plane and local environment at $z10^{10.4-10.6}$ $M_{\odot}$), across all environments. At high redshift ( $0.5, Comment: 14 pages, 14 figures, accepted for publication in MNRAS

We use a large sample of isolated dark matter halo pairs drawn from cosmological N-body simulations to identify candidate systems whose kinematics match that of the Local Group of Galaxies (LG). We find, in agreement with the "timing argument" and earlier work, that the separation and approach velocity of the Milky Way (MW) and Andromeda (M31) galaxies favour a total mass for the pair of $\sim 5\times 10^{12} \,M_{\odot}$. A mass this large, however, is difficult to reconcile with the small relative tangential velocity of the pair, as well as with the small deceleration from the Hubble flow observed for the most distant LG members. Halo pairs that match these three criteria have average masses a factor of $\sim 2$ times smaller than suggested by the timing argument, but with large dispersion. Guided by these results, we have selected $12$ halo pairs with total mass in the range $1.6$-$3.6 \times 10^{12}\,M_{\odot}$ for the APOSTLE project (A Project Of Simulating The Local Environment), a suite of hydrodynamical resimulations at various numerical resolution levels (reaching up to $\sim10^{4}\,M_{\odot}$ per gas particle) that use the subgrid physics developed for the EAGLE project. These simulations reproduce, by construction, the main kinematics of the MW-M31 pair, and produce satellite populations whose overall number, luminosities, and kinematics are in good agreement with observations of the MW and M31 companions. The APOSTLE candidate systems thus provide an excellent testbed to confront directly many of the predictions of the $\Lambda$CDM cosmology with observations of our local Universe., Comment: replaced with the accepted version

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., Comment: 19 pages, 14 figures, MNRAS in press. Comments welcome

Multifield models of inflation predict an inequality between the amplitude τNL of the collapsed limit of the four-point correlator of the primordial curvature perturbation and the amplitude fNL of the squeezed limit of its three-point correlator. While a convincing detection of non-Gaussianity through the squeezed limit of the three-point correlator would rule out all single-field models, a robust confirmation or disproval of the inequality between τNL and fNL would provide crucial information about the validity of multifield models of inflation. In this paper, we discuss to which extent future measurements of the scale dependence of halo bias can test multifield inflationary scenarios. The strong degeneracy between the effect of a non-vanishing fNL and τNL on halo bias can be broken by considering multiple tracer populations of the same surveyed volume. If haloes down to 1013 M⊙ h−1 are resolved in a survey of volume 25 (Gpc3 h−3), then testing multifield models of inflation at the 3σ level would require, for instance, a detection of τNL at the level of τNL ∼ 105 given a measurement of a local bispectrum with amplitude fNL ∼ 10. However, we find that disproving multifield models of inflation with measurements of the non-Gaussian bias only will be very challenging, unless |fNL| ≳ 80 and one can achieve a halo mass resolution of ∼1010 M⊙ h−1

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society ©2011 RAS © 2011 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved, The detection of extremely massive clusters at z > 1 such as SPT-CLJ0546-5345, SPT-CLJ2106-5844 and XMMUJ2235.3-2557 has been considered by some authors as a challenge to the standard Λ cold dark matter cosmology. In fact, assuming Gaussian initial conditions, the theoretical expectation of detecting such objects is as low as ≤1 per cent. In this paper we discuss the probability of the existence of such objects in the light of the vector dark energy paradigm, showing by means of a series of N-body simulations that chances of detection are substantially enhanced in this non-standard framework, EC is supported by the MareNostrum project funded by the Spanish Ministerio de Ciencia e Innovación (MICINN) under grant no. AYA2009-13875-C03-02 and MultiDark Consolider project under grant CSD2009-00064. AK acknowledges support by the MICINN’s Ramón y Cajal programme as well as the grants AYA 2009-13875-C03-02, AYA2009-12792- C03-03, CSD2009- 00064, and CAM S2009/ESP-1496. GY acknowledges support from MICINN’s grants AYA2009- 13875-C03-02 and FPA2009- 08958. JBJ is supported by the Ministerio de Educación under the post-doctoral contract EX2009-0305 and also wishes to acknowledge support from the Norwegian Research Council under the YGGDRASIL programme 2009-2010 and the NILS mobility project grant UCM-EEA-ABEL-03-2010. We also acknowledge support from MICINN (Spain) project numbers FIS 2008-01323, FPA 2008-00592 and CAM/UCM 910309

We use the Millennium Simulation, a large, high resolution N-body simulation of the evolution of structure in a LambdaCDM cosmology, to study the properties and fate of substructures within a large sample of dark matter haloes. We find that the subhalo mass function departs significantly from a power law at the high mass end. We also find that the radial and angular distributions of substructures depend on subhalo mass. In particular, high mass subhaloes tend to be less radially concentrated and to have angular distributions closer to the direction perpendicular to the spin of the host halo than their less massive counterparts. We find that mergers between subhaloes occur. These tend to be between substructures that were already dynamically associated before accretion into the main halo. For subhaloes larger than 0.001 times the mass of the host halo, it is more likely that the subhalo will merge with the central or main subhalo than with another subhalo larger than itself. For lower masses, subhalo-subhalo mergers become equally likely to mergers with the main subhalo. Our results have implications for the variation of galaxy properties with environment and for the treatment of mergers in galaxy formation models., 14 pages, 11 figures. Submitted to MNRAS

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society ©: 2014 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved, Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive method for observing the build-up of cosmological structure, which depends both on the expansion rate of the Universe and our theory of gravity. In preparation for analysis of redshiftspace distortions from the Baryon Oscillation Spectroscopic Survey (BOSS) final data release, we compare a number of analytic and phenomenological models, specified in configuration space, to mock catalogues derived in different ways from several N-body simulations. The galaxies in each mock catalogue have properties similar to those of the higher redshift galaxies measured by BOSS but differ in the details of how small-scale velocities and halo occupancy are determined. We find that all of the analytic models fit the simulations over a limited range of scales while failing at small scales.We discuss which models are most robust and on which scales they return reliable estimates of the rate of growth of structure: we find that models based on some form of resummation can fit our N_body data for BOSS-like galaxies above 30 h-1 Mpc well enough to return unbiased parameter estimates, The MultiDark Database used in this paper and the web application providing online access to it were constructed as part of the activities of the German Astrophysical Virtua lObservatory as result of a collaboration between the Leibniz-Institute forAstrophysics Potsdam (AIP) and the Spanish MultiDark Consolider Project CSD2009-00064

© 2015 The Authors. We investigate the effect of self-interacting dark matter (SIDM) on the density profiles of Vmax ≃ 40km s-1 isolated dwarf dark matter haloes-the scale of relevance for the too big to fail problem (TBTF)-using very high resolution cosmological zoom simulations. Each halo has millions of particles within its virial radius. We find that SIDM models with cross-sections per unitmass spanning the range s/m=0.5-50 cm2 g-1 alleviate TBTF and produce constantdensity cores of size 300-1000 pc, comparable to the half-light radii of M* ~ 105-7M⊙ dwarfs. The largest, lowest density cores develop for cross-sections in the middle of this range, σ/m ~ 5-10 cm2 g-1. Our largest SIDM cross-section run (σ/m = 50 cm2 g-1) develops a slightly denser core owing to mild core-collapse behaviour, but it remains less dense than the cold dark matter case and retains a constant-density core profile. Our work suggests that SIDM cross-sections as large or larger than 50 cm2 g-1 remain viable on velocity scales of dwarf galaxies (vrms ~ 40 km s-1). The range of SIDM cross-sections that alleviate TBTF and the cusp/core problem spans at least two orders of magnitude and therefore need not be particularly fine-tuned.

We investigate the properties of dark matter haloes and subhaloes in an $f(R)$ gravity model with $|f_{R0}|=10^{-6}$, using a very high-resolution N-body simulation. The model is a borderline between being cosmologically interesting and yet still consistent with current data. We find that the halo mass function in this model has a maximum 20% enhancement compared with the $\Lambda$CDM predictions between $z=1$ and $z=0$. Because of the chameleon mechanism which screens the deviation from standard gravity in dense environments, haloes more massive than $10^{13}h^{-1}M_\odot$ in this $f(R)$ model have very similar properties to haloes of similar mass in $\Lambda$CDM, while less massive haloes, such as that of the Milky Way, can have steeper inner density profiles and higher velocity dispersions due to their weaker screening. The halo concentration is remarkably enhanced for low-mass haloes in this model due to a deepening of the total gravitational potential. Contrary to the naive expectation, the halo formation time $z_f$ is later for low-mass haloes in this model, a consequence of these haloes growing faster than their counterparts in $\Lambda$CDM at late times and the definition of $z_f$. Subhaloes, especially those less massive than $10^{11}h^{-1}M_\odot$, are substantially more abundant in this $f(R)$ model for host haloes less massive than $10^{13}h^{-1}M_\odot$. We discuss the implications of these results for the Milky Way satellite abundance problem. Although the overall halo and subhalo properties in this borderline $f(R)$ model are close to their $\Lambda$CDM predictions, our results suggest that studies of the Local Group and astrophysical systems, aided by high-resolution simulations, can be valuable for further tests of it., Comment: 14 pages, 9 figures. Accepted to be published in MNRAS after minor corrections

We use weak gravitational lensing to measure mean mass profiles around Locally Brightest Galaxies (LBGs). These are selected from the SDSS/DR7 spectroscopic and photometric catalogues to be brighter than any neighbour projected within 1.0 Mpc and differing in redshift by $ 83\%$) are expected to be the central galaxies of their dark matter halos. Previous stacking analyses have used this LBG sample to measure mean Sunyaev-Zeldovich flux and mean X-ray luminosity as a function of LBG stellar mass. In both cases, a simulation of the formation of the galaxy population was used to estimate effective halo mass for LBGs of given stellar mass, allowing the derivation of scaling relations between the gas properties of halos and their mass. By comparing results from a variety of simulations to our lensing data, we show that this procedure has significant model dependence reflecting: (i) the failure of any given simulation to reproduce observed galaxy abundances exactly; (ii) a dependence on the cosmology underlying the simulation; and (iii) a dependence on the details of how galaxies populate halos. We use our lensing results to recalibrate the scaling relations, eliminating most of this model dependence and explicitly accounting both for residual modelling uncertainties and for observational uncertainties in the lensing results. The resulting scaling relations link the mean gas properties of dark halos to their mass over an unprecedentedly wide range, $10^{12.5}, Comment: 22 pages, 11 figures, accepted by MNRAS

The evolution of the metal content of galaxies and its relations to other global properties [such as total stellar mass (M*), circular velocity, star formation rate (SFR), halo mass, etc.] provides important constraints on models of galaxy formation. Here we examine the evolution of metallicity scaling relations of simulated galaxies in the Galaxies-Intergalactic Medium Interaction Calculation suite of cosmological simulations. We make comparisons to observations of the correlation of gas-phase abundances with M* (the mass-metallicity relation, MZR), as well as with both M* and SFR or gas mass fraction (the so-called 3D fundamental metallicity relations, FMRs). The simulated galaxies follow the observed local MZR and FMRs over an order of magnitude in M*, but overpredict the metallicity of massive galaxies (log M* > 10.5), plausibly due to inefficient feedback in this regime. We discuss the origin of the MZR and FMRs in the context of galactic outflows and gas accretion. We examine the evolution of mass-metallicity relations defined using different elements that probe the three enrichment channels (SNII, SNIa, and AGB stars). Relations based on elements produced mainly by SNII evolve weakly, whereas those based on elements produced preferentially in SNIa/AGB exhibit stronger evolution, due to the longer timescales associated with these channels. Finally, we compare the relations of central and satellite galaxies, finding systematically higher metallicities for satellites, as observed. We show this is due to the removal of the metal poor gas reservoir that normally surrounds galaxies and acts to dilute their gas-phase metallicity (via cooling/accretion onto the disk), but is lost due to ram pressure stripping for satellites., 18 pages, 13 figures. Accepted for publication in MNRAS

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society ©: 2015 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved, Dark matter cores within galaxy haloes can be formed by energy feedback from star-forming regions: an energy balance suggests that the maximum core formation efficiency arises in galaxies with Mstar ~ 108.5M⊙. We show that a model population of galaxies, in which the density profile has been modified by such baryonic feedback, is able to explain the observed galaxy velocity function and Tully-Fisher relations significantly better than a model in which a universal cuspy density profile is assumed. Alternative models, namely warm or self-interacting dark matter, also provide a better match to these observed relations than a universal profile model does, but make different predictions for how halo density profiles vary with mass compared to the baryonic feedback case. We propose that the expected signatures of the mass dependence of core formation generated by baryonic feedback, CB thanks the MICINN (Spain) for the financial support through the MINECO grant AYA2012-31101 and the Ramon y Cajal program. ADC is supported by the DARK independent fellowship program.

We cross-match galaxy cluster candidates selected via their Sunyaev–Zel'dovich effect (SZE) signatures in 129.1 deg[superscript 2] of the South Pole Telescope 2500d SPT-SZ survey with optically identified clusters selected from the Dark Energy Survey science verification data. We identify 25 clusters between 0.1 ≲ z ≲ 0.8 in the union of the SPT-SZ and redMaPPer (RM) samples. RM is an optical cluster finding algorithm that also returns a richness estimate for each cluster. We model the richness λ-mass relation with the following function 〈ln λ|M500〉 ∝ Bλln M500 + Cλln E(z) and use SPT-SZ cluster masses and RM richnesses λ to constrain the parameters. We find Bλ=1.14[superscript +0.21][subscript −0.18 and Cλ=0.73[superscript +0.77][subscript −0.75]. The associated scatter in mass at fixed richness is σlnM|λ=0.18[superscript +0.08][subscript −0.05] at a characteristic richness λ = 70. We demonstrate that our model provides an adequate description of the matched sample, showing that the fraction of SPT-SZ-selected clusters with RM counterparts is consistent with expectations and that the fraction of RM-selected clusters with SPT-SZ counterparts is in mild tension with expectation. We model the optical-SZE cluster positional offset distribution with the sum of two Gaussians, showing that it is consistent with a dominant, centrally peaked population and a subdominant population characterized by larger offsets. We also cross-match the RM catalogue with SPT-SZ candidates below the official catalogue threshold significance ξ = 4.5, using the RM catalogue to provide optical confirmation and redshifts for 15 additional clusters with ξ ∈ [4, 4.5]., Kavli Foundation, Gordon and Betty Moore Foundation (Grant GBMF 947)

The detection of planar structures within the satellite systems of both the Milky Way (MW) and Andromeda (M31) has been reported as being in stark contradiction to the predictions of the standard cosmological model ($\Lambda$CDM). Given the ambiguity in defining a planar configuration, it is unclear how to interpret the low incidence of the MW and M31 planes in $\Lambda$CDM. We investigate the prevalence of satellite planes around galactic mass haloes identified in high resolution cosmological simulations. We find that planar structures are very common, and that ~10% of $\Lambda$CDM haloes have even more prominent planes than those present in the Local Group. While ubiquitous, the planes of satellite galaxies show a large diversity in their properties. This precludes using one or two systems as small scale probes of cosmology, since a large sample of satellite systems is needed to obtain a good measure of the object-to-object variation. This very diversity has been misinterpreted as a discrepancy between the satellite planes observed in the Local Group and $\Lambda$CDM predictions. In fact, ~10% of $\Lambda$CDM galactic haloes have planes of satellites that are as infrequent as the MW and M31 planes. The look-elsewhere effect plays an important role in assessing the detection significance of satellite planes and accounting for it leads to overestimating the significance level by a factor of 30 and 100 for the MW and M31 systems, respectively., Comment: 13 pages, 15 figures, +3 pages appendix, MNRAS accepted version

We present the first cross-correlation measurement between Sloan Digital Sky Survey (SDSS) Type 1 quasars and the cosmic infrared background (CIB) measured by Herschel. The distribution of the quasars at 0.15=1.4) is $11.1^{+1.6}_{-1.4}$, $7.1^{+1.6}_{-1.3}$ and $3.6^{+1.4}_{-1.0}$ mJy at 250, 350 and 500 microns, respectively, while the mean sub-mm flux densities of the DR9 quasars (=2.5) is $5.7^{+0.7}_{-0.6}$, $5.0^{+0.8}_{-0.7}$ and $1.8^{+0.5}_{-0.4}$ mJy. We find that the correlated sub-mm emission includes both the emission from satellite DSFGs in the same halo as the central quasar and the emission from DSFGs in separate halos (correlated with the quasar-hosting halo). The amplitude of the one-halo term is ~10 times smaller than the sub-mm emission of the quasars, implying the the satellites have a lower star-formation rate than the quasars. The satellite fraction for the DR7 quasars is $0.008^{+0.008}_{-0.005}$ and the host halo mass scale for the central and satellite quasars is $10^{12.36\pm0.87}$ M$_{\odot}$ and $10^{13.60\pm0.38}$ M$_{\odot}$, respectively. The satellite fraction of the DR9 quasars is $0.065^{+0.021}_{-0.031}$ and the host halo mass scale for the central and satellite quasars is $10^{12.29\pm0.62}$ M$_{\odot}$ and $10^{12.82\pm0.39}$ M$_{\odot}$, respectively. Thus, the typical halo environment of the SDSS Type 1 quasars is found to be similar to that of DSFGs, which supports the generally accepted view that dusty starburst and quasar activity are evolutionarily linked., Comment: 20 pages, 14 pages, submitted to MNRAS

We estimate the rate of dark matter scattering in collapsed structures throughout the history of the Universe. If the scattering cross-section is velocity-independent, then the canonical picture is correct that scatterings occur mainly at late times. The scattering rate peaks slightly at redshift z~6, and remains significant today. Half the scatterings occur after z~1, in structures more massive than 10^12 M_sun. Within a factor of two, these numbers are robust to changes in the assumed astrophysics, and the scatterings would be captured in cosmological simulations. However, for particle physics models with a velocity-dependent cross-section (as for Yukawa potential interactions via a massive mediator), the scattering rate peaks before z~20, in objects with mass less than 10^4 M_sun. These precise values are sensitive to the redshift-dependent mass-concentration relation and the small-scale cutoff in the matter power spectrum. In extreme cases, the qualitative effect of early interactions may be reminiscent of warm dark matter and strongly affect the subsequent growth of structure. However, these scatterings are being missed in existing cosmological simulations with limited mass resolution., v2 matches accepted MNRAS version

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society ©: 2015 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved, We present a simple method of generating approximate synthetic halo catalogues: HALOGEN. This method uses a combination of second-order Lagrangian Perturbation Theory (2LPT) in order to generate the large-scale matter distribution, analytical mass functions to generate halo masses, and a single-parameter stochastic model for halo bias to position haloes. HALOGEN represents a simplification of similar recently published methods. Our method is constrained to recover the two-point function at intermediate (10 h-1 Mpc, SA and JGB acknowledge financial support from the Spanish MINECO under grant FPA39684-C03-02 and Consolider-Ingenio “Physics of the Accelerating Universe (PAU)” (CSD2007-00060). They also acknowledge the support from the Spanish MINECO’s “Centro de Excelencia Severo Ochoa” Programme under Grant No. SEV-2012-0249. SA is also supported by a PhD FPI-fellowship from the Universidad Autónoma de Madrid. AK is supported by the Ministerio de Economía y Competitividad (MINECO) in Spain through grant AYA2012-31101 as well as the Consolider-Ingenio 2010 Programme of the Spanish Ministerio de Ciencia e Innovación (MICINN) under grant MultiDark CSD2009-00064. He also acknowledges support from the Australian Research Council (ARC) grants DP130100117 and DP140100198

We report the alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations. The combination of these state-of-the-art hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($11 < log_{10}(M_{200}/ [h^{-1}M_\odot]) < 15$), a wide radial range ($-2.3 < log_{10}(r/[h^{-1}Mpc ]) < 1.3$) and redshifts $0 < z < 1$. The shape parameters of the dark matter, stellar and hot gas distributions follow qualitatively similar trends: they become more aspherical (and triaxial) with increasing halo mass, radius and redshift. We measure the misalignment of the baryonic components (hot gas and stars) of galaxies with their host halo as a function of halo mass, radius, redshift, and galaxy type (centrals vs satellites and early- vs late-type). Overall, galaxies align well with the local distribution of the total (mostly dark) matter. However, the stellar distributions on galactic scales exhibit a median misalignment of about 45-50 degrees with respect to their host haloes. This misalignment is reduced to 25-30 degrees in the most massive haloes ($13 < log_{10}(M_{200}/ [h^{-1}M_\odot ]) < 15$). Half of the disc galaxies in the EAGLE simulations have a misalignment angle with respect to their host haloes larger than 40 degrees. We present fitting functions and tabulated values for the probability distribution of galaxy-halo misalignment to enable a straightforward inclusion of our results into models of galaxy formations based on purely collisionless N-body simulations., 19 pages, 16 figures. Accepted for publication in MNRAS. New figure added, e_rms values corrected, no changes to the conclusions

We present distributions of orbital parameters of infalling satellites of $\Lambda$CDM haloes in the mass range $10^{12}-10^{14}$M$_\odot$, which represent the initial conditions for the subsequent evolution of substructures within the host halo. We use merger trees constructed in a high resolution cosmological N-body simulation to trace satellite haloes, and identify the time of infall. We find signficant trends in the distribution of orbital parameters with both the host halo mass and the ratio of satellite-to-host halo masses. For all host halo masses, satellites whose infall mass is a larger fraction of the host halo mass have more eccentric, radially biased orbits. At fixed satellite-to-host halo mass ratio, high mass haloes are biased towards accreting satellites on slightly more radial orbits. To charactise the orbital distributions fully requires fitting the correlated bivariate distribution of two chosen orbital parameters (e.g. radial and tangential velocity or energy and angular momentum). We provide simple fits to one choice of the bivariate distributions, which when transformed faithfully, captures the behaviour of any of the projected one-dimensional distributions., Comment: 14 pages, 13 figures, accepted for publication in MNRAS

We present one thousand mock galaxy catalogues for the analysis of the Low Redshift Sample (LOWZ, effective redshift z ~ 10.32) of the Baryon Oscillation Spectroscopic Survey Data Releases 10 and 11. These mocks have been created following the PTHalos method of Manera13 et al. (2013) revised to include new developments. The main improvement is the introduction of a redshift dependence in the Halo Occupation Distribution in order to account for the change of the galaxy number density with redshift. These mock catalogues are used in the analyses of the LOWZ galaxy clustering by the BOSS collaboration., 10 pages, 8 figures

While wide-field surveys of M31 have revealed much substructure at large radii, understanding the nature and origin of this material is not straightforward from morphology alone. Using deep HST/ACS data, we have derived further constraints in the form of quantitative star formation histories (SFHs) for 14 inner halo fields which sample diverse substructures. In agreement with our previous analysis of colour-magnitude diagram morphologies, we find the resultant behaviours can be broadly separated into two categories. The SFHs of 'disc-like' fields indicate that most of their mass has formed since z~1, with one quarter of the mass formed in the last 5 Gyr. We find 'stream-like' fields to be on average 1.5 Gyr older, with, 14 pages, 10 figures (+12 figures in appendix). MNRAS, in press

Understanding the relationship between galaxies hosting active galactic nuclei (AGN) and the dark matter halos in which they reside is key to constraining how black-hole fueling is triggered and regulated. Previous efforts have relied on simple halo mass estimates inferred from clustering, weak gravitational lensing, or halo occupation distribution modeling. In practice, these approaches remain uncertain because AGN, no matter how they are identified, potentially live a wide range of halo masses with an occupation function whose general shape and normalization are poorly known. In this work, we show that better constraints can be achieved through a rigorous comparison of the clustering, lensing, and cross-correlation signals of AGN hosts to a fiducial stellar-to-halo mass relation (SHMR) derived for all galaxies. Our technique exploits the fact that the global SHMR can be measured with much higher accuracy than any statistic derived from AGN samples alone. Using 382 moderate luminosity X-ray AGN at z, 15 pages, accepted to MNRAS

We use the "Evolution and assembly of galaxies and their environments" (EAGLE) cosmological simulation to investigate the effect of baryons on the density profiles of rich galaxy clusters. We focus on EAGLE clusters with $M_{200}>10^{14}~M_\odot$ of which we have six examples. The central brightest cluster galaxies (BCGs) in the simulation have steep stellar density profiles, $��_*(r) \propto r^{-3}$. Stars dominate the mass density for $r, 13 pages, 8 figures, accepted for publication in MNRAS

© 2014 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society. We present the measurements and modelling of the projected and redshift-space clustering of CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Data Release 11. For a volume-limited luminous red galaxy sample in the redshift range of 0.48 < z < 0.55, we perform halo occupation distribution modelling of the small- and intermediate-scale (0.1-60 h-1Mpc) projected and redshift-space two-point correlation functions, with an accurate model built on high-resolution N-body simulations. To interpret the measured redshift-space distortions, the distribution of galaxy velocities must differ from that of the dark matter inside haloes of ~1013-1014h-1M⊙, i.e. the data require the existence of galaxy velocity bias. Most notably, central galaxies on average are not at rest with respect to the core of their host haloes, but rather move around it with a 1D velocity dispersion of 0.22-0.04+0.03times that of the dark matter, implying a spatial offset from the centre at the level of ≲1 per cent of the halo virial radius. The luminous satellite galaxies move more slowly than the dark matter, with velocities 0.86-0.03+0.08times those of the dark matter, which suggests that the velocity and spatial distributions of these satellites cannot both be unbiased. The constraints mainly arise from the Fingers-of-God effect at non-linear scales and the smoothing to the Kaiser effect in the translinear regime; the robustness of the results is demonstrated by a variety of tests. We discuss the implications of the existence of galaxy velocity bias for investigations of galaxy formation and cosmology.

We report results for the alignments of galaxies in the EAGLE and cosmo-OWLS simulations as a function of galaxy separation and halo mass. The combination of these hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($10.7, Comment: 13 pages, 12 figures. Accepted for publication in MNRAS