Source author record

John A. Peacock

John A. Peacock 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

48works
9topics
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

48 published item(s)

preprint2022arXiv

Accelerating astronomical and cosmological inference with Preconditioned Monte Carlo

We introduce Preconditioned Monte Carlo (PMC), a novel Monte Carlo method for Bayesian inference that facilitates efficient sampling of probability distributions with non-trivial geometry. PMC utilises a Normalising Flow (NF) in order to decorrelate the parameters of the distribution and then proceeds by sampling from the preconditioned target distribution using an adaptive Sequential Monte Carlo (SMC) scheme. The results produced by PMC include samples from the posterior distribution and an estimate of the model evidence that can be used for parameter inference and model comparison respectively. The aforementioned framework has been thoroughly tested in a variety of challenging target distributions achieving state-of-the-art sampling performance. In the cases of primordial feature analysis and gravitational wave inference, PMC is approximately 50 and 25 times faster respectively than Nested Sampling (NS). We found that in higher dimensional applications the acceleration is even greater. Finally, PMC is directly parallelisable, manifesting linear scaling up to thousands of CPUs.

preprint2022arXiv

An empirical measurement of the Halo Mass Function from the combination of GAMA DR4, SDSS DR12, and REFLEX II data

We construct the halo mass function (HMF) from the GAMA galaxy group catalogue over the mass range 10^12.7M_sol to 10^15.5M_sol, and find good agreement with the expectation from LambdaCDM. In comparison to previous studies, this result extends the mass range over which the HMF has now been measured over by an order of magnitude. We combine the GAMA DR4 HMF with similar data from the SDSS DR12 and REFLEX II surveys, and fit a four-parameter Murray-Robotham-Power (MRP) function, valid at z~0.1, yielding: a density normalisation of: log10 (phi Mpc^3)=-3.96[+0.55,-0.82], a high mass turn-over of: log10(M/M_sol)=14.13[+0.43,-0.40], a low mass power law slope of: alpha=-1.68[+0.21,-0.24] , and a high mass softening parameter of: beta= 0.63[+0.25,-0.11]. If we fold in the constraint on Omega_M from Planck 2018 Cosmology, we are able to reduce these uncertainties further, but this relies on the assumption that the power-law trend can be extrapolated from 10^12.7M_sol to zero mass. Throughout, we highlight the effort needed to improve on our HMF measurement: improved halo mass estimates that do not rely on calibration to simulations; reduced halo mass uncertainties needed to mitigate the strong Eddington Bias that arises from the steepness of the HMF low mass slope; and deeper wider area spectroscopic surveys. To our halo mass limit of 10^12.7 M_sol, we are directly resolving (`seeing') 41+/-5 per cent of the total mass density, i.e. Omega_[M>12.7]=0.128+/-0.016, opening the door for the direct construction of 3D dark matter mass maps at Mpc resolution.

preprint2022arXiv

Extended Hernquist-Springel formalism for cosmic star formation

We present a revised and extended version of the analytic model for cosmic star formation originally given by Hernquist & Springel in 2003. The key assumption of this formalism is that star formation proceeds from cold gas, at a rate that is limited by an internal consumption timescale at early times, or by the rate of generation of gas via cooling at late times. These processes are analysed as a function of the mass of dark matter haloes and integrated over the halo population. We modify this approach in two main ways to make it more general: (1) halo collapse times are included explicitly, so that the behaviour is physically reasonable at late times; (2) allowance is made for a mass-dependent baryon fraction in haloes, which incorporates feedback effects. This model reproduces the main features of the observed baryonic Tully-Fisher relationship, and is consistent with observational estimates of the baryon mass fraction in the intergalactic medium. With minimal adjustment of parameters, our approach reproduces the observed history of cosmic star formation within a factor of two over the redshift range $0 < z < 10$. This level of agreement is comparable to that achieved by state-of-the-art cosmological simulations. Our simplified apparatus has pedagogical value in illuminating the results of such detailed calculations, and also serves as a means for rapid approximate exploration of non-standard cosmological models.

preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Data Release 4 and the z < 0.1 total and z < 0.08 morphological galaxy stellar mass functions

In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z<0.25) galaxy population. DR4 includes 32 value-added tables or Data Management Units (DMUs) that provide a number of measured and derived data products including GALEX, ESO KiDS, ESO VIKING, WISE and Herschel Space Observatory imaging. Within this release, we provide visual morphologies for 15330 galaxies to z<0.08, photometric redshift estimates for all 18million objects to r_KIDS~25mag, and stellar velocity dispersions for 111830 galaxies. We conclude by deriving the total galaxy stellar mass function (GSMF) and its sub-division by morphological class (elliptical, compact-bulge and disc, diffuse-bulge and disc, and disc only). This extends our previous measurement of the total GSMF down to 10^6.75 M_sol h^-2_70 and we find a total stellar mass density of rho_*=(2.97+/-0.04)x10^8 M_sol h_70 Mpc^-3 or Omega_*=(2.17+/-0.03)x10^-3 h^-1_70. We conclude that at z<0.1, the Universe has converted 4.9+/-0.1 per cent of the baryonic mass implied by Big Bang Nucleosynthesis into stars that are gravitationally bound within the galaxy population.

preprint2022arXiv

pocoMC: A Python package for accelerated Bayesian inference in astronomy and cosmology

pocoMC is a Python package for accelerated Bayesian inference in astronomy and cosmology. The code is designed to sample efficiently from posterior distributions with non-trivial geometry, including strong multimodality and non-linearity. To this end, pocoMC relies on the Preconditioned Monte Carlo algorithm which utilises a Normalising Flow in order to decorrelate the parameters of the posterior. It facilitates both tasks of parameter estimation and model comparison, focusing especially on computationally expensive applications. It allows fitting arbitrary models defined as a log-likelihood function and a log-prior probability density function in Python. Compared to popular alternatives (e.g. nested sampling) pocoMC can speed up the sampling procedure by orders of magnitude, cutting down the computational cost substantially. Finally, parallelisation to computing clusters manifests linear scaling.

preprint2022arXiv

Understanding the relation between thermal Sunyaev-Zeldovich decrement and halo mass using the SIMBA and TNG simulations

The relation between the integrated thermal Sunyaev-Zeldovich (tSZ) $y$-decrement versus halo mass ($Y$--$M$) can potentially constrain galaxy formation models, if theoretical and observational systematics can be properly assessed. We investigate the $Y$--$M$ relation in the SIMBA and IllustrisTNG-100 cosmological hydrodynamic simulations, quantifying the effects of feedback, line-of-sight projection, and beam convolution. We find that SIMBA's AGN jet feedback generates strong deviations from self-similar expectations for the $Y$--$M$ relation, especially at $M_{500}<10^{13}M_{\odot}$. In SIMBA, this is driven by suppressed in-halo $y$ contributions owing to lowered halo baryon fractions. IllustrisTNG results more closely resemble SIMBA without jets. Projections of line-of-sight structures weaken these model differences slightly, but they remain significant -- particularly at group and lower halo masses. In contrast, beam smearing at $\textit{Planck}$ resolution makes the models indistinguishable, and both models appear to agree well with $\textit{Planck}$ data down to the lowest masses probed. We show that the arcminute resolution expected from forthcoming facilities would retain the differences between model predictions, and thereby provide strong constraints on AGN feedback.

preprint2020arXiv

Calibration of a star formation and feedback model for cosmological simulations with Enzo

We present results from seventy-one zoom simulations of a Milky Way-sized (MW) halo, exploring the parameter space for a widely-used star formation and feedback model in the {\tt Enzo} simulation code. We propose a novel way to match observations, using functional fits to the observed baryon makeup over a wide range of halo masses. The model MW galaxy is calibrated using three parameters: the star formation efficiency $\left(f_*\right)$, the efficiency of thermal energy from stellar feedback $\left(ε\right)$ and the region into which feedback is injected $\left(r\ {\rm and}\ s\right)$. We find that changing the amount of feedback energy affects the baryon content most significantly. We then identify two sets of feedback parameter values that are both able to reproduce the baryonic properties for haloes between $10^{10}\,\mathrm{M_\odot}$ and $10^{12}\,\mathrm{M_\odot}$. We can potentially improve the agreement by incorporating more parameters or physics. If we choose to focus on one property at a time, we can obtain a more realistic halo baryon makeup. We show that the employed feedback prescription is insensitive to dark matter mass resolution between $10^5\,{\rm M_\odot}$ and $10^7\,{\rm M_\odot}$. Contrasting both star formation criteria and the corresponding combination of optimal feedback parameters, we also highlight that feedback is self-consistent: to match the same baryonic properties, with a relatively higher gas to stars conversion efficiency, the feedback strength required is lower, and vice versa. Lastly, we demonstrate that chaotic variance in the code can cause deviations of approximately 10\% and 25\% in the stellar and baryon mass in simulations evolved from identical initial conditions.

preprint2020arXiv

Multi-tracer extension of the halo model: probing quenching and conformity in eBOSS

We develop a new Multi-Tracer Halo Occupation Distribution (\texttt{MTHOD}) framework for the galaxy distribution and apply it to the extended Baryon Oscillation Spectroscopic Survey (eBOSS) final data between $z=0.7-1.1$. We obtain a best fit \mthod\, for each tracer and describe the host halo properties of these galaxies. The mean halo masses for LRGs, ELGs and QSOs are found to be $1.9 \times 10^{13} \msolaroh$, $1.1 \times 10^{12} \msolaroh$ and $5 \times 10^{12} \msolaroh$ respectively in the eBOSS data. We use the \texttt{MTHOD} framework to create mock galaxy catalogues and predict auto- and cross-correlation functions for all the tracers. Comparing these results with data, we investigate galactic conformity, the phenomenon whereby the properties of neighbouring galaxies are mutually correlated in a manner that is not captured by the basic halo model. We detect \textsl{1-halo} conformity at more than 3$σ$ statistical significance, while obtaining upper limit on \textsl{2-halo} conformity. We also look at the environmental dependence of the galaxy quenching efficiency and find that halo mass driven quenching successfully explains the behaviour in high density regions, but it fails to describe the quenching efficiency in low density regions. In particular, we show that the quenching efficiency in low density filaments is higher in the observed data, as compared to the prediction of the \mthod\ with halo mass driven quenching. The mock galaxy catalogue constructed in this paper is publicly available on https://www.roe.ac.uk/~salam/MTHOD/ .

preprint2020arXiv

Primordial black hole merger rates: distributions for multiple LIGO observables

We have calculated the detectable merger rate of primordial black holes, as a function of the redshift, as well as the binary's mass ratio, total mass and chirp mass (observables that have not previously been explored in great detail for PBHs). We consider both the current and design sensitivity of LIGO and five different primordial black hole mass functions, as well as showing a comparison to a predicted astrophysical black hole merger rate. We show that the empirical preference for nearly equal-mass binaries in current LIGO/Virgo data can be consistent with a PBH hypothesis once observational selection effects are taken into account. However, current data do exclude some PBH mass distributions, and future data may be able to rule out the possibility that all observed BH mergers had a primordial origin.

preprint2020arXiv

Protoclusters at z=5.7: A view from the MultiDark galaxies

Protoclusters, which will yield galaxy clusters at lower redshift, can provide valuable information on the formation of galaxy clusters. However, identifying progenitors of galaxy clusters in observations is not an easy task, especially at high redshift. Different priors have been used to estimate the overdense regions that are thought to mark the locations of protoclusters. In this paper, we use mimicked Ly$α$-emitting galaxies at $z=5.7$ to identify protoclusters in the MultiDark galaxies, which are populated by applying three different semi-analytic models to the 1 $Gpc h^{-1}$ MultiDark Planck2 simulation. To compare with observational results, we extend the criterion 1 (a Ly$α$ luminosity limited sample), to criterion 2 (a match to the observed mean galaxy number density). To further statistically study the finding efficiency of this method, we enlarge the identified protocluster sample (criterion 3) to about 3500 at $z=5.7$ and study their final mass distribution. The number of overdense regions and their selection probability depends on the semi-analytic models and strongly on the three selection criteria (partly by design). The protoclusters identified with criterion 1 are associated with a typical final cluster mass of $2.82\pm0.92 \times 10^{15} M_\odot$ which is in agreement with the prediction (within $\pm 1 σ$) of an observed massive protocluster at $z=5.7$. Identifying more protoclusters allows us to investigate the efficiency of this method, which is more suitable for identifying the most massive clusters: completeness ($\mathbb{C}$) drops rapidly with decreasing halo mass. We further find that it is hard to have a high purity ($\mathbb{P}$) and completeness simultaneously.

preprint2020arXiv

The Three Hundred Project: the stellar and gas profiles

Using the catalogues of galaxy clusters from The Three Hundred project, modelled with both hydrodynamic simulations, (Gadget-X and Gadget-MUSIC), and semi-analytic models (SAMs), we study the scatter and self-similarity of the profiles and distributions of the baryonic components of the clusters: the stellar and gas mass, metallicity, the stellar age, gas temperature, and the (specific) star formation rate. Through comparisons with observational results, we find that the shape and the scatter of the gas density profiles matches well the observed trends including the reduced scatter at large radii which is a signature of self-similarity suggested in previous studies. One of our simulated sets, Gadget-X, reproduces well the shape of the observed temperature profile, while Gadget-MUSIC has a higher and flatter profile in the cluster centre and a lower and steeper profile at large radii. The gas metallicity profiles from both simulation sets, despite following the observed trend, have a relatively lower normalisation. The cumulative stellar density profiles from SAMs are in better agreement with the observed result than both hydrodynamic simulations which show relatively higher profiles. The scatter in these physical profiles, especially in the cluster centre region, shows a dependence on the cluster dynamical state and on the cool-core/non-cool-core dichotomy. The stellar age, metallicity and (s)SFR show very large scatter, which are then presented in 2D maps. We also do not find any clear radial dependence of these properties. However, the brightest central galaxies have distinguishable features compared to the properties of the satellite galaxies.

preprint2020arXiv

Tomographic measurement of the intergalactic gas pressure through galaxy-tSZ cross-correlations

We cross-correlate maps of the thermal Sunyaev-Zeldovich (tSZ) Compton-$y$ parameter published by Planck with the projected distribution of galaxies in a set of low-redshift tomographic bins. We use the nearly full-sky 2MASS Photometric Redshift and WISE $\times$ SuperCOSMOS public catalogues, covering the redshift range $z\lesssim0.4$. Our measurements allow us to place constraints on the redshift dependence of the mass-observable relation for tSZ cluster count analyses in terms of the so-called 'hydrostatic mass bias' parameter $1-b_{\rm H}$. These results can also be interpreted as measurements of the bias-weighted average gas pressure $\langle bP_e\rangle$ as a function of redshift, a quantity that can be related to the thermodynamics of gas inside haloes and used to constrain energy injection processes. We measure $1-b_{\rm H}$ with $\sim6\%$ precision in 6 equispaced redshift bins, and find no evidence for a redshift-dependent mass bias parameter, in agreement with previous analyses. Our mean value of $1-b_{\rm H} = 0.75\pm0.03$ is also in good agreement with the one estimated by the joint analysis of Planck cluster counts and CMB anisotropies calibrated with CMB lensing. Our measurements of $\langle bP_e\rangle$, at the level of $\sim10\%$ in each bin, are the most stringent constraints on the redshift dependence of this parameter to date, and agree well both with previous measurements and with theoretical expectations from shock-heating models.

preprint2016arXiv

Dependence of GAMA galaxy halo masses on the cosmic web environment from 100 square degrees of KiDS weak lensing data

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 < 0.263, from the spectroscopic Galaxy And Mass Assembly (GAMA) survey, using ~100 square degrees of overlapping data from the Kilo-Degree Survey (KiDS). In each of the four cosmic environments we model the contributions from group centrals, satellites and neighbouring groups to the stacked galaxy-galaxy lensing profiles. After correcting the lens samples for differences in the stellar mass distribution, we find no dependence of the average halo mass of central galaxies on their cosmic environment. We do find a significant increase in the average contribution of neighbouring groups to the lensing profile in increasingly dense cosmic environments. We show, however, that the observed effect can be entirely attributed to the galaxy density at much smaller scales (within 4 Mpc/h), which is correlated with the density of the cosmic environments. Within our current uncertainties we find no direct dependence of galaxy halo mass on their cosmic environment.

preprint2016arXiv

Galaxy and Mass Assembly (GAMA): Redshift Space Distortions from the Clipped Galaxy Field

We present the first cosmological measurement derived from a galaxy density field subject to a `clipping' transformation. By enforcing an upper bound on the galaxy number density field in the Galaxy and Mass Assembly survey (GAMA), contributions from the nonlinear processes of virialisation and galaxy bias are greatly reduced. This leads to a galaxy power spectrum which is easier to model, without calibration from numerical simulations. We develop a theoretical model for the power spectrum of a clipped field in redshift space, which is exact for the case of anisotropic Gaussian fields. Clipping is found to extend the applicability of the conventional Kaiser prescription by more than a factor of three in wavenumber, or a factor of thirty in terms of the number of Fourier modes. By modelling the galaxy power spectrum on scales k < 0.3 h/Mpc and density fluctuations $δ_g < 4$ we measure the normalised growth rate $fσ_8(z = 0.18) = 0.29 \pm 0.10$.

preprint2016arXiv

Group-galaxy correlations in redshift space as a probe of the growth of structure

We investigate the use of the cross-correlation between galaxies and galaxy groups to measure redshift-space distortions (RSD) and thus probe the growth rate of cosmological structure. This is compared to the classical approach based on using galaxy auto-correlation. We make use of realistic simulated galaxy catalogues that have been constructed by populating simulated dark matter haloes with galaxies through halo occupation prescriptions. We adapt the classical RSD dispersion model to the case of the group-galaxy cross-correlation function and estimate the RSD parameter $β$ by fitting both the full anisotropic correlation function $ξ(r_p,π)$ and its multipole moments. In addition, we define a modified version of the latter statistics by truncating the multipole moments to exclude strongly non-linear distortions at small transverse scales. We fit these three observable quantities in our set of simulated galaxy catalogues and estimate statistical and systematic errors on $β$ for the case of galaxy-galaxy, group-group, and group-galaxy correlation functions. When ignoring off-diagonal elements of the covariance matrix in the fitting, the truncated multipole moments of the group-galaxy cross-correlation function provide the most accurate estimate, with systematic errors below 3% when fitting transverse scales larger than $10Mpc/h$. Including the full data covariance enlarges statistical errors but keep unchanged the level of systematic error. Although statistical errors are generally larger for groups, the use of group-galaxy cross-correlation can potentially allow the reduction of systematics while using simple linear or dispersion models.

preprint2016arXiv

Redshift-space distortions around voids

We have derived estimators for the linear growth rate of density fluctuations using the cross-correlation function of voids and haloes in redshift space, both directly and in Fourier form. In linear theory, this cross-correlation contains only monopole and quadrupole terms. At scales greater than the void radius, linear theory is a good match to voids traced out by haloes in N-body simulations; small-scale random velocities are unimportant at these radii, only tending to cause small and often negligible elongation of the redshift-space cross-correlation function near its origin. By extracting the monopole and quadrupole from the cross-correlation function, we measure the linear growth rate without prior knowledge of the void profile or velocity dispersion. We recover the linear growth parameter $β$ to 9% precision from an effective volume of 3(Gpc/h)^3 using voids with radius greater than 25Mpc/h. Smaller voids are predominantly sub-voids, which may be more sensitive to the random velocity dispersion; they introduce noise and do not help to improve the measurement. Adding velocity dispersion as a free parameter allows us to use information at radii as small as half of the void radius. The precision on $β$ is reduced to approximately 5%. Contrary to the simple redshift-space distortion pattern in overdensities, voids show diverse shapes in redshift space, and can appear either elongated or flattened along the line of sight. This can be explained by the competing amplitudes of the local density contrast, plus the radial velocity profile and its gradient, with the latter two factors being determined by the cumulative density profile of voids. The distortion pattern is therefore determined solely by the void profile and is different for void-in-cloud and void-in-void. This diversity of redshift-space void morphology complicates measurements of the Alcock-Paczynski effect using voids.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2016arXiv

The lensing and temperature imprints of voids on the Cosmic Microwave Background

We have searched for the signature of cosmic voids in the CMB, in both the Planck temperature and lensing-convergence maps; voids should give decrements in both. We use zobov voids from the DR12 SDSS CMASS galaxy sample. We base our analysis on N-body simulations, to avoid a posteriori bias. For the first time, we detect the signature of voids in CMB lensing: the significance is $3.2σ$, close to $Λ$CDM in both amplitude and projected density-profile shape. A temperature dip is also seen, at modest significance ($2.3σ$), with amplitude about 6 times the prediction. This temperature signal is induced mostly by voids with radius between 100 and 150 Mpc/h, while the lensing signal is mostly contributed by smaller voids -- as expected; lensing relates directly to density, while ISW depends on gravitational potential. The void abundance in observations and simulations agree, as well. We also repeated the analysis excluding lower-significance voids: no lensing signal is detected, with an upper limit of about twice the $Λ$CDM prediction. But the mean temperature decrement now becomes non-zero at the $3.7σ$ level (similar to that found by Granett et al.), with amplitude about 20 times the prediction. However, the observed dependence of temperature on void size is in poor agreement with simulations, whereas the lensing results are consistent with $Λ$CDM theory. Thus, the overall tension between theory and observations does not favour non-standard theories of gravity, despite the hints of an enhanced amplitude for the ISW effect from voids.

preprint2016arXiv

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

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

Halo abundances within the cosmic web

We investigate the dependence of the mass function of dark-matter haloes on their environment within the cosmic web of large-scale structure. A dependence of the halo mass function on large-scale mean density is a standard element of cosmological theory, allowing mass-dependent biasing to be understood via the peak-background split. On the assumption of a Gaussian density field, this analysis can be extended to ask how the mass function depends on the geometrical environment: clusters, filaments, sheets and voids, as classified via the tidal tensor (the Hessian matrix of the gravitational potential). In linear theory, the problem can be solved exactly, and the result is attractively simple: the conditional mass function has no explicit dependence on the local tidal field, and is a function only of the local density on the filtering scale used to define the tidal tensor. There is nevertheless a strong implicit predicted dependence on geometrical environment, because the local density couples statistically to the derivatives of the potential. We compute the predictions of this model and study the limits of their validity by comparing them to results deduced empirically from $N$-body simulations. We have verified that, to a good approximation, the abundance of haloes in different environments depends only on their densities, and not on their tidal structure. In this sense we find relative differences between halo abundances in different environments with the same density which are smaller than 13%. Furthermore, for sufficiently large filtering scales, the agreement with the theoretical prediction is good, although there are important deviations from the Gaussian prediction at small, non-linear scales. We discuss how to obtain improved predictions in this regime, using the 'effective-universe' approach.

preprint2015arXiv

On the Bias of the Distance-Redshift Relation from Gravitational Lensing

A long standing question in cosmology is whether gravitational lensing changes the distance-redshift relation $D(z)$ or the mean flux density of sources. Interest in this has been rekindled by recent studies in non-linear relativistic perturbation theory that find biases in both the area of a surface of constant redshift and in the mean distance to this surface, with a fractional bias in both cases on the order of the mean squared convergence $\langle κ^2 \rangle$. Any such area bias could alter CMB cosmology, and the corresponding bias in mean flux density could affect supernova cosmology. Here we show that, in an ensemble averaged sense, the perturbation to the area of a surface of constant redshift is in reality much smaller, being on the order of the cumulative bending angle squared, or roughly a part-in-a-million effect. This validates the arguments of Weinberg (1976) that the mean magnification $μ$ of sources is unity and of Kibble \& Lieu (2005) that the mean direction-averaged inverse magnification is unity. It also validates the conventional treatment of lensing in analysis of CMB anisotropies. But the existence of a scatter in magnification will cause any non-linear function of these conserved quantities to be statistically biased. The distance $D$, for example, is proportional to $μ^{-1/2}$ so lensing will bias $\langle D\rangle$ even if $\langle μ\rangle=1$. The fractional bias in such quantities is generally of order $\langle κ^2 \rangle$, which is orders of magnitude larger than the area perturbation. Claims for large bias in area or flux density of sources appear to have resulted from misinterpretation of such effects: they do not represent a new non-Newtonian effect, nor do they invalidate standard cosmological analyses.

preprint2014arXiv

Galaxy and Mass Assembly (GAMA): The halo mass of galaxy groups from maximum-likelihood weak lensing

We present a maximum-likelihood weak lensing analysis of the mass distribution in optically selected spectroscopic Galaxy Groups (G$^3$Cv5) in the Galaxy And Mass Assembly (GAMA) survey, using background Sloan Digital Sky Survey (SDSS) photometric galaxies. The scaling of halo mass, $M_h$, with various group observables is investigated. Our main results are: 1) the measured relations of halo mass with group luminosity, virial volume and central galaxy stellar mass, $M_\star$, agree very well with predictions from mock group catalogues constructed from a GALFORM semi-analytical galaxy formation model implemented in the Millennium $Λ$CDM N-body simulation; 2) the measured relations of halo mass with velocity dispersion and projected half-abundance radius show weak tension with mock predictions, hinting at problems in the mock galaxy dynamics and their small scale distribution; 3) the median $M_h|M_\star$ measured from weak lensing depends more sensitively on the lognormal dispersion in $M_\star$ at fixed $M_h$ than it does on the median $M_\star|M_h$. Our measurements suggest an intrinsic dispersion of $σ_{\log(M_\star)}\sim 0.15$; 4) Comparing our mass estimates with those in the catalogue, we find that the G$^3$Cv5 mass can give biased results when used to select subsets of the group sample. Of the various new halo mass estimators that we calibrate using our weak lensing measurements, group luminosity is the best single-proxy estimator of group mass.

