• The Pan-STARRS1 Surveys(1612.05560)

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($\Delta ra, \Delta dec$) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.
• This work is on the Physics of the B Factories. Part A of this book contains a brief description of the SLAC and KEK B Factories as well as their detectors, BaBar and Belle, and data taking related issues. Part B discusses tools and methods used by the experiments in order to obtain results. The results themselves can be found in Part C. Please note that version 3 on the archive is the auxiliary version of the Physics of the B Factories book. This uses the notation alpha, beta, gamma for the angles of the Unitarity Triangle. The nominal version uses the notation phi_1, phi_2 and phi_3. Please cite this work as Eur. Phys. J. C74 (2014) 3026.
• Optical Confirmation and Redshift Estimation of the Planck Cluster Candidates overlapping the Pan-STARRS Survey(1407.6001)

May 29, 2015 astro-ph.CO
We report results of a study of Planck Sunyaev-Zel'dovich effect (SZE) selected galaxy cluster candidates using the Panoramic Survey Telescope & Rapid Response System (Pan-STARRS) imaging data. We first examine 150 Planck confirmed galaxy clusters with spectroscopic redshifts to test our algorithm for identifying optical counterparts and measuring their redshifts; our redshifts have a typical accuracy of $\sigma_{z/(1+z)} \sim 0.022$ for this sample. Using 60 random sky locations, we estimate that our chance of contamination through a random superposition is ~ 3 per cent. We then examine an additional 237 Planck galaxy cluster candidates that have no redshift in the source catalogue. Of these 237 unconfirmed cluster candidates we are able to confirm 60 galaxy clusters and measure their redshifts. A further 83 candidates are so heavily contaminated by stars due to their location near the Galactic plane that we do not attempt to identify counterparts. For the remaining 94 candidates we find no optical counterpart but use the depth of the Pan-STARRS1 data to estimate a redshift lower limit $z_{\text{lim}(10^{15})}$ beyond which we would not have expected to detect enough galaxies for confirmation. Scaling from the already published Planck sample, we expect that $\sim$12 of these unconfirmed candidates may be real clusters.
• The Pan-STARRS1 Small Area Survey 2(1310.6368)

Oct. 23, 2013 astro-ph.CO
The Pan-STARRS1 survey is acquiring multi-epoch imaging in 5 bands (grizy) over the entire sky north of declination -30deg (the $3\pi$ survey). In July 2011 a test area of about 70 sq.deg. was observed to the expected final depth of the main survey. In this, the first of a series of papers targetting the galaxy count and clustering properties of the combined multi-epoch test area data, we present a detailed investigation into the depth of the survey and the reliability of the Pan-STARRS1 analysis software. We show that the Pan-STARRS1 reduction software can recover the properties of fake sources, and show good agreement between the magnitudes measured by Pan-STARRS1 and those from Sloan Digital Sky Survey. We also examine the number of false detections apparent in the Pan-STARRS1 data. Our comparisons show that the test area survey is somewhat deeper than the Sloan Digital Sky Survey in all bands, and, in particular, the z band approaches the depth of the stacked Sloan Stripe 82 data.
• The journey of QSO haloes from z=6 to the present(1203.5339)

Sept. 25, 2013 astro-ph.CO, astro-ph.GA
We apply a recently developed scaling technique to the Millennium-XXL, one of the largest cosmological N-body simulations carried out to date 3x10^11 particles within a cube of volume ~70Gpc^3). This allows us to investigate the cosmological parameter dependence of the mass and evolution of haloes in the extreme high-mass tail of the z=6 distribution. We assume these objects to be likely hosts for the population of rare but ultraluminous high-redshift quasars discovered by the Sloan Digital Sky Survey. Haloes with a similar abundance to these quasars have a median mass of 9x10^12 Msun in the currently preferred cosmology, but do not evolve into equally extreme objects at z=0. Rather, their descendants span the full range conventionally assigned to present-day clusters, 6x10^13 to 2.5x10^15 Msun for this same cosmology. The masses both at z=6 and at z=0 shift up or down by factors exceeding two if cosmological parameters are pushed to the boundaries of the range discussed in published interpretations of data from the WMAP satellite. The main factor determining the future growth of a high-mass z=6 halo is the mean overdensity of its environment on scales of 7 to 14 Mpc, and descendant masses can be predicted 6 to 8 times more accurately if this density is known than if it is not. All these features are not unique to extreme high-z haloes, but are generic to hierarchical growth. Finally, we find that extreme haloes at z=6 typically acquired about half of their total mass in the preceding 100 Myr, implying very large recent accretion rates which may be related to the large black hole masses and high luminosities of the SDSS quasars.
• How well can we really estimate the stellar masses of galaxies from broadband photometry?(1303.7228)

