Source author record

Ofer Lahav

Ofer Lahav 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

89works
10topics
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

89 published item(s)

preprint2026arXiv

Kinematic lensing with high-resolution spectroscopic surveys. A unique opportunity for transformative cosmology at high redshifts in the 2040s

We present a science case to perform high-redshift cosmic shear surveys for cosmology with next-generation spectroscopic instruments, such as the proposed MegaMapper and Wide-field Spectroscopic Telescope. We argue that by using the novel technique called 'kinematic lensing' (KL) it will be possible to obtain shear catalogues at redshifts between 2 and 5. We show that the signal-to-noise ratio of KL at such high redshifts is on average twice as much that expected from current weak lensing (WL) surveys such as Euclid or LSST, and several times that of the previous generation of WL surveys like DES and KiDS, even with very conservative assumptions about the fraction of spectroscopically-detected sources for which KL shear estimates will be available. This will allow cosmologists to perform joint galaxy clustering-cosmic shear analyses over unprecedented cosmic volumes and to probe the growth of structures deep in the matter-dominated era and across the onset of dark-energy domination, offering a unique opportunity to unveil the mystery of cosmic acceleration.

preprint2026arXiv

Probing General Relativity on Cosmological Scales in the 2040s

General relativity is exquisitely tested in strong-field regimes, yet its validity on cosmological scales remains largely unexplored. Upcoming wide and deep large-scale structure surveys will access the ultra-large, linear scales where relativistic effects - Doppler terms, gravitational redshift, lensing magnification, and potential evolution - leave significant imprints in the clustering of galaxies. These signatures represent unique probes of spacetime that are inaccessible to standard Newtonian analyses but increasingly important as survey volumes grow. We outline the scientific potential of next-generation facilities, such as those envisioned within ESO's Expanding Horizons programme, to deliver the first robust measurements of relativistic effects in large-scale structure through multi-tracer power spectra and the single-tracer bispectrum of high-redshift Lyman-break galaxies. Detecting these contributions would open a new window on gravity, enabling precision tests of general relativity and its alternatives on cosmological scales in the 2040s.

preprint2025arXiv

Galaxy-Multiplet Clustering from DESI DR2

We present an efficient estimator for higher-order galaxy clustering using small groups of nearby galaxies, or multiplets. Using the Luminous Red Galaxy (LRG) sample from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2, we identify galaxy multiplets as discrete objects and measure their cross-correlations with the general galaxy field. Our results show that the multiplets exhibit stronger clustering bias as they trace more massive dark matter halos than individual galaxies. When comparing the observed clustering statistics with the mock catalogs generated from the N-body simulation AbacusSummit, we find that the mocks underpredict multiplet clustering despite reproducing the galaxy two-point auto-correlation reasonably well. This discrepancy indicates that the standard Halo Occupation Distribution (HOD) model is insufficient to describe the properties of galaxy multiplets, revealing the greater constraining power of this higher-order statistic on galaxy-halo connection and the possibility that multiplets are specific to additional assembly bias. We demonstrate that incorporating secondary biases into the HOD model improves agreement with the observed multiplet statistics, specifically by allowing galaxies to preferentially occupy halos in denser environments. Our results highlight the potential of utilizing multiplet clustering, beyond traditional two-point correlation measurements, to break degeneracies in models describing the galaxy-dark matter connection.

preprint2024arXiv

The Gravitational Lensing Imprints of DES Y3 Superstructures on the CMB: A Matched Filtering Approach

$ $Low density cosmic voids gravitationally lens the cosmic microwave background (CMB), leaving a negative imprint on the CMB convergence $κ$. This effect provides insight into the distribution of matter within voids, and can also be used to study the growth of structure. We measure this lensing imprint by cross-correlating the Planck CMB lensing convergence map with voids identified in the Dark Energy Survey Year 3 data set, covering approximately 4,200 deg$^2$ of the sky. We use two distinct void-finding algorithms: a 2D void-finder which operates on the projected galaxy density field in thin redshift shells, and a new code, Voxel, which operates on the full 3D map of galaxy positions. We employ an optimal matched filtering method for cross-correlation, using the MICE N-body simulation both to establish the template for the matched filter and to calibrate detection significances. Using the DES Y3 photometric luminous red galaxy sample, we measure $A_κ$, the amplitude of the observed lensing signal relative to the simulation template, obtaining $A_κ= 1.03 \pm 0.22$ ($4.6σ$ significance) for Voxel and $A_κ= 1.02 \pm 0.17$ ($5.9σ$ significance) for 2D voids, both consistent with $Λ$CDM expectations. We additionally invert the 2D void-finding process to identify superclusters in the projected density field, for which we measure $A_κ= 0.87 \pm 0.15$ ($5.9σ$ significance). The leading source of noise in our measurements is Planck noise, implying that future data from the Atacama Cosmology Telescope (ACT), South Pole Telescope (SPT) and CMB-S4 will increase sensitivity and allow for more precise measurements.

preprint2022arXiv

A Spectroscopic Road Map for Cosmic Frontier: DESI, DESI-II, Stage-5

In this white paper, we present an experimental road map for spectroscopic experiments beyond DESI. DESI will be a transformative cosmological survey in the 2020s, mapping 40 million galaxies and quasars and capturing a significant fraction of the available linear modes up to z=1.2. DESI-II will pilot observations of galaxies both at much higher densities and extending to higher redshifts. A Stage-5 experiment would build out those high-density and high-redshift observations, mapping hundreds of millions of stars and galaxies in three dimensions, to address the problems of inflation, dark energy, light relativistic species, and dark matter. These spectroscopic data will also complement the next generation of weak lensing, line intensity mapping and CMB experiments and allow them to reach their full potential.

preprint2022arXiv

Cosmology and neutrino mass with the Minimum Spanning Tree

The information content of the minimum spanning tree (MST), used to capture higher-order statistics and information from the cosmic web, is compared to that of the power spectrum for a $νΛ$CDM model. The measurements are made in redshift space using haloes from the Quijote simulation of mass $\geq 3.2\times 10^{13}\,h^{-1}{\rm M}_{\odot}$ in a box of length $L_{\rm box}=1\,h^{-1}{\rm Gpc}$. The power spectrum multipoles (monopole and quadrupole) are computed for Fourier modes in the range $0.006 < k < 0.5\, h{\rm Mpc}^{-1}$. For comparison the MST is measured with a minimum length scale of $l_{\min}\simeq13\,h^{-1}{\rm Mpc}$. Combining the MST and power spectrum allows for many of the individual degeneracies to be broken; on its own the MST provides tighter constraints on the sum of neutrino masses $M_ν$ and cosmological parameters $h$, $n_{\rm s}$, and $Ω_{\rm b}$ but the power spectrum alone provides tighter constraints on $Ω_{\rm m}$ and $σ_{8}$. Combined we find constraints that are a factor of two (or greater) on all parameters with respect to the power spectrum (for $M_ν$ there is a factor of four improvement). These improvements appear to be driven by the MST's sensitivity to small scale clustering, where the effect of neutrino free-streaming becomes relevant, and high-order statistical information in the cosmic web. The MST is shown to be a powerful tool for cosmology and neutrino mass studies, and therefore could play a pivotal role in ongoing and future galaxy redshift surveys (such as DES, DESI, \emph{Euclid}, and Rubin-LSST).

preprint2022arXiv

Explaining deep learning of galaxy morphology with saliency mapping

We successfully demonstrate the use of explainable artificial intelligence (XAI) techniques on astronomical datasets in the context of measuring galactic bar lengths. The method consists of training convolutional neural networks on human classified data from Galaxy Zoo in order to predict general galaxy morphologies, and then using SmoothGrad (a saliency mapping technique) to extract the bar for measurement by a bespoke algorithm. We contrast this to another method of using a convolutional neural network to directly predict galaxy bar lengths. These methods achieved correlation coefficients of 0.76 and 0.59, and root mean squared errors of 1.69 and 2.10 respective to human measurements. We conclude that XAI methods outperform conventional deep learning in this case, which could be reasonably explained by the larger datasets available when training the models. We suggest that our XAI method can be used to extract other galactic features (such as the bulge-to-disk ratio) without needing to collect new datasets or train new models. We also suggest that these techniques can be used to refine deep learning models as well as identify and eliminate bias within training datasets.

preprint2022arXiv

SOAR/Goodman Spectroscopic Assessment of Candidate Counterparts of the LIGO-Virgo Event GW190814

On 2019 August 14 at 21:10:39 UTC, the LIGO/Virgo Collaboration (LVC) detected a possible neutron star-black hole merger (NSBH), the first ever identified. An extensive search for an optical counterpart of this event, designated GW190814, was undertaken using the Dark Energy Camera (DECam) on the 4m Victor M. Blanco Telescope at the Cerro Tololo Inter-American Observatory. Target of Opportunity interrupts were issued on 8 separate nights to observe 11 candidates using the 4.1m Southern Astrophysical Research (SOAR) telescope's Goodman High Throughput Spectrograph in order to assess whether any of these transients was likely to be an optical counterpart of the possible NSBH merger. Here, we describe the process of observing with SOAR, the analysis of our spectra, our spectroscopic typing methodology, and our resultant conclusion that none of the candidates corresponded to the gravitational wave merger event but were all instead other transients. Finally, we describe the lessons learned from this effort. Application of these lessons will be critical for a successful community spectroscopic follow-up program for LVC observing run 4 (O4) and beyond.

preprint2022arXiv