preprint2014arXiv

Mapping the Cosmic Web with the largest all-sky surveys

Our view of the low-redshift Cosmic Web has been revolutionized by galaxy redshift surveys such as 6dFGS, SDSS and 2MRS. However, the trade-off between depth and angular coverage limits a systematic three-dimensional account of the entire sky beyond the Local Volume (z<0.05). In order to reliably map the Universe to cosmologically significant depths over the full celestial sphere, one must draw on multiwavelength datasets and state-of-the-art photometric redshift techniques. We have undertaken a dedicated program of cross-matching the largest photometric all-sky surveys -- 2MASS, WISE and SuperCOSMOS -- to obtain accurate redshift estimates of millions of galaxies. The first outcome of these efforts -- the 2MASS Photometric Redshift catalog (2MPZ, Bilicki et al. 2014a) -- has been publicly released and includes almost 1 million galaxies with a mean redshift of z=0.08. Here we summarize how this catalog was constructed and how using the WISE mid-infrared sample together with SuperCOSMOS optical data allows us to push to redshift shells of z~0.2--0.3 on unprecedented angular scales. Our catalogs, with ~20 million sources in total, provide access to cosmological volumes crucial for studies of local galaxy flows (clustering dipole, bulk flow) and cross-correlations with the cosmic microwave background such as the integrated Sachs-Wolfe effect or lensing studies.