Sept. 4, 2013 astro-ph.CO
The estimated stellar masses of galaxies are widely used to characterize how the galaxy population evolves over cosmic time. If stellar masses can be estimated in a robust manner, free from any bias, global diagnostics such as the stellar mass function can be used to constrain the physics of galaxy formation. We explore how galaxy stellar masses, estimated by fitting broad-band spectral energy distributions (SEDs) with stellar population models, can be biased as a result of commonly adopted assumptions for the star-formation and chemical enrichment histories, recycled fractions and dust attenuation curves of galaxies. We apply the observational technique of broad-band SED fitting to model galaxy SEDs calculated by the theoretical galaxy formation model GALFORM, isolating the effect of each of these assumptions. We find that, averaged over the entire galaxy population, the common assumption of exponentially declining star-formation histories does not adversely affect stellar mass estimation. We show that fixing the metallicity in SED fitting or using sparsely sampled metallicity grids can introduce mass dependent systematics into stellar mass estimates. We find that the common assumption of a star-dust geometry corresponding to a uniform foreground dust screen can cause the stellar masses of dusty model galaxies to be significantly underestimated. Finally, we show that stellar mass functions recovered by applying SED fitting to model galaxies at high redshift can differ significantly in both shape and normalization from the intrinsic mass functions predicted by a given model. Given these differences, our methodology of using stellar masses estimated from model galaxy SEDs offers a new, self-consistent way to compare model predictions with observations.
• Principles of supernova-driven winds(1111.2529)

May 22, 2012 astro-ph.CO, astro-ph.GA
The formation of galaxies is regulated by a balance between the supply of gas and the rate at which it is ejected. Traditional explanations of gas ejection equate the energy required to escape the galaxy or host halo to an estimate for the energy yield from supernovae. This yield is usually assumed to be a constant fraction of the total available from the supernova, or is derived from the assumption of a consistent momentum yield. By applying these ideas in the context of a cold dark matter cosmogony, we derive a 1st-order analytic connection between these working assumptions and the expected relationship between baryon content and galaxy circular velocity, and find that these quick predictions straddle recent observational estimates. To examine the premises behind these theories in more detail, we then explore their applicability to a set of gasdynamical simulations of idealised galaxies. We show that different premises dominate to differing degrees in the simulated outflow, depending on the mass of the system and the resolution with which it is simulated. Using this study to anticipate the emergent behaviour at arbitrarily high resolution, we motivate more comprehensive analytic model which allows for the range of velocities with which the gas may exit the system, and incorporates both momentum and energy-based constraints on the outflow. Using a trial exit velocity distribution, this is shown to be compatible with the observed baryon fractions in intermediate-mass systems, but implies that current estimates for low-mass systems can not be solely accounted for by supernova winds under commonly-held assumptions.
• Substructure in the stellar halos of the Aquarius simulations(1101.2544)

Oct. 31, 2011 astro-ph.GA
We characterize substructure in the simulated stellar halos of Cooper et al. (2010) which were formed by the disruption of satellite galaxies within the cosmological N-body simulations of galactic halos of the Aquarius Project. These stellar halos exhibit a wealth of tidal features: broad overdensities and very narrow faint streams akin to those observed around the Milky Way. The substructures are distributed anisotropically on the sky, a characteristic that should become apparent in the next generation of photometric surveys. The normalized RMS of the density of stars on the sky appears to be systematically larger for our halos compared to the value estimated for the Milky Way from main sequence turn-off stars in the Sloan Digital Sky Survey. We show that this is likely to be due in part to contamination by faint QSOs and redder main sequence stars, and might suggest that ~10% of the Milky Way halo stars have formed in-situ.
• A Two-Point Correlation Function For Galactic Halo Stars(1011.1926)