The MegaMapper: A Stage-5 Spectroscopic Instrument Concept for the Study of Inflation and Dark Energy

In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2<z<5$. In order to achieve path-breaking results with a mid-scale investment, the MegaMapper combines existing technologies for critical path elements and pushes innovative development in other design areas. To this aim, we envision a 6.5-m Magellan-like telescope, with a newly designed wide field, coupled with DESI spectrographs, and small-pitch robots to achieve multiplexing of at least 26,000. This will match the expected achievable target density in the redshift range of interest and provide a 10x capability over the existing state-of the art, without a 10x increase in project budget.

preprint2022arXiv

Weak lensing magnification of Type Ia Supernovae from the Pantheon sample

Using data from the Pantheon SN Ia compilation and the Sloan Digital Sky Survey (SDSS), we propose an estimator for weak lensing convergence incorporating positional and photometric data of foreground galaxies. The correlation between this and the Hubble diagram residuals of the supernovae has $3.6σ$ significance, and is consistent with weak lensing magnification due to dark matter halos centered on galaxies. We additionally constrain the properties of the galactic haloes, such as the mass-to-light ratio $Γ$ and radial profile of the halo matter density $ρ(r)$. We derive a new relationship for the additional r.m.s. scatter in magnitudes caused by lensing, finding $σ_{\rm lens} = (0.06 \pm 0.017) (d_{\rm C}(z)/ d_{\rm C}(z=1))^{3/2}$ where $d_{\rm C}(z)$ is the comoving distance to redshift $z$. Hence the scatter in apparent magnitudes due lensing will be of the same size as the intrinsic scatter of SN Ia by $z \sim 1.2$. We propose a modification of the distance modulus estimator for SN Ia to incorporate lensing, which can be easily calculated from observational data. We anticipate this will improve the accuracy of cosmological parameter estimation for high-redshift SN Ia data.

preprint2021arXiv

Investigating Deep Learning Methods for Obtaining Photometric Redshift Estimations from Images

Knowing the redshift of galaxies is one of the first requirements of many cosmological experiments, and as it's impossible to perform spectroscopy for every galaxy being observed, photometric redshift (photo-z) estimations are still of particular interest. Here, we investigate different deep learning methods for obtaining photo-z estimates directly from images, comparing these with traditional machine learning algorithms which make use of magnitudes retrieved through photometry. As well as testing a convolutional neural network (CNN) and inception-module CNN, we introduce a novel mixed-input model which allows for both images and magnitude data to be used in the same model as a way of further improving the estimated redshifts. We also perform benchmarking as a way of demonstrating the performance and scalability of the different algorithms. The data used in the study comes entirely from the Sloan Digital Sky Survey (SDSS) from which 1 million galaxies were used, each having 5-filter (ugriz) images with complete photometry and a spectroscopic redshift which was taken as the ground truth. The mixed-input inception CNN achieved a mean squared error (MSE)=0.009, which was a significant improvement (30%) over the traditional Random Forest (RF), and the model performed even better at lower redshifts achieving a MSE=0.0007 (a 50% improvement over the RF) in the range of z<0.3. This method could be hugely beneficial to upcoming surveys such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) which will require vast numbers of photo-z estimates produced as quickly and accurately as possible.

preprint2020arXiv