preprint2013arXiv

2MASS Photometric Redshift catalog: a comprehensive three-dimensional census of the whole sky

Key cosmological applications require the three-dimensional galaxy distribution on the entire celestial sphere. These include measuring the gravitational pull on the Local Group, estimating the large-scale bulk flow and testing the Copernican principle. However, the largest all-sky redshift surveys -- the 2MRS and IRAS PSCz -- have median redshifts of only z=0.03 and sample the very local Universe. There exist all-sky galaxy catalogs reaching much deeper -- SuperCOSMOS in the optical, 2MASS in the near-IR and WISE in the mid-IR -- but these lack complete redshift information. At present, the only rapid way towards larger 3D catalogs covering the whole sky is through photometric redshift techniques. In this paper we present the 2MASS Photometric Redshift catalog (2MPZ) containing 1 million galaxies, constructed by cross-matching 2MASS XSC, WISE and SuperCOSMOS all-sky samples and employing the artificial neural network approach (the ANNz algorithm), trained on such redshift surveys as SDSS, 6dFGS and 2dFGRS. The derived photometric redshifts have errors nearly independent of distance, with an all-sky accuracy of sigma_z=0.015 and a very small percentage of outliers. In this way, we obtain redshift estimates with a typical precision of 12% for all the 2MASS XSC galaxies that lack spectroscopy. In addition, we have made an early effort towards probing the entire 3D sky beyond 2MASS, by pairing up WISE with SuperCOSMOS and training the ANNz on GAMA redshift data reaching currently to z_med~0.2. This has yielded photo-z accuracies comparable to those in the 2MPZ. These all-sky photo-z catalogs, with a median z~0.1 for the 2MPZ, and significantly deeper for future WISE-based samples, will be the largest and most complete of their kind for the foreseeable future.