Sept. 6, 2011 astro-ph.GA
We describe a correlation function statistic that quantifies the amount of spatial and kinematic substructure in the stellar halo. We test this statistic using model stellar halo realizations constructed from the Aquarius suite of six high-resolution N-body simulations in combination with the Galform semi-analytic galaxy formation model. These simulations show considerable scatter in the properties of stellar haloes. We find that our statistic can distinguish between these plausible alternatives for the global structure of the Milky Way stellar halo. We compare with observational data and show that pencil beam surveys of ~100 tracer stars (such as the Spaghetti Survey) are not sufficient to constrain the degree of structure in the Milky Way halo with this statistic. Larger area surveys with >1000 tracer stars (such as BHB stars in the Sloan Digital Sky Survey) provide much tighter constraints on comparisons between models and data. In our simulations, we find examples of haloes with spatial and kinematic substructure consistent with the available Milky Way data.
• The evolution of AGN across cosmic time: what is downsizing?(1011.5222)

Nov. 23, 2010 astro-ph.CO
We use a coupled model of the formation and evolution of galaxies and black holes (BH) to study the evolution of active galactic nuclei (AGN) in a cold dark matter universe. The model predicts the BH mass, spin and mass accretion history. BH mass grows via accretion triggered by discs becoming dynamically unstable or galaxy mergers (called the starburst mode) and accretion from quasi-hydrostatic hot gas haloes (called the hot-halo mode). By taking into account AGN obscuration, we obtain a very good fit to the observed luminosity functions (LF) of AGN (optical, soft and hard X-ray, and bolometric) for a wide range of redshifts (0<z<6). The model predicts a hierarchical build up of BH mass, with the typical mass of actively growing BHs increasing with decreasing redshift. Remarkably, despite this, we find downsizing in the AGN population, in terms of the differential growth with redshift of the space density of faint and bright AGN. This arises naturally from the interplay between the starburst and hot-halo accretion modes. The faint end of the LF is dominated by massive BHs experiencing quiescent accretion via a thick disc, primarily during the hot-halo mode. The bright end of the LF, on the other hand, is dominated by AGN which host BHs accreting close to or in excess of the Eddington limit during the starburst mode. The model predicts that the comoving space density of AGN peaks at z~3, similar to the star formation history. However, when taking into account obscuration, the space density of faint AGN peaks at lower redshift (z<2) than that of bright AGN (z~2-3). This implies that the cosmic evolution of AGN is shaped in part by obscuration.
• Grand unification of AGN activity in the LambdaCDM cosmology(0911.1128)

Aug. 31, 2010 astro-ph.CO
We track the co-evolution of supermassive black holes (SMBHs) and their host galaxies. The calculation is embedded in the GALFORM semi-analytical model which simulates the formation and evolution of galaxies in a cold dark matter (CDM) universe. During the evolution of the host galaxy, hot and cold gas are added to the SMBH by flows triggered by halo gas cooling, disc instabilities and galaxy mergers. This builds up the mass and spin of the BH, and the resulting accretion power regulates the gas cooling and subsequent star formation. The accretion flow is assumed to form a geometrically thin cool disc when the accretion rate exceeds 0.01\dot{M}_Edd, and a geometrically thick, radiatively inefficient hot flow when the accretion rate falls below this value. The resulting quasar optical luminosity function matches observations very well, and the mass of the SMBH correlates with the mass of the galaxy bulge as observed. The BH spin distribution depends strongly on whether the gas in any given accretion episode remains in the same plane (prolonged accretion) or whether, due to self-gravity, it fragments into multiple, randomly aligned accretion episodes (chaotic accretion). In the chaotic accretion model there is a clear correlation of spin with SMBH mass. Massive BHs (M>5\times10^8\Msun) are hosted by giant elliptical galaxies and are rapidly spinning, while lower mass BHs are hosted in spiral galaxies and have much lower spin. Using the Blandford-Znajek mechanism for jet production to calculate the jet power, our model is able to reproduce the radio loudness of radio galaxies, LINERS and Seyferts. This is the first confirmation that a CDM galaxy formation model can reproduce the observed phenomenology of AGN.
• Galactic stellar haloes in the CDM model(0910.3211)