Dark Energy Survey Year 1 Results: Cosmological Constraints from Cluster Abundances and Weak Lensing

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =σ_8(Ω_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $Ω_{\rm m}=0.179^{+0.031}_{-0.038}$, with $σ_8-Ω_{\rm m}$ posteriors in $2.4σ$ tension with the DES Y1 3x2pt results, and in $5.6σ$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($λ\ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

preprint2020arXiv

Deep learning dark matter map reconstructions from DES SV weak lensing data

We present the first reconstruction of dark matter maps from weak lensing observational data using deep learning. We train a convolution neural network (CNN) with a Unet based architecture on over $3.6\times10^5$ simulated data realizations with non-Gaussian shape noise and with cosmological parameters varying over a broad prior distribution. We interpret our newly created DES SV map as an approximation of the posterior mean $P(κ| γ)$ of the convergence given observed shear. Our DeepMass method is substantially more accurate than existing mass-mapping methods. With a validation set of 8000 simulated DES SV data realizations, compared to Wiener filtering with a fixed power spectrum, the DeepMass method improved the mean-square-error (MSE) by 11 per cent. With N-body simulated MICE mock data, we show that Wiener filtering with the optimal known power spectrum still gives a worse MSE than our generalized method with no input cosmological parameters; we show that the improvement is driven by the non-linear structures in the convergence. With higher galaxy density in future weak lensing data unveiling more non-linear scales, it is likely that deep learning will be a leading approach for mass mapping with Euclid and LSST.

preprint2020arXiv

GEOMAX: beyond linear compression for 3pt galaxy clustering statistics

We present the GEOMAX algorithm and its Python implementation for a two-step compression of bispectrum measurements. The first step groups bispectra by the geometric properties of their arguments; the second step then maximises the Fisher information with respect to a chosen set of model parameters in each group. The algorithm only requires the derivatives of the data vector with respect to the parameters and a small number of mock data, producing an effective, non-linear compression. By applying GEOMAX to bispectrum monopole measurements from BOSS DR12 CMASS redshift-space galaxy clustering data, we reduce the $68\%$ credible intervals for the inferred parameters $\left(b_1,b_2,f,σ_8\right)$ by $\left(50.4\%,56.1\%,33.2\%,38.3\%\right)$ with respect to standard MCMC on the full data vector. We run the analysis and comparison between compression methods over one hundred galaxy mocks to test the statistical significance of the improvements. On average GEOMAX performs $\sim15\%$ better than geometrical or maximal linear compression alone and is consistent with being lossless. Given its flexibility, the GEOMAX approach has the potential to optimally exploit three-point statistics of various cosmological probes like weak lensing or line-intensity maps from current and future cosmological data-sets such as DESI, Euclid, PFS and SKA.

preprint2020arXiv

The Impact of Peculiar Velocities on the Estimation of the Hubble Constant from Gravitational Wave Standard Sirens

In this work we investigate the systematic uncertainties that arise from the calculation of the peculiar velocity when estimating the Hubble constant ($H_0$) from gravitational wave standard sirens. We study the GW170817 event and the estimation of the peculiar velocity of its host galaxy, NGC 4993, when using Gaussian smoothing over nearby galaxies. NGC 4993 being a relatively nearby galaxy, at $\sim 40 \ {\rm Mpc}$ away, is subject to a significant effect of peculiar velocities. We demonstrate a direct dependence of the estimated peculiar velocity value on the choice of smoothing scale. We show that when not accounting for this systematic, a bias of $\sim 200 \ {\rm km \ s ^{-1}}$ in the peculiar velocity incurs a bias of $\sim 4 \ {\rm km \ s ^{-1} \ Mpc^{-1}}$ on the Hubble constant. We formulate a Bayesian model that accounts for the dependence of the peculiar velocity on the smoothing scale and by marginalising over this parameter we remove the need for a choice of smoothing scale. The proposed model yields $H_0 = 68.6 ^{+14.0}_{-8.5}~{\rm km\ s^{-1}\ Mpc^{-1}}$. We demonstrate that under this model a more robust unbiased estimate of the Hubble constant from nearby GW sources is obtained.

preprint2020arXiv

The sum of the masses of the Milky Way and M31: a likelihood-free inference approach

We use Density Estimation Likelihood-Free Inference, $Λ$ Cold Dark Matter simulations of $\sim 2M$ galaxy pairs, and data from Gaia and the Hubble Space Telescope to infer the sum of the masses of the Milky Way and Andromeda (M31) galaxies, the two main components of the Local Group. This method overcomes most of the approximations of the traditional timing argument, makes the writing of a theoretical likelihood unnecessary, and allows the non-linear modelling of observational errors that take into account correlations in the data and non-Gaussian distributions. We obtain an $M_{200}$ mass estimate $M_{\rm MW+M31} = 4.6^{+2.3}_{-1.8} \times 10^{12} M_{\odot}$ ($68 \%$ C.L.), in agreement with previous estimates both for the sum of the two masses and for the individual masses. This result is not only one of the most reliable estimates of the sum of the two masses to date, but is also an illustration of likelihood-free inference in a problem with only one parameter and only three data points.

preprint2019arXiv

Quantifying Suspiciousness Within Correlated Data Sets

We propose a principled Bayesian method for quantifying tension between correlated datasets with wide uninformative parameter priors. This is achieved by extending the Suspiciousness statistic, which is insensitive to priors. Our method uses global summary statistics, and as such it can be used as a diagnostic for internal consistency. We show how our approach can be combined with methods that use parameter space and data space to identify the existing internal discrepancies. As an example, we use it to test the internal consistency of the KiDS-450 data in 4 photometric redshift bins, and to recover controlled internal discrepancies in simulated KiDS data. We propose this as a diagnostic of internal consistency for present and future cosmological surveys, and as a tension metric for data sets that have non-negligible correlation, such as LSST and Euclid.

preprint2016arXiv

ANNz2 - photometric redshift and probability distribution function estimation using machine learning

We present ANNz2, a new implementation of the public software for photometric redshift (photo-z) estimation of Collister and Lahav (2004), which now includes generation of full probability distribution functions (PDFs). ANNz2 utilizes multiple machine learning methods, such as artificial neural networks and boosted decision/regression trees. The objective of the algorithm is to optimize the performance of the photo-z estimation, to properly derive the associated uncertainties, and to produce both single-value solutions and PDFs. In addition, estimators are made available, which mitigate possible problems of non-representative or incomplete spectroscopic training samples. ANNz2 has already been used as part of the first weak lensing analysis of the Dark Energy Survey, and is included in the experiment's first public data release. Here we illustrate the functionality of the code using data from the tenth data release of the Sloan Digital Sky Survey and the Baryon Oscillation Spectroscopic Survey. The code is available for download at https://github.com/IftachSadeh/ANNZ .

preprint2016arXiv

Cosmology from large scale galaxy clustering and galaxy-galaxy lensing with Dark Energy Survey Science Verification data

We present cosmological constraints from the Dark Energy Survey (DES) using a combined analysis of angular clustering of red galaxies and their cross-correlation with weak gravitational lensing of background galaxies. We use a 139 square degree contiguous patch of DES data from the Science Verification (SV) period of observations. Using large scale measurements, we constrain the matter density of the Universe as Omega_m = 0.31 +/- 0.09 and the clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/- 0.13 after marginalizing over seven nuisance parameters and three additional cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16} = 0.74 +/- 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 = 0.78 +/- 0.09 using two bins over the range 0.2 <z< 0.5. We study the robustness of the results under changes in the data vectors, modelling and systematics treatment, including photometric redshift and shear calibration uncertainties, and find consistency in the derived cosmological parameters. We show that our results are consistent with previous cosmological analyses from DES and other data sets and conclude with a joint analysis of DES angular clustering and galaxy-galaxy lensing with Planck CMB data, Baryon Accoustic Oscillations and Supernova type Ia measurements.

preprint2016arXiv

Could multiple voids explain the Cosmic Microwave Background Cold Spot anomaly?

Understanding the observed Cold Spot (CS) (temperature of ~ -150 mu K at its centre) on the Cosmic Microwave Background (CMB) is an outstanding problem. Explanations vary from assuming it is just a > 3 sigma primordial Gaussian fluctuation to the imprint of a supervoid via the Integrated Sachs-Wolfe and Rees-Sciama (ISW+RS) effects. Since single spherical supervoids cannot account for the full profile, the ISW+RS of multiple line-of-sight voids is studied here to mimic the structure of the cosmic web. Two structure configurations are considered. The first, through simulations of 20 voids, produces a central mean temperature of ~-50 mu K. In this model the central CS temperature lies at ~ 2 sigma but fails to explain the CS hot ring. An alternative multi-void model (using more pronounced compensated voids) produces much smaller temperature profiles, but contains a prominent hot ring. Arrangements containing closely placed voids at low redshift are found to be particularly well suited to produce CS-like profiles. We then measure the significance of the CS if CS-like profiles (which are fitted to the ISW+RS of multi-void scenarios) are removed. The CS tension with the LCDM model can be reduced dramatically for an array of temperature profiles smaller than the CS itself.

preprint2016arXiv

Halo detection via large-scale Bayesian inference

We present a proof-of-concept of a novel and fully Bayesian methodology designed to detect halos of different masses in cosmological observations subject to noise and systematic uncertainties. Our methodology combines the previously published Bayesian large-scale structure inference algorithm, HADES, and a Bayesian chain rule (the Blackwell-Rao Estimator), which we use to connect the inferred density field to the properties of dark matter halos. To demonstrate the capability of our approach we construct a realistic galaxy mock catalogue emulating the wide-area 6-degree Field Galaxy Survey, which has a median redshift of approximately 0.05. Application of HADES to the catalogue provides us with accurately inferred three-dimensional density fields and corresponding quantification of uncertainties inherent to any cosmological observation. We then use a cosmological simulation to relate the amplitude of the density field to the probability of detecting a halo with mass above a specified threshold. With this information we can sum over the HADES density field realisations to construct maps of detection probabilities and demonstrate the validity of this approach within our mock scenario. We find that the probability of successful of detection of halos in the mock catalogue increases as a function of the signal-to-noise of the local galaxy observations. Our proposed methodology can easily be extended to account for more complex scientific questions and is a promising novel tool to analyse the cosmic large-scale structure in observations.

preprint2016arXiv

Photometric Supernova Classification With Machine Learning

Automated photometric supernova classification has become an active area of research in recent years in light of current and upcoming imaging surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope, given that spectroscopic confirmation of type for all supernovae discovered will be impossible. Here, we develop a multi-faceted classification pipeline, combining existing and new approaches. Our pipeline consists of two stages: extracting descriptive features from the light curves and classification using a machine learning algorithm. Our feature extraction methods vary from model-dependent techniques, namely SALT2 fits, to more independent techniques fitting parametric models to curves, to a completely model-independent wavelet approach. We cover a range of representative machine learning algorithms, including naive Bayes, k-nearest neighbors, support vector machines, artificial neural networks and boosted decision trees (BDTs). We test the pipeline on simulated multi-band DES light curves from the Supernova Photometric Classification Challenge. Using the commonly used area under the curve (AUC) of the Receiver Operating Characteristic as a metric, we find that the SALT2 fits and the wavelet approach, with the BDTs algorithm, each achieves an AUC of 0.98, where 1 represents perfect classification. We find that a representative training set is essential for good classification, whatever the feature set or algorithm, with implications for spectroscopic follow-up. Importantly, we find that by using either the SALT2 or the wavelet feature sets with a BDT algorithm, accurate classification is possible purely from light curve data, without the need for any redshift information.

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.

preprint2015arXiv

CLASH-X: A Comparison of Lensing and X-ray Techniques for Measuring the Mass Profiles of Galaxy Clusters

We present profiles of temperature (Tx), gas mass, and hydrostatic mass estimated from new and archival X-ray observations of CLASH clusters. We compare measurements derived from XMM and Chandra observations with one another and compare both to gravitational lensing mass profiles derived with CLASH HST and ground-based lensing data. Radial profiles of Chandra and XMM electron density and enclosed gas mass are nearly identical, indicating that differences in hydrostatic masses inferred from X-ray observations arise from differences in Tx measurements. Encouragingly, cluster Txs are consistent with one another at ~100-200 kpc radii but XMM Tx systematically decline relative to Chandra Tx at larger radii. The angular dependence of the discrepancy suggests additional investigation on systematics such as the XMM point spread function correction, vignetting and off-axis responses. We present the CLASH-X mass-profile comparisons in the form of cosmology-independent and redshift-independent circular-velocity profiles. Ratios of Chandra HSE mass profiles to CLASH lensing profiles show no obvious radial dependence in the 0.3-0.8 Mpc range. However, the mean mass biases inferred from the WL and SaWLens data are different. e.g., the weighted-mean value at 0.5 Mpc is <b> = 0.12 for the WL comparison and <b> = -0.11 for the SaWLens comparison. The ratios of XMM HSE mass profiles to CLASH lensing profiles show a pronounced radial dependence in the 0.3-1.0 Mpc range, with a weighted mean mass bias of value rising to <b>~0.3 at ~1 Mpc for the WL comparison and <b> of 0.25 for SaWLens comparison. The enclosed gas mass profiles from both Chandra and XMM rise to a value 1/8 times the total-mass profiles inferred from lensing at 0.5 Mpc and remain constant outside of that radius, suggesting that [8xMgas] profiles may be an excellent proxy for total-mass profiles at >0.5 Mpc in massive galaxy clusters.

preprint2015arXiv

Discovery of two gravitationally lensed quasars in the Dark Energy Survey

We present spectroscopic confirmation of two new lensed quasars via data obtained at the 6.5m Magellan/Baade Telescope. The lens candidates have been selected from the Dark Energy Survey (DES) and WISE based on their multi-band photometry and extended morphology in DES images. Images of DES J0115-5244 show two blue point sources at either side of a red galaxy. Our long-slit data confirm that both point sources are images of the same quasar at $z_{s}=1.64.$ The Einstein Radius estimated from the DES images is $0.51$". DES J2200+0110 is in the area of overlap between DES and the Sloan Digital Sky Survey (SDSS). Two blue components are visible in the DES and SDSS images. The SDSS fiber spectrum shows a quasar component at $z_{s}=2.38$ and absorption compatible with Mg II and Fe II at $z_{l}=0.799$, which we tentatively associate with the foreground lens galaxy. The long-slit Magellan spectra show that the blue components are resolved images of the same quasar. The Einstein Radius is $0.68$" corresponding to an enclosed mass of $1.6\times10^{11}\,M_{\odot}.$ Three other candidates were observed and rejected, two being low-redshift pairs of starburst galaxies, and one being a quasar behind a blue star. These first confirmation results provide an important empirical validation of the data-mining and model-based selection that is being applied to the entire DES dataset.

preprint2015arXiv

Gravitational Redshift of Galaxies in Clusters from the Sloan Digital Sky Survey and the Baryon Oscillation Spectroscopic Survey

The gravitational redshift effect allows one to directly probe the gravitational potential in clusters of galaxies. Following up on Wojtak et al. [Nature (London) 477, 567 (2011)], we present a new measurement. We take advantage of new data from the tenth data release of the Sloan Digital Sky Survey and the Baryon Oscillation Spectroscopic Survey. We compare the spectroscopic redshift of the brightest cluster galaxies (BCGs) with that of galaxies at the outskirts of clusters, using a sample with an average cluster mass of $10^{14} M_{\odot}$. We find that these galaxies have an average relative redshift of -11 km/s compared with that of BCGs, with a standard deviation of +7 and -5 km/s. Our measurement is consistent with that of Wojtak et al. However, our derived standard deviation is larger, as we take into account various systematic effects, beyond the size of the dataset. The result is in good agreement with the predictions from general relativity.

preprint2014arXiv

CLASH: Extreme Emission Line Galaxies and Their Implication on Selection of High-Redshift Galaxies

We utilize the CLASH (Cluster Lensing And Supernova survey with Hubble) observations of 25 clusters to search for extreme emission-line galaxies (EELGs). The selections are carried out in two central bands: F105W (Y105) and F125W (J125), as the flux of the central bands could be enhanced by the presence of [O III] 4959, 5007 at redshift of about 0.93-1.14 and 1.57-1.79, respectively. The multi-band observations help to constrain the equivalent widths of emission lines. Thanks to cluster lensing, we are able to identify 52 candidates down to an intrinsic limiting magnitude of 28.5 and to a rest-frame [O III] 4959,5007 equivalent width of about 3737 angstrom. Our samples include a number of EELGs at lower luminosities that are missed in other surveys, and the extremely high equivalent width can be only found in such faint galaxies. These EELGs can mimic the dropout feature similar to that of high redshift galaxies and contaminate the color-color selection of high redshift galaxies when the S/N ratio is limited or the band coverage is incomplete. We predict that the fraction of EELGs in the future high redshift galaxy selections cannot be neglected.

preprint2014arXiv

Dark Energy, Paradigm Shifts, and the Role of Evidence

We comment on cases in the history of Astronomy, which may shed some light on the current established but enigmatic concordance model of Cosmology. Should the model be understood by adding new entities such as Dark Matter and Dark Energy, or by modifying the underlying theory? For example, the prediction and discovery of planet Neptune can be regarded as analogous to finding a dark component; while explaining the anomalous perihelion precession of Mercury by General Relativity can be taken as analogous to the possibility that modified gravity is an alternative to dark components of the universe. In this paper, we revise this analogy coming from the history of astronomy with an eye to illustrating some of the similarities and differences between the two cases.

preprint2014arXiv

Three Gravitationally Lensed Supernovae Behind CLASH Galaxy Clusters

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CLO12Car (z = 1.28), SN CLN12Did (z = 0.85), and SN CLA11Tib (z = 1.14), are located behind three different clusters, MACSJ1720.2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and Abell 383 (z = 0.187), respectively. Each SN was detected in Hubble Space Telescope (HST) optical and infrared images. Based on photometric classification, we find that SNe CLO12Car and CLN12Did are likely to be Type Ia supernovae (SNe Ia), while the classification of SN CLA11Tib is inconclusive. Using multi-color light-curve fits to determine a standardized SN Ia luminosity distance, we infer that SN CLO12Car was approximately 1.0 +/- 0.2 mag brighter than field SNe Ia at a similar redshift and ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is approximately 0.2 +/- 0.2 mag brighter than field SNe Ia. We derive independent estimates of the predicted magnification from CLASH strong+weak lensing maps of the clusters: 0.83 +/- 0.16 mag for SN CLO12Car, 0.28 +/- 0.08 mag for SN CLN12Did, and 0.43 +/- 0.11 mag for SN CLA11Tib. The two SNe Ia provide a new test of the cluster lens model predictions: we find that the magnifications based on the SN Ia brightness and those predicted by the lens maps are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe Ia (from the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct determination of absolute magnifications.

preprint2013arXiv

A Search for Correlation of Ultra-High Energy Cosmic Rays with IRAS-PSCz and 2MASS-6dF Galaxies

We study the arrival directions of 69 ultra-high energy cosmic rays (UHECRs) observed at the Pierre Auger Observatory (PAO) with energies exceeding 55 EeV. We investigate whether the UHECRs exhibit the anisotropy signal expected if the primary particles are protons that originate in galaxies in the local universe, or in sources correlated with these galaxies. We cross-correlate the UHECR arrival directions with the positions of IRAS-PSCz and 2MASS-6dF galaxies taking into account particle energy losses during propagation. This is the first time that the 6dF survey is used in a search for the sources of UHECRs and the first time that the PSCz survey is used with the full 69 PAO events. The observed cross-correlation signal is larger for the PAO UHECRs than for 94% (98%) of realisations from an isotropic distribution when cross-correlated with the PSCz (6dF). On the other hand the observed cross-correlation signal is lower than that expected from 85% of realisations, had the UHECRs originated in galaxies in either survey. The observed cross-correlation signal does exceed that expected by 50% of the realisations if the UHECRs are randomly deflected by intervening magnetic fields by 5 degrees or more. We propose a new method of analysing the expected anisotropy signal, by dividing the predicted UHECR source distribution into equal predicted flux radial shells, which can help localise and constrain the properties of UHECR sources. We find that the 69 PAO events are consistent with isotropy in the nearest of three shells we define, whereas there is weak evidence for correlation with the predicted source distribution in the two more distant shells in which the galaxy distribution is less anisotropic.

preprint2013arXiv

CLASH: Complete Lensing Analysis of the Largest Cosmic Lens MACS J0717.5+3745 and Surrounding Structures

The galaxy cluster MACS J0717.5+3745 (z=0.55) is the largest known cosmic lens, with complex internal structures seen in deep X-ray, Sunyaev-Zel'dovich effect and dynamical observations. We perform a combined weak and strong lensing analysis with wide-field BVRi'z' Subaru/Suprime-Cam observations and 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). We find consistent weak distortion and magnification measurements of background galaxies, and combine these signals to construct an optimally estimated radial mass profile of the cluster and its surrounding large-scale structure out to 5 Mpc/h. We find consistency between strong-lensing and weak-lensing in the region where these independent data overlap, <500 kpc/h. The two-dimensional weak-lensing map reveals a clear filamentary structure traced by distinct mass halos. We model the lensing shear field with 9 halos, including the main cluster, corresponding to mass peaks detected above 2.5σ_κ. The total mass of the cluster as determined by the different methods is M_{vir}=(2.8\pm0.4) \times 10^15 M_sun. Although this is the most massive cluster known at z>0.5, in terms of extreme value statistics we conclude that the mass of MACS J0717.5+3745 by itself is not in serious tension with LambdaCDM, representing only a ~2σ departure above the maximum simulated halo mass at this redshift.

preprint2013arXiv

Galaxy halo truncation and Giant Arc Surface Brightness Reconstruction in the Cluster MACSJ1206.2-0847

In this work we analyze the mass distribution of MACSJ1206.2-0847, especially focusing on the halo properties of its cluster members. The cluster appears relaxed in its X-ray emission, but has significant amounts of intracluster light which is not centrally concentrated, suggesting that galaxy-scale interactions are still ongoing despite the overall relaxed state. The cluster lenses 12 background galaxies into multiple images and one galaxy at $z=1.033$ into a giant arc and its counterimage. The multiple image positions and the surface brightness distribution (SFB) of the arc which is bent around several cluster members are sensitive to the cluster galaxy halo properties. We model the cluster mass distribution with a NFW profile and the galaxy halos with two parameters for the mass normalization and extent of a reference halo assuming scalings with their observed NIR--light. We match the multiple image positions at an r.m.s. level of $0.85\arcsec$ and can reconstruct the SFB distribution of the arc in several filters to a remarkable accuracy based on this cluster model. The length scale where the enclosed galaxy halo mass is best constrained is about 5 effective radii -- a scale in between those accessible to dynamical and field strong lensing mass estimates on one hand and galaxy--galaxy weak lensing results on the other hand. The velocity dispersion and halo size of a galaxy with $m_{\rm 160W,AB}=19.2$ or $M_{\rm B,Vega}=-20.7$ are $σ=150 \rm kms^{-1}$ and $r\approx 26\pm 6 \rm kpc$, indicating that the halos of the cluster galaxies are tidally stripped. We also reconstruct the unlensed source (which is smaller by a factor of $\sim5.8$ in area), demonstrating the increase of morphological information due to lensing and conclude that this galaxy has likely star--forming spiral arms with a red (older) central component.

preprint2013arXiv

PkANN - II. A non-linear matter power spectrum interpolator developed using artificial neural networks

In this paper we introduce PkANN, a freely available software package for interpolating the non-linear matter power spectrum, constructed using Artificial Neural Networks (ANNs). Previously, using Halofit to calculate matter power spectrum, we demonstrated that ANNs can make extremely quick and accurate predictions of the power spectrum. Now, using a suite of 6380 N-body simulations spanning 580 cosmologies, we train ANNs to predict the power spectrum over the cosmological parameter space spanning $3σ$ confidence level (CL) around the concordance cosmology. When presented with a set of cosmological parameters ($Ω_{\rm m} h^2, Ω_{\rm b} h^2, n_s, w, σ_8, \sum m_ν$ and redshift $z$), the trained ANN interpolates the power spectrum for $z\leq2$ at sub-per cent accuracy for modes up to $k\leq0.9\,h\textrm{Mpc}^{-1}$. PkANN is faster than computationally expensive N-body simulations, yet provides a worst-case error $<1$ per cent fit to the non-linear matter power spectrum deduced through N-body simulations. The overall precision of PkANN is set by the accuracy of our N-body simulations, at 5 per cent level for cosmological models with $\sum m_ν<0.5$ eV for all redshifts $z\leq2$. For models with $\sum m_ν>0.5$ eV, predictions are expected to be at 5 (10) per cent level for redshifts $z>1$ ($z\leq1$). The PkANN interpolator may be freely downloaded from http://zuserver2.star.ucl.ac.uk/~fba/PkANN

preprint2013arXiv

The Contribution of Halos with Different Mass Ratios to the Overall Growth of Cluster-Sized Halos

We provide a new observational test for a key prediction of the ΛCDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $σ$ agreement with ΛCDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated.

preprint2013arXiv

The DESI Experiment, a whitepaper for Snowmass 2013

The Dark Energy Spectroscopic Instrument (DESI) is a massively multiplexed fiber-fed spectrograph that will make the next major advance in dark energy in the timeframe 2018-2022. On the Mayall telescope, DESI will obtain spectra and redshifts for at least 18 million emission-line galaxies, 4 million luminous red galaxies and 3 million quasi-stellar objects, in order to: probe the effects of dark energy on the expansion history using baryon acoustic oscillations (BAO), measure the gravitational growth history through redshift-space distortions, measure the sum of neutrino masses, and investigate the signatures of primordial inflation. The resulting 3-D galaxy maps at z<2 and Lyman-alpha forest at z>2 will make 1%-level measurements of the distance scale in 35 redshift bins, thus providing unprecedented constraints on cosmological models.

preprint2013arXiv

Weighing the Local Group in the Presence of Dark Energy

We revise the mass estimate of the Local Group (LG) when Dark Energy (in the form of the Cosmological Constant) is incorporated into the Timing Argument (TA) mass estimator for the Local Group (LG). Assuming the age of the Universe and the Cosmological Constant according to the recent values from the Planck CMB experiment, we find the mass of the LG to be M_TAL = (4.73 +- 1.03) x 10^{12} M_sun, which is 13% higher than the classical TA mass estimate. This partly explains the discrepancy between earlier results from LCDM simulations and the classical TA. When a similar analysis is performed on 16 LG-like galaxy pairs from the CLUES simulations, we find that the scatter in the ratio of the virial to the TA estimated mass is given by M_vir/M_TAL = 1.04 +-0.16. Applying it to the LG mass estimation we find a calibrated M_vir = (4.92 +- 1.08 (obs) +- 0.79 (sys)) x 10^{12} M_sun.

preprint2012arXiv

A highly magnified candidate for a young galaxy seen when the Universe was 500 Myrs old

The early Universe at redshift z\sim6-11 marks the reionization of the intergalactic medium, following the formation of the first generation of stars. However, those young galaxies at a cosmic age of \lesssim 500 million years (Myr, at z \gtrsim 10) remain largely unexplored as they are at or beyond the sensitivity limits of current large telescopes. Gravitational lensing by galaxy clusters enables the detection of high-redshift galaxies that are fainter than what otherwise could be found in the deepest images of the sky. We report the discovery of an object found in the multi-band observations of the cluster MACS1149+22 that has a high probability of being a gravitationally magnified object from the early universe. The object is firmly detected (12 sigma) in the two reddest bands of HST/WFC3, and not detected below 1.2 μm, matching the characteristics of z\sim9 objects. We derive a robust photometric redshift of z = 9.6 \pm 0.2, corresponding to a cosmic age of 490 \pm 15Myr (i.e., 3.6% of the age of the Universe). The large number of bands used to derive the redshift estimate make it one of the most accurate estimates ever obtained for such a distant object. The significant magnification by cluster lensing (a factor of \sim15) allows us to analyze the object's ultra-violet and optical luminosity in its rest-frame, thus enabling us to constrain on its stellar mass, star-formation rate and age. If the galaxy is indeed at such a large redshift, then its age is less than 200 Myr (at the 95% confidence level), implying a formation redshift of zf \lesssim 14. The object is the first z>9 candidate that is bright enough for detailed spectroscopic studies with JWST, demonstrating the unique potential of galaxy cluster fields for finding highly magnified, intrinsically faint galaxies at the highest redshifts.

preprint2012arXiv

A Uniformly Derived Catalogue of Exoplanets from Radial Velocities

A new catalogue of extrasolar planets is presented by re-analysing a selection of published radial velocity data sets using EXOFIT (Balan & Lahav 2009). All objects are treated on an equal footing within a Bayesian framework, to give orbital parameters for 94 exoplanetary systems. Model selection (between one- and two-planet solutions) is then performed, using both a visual flagging method and a standard chi-square analysis, with agreement between the two methods for 99% of the systems considered. The catalogue is to be made available online, and this 'proof of concept' study may be maintained and extended in the future to incorporate all systems with publicly available radial velocity data, as well as transit and microlensing data.

preprint2012arXiv

CLASH: Mass Distribution in and around MACS J1206.2-0847 from a Full Cluster Lensing Analysis

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dlnΣ/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

preprint2012arXiv

CLASH: Precise New Constraints on the Mass Profile of Abell 2261

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with ΛCDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent ΛCDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

preprint2012arXiv

Cluster Lensing And Supernova survey with Hubble (CLASH): An Overview

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated.

preprint2012arXiv

On Combining Galaxy Clustering and Weak Lensing to Unveil Galaxy Biasing via the Halo Model

We formulate the concept of non-linear and stochastic galaxy biasing in the framework of halo occupation statistics. Using two-point statistics in projection, we define the galaxy bias function, b_g(r_p), and the galaxy-dark matter cross-correlation function, R_{gm}(r_p), where r_p is the projected distance. We use the analytical halo model to predict how the scale dependence of b_g and R_{gm}, over the range 0.1 Mpc/h < r_p < 30 Mpc/h, depends on the non-linearity and stochasticity in halo occupation models. In particular we quantify the effect due to the presence of central galaxies, the assumption for the radial distribution of satellite galaxies, the richness of the halo, and the Poisson character of the probability to have a certain number of satellite galaxies in a halo of a certain mass. Overall, brighter galaxies reveal a stronger scale dependence, and out to a larger radius. In real-space, we find that galaxy bias becomes scale independent, with R_{gm} = 1, for radii r > 1 - 5 Mpc/h, depending on luminosity. However, galaxy bias is scale-dependent out to much larger radii when one uses the projected quantities defined in this paper. These projected bias functions have the advantage that they are more easily accessible observationally and that their scale dependence carries a wealth of information regarding the properties of galaxy biasing. To observationally constrain the parameters of the halo occupation statistics and to unveil the origin of galaxy biasing we propose the use of the bias function Gamma_{gm}(r_p)=b_g(r_p)/R_{gm}(r_p). This function is obtained via a combination of weak gravitational lensing and galaxy clustering, and it can be measured using existing and forthcoming imaging and spectroscopic galaxy surveys.

preprint2012arXiv

PkANN - I. Non-linear matter power spectrum interpolation through artificial neural networks

We investigate the interpolation of power spectra of matter fluctuations using Artificial Neural Network (PkANN). We present a new approach to confront small-scale non-linearities in the power spectrum of matter fluctuations. This ever-present and pernicious uncertainty is often the Achilles' heel in cosmological studies and must be reduced if we are to see the advent of precision cosmology in the late-time Universe. We show that an optimally trained artificial neural network (ANN), when presented with a set of cosmological parameters (Omega_m h^2, Omega_b h^2, n_s, w_0, sigma_8, m_nu and redshift z), can provide a worst-case error <=1 per cent (for z<=2) fit to the non-linear matter power spectrum deduced through N-body simulations, for modes up to k<=0.7 h/Mpc. Our power spectrum interpolator is accurate over the entire parameter space. This is a significant improvement over some of the current matter power spectrum calculators. In this paper, we detail how an accurate interpolation of the matter power spectrum is achievable with only a sparsely sampled grid of cosmological parameters. Unlike large-scale N-body simulations which are computationally expensive and/or infeasible, a well-trained ANN can be an extremely quick and reliable tool in interpreting cosmological observations and parameter estimation. This paper is the first in a series. In this method paper, we generate the non-linear matter power spectra using HaloFit and use them as mock observations to train the ANN. This work sets the foundation for Paper II, where a suite of N-body simulations will be used to compute the non-linear matter power spectra at sub-per cent accuracy, in the quasi-non-linear regime 0.1 h/Mpc <= k <= 0.9 h/Mpc. A trained ANN based on this N-body suite will be released for the scientific community.

preprint2012arXiv

The impact of camera optical alignments on weak lensing measures for the Dark Energy Survey

Telescope Point Spread Function (PSF) quality is critical for realising the potential of cosmic weak lensing observations to constrain dark energy and test General Relativity. In this paper we use quantitative weak gravitational lensing measures to inform the precision of lens optical alignment, with specific reference to the Dark Energy Survey (DES). We compute optics spot diagrams and calculate the shear and flexion of the PSF as a function of position on the focal plane. For perfect optical alignment we verify the high quality of the DES optical design, finding a maximum PSF contribution to the weak lensing shear of 0.04 near the edge of the focal plane. However this can be increased by a factor of approximately three if the lenses are only just aligned within their maximum specified tolerances. We calculate the E and B-mode shear and flexion variance as a function of de-centre or tilt of each lens in turn. We find tilt accuracy to be a few times more important than de-centre, depending on the lens considered. Finally we consider the compound effect of de-centre and tilt of multiple lenses simultaneously, by sampling from a plausible range of values of each parameter. We find that the compound effect can be around twice as detrimental as when considering any one lens alone. Furthermore, this combined effect changes the conclusions about which lens is most important to align accurately. For DES, the tilt of the first two lenses is the most important.

preprint2011arXiv

The 2MASS Redshift Survey - Description and Data Release

We present the results of the 2MASS Redshift Survey (2MRS), a ten-year project to map the full three-dimensional distribution of galaxies in the nearby Universe. The 2 Micron All-Sky Survey (2MASS) was completed in 2003 and its final data products, including an extended source catalog (XSC), are available on-line. The 2MASS XSC contains nearly a million galaxies with Ks <= 13.5 mag and is essentially complete and mostly unaffected by interstellar extinction and stellar confusion down to a galactic latitude of |b|=5 deg for bright galaxies. Near-infrared wavelengths are sensitive to the old stellar populations that dominate galaxy masses, making 2MASS an excellent starting point to study the distribution of matter in the nearby Universe. We selected a sample of 44,599 2MASS galaxies with Ks <= 11.75 mag and |b|>=5 deg (>= 8 deg towards the Galactic bulge) as the input catalog for our survey. We obtained spectroscopic observations for 11,000 galaxies and used previously-obtained velocities for the remainder of the sample to generate a redshift catalog that is 97.6% complete to well-defined limits and covers 91% of the sky. This provides an unprecedented census of galaxy (baryonic mass) concentrations within 300 Mpc. Earlier versions of our survey have been used in a number of publications that have studied the bulk motion of the Local Group, mapped the density and peculiar velocity fields out to 50 Mpc, detected galaxy groups, and estimated the values of several cosmological parameters. Additionally, we present morphological types for a nearly-complete sub-sample of 20,860 galaxies with Ks < 11.25 mag and |b|>= 10 deg.

preprint2010arXiv

Excess Clustering on Large Scales in the MegaZ DR7 Photometric Redshift Survey

We observe a large excess of power in the statistical clustering of Luminous Red Galaxies in the photometric SDSS galaxy sample called MegaZ DR7. This is seen over the lowest multipoles in the angular power spectra C_{\ell} in four equally spaced redshift bins between 0.45 < z < 0.65. However, it is most prominent in the highest redshift band at ~ 4 sigma and it emerges at an effective scale k ~ 0.01 h Mpc^{-1}. Given that MegaZ DR7 is the largest cosmic volume galaxy survey to date (3.3 (Gpc h^{-1})^3) this implies an anomaly on the largest physical scales probed by galaxies. Alternatively, this signature could be a consequence of it appearing at the most systematically susceptible redshift. There are several explanations for this excess power that range from systematics to new physics. This could have important consequences for the next generation of galaxy surveys or the LCDM model. We test the survey, data and excess power, as well as possible origins.

preprint2010arXiv

Galaxy Zoo: Reproducing Galaxy Morphologies Via Machine Learning

We present morphological classifications obtained using machine learning for objects in SDSS DR6 that have been classified by Galaxy Zoo into three classes, namely early types, spirals and point sources/artifacts. An artificial neural network is trained on a subset of objects classified by the human eye and we test whether the machine learning algorithm can reproduce the human classifications for the rest of the sample. We find that the success of the neural network in matching the human classifications depends crucially on the set of input parameters chosen for the machine-learning algorithm. The colours and parameters associated with profile-fitting are reasonable in separating the objects into three classes. However, these results are considerably improved when adding adaptive shape parameters as well as concentration and texture. The adaptive moments, concentration and texture parameters alone cannot distinguish between early type galaxies and the point sources/artifacts. Using a set of twelve parameters, the neural network is able to reproduce the human classifications to better than 90% for all three morphological classes. We find that using a training set that is incomplete in magnitude does not degrade our results given our particular choice of the input parameters to the network. We conclude that it is promising to use machine- learning algorithms to perform morphological classification for the next generation of wide-field imaging surveys and that the Galaxy Zoo catalogue provides an invaluable training set for such purposes.

preprint2010arXiv

Neutrino masses from clustering of red and blue galaxies: a test of astrophysical uncertainties

Combining measurements of the galaxy power spectrum and the cosmic microwave background (CMB) is a powerful means of constraining the summed mass of neutrino species sum(m_nu), but is subject to systematic uncertainties due to non-linear structure formation, redshift-space distortions and galaxy bias. We empirically test the robustness of neutrino mass results to these effects by separately analyzing power spectra of red and blue galaxies from the Sloan Digital Sky Survey (SDSS-II) Data Release 7 (DR7), combined with the CMB five-year Wilkinson Microwave Anisotropy Probe (WMAP5) data. We consider fitting for a range of maximum wavenumber k using twelve different galaxy bias models. For example, using a new model based on perturbation theory and including redshift space distortions (Saito et al. 2009), the all-galaxy power spectrum combined with WMAP5 for a wavenumber range of k<0.2 Mpc/h yields 95% CL sum(m_nu)<0.46 eV. The red and blue galaxy power spectra give 0.41 and 0.63 eV respectively for this model. Using mock catalogues, we find the expected difference in these limits assuming a true neutrino mass of zero is 0.10 + or - 0.14 eV. Thus the difference of 0.22 eV between upper limits on neutrino mass for red and blue galaxies is approximately 1 sigma from the expected value. We find similar results for the other models and k ranges tested. This indicates good agreement for current data but hints at possible issues for next-generation surveys. Being able to perform such systematic tests is advantageous, and future surveys would benefit by including broad galaxy populations and luminosities that enable such a decomposition.

preprint2010arXiv

The Angular Power Spectra of Photometric SDSS LRGs

We construct new galaxy angular power spectra based on the extended, updated and final SDSS II Luminous Red Galaxy (LRG) photometric redshift survey: MegaZ DR7. Encapsulating 7746 deg^{2} we utilise 723,556 photometrically determined LRGs between 0.45 < z < 0.65 in a 3.3 (Gpc h^{-1})^3 spherical harmonic analysis of the galaxy distribution. By combining four photometric redshift bins we find preliminary parameter constraints of f_{b} = Ω_{b}/Ω_{m} = 0.173 +/- 0.046 and Ω_{m} = 0.260 +/- 0.035 assuming H_{0} = 75 km s^{-1} Mpc^{-1}, n_{s}=1 and Ω_{k} = 0. These limits are consistent with the CMB and the previous data release (DR4). The C_{\ell} are sensitive to redshift space distortions and therefore we also recast our constraints into a measurement of β~ Ω_{m}^{0.55}/b in different redshift shells. The robustness of these power spectra with respect to a number of potential systematics such as extinction, photometric redshift and ANNz training set extrapolation are examined. The latter includes a cosmological comparison of available photometric redshift estimation codes where we find excellent agreement between template and empirical estimation methods. MegaZ DR7 represents a methodological prototype to next generation surveys such as the Dark Energy Survey (DES) and, furthermore, is a photometric precursor to the spectroscopic BOSS survey. Our galaxy catalogue and all power spectra data can be found at http://zuserver2.star.ucl.ac.uk/~sat/MegaZ/MegaZDR7.tar.gz.

preprint2010arXiv

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

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

preprint2010arXiv

The Role of Environment on the Formation of Early-Type Galaxies

(Abridged) We present a detailed study of the stellar populations of a volume-limited sample of early-type galaxies from SDSS, across a range of environments -- defined as the mass of the host dark matter halo. The stellar populations are explored through the SDSS spectra, via projection onto a set of two spectral vectors determined from Principal Component Analysis. We find the velocity dispersion of the galaxy to be the main driver behind the different star formation histories of early-type galaxies. However, environmental effects are seen to play a role (although minor). Galaxies populating the lowest mass halos have stellar populations on average ~1Gyr younger than the rest of the sample. The fraction of galaxies with small amounts of recent star formation is also seen to be truncated when occupying halos more massive than 3E13Msun. The sample is split into satellite and central galaxies for a further analysis of environment. Satellites are younger than central galaxies of the same stellar mass. The younger satellite galaxies in 6E12Msun halos have stellar populations consistent with the central galaxies found in the lowest mass halos of our sample (i.e. 1E12Msun). This result is indicative of galaxies in lower mass halos being accreted into larger halos.

preprint2010arXiv

Upper Bound of 0.28 eV on the Neutrino Masses from the Largest Photometric Redshift Survey

We present a new upper limit of sum m_ν eV < 0.28 (95% CL) on the sum of the neutrino masses assuming a flat LCDM cosmology. This relaxes slightly to sum m_ν < 0.34 and sum m_ν < 0.47 when quasi non-linear scales are removed and w is not equal to -1, respectively. These bounds are derived from a new photometric redshift catalogue of over 700,000 Luminous Red Galaxies (MegaZ DR7) with a volume of 3.3 (Gpc h^-1)^3, extending over the redshift range 0.45 < z < 0.65 and up to angular scales of l_max = 300. The data are combined with WMAP 5-year CMB fluctuations, Baryon Acoustic Oscillations (BAO), type 1a Supernovae (SNe) and an HST prior on the Hubble parameter. This is the first combined constraint from a photometric redshift catalogue with other cosmological probes. When combined with WMAP this data set proves to be as constraining as the addition of all SNe and BAO data available to date. The upper limit is one of the tightest and `cleanest' constraints on the neutrino mass from cosmology or particle physics. Furthermore, if the aforementioned bounds hold, they all predict that current-to-next generation neutrino experiments, such as KATRIN, are unlikely to obtain a detection.

preprint2009arXiv

Exploring the Luminosity Evolution and Stellar Mass Assembly of 2SLAQ Luminous Red Galaxies Between Redshift 0.4 and 0.8

We present an analysis of the evolution of 8625 Luminous Red Galaxies between z = 0.4 and z = 0.8 in the 2dF and SDSS LRG and QSO (2SLAQ) survey. The LRGs are split into redshift bins and the evolution of both the luminosity and stellar mass function with redshift is considered and compared to the assumptions of a passive evolution scenario. We draw attention to several sources of systematic error that could bias the evolutionary predictions made in this paper. While the inferred evolution is found to be relatively unaffected by the exact choice of spectral evolution model used to compute K+e corrections, we conclude that photometric errors could be a source of significant bias in colour-selected samples such as this, in particular when using parametric maximum likelihood based estimators. We find that the evolution of the most massive LRGs is consistent with the assumptions of passive evolution and that the stellar mass assembly of the LRGs is largely complete by z ~ 0.8. Our findings suggest that massive galaxies with stellar masses above 10^11 solar masses must have undergone merging and star formation processes at a very early stage (z > 1). This supports the emerging picture of downsizing in both the star formation as well as the mass assembly of early type galaxies. Given that our spectroscopic sample covers an unprecedentedly large volume and probes the most massive end of the galaxy mass function, we find that these observational results present a significant challenge for many current models of galaxy formation

preprint2009arXiv

Forecasting neutrino masses from galaxy clustering in the Dark Energy Survey combined with the Planck Measurements

We study the prospects for detecting neutrino masses from the galaxy angular power spectrum in photometric redshift shells of the Dark Energy Survey (DES) over a volume of 20 (Gpc/h)^3 combined with the Cosmic Microwave Background (CMB) angular fluctuations expected to be measured from the Planck satellite. We find that for a Lambda-CDM concordance model with 7 free parameters in addition to a fiducial neutrino mass of M_nu = 0.24 eV, we recover from DES &Planck the correct value with uncertainty of +- 0.12 eV (95 % CL), assuming perfect knowledge of the galaxy biasing. If the fiducial total mass is close to zero, then the upper limit is 0.11 eV (95 % CL). This upper limit from DES &Planck is over 3 times tighter than using Planck alone, as DES breaks the parameter degeneracies in a CMB-only analysis. The analysis utlilizes spherical harmonics up to 300, averaged in bin of 10 to mimic the DES sky coverage. The results are similar if we supplement DES bands (grizY) with the VISTA Hemisphere Survey (VHS) near infrared band (JHK). The result is robust to uncertainties in non-linear fluctuations and redshift distortions. However, the result is sensitive to the assumed galaxy biasing schemes and it requires accurate prior knowledge of the biasing. To summarize, if the total neutrino mass in nature greater than 0.1eV, we should be able to detect it with DES &Planck, a result with great importance to fundamental Physics.

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

The 6dF Galaxy Survey: Final Redshift Release (DR3) and Southern Large-Scale Structures

We report the final redshift release of the 6dF Galaxy Survey, a combined redshift and peculiar velocity survey over the southern sky (|b|>10 deg). Its 136,304 spectra have yielded 110,256 new extragalactic redshifts and a new catalogue of 125,071 galaxies making near-complete samples with (K, H, J, r_F, b_J) <= (12.65, 12.95, 13.75, 15.60, 16.75). The median redshift of the survey is 0.053. Survey data, including images, spectra, photometry and redshifts, are available through an online database. We describe changes to the information in the database since earlier interim data releases. Future releases will include velocity dispersions, distances and peculiar velocities for the brightest early-type galaxies, comprising about 10% of the sample. Here we provide redshift maps of the southern local universe with z<=0.1, showing nearby large-scale structures in hitherto unseen detail. A number of regions known previously to have a paucity of galaxies are confirmed as significantly underdense regions. The URL of the 6dFGS database is http://www-wfau.roe.ac.uk/6dFGS

preprint2008arXiv

A Comparison of Six Photometric Redshift Methods Applied to 1.5 Million Luminous Red Galaxies

We present an updated version of MegaZ-LRG (Collister et al.,(2007)) with photometric redshifts derived with the neural network method, ANNz as well as five other publicly available photo-z codes (HyperZ, SDSS, Le PHARE, BPZ and ZEBRA) for ~1.5 million Luminous Red Galaxies (LRGs) in SDSS DR6. This allows us to identify how reliable codes are relative to each other if used as described in their public release. We compare and contrast the relative merits of each code using ~13000 spectroscopic redshifts from the 2SLAQ sample. We find that the performance of each code depends on the figure of merit used to assess it. As expected, the availability of a complete training set means that the training method performs best in the intermediate redshift bins where there are plenty of training objects. Codes such as Le PHARE, which use new observed templates perform best in the lower redshift bins. All codes produce reasonable photometric redshifts, the 1-sigma scatters ranging from 0.057 to 0.097 if averaged over the entire redshift range. We also perform tests to check whether a training set from a small region of the sky such as 2SLAQ produces biases if used to train over a larger area of the sky. We conclude that this is not likely to be a problem for future wide-field surveys. The complete photometric redshift catalogue including redshift estimates and errors on these from all six methods can be found at www.star.ucl.ac.uk/~mbanerji/MegaZLRGDR6/megaz.html

preprint2008arXiv

Halo-model signatures from 380,000 SDSS Luminous Red Galaxies with photometric redshifts

We analyze the small-scale clustering in "MegaZ-LRG", a large photometric-redshift catalogue of Luminous Red Galaxies extracted from the imaging dataset of the Sloan Digital Sky Survey. MegaZ-LRG, presented in a companion paper, spans the redshift range 0.4 < z < 0.7 with an r.m.s. redshift error dz ~ 0.03(1+z), covering 5,914 deg^2 to map out a total cosmic volume 2.5 h^-3 Gpc^3. In this study we use 380,000 photometric redshifts to measure significant deviations from the canonical power-law fit to the angular correlation function in a series of narrow redshift slices, in which we construct volume-limited samples. These deviations are direct signatures of the manner in which these galaxies populate the underlying network of dark matter haloes. We cleanly delineate the separate contributions of the "1-halo" and "2-halo" clustering terms and fit our measurements by parameterizing the halo occupation distribution N(M) of the galaxies. Our results are successfully fit by a "central" galaxy contribution with a "soft" transition from zero to one galaxies, combined with a power-law "satellite" galaxy component, the slope of which is a strong function of galaxy luminosity. The large majority of galaxies are classified as central objects of their host dark matter haloes rather than satellites in more massive systems. The effective halo mass of MegaZ-LRG galaxies lies in the range log_10 (M_eff / h^-1 M_sol) = 13.61 - 13.8 (increasing with redshift, assuming large-scale normalization sigma_8 = 0.8) for corresponding number densities in the range n_g = 5.03 - 0.56 x 10^-4 h^3 Mpc^-3. Our results confirm the usefulness of the halo model for gaining physical insight into the patterns of galaxy clustering.

preprint2008arXiv

Photometric Redshifts for the Dark Energy Survey and VISTA and Implications for Large Scale Structure

We conduct a detailed analysis of the photometric redshift requirements for the proposed Dark Energy Survey (DES) using two sets of mock galaxy simulations and an artificial neural network code - ANNz. In particular, we examine how optical photometry in the DES grizY bands can be complemented with near infra-red photometry from the planned VISTA Hemisphere Survey (VHS) in the JHK_s bands. We find that the rms scatter on the photometric redshift estimate over 1<z<2 is sigma_z=0.2 from DES alone and sigma_z=0.15 from DES+VISTA, i.e. an improvement of more than 30%. We draw attention to the effects of galaxy formation scenarios such as reddening on the photo-z estimate and using our neural network code, calculate the extinction, A_v for these reddened galaxies. We also look at the impact of using different training sets when calculating photometric redshifts. In particular, we find that using the ongoing DEEP2 and VVDS-Deep spectroscopic surveys to calibrate photometric redshifts for DES, will prove effective. However we need to be aware of uncertainties in the photometric redshift bias that arise when using different training sets as these will translate into errors in the dark energy equation of state parameter, w. Furthermore, we show that the neural network error estimate on the photometric redshift may be used to remove outliers from our samples before any kind of cosmological analysis, in particular for large-scale structure experiments. By removing all galaxies with a neural network photo-z error estimate of greater than 0.1 from our DES+VHS sample, we can constrain the galaxy power spectrum out to a redshift of 2 and reduce the fractional error on this power spectrum by ~15-20% compared to using the entire catalogue.

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

Measuring our universe from galaxy redshift surveys

Galaxy redshift surveys have achieved significant progress over the last couple of decades. Those surveys tell us in the most straightforward way what our local universe looks like. While the galaxy distribution traces the bright side of the universe, detailed quantitative analyses of the data have even revealed the dark side of the universe dominated by non-baryonic dark matter as well as more mysterious dark energy (or Einstein's cosmological constant). We describe several methodologies of using galaxy redshift surveys as cosmological probes, and then summarize the recent results from the existing surveys. Finally we present our views on the future of redshift surveys in the era of Precision Cosmology.

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.

preprint2004arXiv

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

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

preprint2003arXiv

Upper limits on neutrino masses from the 2dFGRS and WMAP: the role of priors

Solar, atmospheric, and reactor neutrino experiments have confirmed neutrino oscillations, implying that neutrinos have non-zero mass, but without pinning down their absolute masses. While it is established that the effect of neutrinos on the evolution of cosmic structure is small, the upper limits derived from large-scale structure data could help significantly to constrain the absolute scale of the neutrino masses. In a recent paper the 2dF Galaxy Redshift Survey (2dFGRS) team provided an upper limit m_nu,tot < 2.2 eV, i.e. approximately 0.7 eV for each of the three neutrino flavours, or phrased in terms of their contributioin to the matter density, Omega_nu/Omega_m < 0.16. Here we discuss this analysis in greater detail, considering issues of assumed 'priors' like the matter density Omega_m and the bias of the galaxy distribution with respect the dark matter distribution. As the suppression of the power spectrum depends on the ratio Omega_nu/Omega_m, we find that the out-of- fashion Mixed Dark Matter Model, with Omega_nu=0.2, Omega_m=1 and no cosmological constant, fits the 2dFGRS power spectrum and the CMB data reasonably well, but only for a Hubble constant H_0<50 km/s/Mpc. As a consequence, excluding low values of the Hubble constant, e.g. with the HST Key Project is important in order to get a strong constraint on the neutrino masses. We also comment on the improved limit by the WMAP team, and point out that the main neutrino signature comes from the 2dFGRS and the Lyman alpha forest.

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: A targeted study of catalogued clusters of galaxies

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

preprint2002arXiv

The 2dF Galaxy Redshift Survey: The 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 environmental dependence of galaxy star formation rates near clusters

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

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: Near Infrared Galaxy Luminosity Functions

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

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.

preprint2000arXiv

Constraints on Cosmological Anisotropy out to z=1 from Supernovae Ia

A combined sample of 79 high and low redshift supernovae Ia (SNe) is used to set constraints on the degree of anisotropy in the Universe out to $z\simeq1$. First we derive the global most probable values of matter density $Ω_M $, the cosmological constant $Ω_Λ$, and the Hubble constant $H_0$, and find them to be consistent with the published results from the two data sets of Riess et al. 1998 (R98) and Perlmutter et al. 1999 (P99). We then examine the Hubble diagram (HD, i.e., the luminosity-redshift relation) in different directions on the sky by utilising spherical harmonic expansion. In particular, via the analysis of the dipole anisotropy, we divide the sky into the two hemispheres that yield the most discrepant of the three cosmological parameters, and the scatter $χ^2_{\rm HD}$ in each case. The most discrepant values roughly move along the locus $-4Ω_M +3 Ω_Λ = 1$ (cf. P99), but by no more than $Δ\approx 2.5$ along this line. For a perfect FRW universe, Monte Carlo realizations that mimic the current set of SNe yield values higher than the measured $Δ$ in $\sim 1/5$ of the cases. We discuss implications for the validity of the Cosmological Principle, and possible calibration problems in the SNe data sets.

preprint1999arXiv

Testing Homogeneity on Large Scales

We review observational tests for the homogeneity of the Universe on large scales. Redshift and peculiar velocity surveys, radio sources, the X-Ray Background, the Lyman-alpha forest and the Cosmic Microwave Background are used to set constraints on inhomogeneous models and in particular on fractal-like models. Assuming the Cosmological Principle and the FRW metric, we estimate cosmological parameters by joint analysis of peculiar velocities, the CMB, cluster abundance, IRAS and Supernovae. Under certain assumptions the best fit density parameter is Omega_m = 1 - lambda = 0.3-0.5.

preprint1996arXiv

The two-point correlation function and morphological segregation in the Optical Redshift Survey

We study the clustering of galaxies in real and redshift space using the Optical Redshift Survey (ORS). We estimate the two point correlation function in redshift space, $ξ(s)$, for several subsamples of ORS, spanning nearly a factor of 30 in volume and detect significant variations in $ξ(s)$ among the subsamples covering small volumes. For volumes \gtsima $(75 h^{-1} {\rm Mpc})^{3}$ the ORS subsamples present very similar clustering patterns. Powerlaw fits to $ξ(s)$ give best-fit values in the range $1.5 \leq γ_{s} \leq 1.7 $ and $6.5 \leq s_{0} \leq 8.8 h^{-1}$ Mpc for several samples extending to redshifts of 8000 km s$^{-1}$. We find that $ξ(s)$ is larger for the magnitude-limited sample than for diameter-limited one within a radius of 4000 km s$^{-1}$. We interpret this as an indirect result of morphological segregation coupled with differences in morphological mix. We split ORS into morphological subsamples and confirm the existence of morphological segregation of galaxies out to scales of $s \sim 10 h^{-1}$ Mpc. Our results indicate that the relative bias factor between early type galaxies and late-types may be weakly dependent on scale. If real, this would suggest non-linear biasing. We also compute correlations as a function of radial and projected separations, $ξ(r_p, π)$ and derive the real space correlation function, $ξ(r)$. The results obtained in real space confirm those found using $ξ(s)$.

preprint1996arXiv

The X-Ray Background as a Probe of Density Fluctuations at High Redshift

The X-Ray Background (XRB) probes structure on scales intermediate between those explored by local galaxy redshift surveys and by the COBE Microwave Background measurements. We predict the large scale angular fluctuations in the XRB, expressed in terms of spherical harmonics for a range of assumed power-spectra and evolution scenarios. The dipole is due to large scale structure as well as to the observer's motion (the Compton-Getting effect). For a typical observer the two effects turn out to be comparable in amplitude. The coupling of the two effects makes it difficult to use the XRB for independent confirmation of the CMB dipole being due to the observer's motion. The large scale structure dipole (rms per component) relative to the monopole is in the range $a_{1m}/a_{00} \sim (0.5-9.0) \times 10^{-3} $. The spread is mainly due to the assumed redshift evolution scenarios of the X-ray volume emissivity $ρ_x(z)$. The dipole's prediction is consistent with a measured dipole in the HEAO1 XRB map. Typically, the harmonic spectrum drops with $l$ like $a_{lm} \sim l^{-0.4}$. This behaviour allows us to discriminate a true clustering signal against the flux shot noise, which is constant with $l$, and may dominate the signal unless bright resolved sources are removed from the XRB map. We also show that Sachs-Wolfe and Doppler (due to the motion of the sources) effects in the XRB are negligible. Although our analysis focuses on the XRB, the formalism is general and can be easily applied to other cosmological backgrounds.

preprint1996arXiv

Wiener Reconstruction of the IRAS 1.2Jy Galaxy Redshift Survey: Cosmographical Implications

We utilise the method of Wiener reconstruction with spherical harmonics and Bessel functions to recover the density and velocity fields from the IRAS 1.2Jy redshift survey. The results are robust to changes in prior parameters and the number of expansion coefficients. Maps of these fields are presented in a variety of projections. Many known structures are observed, including clear confirmations of the clusters N1600 and A3627. The Perseus-Pisces supercluster appears to extend out to roughly 9000km/s, and the reconstruction shows `backside infall' to the Centaurus/Great Attractor region. A qualitative comparison of the reconstructed IRAS gravity field with that from Tully-Fisher peculiar velocity measurements (Mark III) shows reasonable agreement. The Wiener reconstruction of the density field is also the optimal reconstruction (in the minimum variance sense) of any quantity which is linear in the density contrast. We show reconstructions of three such quantities. The misalignment angle between the IRAS and CMB Local Group dipoles is only 13 degrees out to 5000km/s, but increases to 25 degrees out to 20,000km/s. The reconstructed IRAS bulk flow out to 5000km/s is roughly 300km/s, which agrees in amplitude with that derived from the Mark III peculiar velocities (370km/s). However, the two bulk flow vectors deviate by 70 degrees. Finally, moment of inertia analysis shows that the Wiener reconstructed Supergalactic Plane is aligned within 30 degrees of that defined by de Vaucouleurs.

preprint1994arXiv

Large-Scale Mass Distribution Behind the Galactic Plane

We map the smoothed mass-density distribution in the Galactic zone of avoidance (zoa), within $6000\kms$ of the Local Group, using potent reconstruction from peculiar velocities of galaxies. The interpolation into the zoa is based on the assumed gravitational nature of the velocity field implying a potential flow. The main dynamical features at a distance r=4000 km/s are (a) the peak of the Great Attractor connecting Centaurus and Pavo at l=330, (b) a moderate bridge connecting Perseus-Pisces and Cepheus at l= 140, and (c) an extension of a large void from the southern Galactic hemisphere into the zoa near the direction of Puppis, l=220-270. We find a strong correlation between the mass density and the IRAS and optical galaxy density at $b=+-20, which indicates that the main dynamical features in the zoa should also be seen in galaxy surveys through the Galactic plane. The gravitational acceleration at the Local Group, based on the mass distribution out to 6000 km/s, is strongly affected by the mass distribution in the zoa: its direction changes by 31 when the |b| < 20 zoa is included, bringing it to within 4 +- 19 of the CMB dipole.

preprint1994arXiv

Wiener Reconstruction of Density, Velocity, and Potential Fields from All-Sky Galaxy Redshift Surveys

We present a new method for recovering the cosmological density, velocity, and potential fields from all-sky redshift catalogues. The method is based on an expansion of the fields in orthogonal radial (Bessel) and angular (spherical harmonic) functions. In this coordinate system, peculiar velocities introduce a coupling of the radial harmonics describing the density field in redshift space but leave the angular modes unaffected. In the harmonic transform space, this radial coupling is described by a distortion matrix which can be computed analytically within the context of linear theory; the redshift space harmonics can then be converted to their real space values by inversion of this matrix. Statistical noise is mitigated by regularizing the matrix inversion with a Wiener filter. The method yields a minimum variance estimate of the density field in real space. In this coordinate system, the minimum variance harmonics of the peculiar velocity and potential fields are related to those of the density field by simple linear transformations. Tests of the method with simulations of a CDM universe and comparison with previously proposed methods demonstrate it to be a very promising new reconstruction method for the local density and velocity field. A first application to the 1.2 Jy IRAS redshift survey is presented.