preprint2013arXiv

Galaxy And Mass Assembly (GAMA): improved cosmic growth measurements using multiple tracers of large-scale structure

We present the first application of a "multiple-tracer" redshift-space distortion (RSD) analysis to an observational galaxy sample, using data from the Galaxy and Mass Assembly survey (GAMA). Our dataset is an r < 19.8 magnitude-limited sample of 178,579 galaxies covering redshift interval z < 0.5 and area 180 deg^2. We obtain improvements of 10-20% in measurements of the gravitational growth rate compared to a single-tracer analysis, deriving from the correlated sample variance imprinted in the distributions of the overlapping galaxy populations. We present new expressions for the covariances between the auto-power and cross-power spectra of galaxy samples that are valid for a general survey selection function and weighting scheme. We find no evidence for a systematic dependence of the measured growth rate on the galaxy tracer used, justifying the RSD modelling assumptions, and validate our results using mock catalogues from N-body simulations. For multiple tracers selected by galaxy colour, we measure normalized growth rates in two independent redshift bins f*sigma_8(z=0.18) = 0.36 +/- 0.09 and f*sigma_8(z=0.38) = 0.44 +/- 0.06, in agreement with standard GR gravity and other galaxy surveys at similar redshifts.

preprint2013arXiv

Reconstructing the distribution of haloes and mock galaxies below the resolution limit in cosmological simulations