April 14, 2010 astro-ph.CO, astro-ph.GA
We present six simulations of galactic stellar haloes formed by the tidal disruption of accreted dwarf galaxies in a fully cosmological setting. Our model is based on the Aquarius project, a suite of high resolution N-body simulations of individual dark matter haloes. We tag subsets of particles in these simulations with stellar populations predicted by the Galform semi-analytic model. Our method self-consistently tracks the dynamical evolution and disruption of satellites from high redshift. The luminosity function and structural properties of surviving satellites, which agree well with observations, suggest that this technique is appropriate. We find that accreted stellar haloes are assembled between 1<z<7 from less than 5 significant progenitors. These progenitors are old, metal-rich satellites with stellar masses similar to the brightest Milky Way dwarf spheroidals. In contrast to previous stellar halo simulations, we find that several of these major contributors survive as self-bound systems to the present day. Both the number of these significant progenitors and their infall times are inherently stochastic. This results in great diversity among our stellar haloes, which amplifies small differences between the formation histories of their dark halo hosts. The masses and density/surface-brightness profiles of the stellar haloes are consistent with expectations from the Milky Way and M31. Each halo has a complex structure, consisting of well-mixed components, tidal streams, shells and other subcomponents. This structure is not adequately described by smooth models. We find one example of an accreted thick disk. Contrasts in age and metallicity between halo stars and those in surviving satellites are in broad agreement with recent observations. [Abridged]
• Cosmological parameters from CMB measurements and the final 2dFGRS power spectrum(astro-ph/0507583)

March 24, 2006 astro-ph, hep-ph
We derive constraints on cosmological parameters using the power spectrum of galaxy clustering measured from the final two-degree field galaxy redshift survey (2dFGRS) and a compilation of measurements of the temperature power spectrum and temperature-polarization cross-correlation of the cosmic microwave background radiation. We analyse a range of parameter sets and priors, allowing for massive neutrinos, curvature, tensors and general dark energy models. In all cases, the combination of datasets tightens the constraints, with the most dramatic improvements found for the density of dark matter and the energy-density of dark energy. If we assume a flat universe, we find a matter density parameter of $\Omega_{\rm m}=0.237 \pm 0.020$, a baryon density parameter of $\Omega_{\rm b} = 0.041 \pm 0.002$, a Hubble constant of $H_{0}=74\pm2 {\rm kms}^{-1}{\rm Mpc}^{-1}$, a linear theory matter fluctuation amplitude of $\sigma_{8}=0.77\pm0.05$ and a scalar spectral index of $n_{\rm s}=0.954 \pm 0.023$ (all errors show the 68% interval). Our estimate of $n_{\rm s}$ is only marginally consistent with the scale invariant value $n_{\rm s}=1$; this spectrum is formally excluded at the 95% confidence level. However, the detection of a tilt in the spectrum is sensitive to the choice of parameter space. If we allow the equation of state of the dark energy to float, we find $w_{\rm DE}= -0.85_{-0.17}^{+0.18}$, consistent with a cosmological constant. We also place new limits on the mass fraction of massive neutrinos: $f_{\nu} < 0.105$ at the 95% level, corresponding to $\sum m_{\nu} < 1.2$ eV.
• The Effect of Gravitational Recoil on Black Holes Forming in a Hierarchical Universe(astro-ph/0512073)