We present a method for populating dark matter simulations with haloes of mass below the resolution limit. It is based on stochastically sampling a field derived from the density field of the halo catalogue, using constraints from the conditional halo mass function n(m|δ). We test the accuracy of the method and show its application in the context of building mock galaxy samples. We find that this technique allows precise reproduction of the two-point statistics of galaxies in mock samples constructed with this method. Our results demonstrate that the full information content of a simulation can be communicated efficiently using only a catalogue of the more massive haloes.

preprint2013arXiv

Slipher, galaxies, and cosmological velocity fields

By 1917, V.M. Slipher had singlehandedly established a tendency for 'spiral nebulae' to be redshifted (21 out of 25 cases). From a modern perspective, it could seem surprising that the expansion of the universe was not announced at this point. Examination of Slipher's papers shows that he reached a more subtle conclusion: the identification of cosmological peculiar velocities, including the bulk motion of the Milky Way, leading to a beautiful argument in favour of nebulae as distant stellar systems. Nevertheless, Slipher's data actually contain evidence at >8sigma for a positive mean velocity, even after subtracting the dipole owing to the motion of the observer. In 1929, Hubble estimated distances for a sample of no greater depth, using redshifts due almost entirely to Slipher. Hubble's distances were flawed in two distinct ways: in addition to an incorrect absolute calibration, the largest distances were systematically under-estimated. Nevertheless, he claimed the detection of a linear distance-redshift relation. Statistically, the evidence for such a correlation is less strong than the simple evidence for a positive mean velocity in Hubble's sample. Comparison with modern data shows that a sample of more than twice Hubble's depth would generally be required in order to reveal clearly the global linear expansion in the face of the 'noise' from peculiar velocities. When the theoretical context of the time is examined, the role of the de Sitter model and its prediction of a linear distance-redshift relation looms large. A number of searches for this relation were performed prior to Hubble over the period 1924-1928, with a similar degree of success. All were based on the velocities measured by Slipher, whose work from a Century ago stands out both for the precision of his measurements and for the subtle clarity of the arguments he employed to draw correct conclusions from them.

preprint2012arXiv

Galaxy And Mass Assembly: Estimating galaxy group masses via caustic analysis

We have generated complementary halo mass estimates for all groups in the Galaxy And Mass Assembly Galaxy Group Catalogue (GAMA G3Cv1) using a modified caustic mass estimation algorithm, originally developed by Diaferio & Geller (1997). We calibrate the algorithm by applying it on a series of 9 GAMA mock galaxy light cones and investigate the effects of using different definitions for group centre and size. We select the set of parameters that provide median-unbiased mass estimates when tested on mocks, and generate mass estimates for the real group catalogue. We find that on average, the caustic mass estimates agree with dynamical mass estimates within a factor of 2 in 90.8 +/- 6.1% groups and compares equally well to velocity dispersion based mass estimates for both high and low multiplicity groups over the full range of masses probed by the G3Cv1.

preprint2010arXiv

Difficulties Distinguishing Dark Energy from Modified Gravity via Redshift Distortions

The bulk motion of galaxies induced by the growth of cosmic structure offers a rare opportunity to test the validity of general relativity across cosmological scales. However, modified gravity can be degenerate in its effect with the unknown values of cosmological parameters. More seriously, even the `observed' value of the RSD (redshift-space distortions) used to measure the fluctuation growth rate depends on the assumed cosmological parameters (the Alcock-Paczynski effect). We give a full analysis of these issues, showing how to combine RSD with BAO (baryon acoustic oscillations) and CMB (Cosmic Microwave Background) data, in order to obtain joint constraints on deviations from general relativity and on the equation of state of dark energy whilst allowing for factors such as non-zero curvature. In particular we note that the evolution of Omega_m(z), along with the Alcock-Paczynski effect, produces a degeneracy between the equation of state w and the modified growth parameter gamma. Typically, the total marginalized error on either of these parameters will be larger by a factor ~ 2 compared to the conditional error where one or other is held fixed. We argue that future missions should be judged by their Figure of Merit as defined in the w_p - gamma plane, and note that the inclusion of spatial curvature can degrade this value by an order of magnitude.

preprint2010arXiv

Galaxy and Mass Assembly (GAMA): FUV, NUV, ugrizYJHK Petrosian, Kron and Sèrsic photometry

In order to generate credible 0.1-2 μm SEDs, the GAMA project requires many Gigabytes of imaging data from a number of instruments to be re-processed into a standard format. In this paper we discuss the software infrastructure we use, and create self-consistent ugrizYJHK photometry for all sources within the GAMA sample. Using UKIDSS and SDSS archive data, we outline the pre-processing necessary to standardise all images to a common zeropoint, the steps taken to correct for seeing bias across the dataset, and the creation of Gigapixel-scale mosaics of the three 4x12 deg GAMA regions in each filter. From these mosaics, we extract source catalogues for the GAMA regions using elliptical Kron and Petrosian matched apertures. We also calculate Sérsic magnitudes for all galaxies within the GAMA sample using SIGMA, a galaxy component modelling wrapper for GALFIT 3. We compare the resultant photometry directly, and also calculate the r band galaxy LF for all photometric datasets to highlight the uncertainty introduced by the photometric method. We find that (1) Changing the object detection threshold has a minor effect on the best-fitting Schechter parameters of the overall population (M* +/- 0.055mag, α +/- 0.014, Φ* +/- 0.0005 h^3 Mpc^{-3}). (2) An offset between datasets that use Kron or Petrosian photometry regardless of the filter. (3) The decision to use circular or elliptical apertures causes an offset in M* of 0.20mag. (4) The best-fitting Schechter parameters from total-magnitude photometric systems (such as SDSS modelmag or Sérsic magnitudes) have a steeper faint-end slope than photometry dependent on Kron or Petrosian magnitudes. (5) Our Universe's total luminosity density, when calculated using Kron or Petrosian r-band photometry, is underestimated by at least 15%.

preprint2010arXiv

Unmodified Gravity

By relaxing the conventional assumption of a purely gravitational interaction between dark energy and dark matter, substantial alterations to the growth of cosmological structure can occur. In this work we focus on the homogeneous transfer of energy from a decaying form of dark energy. We present simple analytic solutions to the modified growth rates of matter fluctuations in these models, and demonstrate that neglecting physics within the dark sector may induce a significant bias in the inferred growth rate, potentially offering a false signature of modified gravity.

preprint2009arXiv

Galaxies-Intergalactic Medium Interaction Calculation --I. Galaxy formation as a function of large-scale environment

[Abridged] We present the first results of hydrodynamical simulations that follow the formation of galaxies to z=0 in spherical regions of radius ~20 Mpc/h drawn from the Millennium Simulation. The regions have overdensities that deviate by (-2, -1, 0, +1, +2)sigma from the cosmic mean, where sigma is the rms mass fluctuation on a scale of ~20Mpc/h at z=1.5. The simulations have mass resolution of up to 10^6 Msun/h, cover the entire range of large-scale environments and allow extrapolation of statistics to the entire 500 (Mpc/h)^3 Millennium volume. They include gas cooling, photoheating from an ionising background, SNe feedback and winds, but no AGN. We find that the specific SFR density at z <~ 10 varies systematically from region to region by up to an order of magnitude, but the global value, averaged over all volumes, reproduces observational data. Massive, compact galaxies, similar to those observed in the GOODS fields, form in the overdense regions as early as z=6, but do not appear in the underdense regions until z~3. These environmental variations are not caused by a dependence of the star formation properties on environment, but rather by a strong variation of the halo mass function from one environment to another, with more massive haloes forming preferentially in the denser regions. At all epochs, stars form most efficiently in haloes of circular velocity ~ 250 km/s. However, the star formation history exhibits a form of "downsizing" (even in the absence of AGN): the stars comprising massive galaxies at z=0 have mostly formed by z=1-2, whilst those comprising smaller galaxies typically form at later times. However, additional feedback is required to limit star formation in massive galaxies at late times.

preprint2009arXiv

GAMA: towards a physical understanding of galaxy formation

The Galaxy And Mass Assembly (GAMA) project is the latest in a tradition of large galaxy redshift surveys, and is now underway on the 3.9m Anglo-Australian Telescope at Siding Spring Observatory. GAMA is designed to map extragalactic structures on scales of 1kpc - 1Mpc in complete detail to a redshift of z~0.2, and to trace the distribution of luminous galaxies out to z~0.5. The principal science aim is to test the standard hierarchical structure formation paradigm of Cold Dark Matter (CDM) on scales of galaxy groups, pairs, discs, bulges and bars. We will measure (1) the Dark Matter Halo Mass Function (as inferred from galaxy group velocity dispersions); (2) baryonic processes, such as star formation and galaxy formation efficiency (as derived from Galaxy Stellar Mass Functions); and (3) the evolution of galaxy merger rates (via galaxy close pairs and galaxy asymmetries). Additionally, GAMA will form the central part of a new galaxy database, which aims to contain 275,000 galaxies with multi-wavelength coverage from coordinated observations with the latest international ground- and space-based facilities: GALEX, VST, VISTA, WISE, HERSCHEL, GMRT and ASKAP. Together, these data will provide increased depth (over 2 magnitudes), doubled spatial resolution (0.7"), and significantly extended wavelength coverage (UV through Far-IR to radio) over the main SDSS spectroscopic survey for five regions, each of around 50 deg^2. This database will permit detailed investigations of the structural, chemical, and dynamical properties of all galaxy types, across all environments, and over a 5 billion year timeline.

preprint2009arXiv

On lensing by a cosmological constant

Several recent papers have suggested that the cosmological constant Lambda directly influences the gravitational deflection of light. We place this problem in a cosmological context, deriving an expression for the linear potentials which control the cosmological bending of light, finding that it has no explicit dependence on the cosmological constant. To explore the physical origins of the apparent Lambda-dependent potential that appears in the static Kottler metric, we highlight the two classical effects which lead to the aberration of light. The first relates to the observer's motion relative to the source, and encapsulates the familiar concept of angular-diameter distance. The second term, which has proved to be the source of debate, arises from cosmic acceleration, but is rarely considered since it vanishes for photons with radial motion. This apparent form of light-bending gives the appearance of curved geodesics even within a flat and homogeneous universe. However this cannot be construed as a real lensing effect, since its value depends on the observer's frame of reference. Our conclusion is thus that standard results for gravitational lensing in a universe containing Lambda do not require modification, with any influence of Lambda being restricted to negligible high-order terms.

preprint2005arXiv

The 2dF Galaxy Redshift Survey: luminosity functions by density environment and galaxy type

We use the 2dF Galaxy Redshift Survey to measure the dependence of the bJ-band galaxy luminosity function on large-scale environment, defined by density contrast in spheres of radius 8h-1Mpc, and on spectral type, determined from principal component analysis. We find that the galaxy populations at both extremes of density differ significantly from that at the mean density. The population in voids is dominated by late types and shows, relative to the mean, a deficit of galaxies that becomes increasingly pronounced at magnitudes brighter than M_bJ-5log10h <-18.5. In contrast, cluster regions have a relative excess of very bright early-type galaxies with M_bJ-5log10h < -21. Differences in the mid to faint-end population between environments are significant: at M_bJ-5log10h=-18 early and late-type cluster galaxies show comparable abundances, whereas in voids the late types dominate by almost an order of magnitude. We find that the luminosity functions measured in all density environments, from voids to clusters, can be approximated by Schechter functions with parameters that vary smoothly with local density, but in a fashion which differs strikingly for early and late-type galaxies. These observed variations, combined with our finding that the faint-end slope of the overall luminosity function depends at most weakly on density environment, may prove to be a significant challenge for models of galaxy formation.

preprint2004arXiv

Galaxy ecology: groups and low-density environments in the SDSS and 2dFGRS

We analyse the observed correlation between galaxy environment and H-alpha emission line strength, using volume-limited samples and group catalogues of 24968 galaxies drawn from the 2dF Galaxy Redshift Survey (Mb<-19.5) and the Sloan Digital Sky Survey (Mr<-20.6). We characterise the environment by 1) Sigma_5, the surface number density of galaxies determined by the projected distance to the 5th nearest neighbour; and 2) rho1.1 and rho5.5, three-dimensional density estimates obtained by convolving the galaxy distribution with Gaussian kernels of dispersion 1.1 Mpc and 5.5 Mpc, respectively. We find that star-forming and quiescent galaxies form two distinct populations, as characterised by their H-alpha equivalent width, EW(Ha). The relative numbers of star-forming and quiescent galaxies varies strongly and continuously with local density. However, the distribution of EW(Ha) amongst the star-forming population is independent of environment. The fraction of star-forming galaxies shows strong sensitivity to the density on large scales, rho5.5, which is likely independent of the trend with local density, rho1.1. We use two differently-selected group catalogues to demonstrate that the correlation with galaxy density is approximately independent of group velocity dispersion, for sigma=200-1000 km/s. Even in the lowest density environments, no more than ~70 per cent of galaxies show significant H-alpha emission. Based on these results, we conclude that the present-day correlation between star formation rate and environment is a result of short-timescale mechanisms that take place preferentially at high redshift, such as starbursts induced by galaxy-galaxy interactions.

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: Clustering properties of radio galaxies

The clustering properties of local, S_{1.4 GHz} > 1 mJy, radio sources are investigated for a sample of 820 objects drawn from the joint use of the FIRST and 2dF Galaxy Redshift surveys. To this aim, we present 271 new bj < 19.45 spectroscopic counterparts of FIRST radio sources to be added to those already introduced in Magliocchetti et al. (2002). The two-point correlation function for the local radio population is found to be entirely consistent with estimates obtained for the whole sample of 2dFGRS galaxies. We estimate the parameters of the real-space correlation function xi(r)=(r/r_0)^{-γ}, r_0=6.7^{+0.9}_{-1.1} Mpc and γ=1.6\pm 0.1, where h=0.7 is assumed. Different results are instead obtained if we only consider sources that present signatures of AGN activity in their spectra. These objects are shown to be very strongly correlated, with r_0=10.9^{+1.0}_{-1.2} Mpc and γ=2\pm 0.1, a steeper slope than has been claimed in other recent works. No difference is found in the clustering properties of radio-AGNs of different radio luminosity. These results show that AGN-fuelled sources reside in dark matter halos more massive than \sim 10^{13.4} M_{\sun}},higher the corresponding figure for radio-quiet QSOs. This value can be converted into a minimum black hole mass associated with radio-loud, AGN-fuelled objects of M_{BH}^{min}\sim 10^9 M_{\sun}. The above results then suggest -at least for relatively faint radio objects -the existence of a threshold black hole mass associated with the onset of significant radio activity such as that of radio-loud AGNs; however, once the activity is triggered, there appears to be no evidence for a connection between black hole mass and level of radio output. (abridged)

preprint2004arXiv

The 2dF Galaxy Redshift Survey: Spherical Harmonics analysis of fluctuations in the final catalogue

We present the result of a decomposition of the 2dFGRS galaxy overdensity field into an orthonormal basis of spherical harmonics and spherical Bessel functions. Galaxies are expected to directly follow the bulk motion of the density field on large scales, so the absolute amplitude of the observed large-scale redshift-space distortions caused by this motion is expected to be independent of galaxy properties. By splitting the overdensity field into radial and angular components, we linearly model the observed distortion and obtain the cosmological constraint Omega_m^{0.6} sigma_8=0.46+/-0.06. The amplitude of the linear redshift-space distortions relative to the galaxy overdensity field is dependent on galaxy properties and, for L_* galaxies at redshift z=0, we measure beta(L_*,0)=0.58+/-0.08, and the amplitude of the overdensity fluctuations b(L_*,0) sigma_8=0.79+/-0.03, marginalising over the power spectrum shape parameters. Assuming a fixed power spectrum shape consistent with the full Fourier analysis produces very similar parameter constraints.

preprint2004arXiv

The 2dF Galaxy Redshift Survey: stochastic relative biasing between galaxy populations

It is well known that the clustering of galaxies depends on galaxy type.Such relative bias complicates the inference of cosmological parameters from galaxy redshift surveys, and is a challenge to theories of galaxy formation and evolution. In this paper we perform a joint counts-in-cells analysis on galaxies in the 2dF Galaxy Redshift Survey, classified by both colour and spectral type, eta, as early or late type galaxies. We fit three different models of relative bias to the joint probability distribution of the cell counts, assuming Poisson sampling of the galaxy density field. We investigate the nonlinearity and stochasticity of the relative bias, with cubical cells of side 10Mpc \leq L \leq 45Mpc (h=0.7). Exact linear bias is ruled out with high significance on all scales. Power law bias gives a better fit, but likelihood ratios prefer a bivariate lognormal distribution, with a non-zero `stochasticity' - i.e. scatter that may result from physical effects on galaxy formation other than those from the local density field. Using this model, we measure a correlation coefficient in log-density space (r_LN) of 0.958 for cells of length L=10Mpc, increasing to 0.970 by L=45Mpc. This corresponds to a stochasticity sigma_b/bhat of 0.44\pm0.02 and 0.27\pm0.05 respectively. For smaller cells, the Poisson sampled lognormal distribution presents an increasingly poor fit to the data, especially with regard to the fraction of completely empty cells. We compare these trends with the predictions of semianalytic galaxy formation models: these match the data well in terms of overall level of stochasticity, variation with scale, and fraction of empty cells.