Dec. 2, 2005 astro-ph
Galactic bulges are known to harbour central black holes whose mass is tightly correlated with the stellar mass and velocity dispersion of the bulge. In a hierarchical universe, mergers of subgalactic units are accompanied by the amalgamation of bulges and the likely coalescence of galactocentric black holes. In these mergers, the beaming of gravitational radiation during the plunge phase of the black hole collision can impart a linear momentum kick or gravitational recoil'' to the remnant. If large enough, this kick will eject the remnant from the galaxy and populate intergalactic space with wandering black holes. Using a semi-analytic model of galaxy formation, we investigate the effect of black hole ejections on the scatter in the relation between black hole and bulge mass. We find that although not the dominant source of the measured scatter, they do make a significant contribution and may be used to set a constraint, v_kick<500 km/s, on the typical kick velocity, in agreement with values found from general relativistic calculations. Even for the more modest kick velocities implied by these calculations, we find that a substantial number of central black holes are ejected from the progenitors of present day galaxies, giving rise to a population of wandering intrahalo and intergalactic black holes whose distribution we investigate in high-resolution N-body simulations of Milk-Way mass halos. We find that intergalactic black holes make up only ~2-3% of the total galactic black hole mass but, within a halo, wandering black holes can contribute up to about half of the total black hole mass orbiting the central galaxy. Intrahalo black holes offer a natural explanation for the compact X-ray sources often seen near the centres of galaxies and for the hyperluminous non-central X-ray source in M82.
• The 2dF Galaxy Redshift Survey: Power-spectrum analysis of the final dataset and cosmological implications(astro-ph/0501174)

Aug. 5, 2005 astro-ph
We present a power spectrum analysis of the final 2dF Galaxy Redshift Survey, employing a direct Fourier method. The sample used comprises 221,414 galaxies with measured redshifts. We investigate in detail the modelling of the sample selection. A new angular mask is derived, based on revisions to the photometric calibration. The redshift selection function is determined by dividing the survey according to rest-frame colour, and deducing a self-consistent treatment of k-corrections and evolution for each population. The covariance matrix for the power-spectrum estimates is determined using two different approaches to the construction of mock surveys which are used to demonstrate that the input cosmological model can be correctly recovered. We are confident that the 2dFGRS power spectrum can be used to infer the matter content of the universe. On large scales, our estimated power spectrum shows evidence for the baryon oscillations' that are predicted in CDM models. Fitting to a CDM model, assuming a primordial $n_{s}=1$ spectrum, $h=0.72$ and negligible neutrino mass, the preferred parameters are $\Omega_{M} h = 0.168 \pm 0.016$ and a baryon fraction $\Omega_{b} /\Omega_{M} = 0.185\pm0.046$ (1$\sigma$ errors). The value of $\Omega_{M} h$ is $1\sigma$ lower than the $0.20 \pm 0.03$ in our 2001 analysis of the partially complete 2dFGRS. This shift is largely due to the signal from the newly-sampled regions of space, rather than the refinements in the treatment of observational selection. This analysis therefore implies a density significantly below the standard $\Omega_{M} =0.3$: in combination with CMB data from WMAP, we infer $\Omega_{M} =0.231\pm 0.021$. (Abridged.)
• The 2dF Galaxy Redshift Survey: luminosity functions by density environment and galaxy type(astro-ph/0407537)

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

Jan. 10, 2005 astro-ph
We present predictions for the abundance of sub-mm galaxies (SMGs) and Lyman-break galaxies (LBGs) in the $\Lambda$CDM cosmology. A key feature of our model is the self-consistent calculation of the absorption and emission of radiation by dust. The new model successfully matches the LBG luminosity function, as well reproducing the properties of the local galaxy population in the optical and IR. The model can also explain the observed galaxy number counts at $850\mum$, but only if we assume a top-heavy IMF for the stars formed in bursts. The predicted redshift distribution of SMGs depends relatively little on their flux over the range 1-$10\mjy$, with a median value of $z\approx 2.0$ at a flux of $5\mjy$, in very good agreement with the recent measurement by Chapman et al The counts of SMGs are predicted to be dominated by ongoing starbursts. However, in the model these bursts are responsible for making only a few per cent of the stellar mass locked up in massive ellipticals at the present day.
• Where are the stars?(astro-ph/0412049)