preprint2002arXiv

Parameter constraints for flat cosmologies from CMB and 2dFGRS power spectra

We constrain flat cosmological models with a joint likelihood analysis of a new compilation of data from the cosmic microwave background (CMB) and from the 2dF Galaxy Redshift Survey (2dFGRS). Fitting the CMB alone yields a known degeneracy between the Hubble constant h and the matter density Omega_m, which arises mainly from preserving the location of the peaks in the angular power spectrum. This `horizon-angle degeneracy' is considered in some detail and shown to follow a simple relation Omega_m h^{3.4} = constant. Adding the 2dFGRS power spectrum constrains Omega_m h and breaks the degeneracy. If tensor anisotropies are assumed to be negligible, we obtain values for the Hubble constant h=0.665 +/- 0.047, the matter density Omega_m=0.313 +/- 0.055, and the physical CDM and baryon densities Omega_c h^2 = 0.115 +/- 0.009, Omega_b h^2 = 0.022 +/- 0.002 (standard rms errors). Including a possible tensor component causes very little change to these figures; we set a upper limit to the tensor-to-scalar ratio of r<0.7 at 95% confidence. We then show how these data can be used to constrain the equation of state of the vacuum, and find w<-0.52 at 95% confidence. The preferred cosmological model is thus very well specified, and we discuss the precision with which future CMB data can be predicted, given the model assumptions. The 2dFGRS power-spectrum data and covariance matrix, and the CMB data compilation used here, are available from http://www.roe.ac.uk/~wjp/

preprint2002arXiv

The 2dF Galaxy Redshift Survey: The amplitudes of fluctuations in the 2dFGRS and the CMB, and implications for galaxy biasing

We compare the amplitudes of fluctuations probed by the 2dF Galaxy Redshift Survey and by the latest measurements of the Cosmic Microwave Background anisotropies. By combining the 2dFGRS and CMB data we find the linear-theory rms mass fluctuations in 8 Mpc/h spheres to be sigma_8 = 0.73 +-0.05 (after marginalization over the matter density parameter Omega_m and three other free parameters). This normalization is lower than the COBE normalization and previous estimates from cluster abundance, but it is in agreement with some revised cluster abundance determinations. We also estimate the scale-independent bias parameter of present-epoch L_s = 1.9L_* APM-selected galaxies to be b(L_s,z=0) = 1.10 +- 0.08 on comoving scales of 0.02 < k < 0.15 h/Mpc. If luminosity segregation operates on these scales, L_* galaxies would be almost un-biased, b(L_*,z=0) = 0.96. These results are derived by assuming a flat Lambda-CDM Universe, and by marginalizing over other free parameters and fixing the spectral index n=1 and the optical depth due to reionization tau=0. We also study the best fit pair (Omega_m,b), and the robustness of the results to varying n and tau. Various modelling corrections can each change the resulting b by 5-15 per cent. The results are compared with other independent measurements from the 2dFGRS itself, and from the SDSS, cluster abundance and cosmic shear.

preprint2002arXiv

The 2dF Galaxy Redshift Survey: The Number and Luminosity Density of Galaxies

We present the bivariate brightness distribution (BBD) for the 2dF Galaxy Redshift Survey (2dFGRS) based on a preliminary subsample of 45,000 galaxies. The BBD is an extension of the galaxy luminosity function incorporating surface brightness information. It allows the measurement of the local luminosity density, j_B, and the galaxy luminosity and surface brightness distributions while accounting for surface brightness selection biases. The recovered 2dFGRS BBD shows a strong surface brightness-luminosity relation (M_B~2.4μ_e). The luminosity-density is dominated by normal galaxies and the luminosity-density peak lies away from the selection boundaries implying that the 2dFGRS is complete and that luminous low surface brightness galaxies are rare. The final value we derive for the local luminosity-density, inclusive of surface brightness corrections, is: j_B=2.49+/-0.20x10^8 h L_solar Mpc^-3. Representative Schechter function parameters are: M*=-19.75+/-0.05, phi*=2.02+/-0.02x10^-2 and alpha=-1.09+/-0.03. Extending the conventional methodology to incorporate surface brightness selection effects has resulted in an increase in the luminosity-density of 37%.

preprint2001arXiv

Evidence for a non-zero Lambda and a low matter density from a combined analysis of the 2dF Galaxy Redshift Survey and Cosmic Microwave Background Anisotropies

We perform a joint likelihood analysis of the power spectra of the 2dF Galaxy Redshift Survey (2dFGRS) and the cosmic microwave background (CMB) anisotropies under the assumptions that the initial fluctuations were adiabatic, Gaussian and well described by power laws with scalar and tensor indices of n_s and n_t. On its own, the 2dFGRS sets tight limits on the parameter combination Omega_m h, but relatively weak limits on the fraction of the cosmic matter density in baryons Omega_b/Omega_m. The CMB anisotropy data alone set poor constraints on the cosmological constant and Hubble constant because of a `geometrical degeneracy' among parameters. Furthermore, if tensor modes are allowed, the CMB data allow a wide range of values for the physical densities in baryons and cold dark matter. Combining the CMB and 2dFGRS data sets helps to break both the geometrical and tensor mode degeneracies. The values of the parameters derived here are consistent with the predictions of the simplest models of inflation, with the baryon density derived from primordial nucleosynthesis and with direct measurements of the Hubble parameter. In particular, we find strong evidence for a positive cosmological constant with a pm 2sigma range of 0.65 < Omega_Lambda < 0.85, completely independently of constraints on Omega_Λderived from Type Ia supernovae.

preprint2001arXiv

The 2dF Galaxy Redshift Survey: The power spectrum and the matter content of the universe

The 2dF Galaxy Redshift Survey has now measured in excess of 160000 galaxy redshifts. This paper presents the power spectrum of the galaxy distribution, calculated using a direct FFT-based technique. We argue that, within the k-space region 0.02<k<0.15 h Mpc^-1, the shape of this spectrum should be close to that of the linear density perturbations convolved with the window function of the survey. This window function and its convolving effect on the power spectrum estimate are analyzed in detail. By convolving model spectra, we are able to fit the power-spectrum data and provide a measure of the matter content of the universe. Our results show that models containing baryon oscillations are mildly preferred over featureless power spectra. Analysis of the data yields 68% confidence limits on the total matter density times the Hubble parameter Ω_m h = 0.20 +/- 0.03, and the baryon fraction Ω_b/Ω_m = 0.15 +/- 0.07, assuming scale-invariant primordial fluctuations.

preprint1993arXiv

Power Spectrum Analysis of Three-Dimensional Redshift Surveys

We develop a general method for power spectrum analysis of three dimensional redshift surveys. We present rigorous analytical estimates for the statistical uncertainty in the power and we are able to derive a rigorous optimal weighting scheme under the reasonable (and largely empirically verified) assumption that the long wavelength Fourier components are Gaussian distributed. We apply the formalism to the updated 1-in-6 QDOT IRAS redshift survey, and compare our results to data from other probes: APM angular correlations; the CfA and the Berkeley 1.2Jy IRAS redshift surveys. Our results bear out and further quantify the impression from e.g.\ counts-in-cells analysis that there is extra power on large scales as compared to the standard CDM model with $Ωh\simeq 0.5$. We apply likelihood analysis using the CDM spectrum with $Ωh$ as a free parameter as a phenomenological family of models; we find the best fitting parameters in redshift space and transform the results to real space. Finally, we calculate the distribution of the estimated long wavelength power. This agrees remarkably well with the exponential distribution expected for Gaussian fluctuations, even out to powers of ten times the mean. Our results thus reveal no trace of periodicity or other non-Gaussian behavior.