Dec. 2, 2004 astro-ph
The 2dFGRS is used in conjunction with the 2MASSXSC to study the near-IR light and stellar mass content of the local Universe. Mock galaxy catalogues, constructed from cosmological N-body simulations and semi-analytical galaxy formation models, are used to gauge the accuracy with which quantities can be recovered. The mean luminosity densities of the Universe are found to be rho_J=(3.57+/-0.11)*10^8 h Lsol/Mpc^3 and rho_KS=(7.04+/-0.23)*10^8 h Lsol/Mpc^3 (statistical uncertainty). Using the 2PIGG catalogue, the group dynamical mass-to-light ratio in the K_S band is found to increase by a factor of ~3 when going from groups with total bJ-band luminosities of 3*10^10 h^-2 Lsol to rich clusters, which have typical values of Upsilon_K~80 h Upsilon_sol. Taking into account the bias introduced by uncertainties in estimating galaxy stellar masses from luminosities, a value of Omega_* h=(0.99+/-0.03)*10^-3 is measured, assuming that a Kennicutt stellar IMF is applicable to all galaxies. The 2PIGGs are then used to study the distribution of the stellar content of the local Universe. The three main conclusions are: (1) a slowly rising stellar M/L_KS is found with the clusters having the largest value of ~0.6 Upsilon_sol, (2) in contrast, the fraction of mass in stars decreases with increasing group size, reaching ~5*10^-3 h for the rich clusters, and (3) in answer to the question posed in the title, most stellar mass is contained in Local Group-sized objects (M~2*10^12 h^-1 Msol) with only ~2% in clusters with M>5*10^14 h^-1 Msol.
• The 2dF Galaxy Redshift Survey: stochastic relative biasing between galaxy populations(astro-ph/0404275)

Nov. 26, 2004 astro-ph
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.
• The 2dF Galaxy Redshift Survey: higher order galaxy correlation functions(astro-ph/0401434)

Aug. 23, 2004 astro-ph
We measure moments of the galaxy count probability distribution function in the two-degree field galaxy redshift survey (2dFGRS). The survey is divided into volume limited subsamples in order to examine the dependence of the higher order clustering on galaxy luminosity. We demonstrate the hierarchical scaling of the averaged p-point galaxy correlation functions, xibar_p, up to p=6. The hierarchical amplitudes, S_p = xibar_p/xibar_2^{p-1}, are approximately independent of the cell radius used to smooth the galaxy distribution on small to medium scales. On larger scales we find the higher order moments can be strongly affected by the presence of rare, massive superstructures in the galaxy distribution. The skewness S_3 has a weak dependence on luminosity, approximated by a linear dependence on log luminosity. We discuss the implications of our results for simple models of linear and non-linear bias that relate the galaxy distribution to the underlying mass.
• The 2dF Galaxy Redshift Survey: voids and hierarchical scaling models(astro-ph/0401406)

Aug. 23, 2004 astro-ph
We measure the redshift space reduced void probability function (VPF) for 2dFGRS volume limited galaxy samples covering the absolute magnitude range M_bJ-5logh=-18 to -22. Theoretically, the VPF connects the distribution of voids to the moments of galaxy clustering of all orders, and can be used to discriminate clustering models in the weakly non-linear regime. The reduced VPF measured from the 2dFGRS is in excellent agreement with the paradigm of hierarchical scaling of the galaxy clustering moments. The accuracy of our measurement is such that we can rule out, at a very high significance, popular models for galaxy clustering, including the lognormal distribution. We demonstrate that the negative binomial model gives a very good approximation to the 2dFGRS data over a wide range of scales, out to at least 20h-1Mpc. Conversely, the reduced VPF for dark matter in a LambdaCDM universe does appear to be lognormal on small scales but deviates significantly beyond \approx 4h-1Mpc. We find little dependence of the 2dFGRS reduced VPF on galaxy luminosity. Our results hold independently in both the north and south Galactic pole survey regions.
• The 2dF Galaxy Redshift Survey: The clustering of galaxy groups(astro-ph/0402577)

July 12, 2004 astro-ph
We measure the clustering of galaxy groups in the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. The 2PIGG sample has 29,000 groups with at least two members. The clustering amplitude of the full 2PIGG catalogue is weaker than that of 2dFGRS galaxies, in agreement with theoretical predictions. We have subdivided the 2PIGG catalogue into samples that span a factor of 25 in median total luminosity. Our correlation function measurements span an unprecedented range of clustering strengths, connecting the regimes probed by groups fainter than L* galaxies and rich clusters. There is a steady increase in clustering strength with group luminosity; the most luminous groups are ten times more strongly clustered than the full 2PIGG catalogue. We demonstrate that the 2PIGG results are in very good agreement with the clustering of groups expected in the LCDM model.
• The 2dF Galaxy Redshift Survey: Hierarchical galaxy clustering(astro-ph/0401405)

June 20, 2004 astro-ph
We use the two-degree field Galaxy Redshift Survey (2dFGRS) to test the hierarchical scaling hypothesis: namely, that the $p$-point galaxy correlation functions can be written in terms of the two point correlation function or variance. This scaling is expected if an initially Gaussian distribution of density fluctuations evolves under the action of gravitational instability. We measure the volume averaged $p$-point correlation functions using a counts in cells technique applied to a volume limited sample of 44,931 $L_*$ galaxies. We demonstrate that $L_{*}$ galaxies display hierarchical clustering up to order $p=6$ in redshift space. The variance measured for $L_{*}$ galaxies is in excellent agreement with the predictions from a $\Lambda$-cold dark matter N-body simulation. This applies to all cell radii considered, $0.3<(R/h^{-1}{\rm Mpc})<30$. However, the higher order correlation functions of $L_*$ galaxies have a significantly smaller amplitude than is predicted for the dark matter for $R<10h^{-1}$Mpc. This disagreement implies that a non-linear bias exists between the dark matter and $L_*$ galaxies on these scales. We also show that the presence of two rare, massive superclusters in the 2dFGRS has an impact on the higher-order clustering moments measured on large scales.
• The nature of the relative bias between galaxies of different spectral type in 2dFGRS(astro-ph/0404276)

April 14, 2004 astro-ph
We present an analysis of the relative bias between early- and late-type galaxies in the Two-degree Field Galaxy Redshift Survey (2dFGRS). Our analysis examines the joint counts in cells between early- and late-type galaxies, using approximately cubical cells with sides ranging from 7h^{-1}Mpc to 42h^{-1}Mpc. We measure the variance of the counts in cells using the method of Efstathiou et al. (1990), which we find requires a correction for a finite volume effect. We fit lognormal models to the one-point density distribution and develop methods of dealing with biases in the recovered variances resulting from this technique. We directly fit deterministic models for the joint density distribution function, f(delta_E,delta_L), to the joint counts in cells using a maximum likelihood technique. Our results are consistent with a scale invariant relative bias factor on all scales studied. Linear bias is ruled out on scales less than l=28h^{-1}Mpc. A power-law bias model is a significantly better fit to the data on all but the largest scales studied; the relative goodness of fit of this model as compared to that of the linear bias model suggests that any nonlinearity is negligible for l>~40h^{-1}Mpc, consistent with the expectation from theory that the bias should become linear on large scales. (abridged)
• Galaxy groups in the 2dFGRS: the luminous content of the groups(astro-ph/0402566)

Feb. 26, 2004 astro-ph
The 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue of ~29000 objects is used to study the luminous content of galaxy systems of various sizes. Mock galaxy catalogues constructed from cosmological simulations are used to gauge the accuracy with which intrinsic group properties can be recovered. A Schechter function is found to be a reasonable fit to the galaxy luminosity functions in groups of different mass in the real data. The characteristic luminosity L* is larger for more massive groups. However, the mock data show that the shape of the recovered luminosity function is expected to differ from the true shape, and this must be allowed for when interpreting the data. The variation of halo mass-to-light ratio with group size is studied in both these wavebands. A robust trend of increasing M/L with increasing group luminosity is found in the 2PIGG data. From groups with L_bj=10^{10}Lsol to those 100 times more luminous, the typical bj-band M/L increases by a factor of 5, whereas the rf-band M/L grows by a factor of 3.5. These trends agree well with the simulations, which also predict a minimum M/L on a scale corresponding to the Local Group. Our data indicate that if such a minimum exists, then it must occur at L<~10^{10}Lsol, below the range accurately probed by the 2PIGG catalogue. According to the mock data, the bj M/Ls of the largest groups are expected to be approximately 1.1 times the global value. Assuming that this correction applies to the real data yields an estimate of Omega_m=0.26+/-0.03 (statistical).