We consider black holes which form from an initially spherically symmetric super-Hubble perturbation of a cosmological background filled by a perfect fluid $p=w \rho$ with $w\in (0,1]$. Previous work has shown that when $w = 1/3$ (radiation), there is a critical threshold for black hole formation ($\delta_c$), which, to a very good approximation, only depends upon the curvature of the compaction function around its peak value. We find that this generalizes to all $w\gtrsim 1/3$; for smaller $w$s the knowledge of the full shape of the compaction function is necessary. We provide analytic approximations for $\delta_c$ which are accurate for $w\in [1/3,1]$.

This is the third paper of a series where we study the stellar population gradients (SP; ages, metallicities, $\alpha$-element abundance ratios and stellar initial mass functions) of early type galaxies (ETGs) at $z\le 0.08$ from the MaNGA-DR15 survey. In this work we focus on the S0 population and quantify how the SP varies across the population as well as with galactocentric distance. We do this by measuring Lick indices and comparing them to stellar population synthesis models. This requires spectra with high signal-to-noise which we achieve by stacking in bins of luminosity (L$_r$) and central velocity dispersion ($\sigma_0$). We find that: 1) There is a bimodality in the S0 population: S0s more massive than $3\times 10^{10}M_\odot$ show stronger velocity dispersion and age gradients (age and $\sigma_r$ decrease outwards) but little or no metallicity gradient, while the less massive ones present relatively flat age and velocity dispersion profiles, but a significant metallicity gradient (i.e. [M/H] decreases outwards). Above $2\times10^{11}M_\odot$ the number of S0s drops sharply. These two mass scales are also where global scaling relations of ETGs change slope. 2) S0s have steeper velocity dispersion profiles than fast rotating elliptical galaxies (E-FRs) of the same luminosity and velocity dispersion. The kinematic profiles and stellar population gradients of E-FRs are both more similar to those of slow rotating ellipticals (E-SRs) than to S0s, suggesting that E-FRs are not simply S0s viewed face-on. 3) At fixed $\sigma_0$, more luminous S0s and E-FRs are younger, more metal rich and less $\alpha$-enhanced. Evidently for these galaxies, the usual statement that 'massive galaxies are older' is not true if $\sigma_0$ is held fixed.

Early-type galaxies -- slow and fast rotating ellipticals (E-SRs and E-FRs) and S0s/lenticulars -- define a Fundamental Plane (FP) in the space of half-light radius $R_e$, enclosed surface brightness $I_e$ and velocity dispersion $\sigma_e$. Since $I_e$ and $\sigma_e$ are distance-independent measurements, the thickness of the FP is often expressed in terms of the accuracy with which $I_e$ and $\sigma_e$ can be used to estimate sizes $R_e$. We show that: 1) The thickness of the FP depends strongly on morphology. If the sample only includes E-SRs, then the observed scatter in $R_e$ is $\sim 16\%$, of which only $\sim 9\%$ is intrinsic. Removing galaxies with $M_*<10^{11}M_\odot$ further reduces the observed scatter to $\sim 13\%$ ($\sim 4\%$ intrinsic). The observed scatter increases to the $\sim 25\%$ usually quoted in the literature if E-FRs and S0s are added. If the FP is defined using the eigenvectors of the covariance matrix of the observables, then the E-SRs again define an exceptionally thin FP, with intrinsic scatter of only $5\%$ orthogonal to the plane. 2) The structure within the FP is most easily understood as arising from the fact that $I_e$ and $\sigma_e$ are nearly independent, whereas the $R_e-I_e$ and $R_e-\sigma_e$ correlations are nearly equal and opposite. 3) If the coefficients of the FP differ from those associated with the virial theorem the plane is said to be `tilted'. If we multiply $I_e$ by the global stellar mass-to-light ratio $M_*/L$ and we account for non-homology across the population by using S\'ersic photometry, then the resulting stellar mass FP is less tilted. Accounting self-consistently for $M_*/L$ gradients will change the tilt. The tilt we currently see suggests that the efficiency of turning baryons into stars increases and/or the dark matter fraction decreases as stellar surface brightness increases.

We develop the non-linear statistics of primordial black holes generated by a gaussian spectrum of primordial curvature perturbations. This is done by employing the compaction function as the main statistical variable under the constraints that: a) the over-density has a high peak at a point $\vec{x}_0$, b) the compaction function has a maximum at a smoothing scale $R$, and finally, c) the compaction function amplitude at its maximum is higher than the threshold necessary to trigger a gravitational collapse into a black hole of the initial over-density. Our calculation allows for the fact that the patches which are destined to form PBHs may have a variety of profile shapes and sizes. The predicted PBH abundances depend on the power spectrum of primordial fluctuations. For a very peaked power spectrum, our non-linear statistics, the one based on the linear over-density and the one based on the use of curvature perturbations, all predict a narrow distribution of PBH masses and comparable abundance. For broader power spectra the linear over-density statistics over-estimate the abundance of primordial black holes while the curvature-based approach under-estimates it. Additionally, for very large smoothing scales, the abundance is no longer dominated by the contribution of a mean over-density but rather by the whole statistical realisations of it.

In this letter, we argue and show numerically that the threshold to form primordial black holes from an initial spherically symmetric perturbation is, to an excellent approximation, universal, whenever given in terms of the compaction function averaged over a sphere of radius $r_m$, where $r_m$ is the scale on which the compaction function is maximum. This can be understood as the requirement that, for a black hole to form, each shell of the averaged compaction function should have an amplitude exceeding the so-called Harada-Yoo-Kohri limit. For a radiation dominated universe we argued, supported by the numerical simulations, that this limit is $\delta_c = 0.40$, which is slightly below the one quoted in the literature. Additionally, we show that the profile dependence of the threshold for the compaction function is only sensitive to its curvature at the maximum. We use these results to provide an analytic formula for the threshold amplitude of the compaction function at its maximum in terms of the normalised compaction function curvature at $r_m$.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

The Baryon Acoustic Oscillations feature (BAO) imprinted in the clustering correlation function is known to furnish us cosmic distance determinations that are independent of the cosmological-background model and the primordial perturbation parameters. These measurements can be accomplished rigorously by means of the Purely Geometric BAO methods. To date two different Purely Geometric BAO approaches have been proposed. The first exploits the linear-point standard ruler. The second, called correlation-function model-fitting, exploits the sound-horizon standard ruler. A key difference between them is that, when estimated from clustering data, the linear point makes use of a cosmological-model-independent procedure to extract the ratio of the ruler to the cosmic distance, while the correlation-function model-fitting relies on a phenomenological cosmological model for the correlation function. Nevertheless the two rulers need to be precisely defined independently of any specific observable. We define the linear point and sound horizon and we characterize and compare the two rulers' cosmological-parameter dependence. We find that they are both geometrical within the required accuracy, and they have the same parameter dependence for a wide range of parameter values. We estimate the rulers' best-fit values and errors given the cosmological constraints obtained by the Planck Satellite team from the CMB measurements. We do this for three different cosmological models encompassed by the Purely Geometric BAO methods. In each case we find that the relative errors of the two rulers coincide and they are insensitive to the assumed cosmological model. Interestingly both the linear point and the sound horizon shift by $0.5\sigma$ when we do not fix the spatial geometry to be flat in LCDM. This points toward a sensitivity of the rulers to different cosmological models when they are estimated from the CMB.

A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demography of supermassive black holes thus preventing accurate model comparisons and progress in this field. Here we show that the large-scale spatial distribution of local active galactic nuclei (AGN), believed to be accreting supermassive black holes, can constrain the shape and normalization of the black hole-stellar mass relation thus bypassing resolution-related observational biases. In turn, our results can set more stringent constraints on the so-called "radiative efficiency", a fundamental parameter describing the inner physics of supermassive black holes that is closely linked to their spin, geometry, and ability to release energy. The mean value of the radiative efficiency can be estimated by comparing the average total luminous output of AGN with the relic mass density locked up in quiescent supermassive black holes at galaxy centres today. For currently accepted values of the AGN obscured fractions and bolometric corrections, our newest estimates of the local supermassive black hole mass density favour mean radiative efficiencies of ~10-20\%, suggesting that the vast majority of supermassive black holes are spinning moderately to rapidly. With large-scale AGN surveys coming online, our novel methodology will enable even tighter constraints on the fundamental parameters that regulate the growth of supermassive black holes.

We present estimates of stellar population (SP) gradients from stacked spectra of slow (SR) and fast (FR) rotator elliptical galaxies from the MaNGA-DR15 survey. We find that: 1) FRs are $\sim 5$ Gyrs younger, more metal rich, less $\alpha$-enhanced and smaller than SRs of the same luminosity $L_r$ and central velocity dispersion $\sigma_0$. This explains why when one combines SRs and FRs, objects which are small for their $L_r$ and $\sigma_0$ tend to be younger. Their SP gradients are also different. 2) Ignoring the FR/SR dichotomy leads one to conclude that compact galaxies are older than their larger counterparts of the same mass, even though almost the opposite is true for FRs and SRs individually. 3) SRs with $\sigma_0\le 250$ km s$^{-1}$ are remarkably homogeneous within $\sim R_e$: they are old, $\alpha$-enhanced and only slightly super-solar in metallicity. These SRs show no gradients in age and $M_*/L_r$, negative gradients in metallicity, and slightly positive gradients in [$\alpha$/Fe] (the latter are model dependent). SRs with $\sigma_0\ge 250$ km $s^{-1}$ are slightly younger and more metal rich, contradicting previous work suggesting that age increases with $\sigma_0$. They also show larger $M_*/L_r$ gradients. 4) Self-consistently accounting for $M_*/L$ gradients yields $M_{\rm dyn}\approx M_*$ because gradients reduce $M_{\rm dyn}$ by $\sim 0.2$ dex while only slightly increasing the $M_*$ inferred using a Kroupa (not Salpeter) IMF. 5) The FR population all but disappears above $M_*\ge 3\times 10^{11}M_\odot$; this is the same scale at which the size-mass correlation and other scaling relations change. Our results support the finding that this is an important mass scale which correlates with the environment and above which mergers matter.

We estimate ages, metallicities, $\alpha$-element abundance ratios and stellar initial mass functions of elliptical (E) and S0 galaxies from the MaNGA-DR15 survey. We stack spectra and use a variety of single stellar population synthesis models to interpret the absorption line strengths in these spectra. We quantify how these properties vary across the population, as well as with galactocentric distance. This paper is the first of a series and is based on a sample of pure elliptical galaxies at z $\le$ 0.08. We show that the properties of the inner regions of Es with the largest luminosity (L$_r$) and central velocity dispersion ($\sigma_0$) are consistent with those associated with the commonly used Salpeter IMF, whereas a Kroupa-like IMF is a better description at $\sim$ 0.8R/Re (assuming [Ti/Fe] variations are limited). For these galaxies the stellar mass-to-light ratio decreases at most by a factor of 2 from the central regions to Re. In contrast, for lower L$_r$ and $\sigma_0$ galaxies, the IMF is shallower and M$_{*}$/L$_r$ in the central regions is similar to the outskirts. Although a factor of 2 is smaller than previous reports based on a handful of galaxies, it is still large enough to matter for dynamical mass estimates. Accounting self-consistently for these gradients when estimating both M$_{*}$ and M$_{dyn}$ brings the two into good agreement: gradients reduce M$_{dyn}$ by $\sim$ 0.2 dex while only slightly increasing the M$_{*}$ inferred using a Kroupa IMF. This is a different resolution of the M$_{*}$-M$_{dyn}$ discrepancy than has been followed in the recent literature where M$_{*}$ of massive galaxies is increased by adopting a Salpeter IMF while leaving Mdyn unchanged. A companion paper discusses how stellar population differences are even more pronounced if one separates slow from fast rotators.

The internal properties of dark matter haloes correlate with the large-scale halo clustering strength at fixed halo mass $-$ an effect known as assembly bias $-$ and are also strongly affected by the local, non-linear cosmic web. Characterising a halo's local web environment by its tidal anisotropy $\alpha$ at scales $\sim4$ x the halo radius, we demonstrate that these multi-scale correlations represent two distinct statistical links: one between the internal property and $\alpha$, and the other between $\alpha$ and large-scale ( $>30h^{-1}$Mpc) halo bias $b_1$. We focus on scalar internal properties of haloes related to formation time (concentration $c_{\rm vir}$), shape (mass ellipsoid asphericity $c/a$), velocity dispersion structure (velocity ellipsoid asphericity $c_v/a_v$ and velocity anisotropy $\beta$) and angular momentum (dimensionless spin $\lambda$) in the mass range $8\times10^{11}< M_{\rm vir}/(h^{-1}M_\odot)<5\times10^{14}$. Using conditional correlation coefficients and other detailed tests, we show that the joint distribution of $\alpha$, $b_1$ and any of the internal properties $c\in\{\beta,c_v/a_v,c/a,c_{\rm vir},\lambda\}$ is consistent with $p(\alpha,b_1,c)\simeq p(\alpha)p(b_1|\alpha)p(c|\alpha)$, at all but the largest masses. $\textit{Thus, the assembly bias trends $c-b_1$ reflect the two fundamental correlations $c-\alpha$ and $b_1-\alpha$.}$ Our results are unaffected by the exclusion of haloes with recent major merger events or splashback objects, although the latter are distinguished by the fact that $\alpha$ does not explain their assembly bias trends. The overarching importance of $\alpha$ provides a new perspective on the nature of assembly bias of distinct haloes, with potential ramifications for incorporating realistic assembly bias effects into mock catalogs of future large-scale structure surveys and for detecting galaxy assembly bias.

The simplest models of dark matter halo formation are based on the heuristic assumption, motivated by spherical collapse, that virialized haloes originate from initial regions that are maxima of the smoothed density field. Here, we replace this notion with the dynamical requirement that protohalo patches be regions where the local gravitational flow converges to a point. For this purpose, we look for spheres whose gravitational acceleration at the boundary -- relative to their center of mass -- points towards their geometric center: that is, spheres with null dipole moment. We show that these configurations are minima of the total energy, i.e. the most energetically bound spheres. For this reason, we study peaks of the energy overdensity field, and argue that the approach shows considerable promise. This change simply requires that one modify the standard top-hat filter, with the added important benefit that, for power spectra of cosmological interest, the resulting model is no longer plagued by divergences. Although the formalism is no more complicated than the overdensity based approach, the model is richer in the sense that it naturally predicts scatter in the overdensities of protohalo patches that are destined to form haloes of the same mass, in qualitative agreement with simulations of halo formation.

We propose an analytical approach to study non-Markov random walks by employing an exact enumeration method. Using the method, we derive an exact expansion for the first-passage time (FPT) distribution for any continuous, differentiable non-Markov random walk with Gaussian or non-Gaussian multivariate distribution. As an example, we study the FPT distribution of a fractional Brownian motion with a Hurst exponent $H\in(1/2,1)$ that describes numerous non-Markov stochastic phenomena in physics, biology and geology, and for which the limit $H=1/2$ represents a Markov process.

An understanding of astrophysical feedback is important for constraining models of galaxy formation and for extracting cosmological information from current and future weak lensing surveys. The thermal Sunyaev-Zel'dovich effect, quantified via the Compton-$y$ parameter, is a powerful tool for studying feedback, because it directly probes the pressure of the hot, ionized gas residing in dark matter halos. Cross-correlations between galaxies and maps of Compton-$y$ obtained from cosmic microwave background surveys are sensitive to the redshift evolution of the gas pressure, and its dependence on halo mass. In this work, we use galaxies identified in year one data from the Dark Energy Survey and Compton-$y$ maps constructed from Planck observations. We find highly significant (roughly $12\sigma$) detections of galaxy-$y$ cross-correlation in multiple redshift bins. By jointly fitting these measurements as well as measurements of galaxy clustering, we constrain the halo bias-weighted, gas pressure of the Universe as a function of redshift between $0.15 \lesssim z \lesssim 0.75$. We compare these measurements to predictions from hydrodynamical simulations, allowing us to constrain the amount of thermal energy in the halo gas relative to that resulting from gravitational collapse.

The central regions of nearby elliptical galaxies are dominated by baryons (stars) and provide interesting laboratories for studying the radial acceleration relation (RAR). We carry out exploratory analyses and discuss the possibility of constraining the RAR in the super-critical acceleration range $(10^{-9.5},\hspace{1ex}10^{-8})$~${\rm m}~{\rm s}^{-2}$ using a sample of nearly round pure-bulge (spheroidal, dispersion-dominated) galaxies including 24 ATLAS$^{\rm 3D}$ galaxies and 4201 SDSS galaxies covering a wide range of masses, sizes and luminosity density profiles. We consider a range of current possibilities for the stellar mass-to-light ratio ($M_\star/L$), its gradient and dark or phantom matter (DM/PM) halo profiles. We obtain the probability density functions (PDFs) of the parameters of the considered models via Bayesian inference based on spherical Jeans Monte Carlo modeling of the observed velocity dispersions. We then constrain the DM/PM-to-baryon acceleration ratio $a_{\rm X}/a_{\rm B}$ from the PDFs. Unless we ignore observed radial gradients in $M_\star/L$, or assume unreasonably strong gradients, marginalization over nuisance factors suggests $a_{\rm X}/a_{\rm B} = 10^{p} (a_{\rm B}/a_{+1})^q$ with $p = -1.00 \pm 0.03$ (stat) $^{+0.11}_{-0.06}$ (sys) and $q=-1.02 \pm 0.09$ (stat) $^{+0.16}_{-0.00}$ (sys) around a super-critical acceleration $a_{+1}\equiv 1.2\times 10^{-9}~{\rm m}~{\rm s}^{-2}$. In the context of the $\Lambda$CDM paradigm, this RAR suggests that the NFW DM halo profile is a reasonable description of galactic halos even after the processes of galaxy formation and evolution. In the context of the MOND paradigm, this RAR favors the Simple interpolating function but is inconsistent with the vast majority of other theoretical proposals and fitting functions motivated mainly from sub-critical acceleration data.

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

We investigate the anisotropy of the stellar velocity dispersions within the effective radius, $R_{\rm e}$, in 24 ATLAS$^{\rm 3D}$ pure-bulge galaxies, 16 of which are kinematic slow rotators (SRs). We allow the spherical anisotropy parameter $\beta$ to be radially varying and allow a radial gradient in the stellar mass-to-light ratio ($M_\star/L$) through the parameter $K$ introduced earlier. The median anisotropy for SRs depends on $K$ as follows: $\langle\beta_{\rm m}\rangle = a + b K$ with $a=0.19\pm 0.05$, $b=-0.13\pm 0.07$ ($\Lambda$CDM) or $a=0.21\pm 0.05$, $b=-0.26\pm 0.08$ (MOND), where $\beta_{\rm m}$ refers to the radially averaged quantity. Under the $\Lambda$CDM paradigm this scaling is tied to a scaling of $\langle f_{\rm DM}\rangle = (0.16\pm 0.03) +(0.31\pm 0.06) K$, where $f_{\rm DM}$ refers to the DM fraction within a sphere of $r=R_{\rm e}$. For $K=0$ (constant $M_\star/L$), we obtain radially biased results with $\langle\beta_{\rm m}\rangle \approx 0.2$ consistent with previous results. However, marginalizing over $0 < K < 1.5$ yields $\langle\beta_{\rm m}\rangle = 0.06 ^{+0.11}_{-0.14}$ with $\langle f_{\rm DM}\rangle = 0.35 \pm 0.08$: isotropy is preferred. This isotropy hides the fact that $\beta_{\rm m}$ is correlated with kinematic features such as counter rotating cores (CRCs), kinematically distinct cores (KDCs), and low-level velocities (LVs): SRs with LVs are likely to be radially biased while SRs with CRCs are likely to be tangentially biased, and SRs with KDCs are intermediate. Existing cosmological simulations allow us to understand these results qualitatively in terms of their dynamical structures and formation histories although there exist quantitative tensions. More realistic cosmological simulations, particularly allowing for $M_\star/L$ gradients, may be required to better understand SRs.

Local samples of quiescent galaxies with dynamically measured black hole masses (Mbh) may suffer from an angular resolution-related selection effect, which could bias the observed scaling relations between Mbh and host galaxy properties away from the intrinsic relations. In particular, previous work has shown that the observed Mbh-Mstar (stellar mass) relation is more strongly biased than the Mbh-sigma (velocity dispersion) relation. Local samples of active galactic nuclei (AGN) do not suffer from this selection effect, as in these samples Mbh is estimated from megamasers and/or reverberation mapping-based techniques. With the exception of megamasers, Mbh-estimates in these AGN samples are proportional to a virial coefficient fvir. Direct modelling of the broad line region suggests that fvir~3.5. However, this results in a Mbh-Mstar relation for AGN which lies below and is steeper than the one observed for quiescent black hole samples. A similar though milder trend is seen for the Mbh-sigma relation. Matching the high-mass end of the Mbh-Mstar and Mbh-sigma relations observed in quiescent samples requires fvir~15 and fvir~7, respectively. On the other hand, fvir~3.5 yields Mbh-sigma and Mbh-Mstar relations for AGN which are remarkably consistent with the expected `intrinsic' correlations for quiescent samples (i.e., once account has been made of the angular resolution-related selection effect), providing additional evidence that the sample of local quiescent black holes is biased. We also show that, as is the case for quiescent black holes, the Mbh-Mstar scaling relation of AGN is driven by velocity dispersion, thus providing additional key constraints to black hole-galaxy co-evolution models.

Leveraging the Baryon Acoustic Oscillations (BAO) feature present in clustering 2-point statistics, we aim to measure cosmological distances independently of the underlying background cosmological model. However this inference is complicated by late-time non-linearities that introduce model and tracer dependencies in the clustering correlation function and power spectrum, which must be properly accounted for. With this in mind, we introduce the "Purely-Geometric-BAO," which provides a rigorous tool to measure cosmological distances without assuming a specific background cosmology. We focus on the 2-point clustering correlation function monopole, and show how to implement such an inference scheme employing two different methodologies: the Linear Point standard ruler (LP) and correlation-function model-fitting (CF-MF). For the first time we demonstrate how, by means of the CF-MF, we can measure very precisely the sound-horizon/isotropic-volume-distance ratio, $r_{d}/D_{V}(\bar{z})$, while correctly propagating all the uncertainties. Using synthetic data, we compare the outcomes of the two methodologies, and find that the LP provides up to $50\%$ more precise measurements than the CF-MF. Finally, we test a procedure widely employed in BAO analyses: fitting the 2-point function while fixing the cosmological and the non-linear-damping parameters at fiducial values. We find that this underestimates the distance errors by nearly a factor of $2$. We thus recommend that this practice be reconsidered, whether for parameter determination or model selection.

We study the evolution of the cross-correlation between voids and the mass density field - i.e. of void profiles. We show that approaches based on the spherical model alone miss an important contribution to the evolution on large scales of most interest to cosmology: they fail to capture the well-known fact that the large-scale bias factor of conserved tracers evolves. We also show that the operations of evolution and averaging do not commute, but this difference is only significant within about two effective radii. We show how to include a term which accounts for the evolution of bias, which is directly related to the fact that voids move. The void motions are approximately independent of void size, so they are more significant for smaller voids that are typically more numerous. This term also contributes to void-matter pairwise velocities: including it is necessary for modeling the typical outflow speeds around voids. It is, therefore, important for void redshift space distortions. Finally, we show that the excursion set peaks/troughs approach provides a useful, but not perfect framework for describing void profiles and their evolution.

We carry out spherical Jeans modeling of nearly round pure-bulge galaxies selected from the ATLAS$^{\rm 3D}$ sample. Our modeling allows for gradients in the stellar mass-to-light ratio ($M_\star/L$) through analytic prescriptions parameterized with a `gradient strength' $K$ introduced to accommodate any viable gradient. We use a generalized Osipkov-Merritt model for the velocity dispersion (VD) anisotropy. We produce Monte Carlo sets of models based on the stellar VD profiles under both the $\Lambda$CDM and MOND paradigms. Here, we describe the galaxy data, the empirical inputs, and the modeling procedures of obtaining the Monte Carlo sets. We then present the projected dynamical stellar mass, $M_{\rm \star e}$, within the effective radius $R_{\rm e}$, and the fundamental mass plane (FMP) as a function of $K$. We find the scaling of the $K$-dependent mass with respect to the ATLAS$^{\rm 3D}$ reported mass as: $\log_{10} \left[M_{\star{\rm e}}(K)/M_{\star{\rm e}}^{\rm A3D} \right] = a' + b' K$ with $a'=-0.019\pm 0.012$ and $b'=-0.18\pm 0.02$ ($\Lambda$CDM), or $a'=-0.023\pm 0.014$ and $b'=-0.23\pm 0.03$ (MOND), for $0\le K < 1.5$. The FMP has coefficients consistent with the virial expectation and only the zero-point scales with $K$. The median value of $K$ for the ATLAS$^{\rm 3D}$ galaxies is $\langle K\rangle =0.53^{+0.05}_{-0.04}$. We perform a similar analysis of the much larger SDSS DR7 spectroscopic sample. In this case, only the VD within a single aperture is available, so we impose the additional requirement that the VD slope be similar to that in the ATLAS$^{\rm 3D}$ galaxies. Our analysis of the SDSS galaxies suggests a positive correlation of $K$ with stellar mass.

The excursion set approach is a framework for estimating how the number density of nonlinear structures in the cosmic web depends on the expansion history of the universe and the nature of gravity. A key part of the approach is the estimation of the first crossing distribution of a suitably chosen barrier by random walks having correlated steps: The shape of the barrier is determined by the physics of nonlinear collapse, and the correlations between steps by the nature of the initial density fluctuation field. We describe analytic and numerical methods for calculating such first up-crossing distributions. While the exact solution can be written formally as an infinite series, we show how to approximate it efficiently using the Stratonovich approximation. We demonstrate its accuracy using Monte-Carlo realizations of the walks, which we generate using a novel Cholesky-decomposition based algorithm, which is significantly faster than the algorithm that is currently in the literature.

Within a galaxy the stellar mass-to-light ratio $\Upsilon_*$ is not constant. Spatially resolved kinematics of nearby early-type galaxies suggest that allowing for a variable initial mass function (IMF) returns significantly larger $\Upsilon_*$ gradients than if the IMF is held fixed. If $\Upsilon_*$ is greater in the central regions, then ignoring the IMF-driven gradient can overestimate $M_*^{\rm dyn}$ by as much as a factor of two for the most massive galaxies, though stellar population estimates $M_*^{\rm SP}$ are also affected. Large $\Upsilon_*$-gradients have four main consequences: First, $M_*^{\rm dyn}$ cannot be estimated independently of stellar population synthesis models. Second, if there is a lower limit to $\Upsilon_*$ and gradients are unknown, then requiring $M_*^{\rm dyn}=M_*^{\rm SP}$ constrains them. Third, if gradients are stronger in more massive galaxies, then $M_*^{\rm dyn}$ and $M_*^{\rm SP}$ can be brought into agreement, not by shifting $M_*^{\rm SP}$ upwards by invoking constant bottom-heavy IMFs, as advocated by a number of recent studies, but by revising $M_*^{\rm dyn}$ estimates in the literature downwards. Fourth, accounting for $\Upsilon_*$ gradients changes the high-mass slope of the stellar mass function $\phi(M_*^{\rm dyn})$, and reduces the associated stellar mass density. These conclusions potentially impact estimates of the need for feedback and adiabatic contraction, so our results highlight the importance of measuring $\Upsilon_*$ gradients in larger samples.

We use a 200 $h^{-1}Mpc$ a side N-body simulation to study the mass accretion history (MAH) of dark matter halos to be accreted by larger halos, which we call infall halos. We define a quantity $a_{\rm nf}\equiv (1+z_{\rm f})/(1+z_{\rm peak})$ to characterize the MAH of infall halos, where $z_{\rm peak}$ and $z_{\rm f}$ are the accretion and formation redshifts, respectively. We find that, at given $z_{\rm peak}$, their MAH is bimodal. Infall halos are dominated by a young population at high redshift and by an old population at low redshift. For the young population, the $a_{\rm nf}$ distribution is narrow and peaks at about $1.2$, independent of $z_{\rm peak}$, while for the old population, the peak position and width of the $a_{\rm nf}$ distribution both increases with decreasing $z_{\rm peak}$ and are both larger than those of the young population. This bimodal distribution is found to be closely connected to the two phases in the MAHs of halos. While members of the young population are still in the fast accretion phase at $z_{\rm peak}$, those of the old population have already entered the slow accretion phase at $z_{\rm peak}$. This bimodal distribution is not found for the whole halo population, nor is it seen in halo merger trees generated with the extended Press-Schechter formalism. The infall halo population at $z_{\rm peak}$ are, on average, younger than the whole halo population of similar masses identified at the same redshift. We discuss the implications of our findings in connection to the bimodal color distribution of observed galaxies and to the link between central and satellite galaxies.

We study the role of the local tidal environment in determining the assembly bias of dark matter haloes. Previous results suggest that the anisotropy of a halo's environment (i.e, whether it lies in a filament or in a more isotropic region) can play a significant role in determining the eventual mass and age of the halo. We statistically isolate this effect using correlations between the large-scale and small-scale environments of simulated haloes at $z=0$ with masses between $10^{11.6}\lesssim (m/h^{-1}M_{\odot})\lesssim10^{14.9}$. We probe the large-scale environment using a novel halo-by-halo estimator of linear bias. For the small-scale environment, we identify a variable $\alpha_R$ that captures the $\textit{tidal anisotropy}$ in a region of radius $R=4R_{\textrm{200b}}$ around the halo and correlates strongly with halo bias at fixed mass. Segregating haloes by $\alpha_R$ reveals two distinct populations. Haloes in highly isotropic local environments ($\alpha_R\lesssim0.2$) behave as expected from the simplest, spherically averaged analytical models of structure formation, showing a $\textit{negative}$ correlation between their concentration and large-scale bias at $\textit{all}$ masses. In contrast, haloes in anisotropic, filament-like environments ($\alpha_R\gtrsim0.5$) tend to show a $\textit{positive}$ correlation between bias and concentration at any mass. Our multi-scale analysis cleanly demonstrates how the overall assembly bias trend across halo mass emerges as an average over these different halo populations, and provides valuable insights towards building analytical models that correctly incorporate assembly bias. We also discuss potential implications for the nature and detectability of galaxy assembly bias.

The influence of the Cosmic Web on galaxy formation and evolution is of great observational and theoretical interest. We investigate whether the Cosmic Web leaves an imprint in the spatial clustering of galaxies in the Sloan Digital Sky Survey (SDSS), using the group catalog of Yang et al. and tidal field estimates at $\sim2h^{-1}$Mpc scales from the Mass-Tides-Velocity data set of Wang et al. We use the $\textit{tidal anisotropy}$ $\alpha$ (Paranjape et al.) to characterise the tidal environment of groups, and measure the redshift-space 2-point correlation function (2pcf) of group positions and the luminosity- and colour-dependent clustering of group galaxies using samples segregated by $\alpha$. We find that all the 2pcf measurements depend strongly on $\alpha$, with factors of $\sim20$ between the large-scale 2pcf of objects in the most and least isotropic environments. To test whether these strong trends imply `beyond halo mass' effects for galaxy evolution, we compare our results with corresponding 2pcf measurements in mock catalogs constructed using a halo occupation distribution that only uses halo mass as an input. We find that this prescription qualitatively reproduces $\textit{all}$ observed trends, and also quantitatively matches many of the observed results. Although there are some statistically significant differences between our `halo mass only' mocks and the data -- in the most and least isotropic environments -- which deserve further investigation, our results suggest that if the tidal environment induces additional effects on galaxy properties other than those inherited from their host halos, then these must be weak.

Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy ($M_*^{\rm JAM}/M_*$) correlates with $\sigma_e$, if $M_*$ is estimated using the same IMF for all galaxies and the stellar M/L ratio within each galaxy is constant. This correlation may indicate that, in fact, the IMF is more dwarf-rich for galaxies with large $\sigma$. We use this correlation to estimate a dynamical or IMF-corrected stellar mass, $M_*^{\rm \alpha_{JAM}}$, from $M_{*}$ and $\sigma_e$ for a sample of $6 \times 10^5$ SDSS galaxies for which spatially resolved kinematics is not available. We also compute the `virial' mass estimate $k(n,R)\,R_e\,\sigma_R^2/G$, where $n$ is the Sersic index, in the SDSS and ATLAS$^{\rm 3D}$ samples. We show that an $n$-dependent correction must be applied to the $k(n,R)$ values provided by Prugniel & Simien (1997). Our analysis also shows that the shape of the velocity dispersion profile in the ATLAS$^{\rm 3D}$ sample varies weakly with $n$: $(\sigma_R/\sigma_e) = (R/R_e)^{-\gamma(n)}$. The resulting stellar mass functions, based on $M_*^{\rm \alpha_{JAM}}$ and the recalibrated virial mass, are in good agreement. If the $M_*^{\rm \alpha_{JAM}}/M_* - \sigma_e$ correlation is indeed due to the IMF, and stellar M/L gradients can be ignored, then our $\phi(M_*^{\rm \alpha_{JAM}})$ is an estimate of the stellar mass function in which $\sigma_e$-dependent variations in the IMF across the population have been accounted for. Using a Fundamental Plane based observational proxy for $\sigma_e$ produces comparable results. By demonstrating that cheaper proxies are sufficiently accurate, our analysis should enable a more reliable census of the mass in stars for large galaxy samples, at a fraction of the cost. Our results are provided in tabular form.

We show how a characteristic length scale imprinted in the galaxy two-point correlation function, dubbed the "linear point", can serve as a comoving cosmological standard ruler. In contrast to the Baryon Acoustic Oscillation peak location, this scale is constant in redshift and is unaffected by non-linear effects to within $0.5$ percent precision. We measure the location of the linear point in the galaxy correlation function of the LOWZ and CMASS samples from the Twelfth Data Release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS) collaboration. We combine our linear-point measurement with cosmic-microwave-background constraints from the Planck satellite to estimate the isotropic-volume distance $D_{V}(z)$, without relying on a model-template or reconstruction method. We find $D_V(0.32)=1264\pm 28$ Mpc and $D_V(0.57)=2056\pm 22$ Mpc respectively, consistent with the quoted values from the BOSS collaboration. This remarkable result suggests that all the distance information contained in the baryon acoustic oscillations can be conveniently compressed into the single length associated with the linear point.

Due to late time non-linearities, the location of the acoustic peak in the two-point galaxy correlation function is a redshift-dependent quantity, thus it cannot be simply employed as a cosmological standard ruler. This has motivated the recent proposal of a novel ruler, also located in the Baryon Acoustic Oscillation range of scales of the correlation function, dubbed the "linear point". Unlike the peak, it is insensitive at the $0.5\%$ level to many of the non-linear effects that distort the clustering correlation function and shift the peak. However, this is not enough to make the linear point a useful standard ruler. In addition, we require a model-independent method to estimate its value from real data, avoiding the need to deploy a poorly known non-linear model of the correlation function. In this manuscript, we precisely validate a procedure for model-independent estimation of the linear point. We also identify the optimal set-up to estimate the linear point from the correlation function using galaxy catalogs. The methodology developed here is of general validity, and can be applied to any galaxy correlation-function data. As a working example, we apply this procedure to the LOWZ and CMASS galaxy samples of the Twelfth Data Release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS), for which the estimates of cosmic distances using the linear point have been presented in Anselmi et al. (2017) [1].

The window function for protohalos in Lagrangian space is often assumed to be a tophat in real space. We measure this profile directly and find that it is more extended than a tophat but less extended than a Gaussian; its shape is well-described by rounding the edges of the tophat by convolution with a Gaussian that has a scale length about 5 times smaller. This effective window $W_{\rm eff}$ is particularly simple in Fourier space, and has an analytic form in real space. Together with the excursion set bias parameters, $W_{\rm eff}$ describes the scale-dependence of the Lagrangian halo-matter cross correlation up to $kR_{\rm Lag} \sim 10 $, where $R_{\rm Lag}$ is the Lagrangian size of the protohalo. Moreover, with this $W_{\rm eff}$, all the spectral moments of the power spectrum are finite, allowing a straightforward estimate of the excursion set peak mass function. This estimate requires a prescription of the critical overdensity enclosed within a protohalo if it is to collapse, which we calibrate from simulations. We find that the resulting estimate of halo abundances is only accurate to about 20%, and we discuss why: A tophat in `infall time' towards the protohalo center need not correspond to a tophat in the initial spatial distribution, so models in which infall rather than smoothed overdensity is the relevant variable may be more accurate.

The simplest analyses of halo bias assume that halo mass alone determines halo clustering. However, if the large scale environment is fixed, then halo clustering is almost entirely determined by environment, and is almost completely independent of halo mass. We show why. Our analysis is useful for studies which use the environmental dependence of clustering to constrain cosmological and galaxy formation models. It also shows why many correlations between galaxy properties and environment are merely consequences of the underlying correlations between halos and their environments, and provides a framework for quantifying such inherited correlations.

Recent work has confirmed that the masses of supermassive black holes, estimated from scaling relations with global properties such as the stellar masses of their host galaxies, may be biased high. Much of this may be caused by the requirement that the gravitational sphere of influence of the black hole must be resolved for the black-hole mass to be reliably estimated. We revisit this issue by using a comprehensive galaxy evolution semi-analytic model, which self-consistently evolves supermassive black holes from high-redshift seeds via gas accretion and mergers, and also includes AGN feedback. Once tuned to reproduce the (mean) correlation of black-hole mass with velocity dispersion, the model is unable to also account for the correlation with stellar mass. This behaviour is independent of the model's parameters, thus suggesting an internal inconsistency in the data. The predicted distributions, especially at the low-mass end, are also much broader than observed. However, if selection effects are included, the model's predictions tend to align with the observations. We also demonstrate that the correlations between the residuals of the local scaling relations are more effective than the scaling relations themselves at constraining AGN feedback models. In fact, we find that our semi-analytic model, while in apparent broad agreement with the scaling relations when accounting for selection biases, yields very weak correlations between their residuals at fixed stellar mass, in stark contrast with observations. This problem persists when changing the AGN feedback strength, and is also present in the $z\sim 0$ outputs of the hydrodynamic cosmological simulation Horizon-AGN, which includes state-of-the-art treatments of AGN feedback. This suggests that current AGN feedback models may be too weak or are simply not capturing the effect of the black hole on the stellar velocity dispersion.

The protohalo patches from which halos form are defined by a number of constraints imposed on the Lagrangian dark matter density field. Each of these constraints contributes to biasing the spatial distribution of the protohalos relative to the matter. We show how measurements of this spatial distribution -- linear combinations of protohalo bias factors -- can be used to make inferences about the physics of halo formation. Our analysis exploits the fact that halo bias factors satisfy consistency relations which encode this physics, and that these relations are the same even for sub-populations in which assembly bias has played a role. We illustrate our methods using a model in which three parameters matter: a density threshold, the local slope and the curvature of the smoothed density field. The latter two are nearly degenerate; our approach naturally allows one to build an accurate effective two-parameter model for which the consistency relations still apply. This, with an accurate description of the smoothing window, allows one to describe the protohalo-matter cross-correlation very well, both in Fourier and configuration space. We then use our determination of the large scale bias parameters together with the consistency relations, to estimate the enclosed density and mean slope on the Lagrangian radius scale of the protohalos. Direct measurements of these quantities, made on smaller scales than those on which the bias parameters are typically measured, are in good agreement.

The Sloan Digital Sky Survey pipeline photometry underestimates the brightnesses of the most luminous galaxies. This is mainly because (i) the SDSS overestimates the sky background and (ii) single or two-component Sersic-based models better fit the surface brightness profile of galaxies, especially at high luminosities, than does the de Vaucouleurs model used by the SDSS pipeline. We use the PyMorph photometric reductions to isolate effect (ii) and show that it is the same in the full sample as in small group environments, and for satellites in the most massive clusters as well. None of these are expected to be significantly affected by intracluster light (ICL). We only see an additional effect for centrals in the most massive halos, but we argue that even this is not dominated by ICL. Hence, for the vast majority of galaxies, the differences between PyMorph and SDSS pipeline photometry cannot be ascribed to the semantics of whether or not one includes the ICL when describing the stellar mass of massive galaxies. Rather, they likely reflect differences in star formation or assembly histories. Failure to account for the SDSS underestimate has significantly biased most previous estimates of the SDSS luminosity and stellar mass functions, and therefore Halo Model estimates of the z ~ 0.1 relation between the mass of a halo and that of the galaxy at its center. We also show that when one studies correlations, at fixed group mass, with a quantity which was not used to define the groups, then selection effects appear. We show why such effects arise, and should not be mistaken for physical effects.

We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and S\'{e}rsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a S\'{e}rsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations.

We propose a new way of looking at the Baryon Acoustic Oscillations in the Large Scale Structure clustering correlation function. We identify a scale s_LP that has two fundamental features: its position is insensitive to non-linear gravity, redshift space distortions, and scale-dependent bias at the 0.5% level; it is geometrical, i.e. independent of the power spectrum of the primordial density fluctuation parameters. These two properties together make s_LP, called the "linear point", an excellent cosmological standard ruler. The linear point is also appealing because it is easily identified irrespectively of how non-linearities distort the correlation function. Finally, the correlation function amplitude at s_LP is similarly insensitive to non-linear corrections to within a few percent. Hence, exploiting the particular Baryon features in the correlation function, we propose three new estimators for growth measurements. A preliminary analysis of s_LP in current data is encouraging.

We quantify the systematic effects on the stellar mass function which arise from assumptions about the stellar population, as well as how one fits the light profiles of the most luminous galaxies at z ~ 0.1. When comparing results from the literature, we are careful to separate out these effects. Our analysis shows that while systematics in the estimated comoving number density which arise from different treatments of the stellar population remain of order < 0.5 dex, systematics in photometry are now about 0.1 dex, despite recent claims in the literature. Compared to these more recent analyses, previous work based on Sloan Digital Sky Survey (SDSS) pipeline photometry leads to underestimates of rho_*(> M_*) by factors of 3-10 in the mass range 10^11 - 10^11.6 M_Sun, but up to a factor of 100 at higher stellar masses. This impacts studies which match massive galaxies to dark matter halos. Although systematics which arise from different treatments of the stellar population remain of order < 0.5 dex, our finding that systematics in photometry now amount to only about 0.1 dex in the stellar mass density is a significant improvement with respect to a decade ago. Our results highlight the importance of using the same stellar population and photometric models whenever low and high redshift samples are compared.

Mass around dark matter halos can be divided into "infalling" material and "collapsed" material that has passed through at least one pericenter. Analytical models and simulations predict a rapid drop in the halo density profile associated with the transition between these two regimes. Using data from SDSS, we explore the evidence for such a feature in the density profiles of galaxy clusters and investigate the connection between this feature and a possible phase space boundary. We first estimate the steepening of the outer galaxy density profile around clusters: the profiles show an abrupt steepening, providing evidence for truncation of the halo profile. Next, we measure the galaxy density profile around clusters using two sets of galaxies selected based on color. We find evidence of an abrupt change in the galaxy colors that coincides with the location of the steepening of the density profile. Since galaxies are likely to be quenched of star formation and turn red inside of clusters, this change in the galaxy color distribution can be interpreted as the transition from an infalling regime to a collapsed regime. We also measure this transition using a model comparison approach which has been used recently in studies of the "splashback" phenomenon, but find that this approach is not a robust way to quantify the significance of detecting a splashback-like feature. Finally, we perform measurements using an independent cluster catalog to test for potential systematic errors associated with cluster selection. We identify several avenues for future work: improved understanding of the small-scale galaxy profile, lensing measurements, identification of proxies for the halo accretion rate, and other tests. With upcoming data from the DES, KiDS and HSC surveys, we can expect significant improvements in the study of halo boundaries.

We extend the comparison between the set of local galaxies having dynamically measured black holes with galaxies in the Sloan Digital Sky Survey (SDSS). We first show that the most up-to-date local black hole samples of early-type galaxies with measurements of effective radii, luminosities, and S\'ersic indices of the bulges of their host galaxies, have dynamical mass and S\'ersic index distributions consistent with those of SDSS early-type galaxies of similar bulge stellar mass. The host galaxies of local black hole samples thus do not appear structurally different from SDSS galaxies, sharing similar dynamical masses, light profiles and light distributions. Analysis of the residuals reveals that velocity dispersion is more fundamental than S\'ersic index n in the scaling relations between black holes and galaxies. Indeed, residuals with S\'ersic index could be ascribed to the (weak) correlation with bulge mass or even velocity dispersion. Finally, targetted Monte Carlo simulations that include the effects of the sphere of influence of the black hole, and tuned to reproduce the observed residuals and scaling relations in terms of velocity dispersion and stellar mass, show that, at least for galaxies with Mbulge > 1e10 Msun and n>5, the observed mean black hole mass at fixed S\'ersic index is biased significantly higher than the intrinsic value.

We compare the set of local galaxies having dynamically measured black holes with a large, unbiased sample of galaxies extracted from the Sloan Digital Sky Survey. We confirm earlier work showing that the majority of black hole hosts have significantly higher velocity dispersions sigma than local galaxies of similar stellar mass. We use Monte-Carlo simulations to illustrate the effect on black hole scaling relations if this bias arises from the requirement that the black hole sphere of influence must be resolved to measure black hole masses with spatially resolved kinematics. We find that this selection effect artificially increases the normalization of the Mbh-sigma relation by a factor of at least ~3; the bias for the Mbh-Mstar relation is even larger. Our Monte Carlo simulations and analysis of the residuals from scaling relations both indicate that sigma is more fundamental than Mstar or effective radius. In particular, the Mbh-Mstar relation is mostly a consequence of the Mbh-sigma and sigma-Mstar relations, and is heavily biased by up to a factor of 50 at small masses. This helps resolve the discrepancy between dynamically-based black hole-galaxy scaling relations versus those of active galaxies. Our simulations also disfavour broad distributions of black hole masses at fixed sigma. Correcting for this bias suggests that the calibration factor used to estimate black hole masses in active galaxies should be reduced to values of fvir~1. Black hole mass densities should also be proportionally smaller, perhaps implying significantly higher radiative efficiencies/black hole spins. Reducing black hole masses also reduces the gravitational wave signal expected from black hole mergers.

Cross-correlations between biased tracers and the dark matter field encode information about the physical variables which characterize these tracers. However, if the physical variables of interest are correlated with one another, then extracting this information is not as straightforward as one might naively have thought. We show how to exploit these correlations so as to estimate scale-independent bias factors of all orders in a model-independent way. We also show that failure to account for this will lead to incorrect conclusions about which variables matter and which do not. Morever, accounting for this allows one to use the scale dependence of bias to constrain the physics of halo formation; to date the argument has been phrased the other way around. We illustrate by showing that the scale dependence of linear and nonlinear bias, measured on nonlinear scales, can be used to provide consistent estimates of how the critical density for halo formation depends on halo mass. Our methods work even when the bias is nonlocal and stochastic, such as when, in addition to the spherically averaged density field and its derivatives, the quadrupolar shear field also matters for halo formation. In such models, the nonlocal bias factors are closely related to the more familiar local nonlinear bias factors, which are much easier to measure. Our analysis emphasizes the fact that biased tracers are biased because they do not sample fields (density, velocity, shear, etc.) at all positions in space in the same way that the dark matter does.

Recent analytical work on the modelling of dark halo abundances and clustering has demonstrated the advantages of combining the excursion set approach with peaks theory. We extend these ideas and introduce a model of excursion set peaks that incorporates the role of initial tidal effects or shear in determining the gravitational collapse of dark haloes. The model -- in which the critical density threshold for collapse depends on the tidal influences acting on protohaloes -- is well motivated from ellipsoidal collapse arguments and is also simple enough to be analytically tractable. We show that the predictions of this model are in very good agreement with measurements of the halo mass function and traditional scale dependent halo bias in N-body simulations across a wide range of masses and redshift. The presence of shear in the collapse threshold means that halo bias is naturally predicted to be nonlocal, and that protohalo densities at fixed mass are naturally predicted to have Lognormal-like distributions. We present the first direct estimate of Lagrangian nonlocal bias in N-body simulations, finding broad agreement with the model prediction. Finally, the simplicity of the model (which has essentially a single free parameter) opens the door to building efficient and accurate non-universal fitting functions of halo abundances and bias for use in precision cosmology.

We use the Halo Model to explore the implications of assuming that galaxy luminosities in groups are randomly drawn from an underlying luminosity function. We show that even the simplest of such order statistics models -- one in which this luminosity function $p(L)$ is universal -- naturally produces a number of features associated with previous analyses based on the `central plus Poisson satellites' hypothesis. These include the monotonic relation of mean central luminosity with halo mass, the Lognormal distribution around this mean, and the tight relation between the central and satellite mass scales. In stark contrast to observations of galaxy clustering, however, this model predicts $\textit{no}$ luminosity dependence of large scale clustering. We then show that an extended version of this model, based on the order statistics of a $\textit{halo mass dependent}$ luminosity function $p(L|m)$, is in much better agreement with the clustering data as well as satellite luminosities, but systematically under-predicts central luminosities. This brings into focus the idea that central galaxies constitute a distinct population that is affected by different physical processes than are the satellites. We model this physical difference as a statistical brightening of the central luminosities, over and above the order statistics prediction. The magnitude gap between the brightest and second brightest group galaxy is predicted as a by-product, and is also in good agreement with observations. We propose that this order statistics framework provides a useful language in which to compare the Halo Model for galaxies with more physically motivated galaxy formation models

Supermassive black hole -- host galaxy relations are key to the computation of the expected gravitational wave background (GWB) in the pulsar timing array (PTA) frequency band. It has been recently pointed out that standard relations adopted in GWB computations are in fact biased-high. We show that when this selection bias is taken into account, the expected GWB in the PTA band is a factor of about three smaller than previously estimated. Compared to other scaling relations recently published in the literature, the median amplitude of the signal at $f=1$yr$^{-1}$ drops from $1.3\times10^{-15}$ to $4\times10^{-16}$. Although this solves any potential tension between theoretical predictions and recent PTA limits without invoking other dynamical effects (such as stalling, eccentricity or strong coupling with the galactic environment), it also makes the GWB detection more challenging.

The abundance of galaxy clusters can constrain both the geometry and growth of structure in our Universe. However, this probe could be significantly complicated by recent claims of nonuniversality -- non-trivial dependences with respect to the cosmological model and redshift. In this work we analyse the dependance of the mass function on the way haloes are identified and establish if this can cause departures from universality. In order to explore this dependance, we use a set of different N-body cosmological simulations (Le SBARBINE simulations), with the latest cosmological parameters from the Planck collaboration; this first suite of simulations is followed by a lower resolution set, carried out with different cosmological parameters. We identify dark matter haloes using a Spherical Overdensity algorithm with varying overdensity thresholds (virial, 2000rho_c, 1000rho_c, 500rho_c, 200rho_c and 200rho_b) at all redshifts. We notice that, when expressed in term of the rescaled variable nu, the mass functionfor virial haloes is a nearly universal as a function of redshift and cosmology, while this is clearly not the case for the other overdensities we considered. We provide fitting functions for the halo mass function parameters as a function of overdensity, that allow to predict, to within a few percent accuracy, the halo mass function for a wide range of halo definitions, redshifts and cosmological models. We then show how the departures from universality associated with other halo definitions can be derived by combining the universality of the virial definition with the expected shape of the density profile of halos.

We describe the luminosity function, based on Sersic fits to the light profiles, of CMASS galaxies at z ~ 0.55. Compared to previous estimates, our Sersic-based reductions imply more luminous, massive galaxies, consistent with the effects of Sersic- rather than Petrosian or de Vaucouleur-based photometry on the Sloan Digital Sky Survey (SDSS) main galaxy sample at z ~ 0.1. This implies a significant revision of the high mass end of the correlation between stellar and halo mass. Inferences about the evolution of the luminosity and stellar mass functions depend strongly on the assumed, and uncertain, k+e corrections. In turn, these depend on the assumed age of the population. Applying k+e corrections taken from fitting the models of Maraston et al. (2009) to the colors of both SDSS and CMASS galaxies, the evolution of the luminosity and stellar mass functions appears impressively passive, provided that the fits are required to return old ages. However, when matched in comoving number- or luminosity-density, the SDSS galaxies are less strongly clustered compared to their counterparts in CMASS. This rules out the passive evolution scenario, and, indeed, any minor merger scenarios which preserve the rank ordering in stellar mass of the population. Potential incompletenesses in the CMASS sample would further enhance this mismatch. Our analysis highlights the virtue of combining clustering measurements with number counts.

The mass and structural evolution of massive galaxies is one of the hottest topics in galaxy formation. This is because it may reveal invaluable insights into the still debated evolutionary processes governing the growth and assembly of spheroids. However, direct comparison between models and observations is usually prevented by the so-called "progenitor bias", i.e., new galaxies entering the observational selection at later epochs, thus eluding a precise study of how pre-existing galaxies actually evolve in size. To limit this effect, we here gather data on high-redshift brightest group and cluster galaxies, evolve their (mean) host halo masses down to z=0 along their main progenitors, and assign as their "descendants" local SDSS central galaxies matched in host halo mass. At face value, the comparison between high redshift and local data suggests a noticeable increase in stellar mass of a factor of >2 since z~1, and of >2.5 in mean effective radius. We then compare the inferred stellar mass and size growth with those predicted by hierarchical models for central galaxies, selected at high redshifts to closely match the halo and stellar mass bins as in the data. Only hierarchical models characterized by very limited satellite stellar stripping and parabolic orbits are capable of broadly reproducing the stellar mass and size increase of a factor ~2-4 observed in cluster galaxies since z ~1. The predicted, average (major) merger rate since z~1 is in good agreement with the latest observational estimates.

The excursion set approach provides a framework for predicting how the abundance of dark matter halos depends on the initial conditions. A key ingredient of this formalism comes from the physics of halo formation: the specification of a critical overdensity threshold (barrier) which protohalos must exceed if they are to form bound virialized halos at a later time. Another ingredient is statistical, as it requires the specification of the appropriate statistical ensemble over which to average when making predictions. The excursion set approach explicitly averages over all initial positions, thus implicitly assuming that the appropriate ensemble is that associated with randomly chosen positions in space, rather than special positions such as peaks of the initial density field. Since halos are known to collapse around special positions, it is not clear that the physical and statistical assumptions which underlie the excursion set approach are self-consistent. We argue that they are at least for low mass halos, and illustrate by comparing our excursion set predictions with numerical data from the DEUS simulations.

The stellar mass-halo mass relation is a key constraint in all semi-analytic, numerical, and semi-empirical models of galaxy formation and evolution. However, its exact shape and redshift dependence remain debated. Several recent works support a relation in the local Universe steeper than previously thought. Based on the comparisons with a variety of data on massive central galaxies, we show that this steepening holds up to z~1, for stellar masses Mstar>2e11 Msun. Specifically, we find significant evidence for a high-mass end slope of \beta>0.35-0.70, instead of the usual \beta~0.20-0.30 reported by a number of previous results. When including the independent constraints from the recent BOSS clustering measurements, the data, independent of any systematic errors in stellar masses, tend to favor a model with a very small scatter (< 0.15 dex) in stellar mass at fixed halo mass, in the redshift range z < 0.8 and for Mstar>3e11 Msun, suggesting a close connection between massive galaxies and host halos even at relatively recent epochs. We discuss the implications of our results with respect to the evolution of the most massive galaxies since z~1.

The excursion set approach uses the statistics of the density field smoothed on a wide range of scales, to gain insight into a number of interesting processes in nonlinear structure formation, such as cluster assembly, merging and clustering. The approach treats the curve defined by the height of the overdensity fluctuation field when changing the smoothing scale as a random walk. The steps of the walks are often assumed to be uncorrelated, so that the walk heights are a Markov process, even though this assumption is known to be inaccurate for physically relevant filters. We develop a model in which the walk steps, rather than heights, are a Markov process, and correlations between steps arise because of nearest neighbour interactions. This model is a particular case of a general class, which we call Markov Velocity models. We show how these can approximate the walks generated by arbitrary power spectra and filters, and, unlike walks with Markov heights, provide a very good approximation to physically relevant models. We define a Markov Velocity Monte Carlo algorithm to generate walks whose first crossing distribution is very similar to that of TopHat-smoothed LCDM walks. Finally, we demonstrate that Markov Velocity walks generically exhibit a simple but realistic form of assembly bias, so we expect them to be useful in the construction of more realistic merger history trees.

The simplest stochastic halo formation models assume that the traceless part of the shear field acts to increase the initial overdensity (or decrease the underdensity) that a protohalo (or protovoid) must have if it is to form by the present time. Equivalently, it is the difference between the overdensity and (the square root of the) shear that must be larger than a threshold value. To estimate the effect this has on halo abundances using the excursion set approach, we must solve for the first crossing distribution of a barrier of constant height by the random walks associated with the difference, which is now (even for Gaussian initial conditions) a non-Gaussian variate. The correlation properties of such non-Gaussian walks are inherited from those of the density and the shear, and, since they are independent processes, the solution is in fact remarkably simple. We show that this provides an easy way to understand why earlier heuristic arguments about the nature of the solution worked so well. In addition to modelling halos and voids, this potentially simplifies models of the abundance and spatial distribution of filaments and sheets in the cosmic web.

We quantify the systematics in the size-luminosity relation of galaxies in the SDSS main sample which arise from fitting different 1- and 2-component model profiles to the images. In objects brighter than L*, fitting a single Sersic profile to what is really a two-component SerExp system leads to biases: the half-light radius is increasingly overestimated as n of the fitted single component increases; it is also overestimated at B/T ~ 0.6. However, the net effect on the R-L relation is small, except for the most luminous tail, where it curves upwards towards larger sizes. We also study how this relation depends on morphological type. Our analysis is one of the first to use Bayesian-classifier derived weights, rather than hard cuts, to define morphology. Crudely, there appear to be only two relations: one for early-types (Es, S0s and Sa's) and another for late-types (Sbs and Scds). However, closer inspection shows that within the early-type sample S0s tend to be 15% smaller than Es of the same luminosity, and, among faint late types, Sbs are more than 25% smaller than Scds. Neither the early- nor the late-type relations are pure power-laws: both show significant curvature, which we quantify. However, the R-L relations of the bulges of early-types are almost pure power laws. Our analysis confirms that two mass scales are special for early-type galaxies: M* = 3e10 and 2e11 Msun. These same mass scales are also special for late types: there is almost no correlation between R and M* below the former, and almost no late-types above the latter. We also show that the intrinsic scatter around the relation decreases at large luminosity and/or stellar mass; this should provide additional constraints on models of how the most massive galaxies formed.

Insight into a number of interesting questions in cosmology can be obtained from the first crossing distributions of physically motivated barriers by random walks with correlated steps. We write the first crossing distribution as a formal series, ordered by the number of times a walk upcrosses the barrier. Since the fraction of walks with many upcrossings is negligible if the walk has not taken many steps, the leading order term in this series is the most relevant for understanding the massive objects of most interest in cosmology. This first term only requires knowledge of the bivariate distribution of the walk height and slope, and provides an excellent approximation to the first crossing distribution for all barriers and smoothing filters of current interest. We show that this simplicity survives when extending the approach to the case of non-Gaussian random fields. For non-Gaussian fields which are obtained by deterministic transformations of a Gaussian, the first crossing distribution is simply related to that for Gaussian walks crossing a suitably rescaled barrier. Our analysis shows that this is a useful way to think of the generic case as well. Although our study is motivated by the possibility that the primordial fluctuation field was non-Gaussian, our results are general. In particular, they do not assume the non-Gaussianity is small, so they may be viewed as the solution to an excursion set analysis of the late-time, nonlinear fluctuation field rather than the initial one. They are also useful for models in which the barrier height is determined by quantities other than the initial density, since most other physically motivated variables (such as the shear) are usually stochastic and non-Gaussian. We use the Lognormal transformation to illustrate some of our arguments.

We compare state-of-the-art semi-analytic models of galaxy formation as well as advanced sub-halo abundance matching models with a large sample of early-type galaxies from SDSS at z < 0.3. We focus our attention on the dependence of median sizes of central galaxies on host halo mass. The data do not show any difference in the structural properties of early-type galaxies with environment, at fixed stellar mass. All hierarchical models considered in this work instead tend to predict a moderate to strong environmental dependence, with the median size increasing by a factor of about 1.5-3 when moving from low to high mass host haloes. At face value the discrepancy with the data is highly significant, especially at the cluster scale, for haloes above log Mhalo > 14. The convolution with (correlated) observational errors reduces some of the tension. Despite the observational uncertainties, the data tend to disfavour hierarchical models characterized by a relevant contribution of disc instabilities to the formation of spheroids, strong gas dissipation in (major) mergers, short dynamical friction timescales, and very short quenching timescales in infalling satellites. We also discuss a variety of additional related issues, such as the slope and scatter in the local size-stellar mass relation, the fraction of gas in local early-type galaxies, and the general predictions on satellite galaxies.

We use a large suite of N-body simulations to study departures from universality in halo abundances and clustering in cosmologies with non-vanishing neutrino masses. To this end, we study how the halo mass function and halo bias factors depend on the scaling variable $\sigma^2(M,z)$, the variance of the initial matter fluctuation field, rather than on halo mass $M$ and redshift $z$ themselves. We show that using the variance of the cold dark matter rather than the total mass field, i.e., $\sigma^2_{cdm}(M,z)$ rather than $\sigma^2_{m}(M,z)$, yields more universal results. Analysis of halo bias yields similar conclusions: When large-scale halo bias is defined with respect to the cold dark matter power spectrum, the result is both more universal, and less scale- or $k$-dependent. These results are used extensively in Papers I and III of this series.

We explore the scale dependence of halo bias using real space cross-correlation measurements in N-body simulations and in Pinocchio, an algorithm based on Lagrangian Perturbation Theory. Recent work has shown how to interpret such real space measurements in terms of k-dependent bias in Fourier space, and how to remove the k-dependence to reconstruct the k-independent peak-background split halo bias parameters. We compare our reconstruction of the linear bias, which requires no free parameters, with previous estimates from N-body simulations which were obtained directly in Fourier space at large scales, and find very good agreement. Our reconstruction of the quadratic bias is similarly parameter-free, although in this case there are no previous Fourier space measurements to compare with. Our analysis of N-body simulations explicitly tests the predictions of the excursion set peaks (ESP) formalism of Paranjape et al. (2013) for the scale dependence of bias; we find that the ESP predictions accurately describe our measurements. In addition, our measurements in Pinocchio serve as a useful, successful consistency check between Pinocchio and N-body simulations that is not accessible to traditional measurements.

In addition to the large systematic differences arising from assumptions about the stellar mass-to-light ratio, the massive end of the stellar mass function is rather sensitive to how one fits the light profiles of the most luminous galaxies. We quantify this by comparing the luminosity and stellar mass functions based on SDSS cmodel magnitudes, and PyMorph single-Sersic and Sersic-Exponential fits to the surface brightness profiles of galaxies in the SDSS. The PyMorph fits return more light, so that the predicted masses are larger than when cmodel magnitudes are used. As a result, the total stellar mass density at z~0.1 is about 1.2x larger than in our previous analysis of the SDSS. The differences are most pronounced at the massive end, where the measured number density of objects having M* > 6 x 10^{11} Msun is ~5x larger. Alternatively, at number densities of 10^{-6} Mpc^{-3}, the limiting stellar mass is 2x larger. The differences with respect to fits by other authors, typically based on Petrosian-like magnitudes, are even more dramatic, although some of these differences are due to sky-subtraction problems, and are sometimes masked by large differences in the assumed $M_*/L$ (even after scaling to the same IMF). Our results impact studies of the growth and assembly of stellar mass in galaxies, and of the relation between stellar and halo mass, so we provide simple analytic fits to these new luminosity and stellar mass functions and quantify how they depend on morphology, as well as the binned counts in electronic format.

We examine the two-point correlation function of local maxima in temperature fluctuations at the last scattering surface when this stochastic field is modified by the additional fluctuations produced by straight cosmic strings via the Kaiser-Stebbins effect. We demonstrate that one can detect the imprint of cosmic strings with tension $G\mu \gtrsim 1.2 \times 10^{-8}$ on noiseless $1^\prime$ resolution cosmic microwave background (CMB) maps at 95% confidence interval. Including the effects of foregrounds and anticipated systematic errors increases the lower bound to $G\mu \gtrsim 9.0\times 10^{-8}$ at $2\sigma$ confidence level. Smearing by beams of order 4' degrades the bound further to $G\mu \gtrsim 1.6 \times 10^{-7}$. Our results indicate that two-point statistics are more powerful than 1-point statistics (e.g. number counts) for identifying the non-Gaussianity in the CMB due to straight cosmic strings.

Recently, we provided a simple but accurate formula which closely approximates the first crossing distribution associated with random walks having correlated steps. The approximation is accurate for the wide range of barrier shapes of current interest and is based on the requirement that, in addition to having the right height, the walk must cross the barrier going upwards. Therefore, it only requires knowledge of the bivariate distribution of the walk height and slope, and is particularly useful for excursion set models of the massive end of the halo mass function. However, it diverges at lower masses. We show how to cure this divergence by using a formulation which requires knowledge of just one other variable. While our analysis is general, we use examples based on Gaussian initial conditions to illustrate our results. Our formulation, which is simple and fast, yields excellent agreement with the considerably more computationally expensive Monte-Carlo solution of the first crossing distribution, for a wide variety of moving barriers, even at very low masses.

We present a new parallel implementation of the PINpointing Orbit Crossing-Collapsed HIerarchical Objects (PINOCCHIO) algorithm, a quick tool, based on Lagrangian Perturbation Theory, for the hierarchical build-up of Dark Matter halos in cosmological volumes. To assess its ability to predict halo correlations on large scales, we compare its results with those of an N-body simulation of a 3 Gpc/h box sampled with 2048^3 particles taken from the MICE suite, matching the same seeds for the initial conditions. Thanks to the FFTW libraries and to the relatively simple design, the code shows very good scaling properties. The CPU time required by PINOCCHIO is a tiny fraction (~1/2000) of that required by the MICE simulation. Varying some of PINOCCHIO numerical parameters allows one to produce a universal mass function that lies in the range allowed by published fits, although it underestimates the MICE mass function of FoF halos in the high mass tail. We compare the matter-halo and the halo-halo power spectra with those of the MICE simulation and find that these 2-point statistics are well recovered on large scales. In particular, when catalogs are matched in number density, agreement within ten per cent is achieved for the halo power spectrum. At scales k>0.1 h/Mpc, the inaccuracy of the Zel'dovich approximation in locating halo positions causes an underestimate of the power spectrum that can be modeled as a Gaussian factor with a damping scale of d=3 Mpc/h at z=0, decreasing at higher redshift. Finally, a remarkable match is obtained for the reduced halo bispectrum, showing a good description of nonlinear halo bias. Our results demonstrate the potential of PINOCCHIO as an accurate and flexible tool for generating large ensembles of mock galaxy surveys, with interesting applications for the analysis of large galaxy redshift surveys.

Halos are biased tracers of the dark matter distribution. It is often assumed that the patches from which halos formed are locally biased with respect to the initial fluctuation field, meaning that the halo-patch fluctuation field can be written as a Taylor series in that of the dark matter. If quantities other than the local density influence halo formation, then this Lagrangian bias will generically be nonlocal; the Taylor series must be performed with respect to these other variables as well. We illustrate the effect with Monte-Carlo simulations of a model in which halo formation depends on the local shear (the quadrupole of perturbation theory), and provide an analytic model which provides a good description of our results. Our model, which extends the excursion set approach to walks in more than one dimension, works both when steps in the walk are uncorrelated, as well as when there are correlations between steps. For walks with correlated steps, our model includes two distinct types of nonlocality: one is due to the fact that the initial density profile around a patch which is destined to form a halo must fall sufficiently steeply around it -- this introduces k-dependence to even the linear bias factor, but otherwise only affects the monopole of the clustering signal. The other is due to the surrounding shear field; this affects the quadratic and higher order bias factors, and introduces an angular dependence to the clustering signal. In both cases, our analysis shows that these nonlocal Lagrangian bias terms can be significant, particularly for massive halos; they must be accounted for in analyses of higher order clustering such as the halo bispectrum in Lagrangian or Eulerian space. Although we illustrate these effects using halos, our analysis and conclusions also apply to the other constituents of the cosmic web -- filaments, sheets and voids.

We describe an algorithm for identifying ellipsoidal haloes in numerical simulations, and quantify how the resulting estimates of halo mass and shape differ with respect to spherical halo finders. Haloes become more prolate when fit with ellipsoids, the difference being most pronounced for the more aspherical objects. Although the ellipsoidal mass is systematically larger, this is less than 10% for most of the haloes. However, even this small difference in mass corresponds to a significant difference in shape. We quantify these effects also on the initial mass and deformation tensors, on which most models of triaxial collapse are based. By studying the properties of protohaloes in the initial conditions, we find that models in which protohaloes are identified in Lagrangian space by three positive eigenvalues of the deformation tensor are tenable only at the masses well-above $M_*$. The overdensity $\delta$ within almost any protohalo is larger than the critical value associated with spherical collapse (increasing as mass decreases); this is in good qualitative agreement with models which identify haloes requiring that collapse have occured along all three principal axes, each axis having turned around from the universal expansion at a different time. The distributions of initial values are in agreement with the simplest predictions associated with ellipsoidal collapse, assuming initially spherical protohaloes, collapsed around random positions which were sufficiently overdense. However, most protohaloes are not spherical and departures from sphericity increase as protohalo mass decreases. [Abridged]

It is often assumed that the halo-patch fluctuation field can be written as a Taylor series in the initial Lagrangian dark matter density fluctuation field. We show that if this Lagrangian bias is local, and the initial conditions are Gaussian, then the two-point cross-correlation between halos and mass should be linearly proportional to the mass-mass auto-correlation function. This statement is exact and valid on all scales; there are no higher order contributions, e.g., from terms proportional to products or convolutions of two-point functions, which one might have thought would appear upon truncating the Taylor series of the halo bias function. In addition, the auto-correlation function of locally biased tracers can be written as a Taylor series in the auto-correlation function of the mass; there are no terms involving, e.g., derivatives or convolutions. Moreover, although the leading order coefficient, the linear bias factor of the auto-correlation function is just the square of that for the cross-correlation, it is the same as that obtained from expanding the mean number of halos as a function of the local density only in the large-scale limit. In principle, these relations allow simple tests of whether or not halo bias is indeed local in Lagrangian space. We discuss why things are more complicated in practice. We also discuss our results in light of recent work on the renormalizability of halo bias, demonstrating that it is better to renormalize than not. We use the Lognormal model to illustrate many of our findings.

We describe how to extend the excursion set peaks framework so that its predictions of dark halo abundances and clustering can be compared directly with simulations. These extensions include: a halo mass definition which uses the TopHat filter in real space; the mean dependence of the critical density for collapse delta_c on halo mass m; and the scatter around this mean value. All three of these are motivated by the physics of triaxial rather than spherical collapse. A comparison of the resulting mass function with N-body results shows that, if one uses delta_c(m) and its scatter as determined from simulations, then all three are necessary ingredients for obtaining ~10% accuracy. E.g., assuming a constant value of delta_c with no scatter, as motivated by the physics of spherical collapse, leads to many more massive halos than seen in simulations. The same model is also in excellent agreement with N-body results for the linear halo bias, especially at the high mass end where the traditional peak-background split argument applied to the mass function fit is known to underpredict the measured bias by ~10%. In the excursion set language, our model is about walks centered on special positions (peaks) in the initial conditions -- we discuss what it implies for the usual calculation in which all walks contribute to the statistics.

We analyze environmental correlations using mark clustering statistics with the mock galaxy catalogue constructed by Muldrew et al. (Paper I). We find that mark correlation functions are able to detect even a small dependence of galaxy properties on the environment, quantified by the overdensity $1+\delta$, while such a small dependence would be difficult to detect by traditional methods. We then show that rank ordering the marks and using the rank as a weight is a simple way of comparing the correlation signals for different marks. With this we quantify to what extent fixed-aperture overdensities are sensitive to large-scale halo environments, nearest-neighbor overdensities are sensitive to small-scale environments within haloes, and colour is a better tracer of overdensity than is luminosity.

We describe a simple fully analytic model of the excursion set approach associated with two Gaussian random walks: the first walk represents the initial overdensity around a protohalo, and the second is a crude way of allowing for other factors which might influence halo formation. This model is richer than that based on a single walk, because it yields a distribution of heights at first crossing. We provide explicit expressions for the unconditional first crossing distribution which is usually used to model the halo mass function, the progenitor distributions, and the conditional distributions from which correlations with environment are usually estimated. These latter exhibit perhaps the simplest form of what is often called nonlocal bias, and which we prefer to call stochastic bias, since the new bias effects arise from `hidden-variables' other than density, but these may still be defined locally. We provide explicit expressions for these new bias factors. We also provide formulae for the distribution of heights at first crossing in the unconditional and conditional cases. In contrast to the first crossing distribution, these are exact, even for moving barriers, and for walks with correlated steps. The conditional distributions yield predictions for the distribution of halo concentrations at fixed mass and formation redshift. They also exhibit assembly bias like effects, even when the steps in the walks themselves are uncorrelated. Finally, we show how the predictions are modified if we add the requirement that halos form around peaks: these depend on whether the peaks constraint is applied to a combination of the overdensity and the other variable, or to the overdensity alone. Our results demonstrate the power of requiring models to reproduce not just halo counts but the distribution of overdensities at fixed protohalo mass as well.

Motivated by the recent detection of an enhanced clustering signal along the major axis of haloes in N-body simulations, we derive a formula for the anisotropic density distribution around haloes and voids on large scales. Our model, which assumes linear theory and that the formation and orientation of nonlinear structures are strongly correlated with the Lagrangian shear, is in good agreement with measurements. We also show that the measured amplitude is inconsistent with a model in which the alignment is produced by the initial inertia rather than shear tensor.

We describe a model of dark matter halo abundances and clustering which combines the two most widely used approaches to this problem: that based on peaks and the other based on excursion sets. Our approach can be thought of as addressing the cloud-in-cloud problem for peaks and/or modifying the excursion set approach so that it averages over a special subset, rather than all possible walks. In this respect, it seeks to account for correlations between steps in the walk as well as correlations between walks. We first show how the excursion set and peaks models can be written in the same formalism, and then use this correspondence to write our combined excursion set peaks model. We then give simple expressions for the mass function and bias, showing that even the linear halo bias factor is predicted to be k-dependent as a consequence of the nonlocality associated with the peak constraint. At large masses, our model has little or no need to rescale the variable delta_c from the value associated with spherical collapse, and suggests a simple explanation for why the linear halo bias factor appears to lie above that based on the peak-background split at high masses when such a rescaling is assumed. Although we have concentrated on peaks, our analysis is more generally applicable to other traditionally single-scale analyses of large-scale structure.

If one accounts for correlations between scales, then nonlocal, k-dependent halo bias is part and parcel of the excursion set approach, and hence of halo model predictions for galaxy bias. We present an analysis that distinguishes between a number of different effects, each one of which contributes to scale-dependent bias in real space. We show how to isolate these effects and remove the scale dependence, order by order, by cross-correlating the halo field with suitably transformed versions of the mass field. These transformations may be thought of as simple one-point, two-scale measurements that allow one to estimate quantities which are usually constrained using n-point statistics. As part of our analysis, we present a simple analytic approximation for the first crossing distribution of walks with correlated steps which are constrained to pass through a specified point, and demonstrate its accuracy. Although we concentrate on nonlinear, nonlocal bias with respect to a Gaussian random field, we show how to generalize our analysis to more general fields.

We present a fast method of producing mock galaxy catalogues that can be used to compute covariance matrices of large-scale clustering measurements and test the methods of analysis. Our method populates a 2nd-order Lagrangian Perturbation Theory (2LPT) matter field, where we calibrate masses of dark matter halos by detailed comparisons with N-body simulations. We demonstrate the clustering of halos is recovered at ~10 per cent accuracy. We populate halos with mock galaxies using a Halo Occupation Distribution (HOD) prescription, which has been calibrated to reproduce the clustering measurements on scales between 30 and 80 Mpc/h. We compare the sample covariance matrix from our mocks with analytic estimates, and discuss differences. We have used this method to make catalogues corresponding to Data Release 9 of the Baryon Oscillation Spectroscopic Survey (BOSS),producing 600 mock catalogues of the "CMASS" galaxy sample. These mocks enabled detailed tests of methods and errors that formed an integral part of companion analyses of these galaxy data.

We show that the distribution of luminosities of Brightest Cluster Galaxies in an SDSS-based group catalog suggests that BCG luminosities are just the statistical extremes of the group galaxy luminosity function. This latter happens to be very well approximated by the all-galaxy luminosity function (restricted to Mr<-19.9), provided one uses a parametrization of this function that is accurate at the bright end. A similar analysis of the luminosity distribution of the Brightest Satellite Galaxies suggests that they are best thought of as being the second brightest pick from the same luminosity distribution of which BCGs are the brightest. I.e., BSGs are not the brightest of some universal satellite luminosity function, in contrast to what Halo Model analyses of the luminosity dependence of clustering suggest. However, we then use mark correlations to provide a novel test of these order statistics, showing that the hypothesis of a universal luminosity function (i.e. no halo mass dependence) from which the BCGs and BSGs are drawn is incompatible with the data, despite the fact that there was no hint of this in the BCG and BSG luminosity distributions themselves. We also discuss why, since extreme value statistics are explicitly a function of the number of draws, the consistency of BCG luminosities with extreme value statistics is most clearly seen if one is careful to perform the test at fixed group richness N. Tests at, e.g., fixed total group luminosity Ltot, will generally be biased and may lead to erroneous conclusions.

We apply CMB lensing techniques to large scale structure and solve for the 3-D cosmic tidal field. We use small scale filamentary structures to solve for the large scale tidal shear and gravitational potential. By comparing this to the redshift space density field, one can measure the gravitational growth factor on large scales without cosmic variance. This potentially enables accurate measurements of neutrino masses and reconstruction of radial modes lost in 21 cm intensity mapping, which are essential for CMB and other cross correlations. We relate the tidal fields to the squeezed limit bispectrum, and present initial results from simulations and data from the SDSS.

The relationship between galaxy and matter overdensities, bias, is most often assumed to be local. This is however unstable under time evolution, we provide proofs under several sets of assumptions. In the simplest model galaxies are created locally and linearly biased at a single time, and subsequently move with the matter (no velocity bias) conserving their comoving number density (no merging). We show that, after this formation time, the bias becomes unavoidably non-local and non-linear at large scales. We identify the non-local gravitationally induced fields in which the galaxy overdensity can be expanded, showing that they can be constructed out of the invariants of the deformation tensor (Galileons). In addition, we show that this result persists if we include an arbitrary evolution of the comoving number density of tracers. We then include velocity bias, and show that new contributions appear, a dipole field being the signature at second order. We test these predictions by studying the dependence of halo overdensities in cells of fixed matter density: measurements in simulations show that departures from the mean bias relation are strongly correlated with the non-local gravitationally induced fields identified by our formalism. The effects on non-local bias seen in the simulations are most important for the most biased halos, as expected from our predictions. The non-locality seen in the simulations is not fully captured by assuming local bias in Lagrangian space. Accounting for these effects when modeling galaxy bias is essential for correctly describing the dependence on triangle shape of the galaxy bispectrum, and hence constraining cosmological parameters and primordial non-Gaussianity. We show that using our formalism we remove an important systematic in the determination of bias parameters from the galaxy bispectrum, particularly for luminous galaxies. (abridged)

The coefficients a and b of the Fundamental Plane relation R ~ Sigma^a I^b depend on whether one minimizes the scatter in the R direction or orthogonal to the Plane. We provide explicit expressions for a and b (and confidence limits) in terms of the covariances between logR, logSigma and logI. Our analysis is more generally applicable to any other correlations between three variables: e.g., the color-magnitude-Sigma relation, the L-Sigma-Mbh relation, or the relation between the X-ray luminosity, Sunyaev-Zeldovich decrement and optical richness of a cluster, so we provide IDL code which implements these ideas, and we show how our analysis generalizes further to correlations between more than three variables. We show how to account for correlated errors and selection effects, and quantify the difference between the direct, inverse and orthogonal fit coefficients. We show that the three vectors associated with the Fundamental Plane can all be written as simple combinations of a and b because the distribution of I is much broader than that of Sigma, and Sigma and I are only weakly correlated. Why this should be so for galaxies is a fundamental open question about the physics of early-type galaxy formation. If luminosity evolution is differential, and Rs and Sigmas do not evolve, then this is just an accident: Sigma and I must have been correlated in the past. On the other hand, if the (lack of) correlation is similar to that at the present time, then differential luminosity evolution must have been accompanied by structural evolution. A model in which the luminosities of low-L galaxies evolve more rapidly than do those of higher-L galaxies is able to produce the observed decrease in a (by a factor of 2 at z~1) while having b decrease by only about 20 percent. In such a model, the Mdyn/L ratio is a steeper function of Mdyn at higher z.

We develop a simple model for estimating the mass growth histories of dark matter halos. The model is based on a fit to the formation time distribution, where formation is defined as the earliest time that the main branch of the merger tree contains a fraction f of the final mass M . Our analysis exploits the fact that the median formation time as a function of f is the same as the median of the main progenitor mass distribution as a function of time. When coupled with previous work showing that the concentration c of the final halo is related to the formation time tf associated with f \sim 0.04, our approach provides a simple algorithm for estimating how the distribution of halo concentrations may be expected to depend on mass, redshift and the expansion history of the background cosmology. We also show that one can predict log10 c with a precision of about 0.13 and 0.09 dex if only its mass, or both mass and tf are known. And, conversely, one can predict log10 tf from mass or c with a precision of 0.12 and 0.09 dex, approximately independent of f. Adding the mass to the c-based estimate does not result in further improvement. These latter results may be useful for studies which seek to compare the age of the stars in the central galaxy in a halo with the time the core was first assembled.

We provide a simple formula that accurately approximates the first crossing distribution of barriers having a wide variety of shapes, by random walks with a wide range of correlations between steps. Special cases of it are useful for estimating halo abundances, evolution, and bias, as well as the nonlinear counts in cells distribution. We discuss how it can be extended to allow for the dependence of the barrier on quantities other than overdensity, to construct an excursion set model for peaks, and to show why assembly and scale dependent bias are generic even at the linear level.

We extend earlier work on the problem of estimating the void-volume function -- the abundance and evolution of large voids which grow gravitationally in an expanding universe -- in two ways. The first removes an ambiguity about how the void-in-cloud process, which erases small voids, should be incorporated into the excursion set approach. The main technical change here is to think of voids within a fully Eulerian, rather than purely Lagrangian, framework. The second accounts for correlations between different spatial scales in the initial conditions. We provide numerical and analytical arguments showing how and why both changes modify the predicted abundances substantially. In particular, we show that the predicted importance of the void-in-cloud process depends strongly on whether or not one accounts for correlations between scales. With our new formulation, the void-in-cloud process dramatically reduces the predicted abundances of voids if such correlations are ignored, but only matters for the smallest voids in the more realistic case in which the spatial correlations are included.

The Excursion Set approach has been used to make predictions for a number of interesting quantities in studies of nonlinear hierarchical clustering. These include the halo mass function, halo merger rates, halo formation times and masses, halo clustering, analogous quantities for voids, and the distribution of dark matter counts in randomly placed cells. The approach assumes that all these quantities can be mapped to problems involving the first crossing distribution of a suitably chosen barrier by random walks. Most analytic expressions for these distributions ignore the fact that, although different k-modes in the initial Gaussian field are uncorrelated, this is not true in real space: the values of the density field at a given spatial position, when smoothed on different real-space scales, are correlated in a nontrivial way. As a result, the problem is to estimate first crossing distribution by random walks having correlated rather than uncorrelated steps. In 1990, Peacock & Heavens presented a simple approximation for the first crossing distribution of a single barrier of constant height by walks with correlated steps. We show that their approximation can be thought of as a correction to the distribution associated with what we call smooth completely correlated walks. We then use this insight to extend their approach to treat moving barriers, as well as walks that are constrained to pass through a certain point before crossing the barrier. For the latter, we show that a simple rescaling, inspired by bivariate Gaussian statistics, of the unconditional first crossing distribution, accurately describes the conditional distribution, independently of the choice of analytical prescription for the former. In all cases, comparison with Monte-Carlo solutions of the problem shows reasonably good agreement. (Abridged)

In the Excursion Set approach, halo abundances and clustering are closely related. This relation is exploited in many modern methods which seek to constrain cosmological parameters on the basis of the observed spatial distribution of clusters. However, to obtain analytic expressions for these quantities, most Excursion Set based predictions ignore the fact that, although different k-modes in the initial Gaussian field are uncorrelated, this is not true in real space: the values of the density field at a given spatial position, when smoothed on different real-space scales, are correlated in a nontrivial way. We show that when the excursion set approach is extended to include such correlations, then one must be careful to account for the fact that the associated prediction for halo bias is explicitly a real-space quantity. Therefore, care must be taken when comparing the predictions of this approach with measurements in simulations, which are typically made in Fourier-space. We show how to correct for this effect, and demonstrate that ignorance of this effect in recent analyses of halo bias has led to incorrect conclusions and biased constraints.

Using analytical models and cosmological N-body simulations, we study the free-free radio emission from ionized gas in clusters and groups of galaxies. The results obtained with the simulations are compared with analytical predictions based on the mass function and scaling relations. Earlier works based on analytical models have shown that the average free-free signal from small haloes (galaxies) during and after the reionization time could be detected with future experiments as a distortion of the CMB spectrum at low frequencies ($\nu <$ 5 GHz). We focus on the period after the reionization time (from redshift $z=0$ up to $z=7$) and on haloes that are more massive than in previous works (groups and clusters). We show how the average signal from haloes with $M > 10^{13} h^{-1} M_{\odot}$ is less than 10% the signal from the more abundant and colder smaller mass haloes. However, the individual signal from the massive haloes could be detected with future experiments opening the door for a new window to study the intracluster medium.

We use extreme value statistics to assess the significance of two of the most dramatic structures in the local Universe: the Shapley supercluster and the Sloan Great Wall. If we assume that Shapley (volume ~ 1.2 x 10^5 (Mpc/h)^3) evolved from an overdense region in the initial Gaussian fluctuation field, with currently popular choices for the background cosmological model and the shape and amplitude sigma8 of the initial power spectrum, we estimate that the total mass of the system is within 20 percent of 1.8 x 10^16 Msun/h. Extreme value statistics show that the existence of this massive concentration is not unexpected if the initial fluctuation field was Gaussian, provided there are no other similar objects within a sphere of radius 200 Mpc/h centred on our Galaxy. However, a similar analysis of the Sloan Great Wall, a more distant (z ~ 0.08) and extended concentration of structures (volume ~ 7.2 x 10^5 (Mpc/h)^3) suggests that it is more unusual. We estimate its total mass to be within 20 percent of 1.2 x 10^17 Msun/h; even if it is the densest such object of its volume within z=0.2, its existence is difficult to reconcile with Gaussian initial conditions if sigma8 < 0.9. This tension can be alleviated if this structure is the densest within the Hubble volume. Finally, we show how extreme value statistics can be used to address the likelihood that an object like Shapley exists in the same volume which contains the Great Wall, finding, again, that Shapley is not particularly unusual. It is straightforward to incorporate other models of the initial fluctuation field into our formalism.

We combine the physics of the ellipsoidal collapse model with the excursion set theory to study the shapes of dark matter halos. In particular, we develop an analytic approximation to the nonlinear evolution that is more accurate than the Zeldovich approximation; we introduce a planar representation of halo axis ratios, which allows a concise and intuitive description of the dynamics of collapsing regions and allows one to relate the final shape of a halo to its initial shape; we provide simple physical explanations for some empirical fitting formulae obtained from numerical studies. Comparison with simulations is challenging, as there is no agreement about how to define a non-spherical gravitationally bound object. Nevertheless, we find that our model matches the conditional minor-to-intermediate axis ratio distribution rather well, although it disagrees with the numerical results in reproducing the minor-to-major axis ratio distribution. In particular, the mass dependence of the minor-to-major axis distribution appears to be the opposite to what is found in many previous numerical studies, where low-mass halos are preferentially more spherical than high-mass halos. In our model, the high-mass halos are predicted to be more spherical, consistent with results based on a more recent and elaborate halo finding algorithm, and with observations of the mass dependence of the shapes of early-type galaxies. We suggest that some of the disagreement with some previous numerical studies may be alleviated if we consider only isolated halos.

For early-type galaxies, the correlations between stellar mass and size, velocity dispersion, surface brightness, color, axis ratio and color-gradient all indicate that two mass scales, M* = 3 x 10^10 Msun and M* = 2 x 10^11 Msun, are special. The smaller scale could mark the transition between wet and dry mergers, or it could be related to the interplay between SN and AGN feedback, although quantitative measures of this transition may be affected by morphological contamination. At the more massive scale, mean axis ratios and color gradients are maximal, and above it, the colors are redder, the sizes larger and the velocity dispersions smaller than expected based on the scaling at lower M*. In contrast, the color-sigma relation, and indeed, most scaling relations with sigma, are not curved: they are well-described by a single power law, or in some cases, are almost completely flat. When major dry mergers change masses, sizes, axis ratios and color gradients, they are expected to change the colors or velocity dispersions much less. Therefore, the fact that scaling relations at sigma > 150 km/s show no features, whereas the size-M*, b/a-M*, color-M* and color gradient-M* relations do, suggests that M* = 2 x 10^11 Msun is the scale above which major dry mergers dominate the assembly histories of early-type galaxies.

We derive a general expression for the large-scale halo bias, in theories with a scale-dependent linear growth, using the excursion set formalism. Such theories include modified gravity models, and models in which the dark energy clustering is non-negligible. A scale dependence is imprinted in both the formation and evolved biases by the scale-dependent growth. Mergers are accounted for in our derivation, which thus extends earlier work which focused on passive evolution. There is a simple analytic form for the bias for those theories in which the nonlinear collapse of perturbations is approximately the same as in general relativity. As an illustration, we apply our results to a simple Yukawa modification of gravity, and use SDSS measurements of the clustering of luminous red galaxies to constrain the theory's parameters.

We derive the weight function w(M) to apply to dark-matter halos that minimizes the stochasticity between the weighted halo distribution and its underlying mass density field. The optimal w(M) depends on the range of masses being used in the estimator. In N-body simulations, the Poisson estimator is up to 15 times noisier than the optimal. Implementation of the optimal weight yields significantly lower stochasticity than weighting halos by their mass, bias or equal. Optimal weighting could make cosmological tests based on the matter power spectrum or cross-correlations much more powerful and/or cost-effective. A volume-limited measurement of the mass power spectrum at k=0.2h/Mpc over the entire z<1 universe could ideally be done using only 6 million redshifts of halos with mass M>6\times10^{13}h^{-1}M_\odot (1\times10^{13}) at z=0 (z=1); this is 5 times fewer than the Poisson model predicts. Using halo occupancy distributions (HOD) we find that uniformly-weighted catalogs of luminous red galaxies require >3 times more redshifts than an optimally-weighted halo catalog to reconstruct the mass to the same accuracy. While the mean HODs of galaxies above a threshold luminosity are similar to the optimal w(M), the stochasticity of the halo occupation degrades the mass estimator. Blue or emission-line galaxies are about 100 times less efficient at reconstructing mass than an optimal weighting scheme. This suggests an efficient observational approach of identifying and weighting halos with a deep photo-z survey before conducting a spectroscopic survey. The optimal w(M) and mass-estimator stochasticity predicted by the standard halo model for M>10^{12}h^{-1}M_\odot are in reasonable agreement with our measurements, with the important exceptions that the halos must be assumed to be linearly biased samples of a "halo field" that is distinct from the mass field. (Abridged)

The color-magnitude relation of early-type galaxies differs slightly but significantly from a pure power-law, curving downwards at low and upwards at large luminosities (Mr>-20.5 and Mr<-22.5). This remains true of the color-size relation, and is even more apparent with stellar mass (M* < 3x10^10 Msun and M* > 2x10^11 Msun). The upwards curvature at the massive end does not appear to be due to stellar population effects. In contrast, the color-sigma relation is well-described by a single power law. Since major dry mergers change neither the colors nor sigma, but they do change masses and sizes, the clear features observed in the scaling relations with M*, but not with sigma > 150 km/s, suggest that M* > 2x10^11 Msun is the scale above which major dry mergers dominate the assembly history. We discuss three models of the merger histories since z ~ 1 which are compatible with our measurements. In all three models, dry mergers are responsible for the flattening of the color-M* relation at M* > 3x10^10 Msun - wet mergers only matter at smaller masses. At M* > 2 x 10^11 Msun, the merger histories in one model are dominated by major rather than minor dry mergers, as suggested by the axis ratio and color gradient trends. In another, although both major and minor mergers occur at the high mass end, the minor mergers contribute primarily to the formation of the ICL, rather than to the mass growth of the central massive galaxy. A final model assumes that the reddest objects were assembled by a mix of major and minor dry mergers.

Models of galaxy and halo clustering commonly assume that the tracers can be treated as a continuous field locally biased with respect to the underlying mass distribution. In the peak model pioneered by BBKS, one considers instead density maxima of the initial, Gaussian mass density field as an approximation to the formation site of virialized objects. In this paper, the peak model is extended in two ways to improve its predictive accuracy. Firstly, we derive the two-point correlation function of initial density peaks up to second order and demonstrate that a peak-background split approach can be applied to obtain the k-independent and k-dependent peak bias factors at all orders. Secondly, we explore the gravitational evolution of the peak correlation function within the Zel'dovich approximation. We show that the local (Lagrangian) bias approach emerges as a special case of the peak model, in which all bias parameters are scale-independent and there is no statistical velocity bias. We apply our formulae to study how the Lagrangian peak biasing, the diffusion due to large scale flows and the mode-coupling due to nonlocal interactions affect the scale dependence of bias from small separations up to the baryon acoustic oscillation (BAO) scale. For 2-sigma density peaks collapsing at z=0.3, our model predicts a ~ 5% residual scale-dependent bias around the acoustic scale that arises mostly from first-order Lagrangian peak biasing (as opposed to second-order gravity mode-coupling). We also search for a scale dependence of bias in the large scale auto-correlation of massive halos extracted from a very large N-body simulation provided by the MICE collaboration. For halos with mass M>10^{14}Msun/h, our measurements demonstrate a scale-dependent bias across the BAO feature which is very well reproduced by a prediction based on the peak model.

We use the ellipsoidal collapse approximation to investigate the nonlinear redshift space evolution of the density field with primordial non-Gaussianity of the local f_{nl}-type. We utilize the joint distribution of eigenvalues of the initial non-Gaussian shear field and evaluate the evolved redshift space probability distribution function (PDF). It is shown that, similar to the real space analysis, the underdense tail of the nonlinear redshift space PDF differs significantly from that for Gaussian initial conditions. We also derive the lowest order correction of the Kaiser's formulain the presence of a non-zero f_{nl}.

Accurate knowledge of the non-linear dark-matter power spectrum is important for understanding the large-scale structure of the Universe, the statistics of dark-matter haloes and their evolution, and cosmological gravitational lensing. We analytically model the dark-matter power spectrum and its cross-power spectrum with dark-matter haloes. Our model extends the halo-model formalism, including realistic substructure population within individual dark-matter haloes and the scatter of the concentration parameter at fixed halo mass. We consider three prescriptions for the mass-concentration relation and two for the substructure distribution in dark-matter haloes. We show that this extension of the halo model mainly increases the predicted power on the small scales, and is crucial for proper modeling the cosmological weak-lensing signal due to low-mass haloes. Our extended formalism shows how the halo model approach can be improved in accuracy as one increases the number of ingredients that are calibrated from n-body simulations.

We provide fits to the distribution of galaxy luminosity, size, velocity dispersion and stellar mass as a function of concentration index C_r and morphological type in the SDSS. We also quantify how estimates of the fraction of `early' or `late' type galaxies depend on whether the samples were cut in color, concentration or light profile shape, and compare with similar estimates based on morphology. Our fits show that Es account for about 20% of the r-band luminosity density, rho_Lr, and 25% of the stellar mass density, rho_*; including S0s and Sas increases these numbers to 33% and 40%, and 50% and 60%, respectively. Summed over all galaxy types, we find rho_* ~ 3 * 10^8 M_Sun Mpc^{-3} at z ~ 0. This is in good agreement with expectations based on integrating the star formation history. However, compared to most previous work, we find an excess of objects at large masses, up to a factor of ~ 10 at M_* ~ 5*10^{11} M_Sun. The stellar mass density further increases at large masses if we assume different IMFs for Es and spiral galaxies, as suggested by some recent chemical evolution models, and results in a better agreement with the dynamical mass function. We also show that the trend for ellipticity to decrease with luminosity is primarily because the E/S0 ratio increases at large L. However, the most massive galaxies, M_* > 5 * 10^{11} M_Sun, are less concentrated and not as round as expected if one extrapolates from lower L, and they are not well-fit by pure deVaucouleur laws. This suggests formation histories with recent radial mergers. Finally, we show that the age-size relation is flat for Es of fixed dynamical mass, but, at fixed M_dyn, S0s and Sas with large sizes tend to be younger. Explaining this difference between E and S0 formation is a new challenge for models of early-type galaxy formation.

We calculate the redshift space correlation function and the power spectrum of density peaks of a Gaussian random field. In the linear regime k < 0.1 h/Mpc, the redshift space power spectrum is P^s_{pk}(k,u) = exp(-f^2 s_{vel}^2 k^2 u^2) * [b_{pk}(k) + b_{vel}(k) f u^2]^2 * P_m(k), where u is the angle with respect to the line of sight, s_{vel} is the one-dimensional velocity dispersion, f is the growth rate, and b_{pk}(k) and b_{vel}(k) are k-dependent linear spatial and velocity bias factors. For peaks, the value of s_{vel} depends upon the functional form of b_{vel}. The peaks model is remarkable because it has unbiased velocities -- peak motions are driven by dark matter flows -- but, in order to achieve this, b_{vel} is k-dependent. We speculate that this is true in general: k-dependence of the spatial bias will lead to k-dependence of b_{vel} even if the biased tracers flow with the dark matter. Because of the k-dependence of the linear bias parameters, standard manipulations applied to the peak model will lead to k-dependent estimates of the growth factor that could erroneously be interpreted as a signature of modified dark energy or gravity. We use the Fisher formalism to show that the constraint on the growth rate f is degraded by a factor of two if one allows for a k-dependent velocity bias of the peak type. We discuss a simple estimate of nonlinear evolution and illustrate the effect of the peak bias on the redshift space multipoles. For k < 0.1 h/Mpc, the peak bias is deterministic but k-dependent, so the configuration space bias is stochastic and scale dependent, both in real and redshift space. We provide expressions for this stochasticity and its evolution (abridged).

We present a new algorithm for identifying the substructure within simulated dark matter haloes. The method is an extension of that proposed by Tormen et al. (2004) and Giocoli et al. (2008a), which identifies a subhalo as a group of self-bound particles that prior to being accreted by the main progenitor of the host halo belonged to one and the same progenitor halo (hereafter satellite). However, this definition does not account for the fact that these satellite haloes themselves may also have substructure, which thus gives rise to sub-subhaloes, etc. Our new algorithm identifies substructures at all levels of this hierarchy, and we use it to determine the mass function of all substructure (counting sub-haloes, sub-subhaloes, etc.). On average, haloes which formed more recently tend to have a larger mass fraction in substructure and to be less concentrated than average haloes of the same mass. We provide quantitative fits to these correlations. Even though our algorithm is very different from that of Gao et al. (2004), we too find that the subhalo mass function per unit mass at redshift z = 0 is universal. This universality extends to any redshift only if one accounts for the fact that host haloes of a given mass are less concentrated at higher redshifts, and concentration and substructure abundance are anti-correlated. This universality allows a simple parametrization of the subhalo mass function integrated over all host halo masses, at any given time. We provide analytic fits to this function which should be useful in halo model analyses which equate galaxies with halo substructure when interpreting clustering in large sky surveys. Finally, we discuss systematic differences in the subhalo mass function that arise from different definitions of (host) halo mass.

In addition to the maximum likelihood approach, there are two other methods which are commonly used to reconstruct the true redshift distribution from photometric redshift datasets: one uses a deconvolution method, and the other a convolution. We show how these two techniques are related, and how this relationship can be extended to include the study of galaxy scaling relations in photometric datasets. We then show what additional information photometric redshift algorithms must output so that they too can be used to study galaxy scaling relations, rather than just redshift distributions. We also argue that the convolution based approach may permit a more efficient selection of the objects for which calibration spectra are required.

In a sample of about 45,700 early-type galaxies extracted from SDSS, we find that the shape, normalization, and dispersion around the mean size-stellar mass relation is the same for young and old systems, provided the stellar mass is greater than 3*10^10 Msun. This is difficult to reproduce in pure passive evolution models, which generically predict older galaxies to be much more compact than younger ones of the same stellar mass. However, this aspect of our measurements is well reproduced by hierarchical models of galaxy formation. Whereas the models predict more compact galaxies at high redshifts, subsequent minor, dry mergers increase the sizes of the more massive objects, resulting in a flat size-age relation at the present time. At lower masses, the models predict that mergers are less frequent, so that the expected anti-correlation between age and size is not completely erased. This is in good agreement with our data: below 3*10^10 Msun, the effective radius R_e is a factor of ~2 lower for older galaxies. These successes of the models are offset by the fact that the predicted sizes have other serious problems, which we discuss.

We study the formation of voids in a modified gravity model in which gravity is generically stronger or weaker on large scales. We show that void abundances provide complementary information to halo abundances: if normalized to the CMB, models with weaker large-scale gravity have smaller large scale power, fewer massive halos and fewer large voids, although the scalings are not completely degenerate with $\sigma_8$. Our results suggest that, in addition to their abundances, halo and void density profiles may also provide interesting constraints on such models: stronger large scale gravity produces more concentrated halos, and thinner void walls. This potentially affects the scaling relations commonly assumed to translate cluster observables to halo masses, potentially making these too, useful probes of gravity.

The peak-background split argument is commonly used to relate the abundance of dark matter halos to their spatial clustering. Testing this argument requires an accurate determination of the halo mass function. We present a Maximum Likelihood method for fitting parametric functional forms to halo abundances which differs from previous work because it does not require binned counts. Our conclusions do not depend on whether we use our method or more conventional ones. In addition, halo abundances depend on how halos are defined. Our conclusions do not depend on the choice of link length associated with the friends-of-friends halo-finder, nor do they change if we identify halos using a spherical overdensity algorithm instead. The large scale halo bias measured from the matter-halo cross spectrum b_x and the halo autocorrelation function b_xi (on scales k~0.03h/Mpc and r ~50 Mpc/h) can differ by as much as 5% for halos that are significantly more massive than the characteristic mass M*. At these large masses, the peak background split estimate of the linear bias factor b1 is 3-5% smaller than b_xi, which is 5% smaller than b_x. We discuss the origin of these discrepancies: deterministic nonlinear local bias, with parameters determined by the peak-background split argument, is unable to account for the discrepancies we see. A simple linear but nonlocal bias model, motivated by peaks theory, may also be difficult to reconcile with our measurements. More work on such nonlocal bias models may be needed to understand the nature of halo bias at this level of precision.

Noisy distance estimates associated with photometric rather than spectroscopic redshifts lead to a mis-estimate of the luminosities, and produce a correlated mis-estimate of the sizes. We consider a sample of early-type galaxies from the SDSS DR6 for which both spectroscopic and photometric information is available, and apply the generalization of the V_max method to correct for these biases. We show that our technique recovers the true redshift, magnitude and size distributions, as well as the true size-luminosity relation. We find that using only 10% of the spectroscopic information randomly spaced in our catalog is sufficient for the reconstructions to be accurate within about 3%, when the photometric redshift error is dz = 0.038. We then address the problem of extending our method to deep redshift catalogs, where only photometric information is available. In addition to the specific applications outlined here, our technique impacts a broader range of studies, when at least one distance-dependent quantity is involved. It is particularly relevant for the next generation of surveys, some of which will only have photometric information.

We show how the excursion set moving barrier model for halo abundances may be generalized to the local non-Gaussian f_{nl} model. Our estimate assumes that the distribution of step sizes depends on f_{nl}, but that they are otherwise uncorrelated. Our analysis is consistent with previous results for the case of a constant barrier, and highlights some implicit assumptions. It also clarifies the basis of an approximate analytic solution to the moving barrier problem in the Gaussian case, and shows how it might be improved.

We generalize Doroshkevich's celebrated formulae for the eigenvalues of the initial shear field associated with Gaussian statistics to the local non-Gaussian f_{nl} model. This is possible because, to at least second order in f_{nl}, distributions at fixed overdensity are unchanged from the case f_{nl}=0. We use this generalization to estimate the effect of f_{nl}\ne 0 on the abundance of virialized halos. Halo abundances are expected to be related to the probability that a certain quantity in the initial fluctuation field exceeds a threshold value, and we study two choices for this variable: it can either be the sum of the eigenvalues of the initial deformation tensor (the initial overdensity), or its smallest eigenvalue. The approach based on a critical overdensity yields results which are in excellent agreement with numerical measurements. We then use these same methods to develop approximations describing the sensitivity of void abundances on f_{nl}. While a positive f_{nl} produces more extremely massive halos, it makes fewer extremely large voids. Its effect thus is qualitatively different from a simple rescaling of the normalisation of the density fluctuation field \sigma_8. Therefore, void abundances furnish complementary information to cluster abundances, and a joint comparison of both might provide interesting constraints on primordial non-Gaussianity.

In hierarchical models, the time derivative of the halo mass function may be thought of as the difference of two terms - a creation term, which describes the increase in the number of haloes of mass m from mergers of less massive objects, and a destruction term, which describes the decrease in the number of m-haloes as these merge with other haloes, creating more massive haloes as a result. In models where haloes form from a spherical collapse, the distribution of halo creation times can be estimated from the same formalism which is used to estimate halo abundances: the constant-barrier excursion-set approach. In the excursion set approach, moving, rather than constant-barriers, are necessary for estimating halo abundances when the collapse is triaxial. We show how the corresponding estimates are modified by providing analytic estimates of the creation rate for a wide range of halo masses. We then show that these moving-barrier based predictions are in better agreement with measurements in numerical simulations than are the corresponding predictions of the spherical collapse model. These results should be useful for studies of merger-driven star-formation rates and AGN activity. We also present a similar study of creation of haloes conditioned on belonging to an object of a certain mass today, and reach similar conclusions - the moving barrier based estimates are in substantially better agreement with the simulations. This part of the study may be useful for understanding the tendency for the oldest stars to exist in the most massive objects.

We present measurements of the clustering of hot and cold patches in the microwave background sky as measured from the Wilkinson Microwave Anisotropy Probe (WMAP) five-year data. These measurements are compared with theoretical predictions which assume that the cosmological signal obeys Gaussian statistics. We find significant differences from the simplest Gaussian-based prediction. However, the measurements are sensitive to the fact that the noise is spatially inhomogeneous (e.g., because different parts of the sky were observed for different lengths of time). We show how to account for this spatial inhomogeneity when making predictions. Differences from the Gaussian-based expectation remain even after this more careful accounting of the noise. In particular, we note that hot and cold pixels cluster differently within the same temperature thresholds at few-degree scales. While these findings may indicate primordial non-Gaussianity, we discuss other plausible explanations for these discrepancies. In addition, we find some deviations from Gaussianity at sub-degree scales, especially in the W band, whose origin may be associated with extragalactic dust emission.

Modifications to the gravitational potential affect the nonlinear gravitational evolution of large scale structures in the Universe. To illustrate some generic features of such changes, we study the evolution of spherically symmetric perturbations when the modification is of Yukawa type; this is non-trivial, because we should not and do not assume that Birkhoff's theorem applies. We then show how to estimate the abundance of virialized objects in such models. Comparison with numerical simulations shows reasonable agreement: When normalized to have the same fluctuations at early times, weaker large scale gravity produces fewer massive halos. However, the opposite can be true for models that are normalized to have the same linear theory power spectrum today, so the abundance of rich clusters potentially places interesting constraints on such models. Our analysis also indicates that the formation histories and abundances of sufficiently low mass objects are unchanged from standard gravity. This explains why simulations have found that the nonlinear power-spectrum at large k is unaffected by such modifications to the gravitational potential. In addition, the most massive objects in CMB-normalized models with weaker gravity are expected to be similar to the high-redshift progenitors of the most massive objects in models with stronger gravity. Thus, the difference between the cluster and field galaxy populations is expected to be larger in models with stronger large-scale gravity.

In studies of the environmental dependence of structure formation, the large scale environment is often thought of as providing an effective background cosmology: e.g. the formation of structure in voids is expected to be just like that in a less dense universe with appropriately modified Hubble and cosmological constants. However, in the excursion set description of structure formation which is commonly used to model this effect, no explicit mention is made of the effective cosmology. Rather, this approach uses the spherical evolution model to compute an effective linear theory growth factor, which is then used to predict the growth and evolution of nonlinear structures. We show that these approaches are, in fact, equivalent: a consequence of Birkhoff's theorem. We speculate that this equivalence will not survive in models where the gravitational force law is modified from an inverse square, potentially making the environmental dependence of clustering a good test of such models.

We present an algorithm for generating merger histories of dark matter haloes. The algorithm is based on the excursion set approach with moving barriers whose shape is motivated by the ellipsoidal collapse model of halo formation. In contrast to most other merger-tree algorithms, ours takes discrete steps in mass rather than time. This allows us to quantify effects which arise from the fact that outputs from numerical simulations are usually in discrete time bins. In addition, it suggests a natural set of scaling variables for describing the abundance of halo progenitors; this scaling is not as general as that associated with a spherical collapse. We test our algorithm by comparing its predictions with measurements in numerical simulations. The progenitor mass fractions and mass functions are in good agreement, as is the predicted scaling law. We also test the formation-redshift distribution, the mass distribution at formation, and the redshift distribution of the most recent major merger; all are in reasonable agreement with N-body simulation data, over a broad range of masses and redshifts. Finally, we study the effects of sampling in discrete time snapshots. In all cases, the improvement over algorithms based on the spherical collapse assumption is significant.

We use the spherical evolution approximation to investigate nonlinear evolution from the non-Gaussian initial conditions characteristic of the local f_nl model. We provide an analytic formula for the nonlinearly evolved probability distribution function of the dark matter which shows that the under-dense tail of the nonlinear PDF in the f_nl model should differ significantly from that for Gaussian initial conditions. Measurements of the under-dense tail in numerical simulations may be affected by discreteness effects, and we use a Poisson counting model to describe this effect. Once this has been accounted for, our model is in good quantitative agreement with the simulations.

Linear theory provides a reasonable description of the velocity correlations of biased tracers both perpendicular and parallel to the line of separation, provided one accounts for the fact that the measurement is almost always made using pair-weighted statistics. This introduces an additional term which, for sufficiently biased tracers, may be large. Previous work suggesting that linear theory was grossly in error for the components parallel to the line of separation ignored this term.

We estimate the absolute magnitude distribution of galaxies which lie within about a Mpc of Mg II absorption systems. The absorption systems themselves lie along 1880 lines of sight to QSOs from the Sloan Digital Sky Survey Data Release 3, have rest equivalent widths greater than 0.88 Angstroms, and redshifts between 0.37 < z < 0.82. Our measurement is based on all galaxies which lie within a projected distance of about 900 kpc/h of each QSO demonstrating absorption. The redshifts of these projected neighbors are not available, so we use a background subtraction technique to estimate the absolute magnitude distribution of true neighbors. (Our method exploits the fact that, although we do not know the redshifts of the neighbors, we do know the redshift of the absorbers.) The absolute magnitude distribution we find is well described by a bell-shaped curve peaking at about rest-frame M_B = -20, corresponding to L/L* = 1.4. A comparison of this observed distribution to ones in the literature suggests that it is unlikely to be drawn from a population dominated by late-type galaxies. However, the strong equivalent width systems may be associated with later galaxy types. Finally we use the absolute magnitude distribution, along with the observed covering fraction of about 8 percent, to estimate the extent of the MgII absorbing gas around a galaxy. For an L* galaxy, this scale is about 70 kpc/h. We provide an analytic description of our method, which is generally applicable to any dataset in which redshifts are only available for a small sub-sample. Hence, we expect it to aid in the analysis of galaxy scaling relations from photometric redshift datasets.

Successful halo-model descriptions of the luminosity dependence of clustering distinguish between the central galaxy in a halo and all the others (satellites). To include colors, we provide a prescription for how the color-magnitude relation of centrals and satellites depends on halo mass. This follows from two assumptions: (i) the bimodality of the color distribution at fixed luminosity is independent of halo mass, and (ii) the fraction of satellite galaxies which populate the red sequence increases with luminosity. We show that these two assumptions allow one to build a model of how galaxy clustering depends on color without any additional free parameters than those required to model the luminosity dependence of galaxy clustering. We then show that the resulting model is in good agreement with the distribution and clustering of colors in the SDSS, both by comparing the predicted correlation functions of red and blue galaxies with measurements, and by comparing the predicted color mark correlation function with the measured one. Mark correlation functions are powerful tools for identifying and quantifying correlations between galaxy properties and their environments: our results indicate that the correlation between halo mass and environment is the primary driver for correlations between galaxy colors and the environment; additional correlations associated with halo `assembly bias' are relatively small. Our approach shows explicitly how to construct mock catalogs which include both luminosities {\em and} colors -- thus providing realistic training sets for, e.g., galaxy cluster finding algorithms. Our prescription is the first step towards incorporating the entire spectral energy distribution into the halo model approach.

We used the Advanced Camera for Surveys on board the Hubble Space Telescope to obtain high resolution i-band images of the centers of 23 single galaxies, which were selected because they have SDSS velocity dispersions larger than 350 km/s. The surface brightness profiles of the most luminous of these objects (M_i<-24) have well-resolved `cores' on scales of 150-1000 pc, and share similar properties to BCGs. The total luminosity of the galaxy is a better predictor of the core size than is the velocity dispersion. The correlations of luminosity and velocity dispersion with core size agree with those seen in previous studies of galaxy cores. Because of high velocity dispersions, our sample of galaxies can be expected to harbor the most massive black holes, and thus have large cores with large amounts of mass ejection. The mass-deficits inferred from core-Sersic fits to the surface-brightness profiles are approximately double the black-hole masses inferred from the M_bh-sigma relation and the same as those inferred from the M_bh-L relation. The less luminous galaxies (M_i>-23) tend to have steeper `power-law' inner profiles, higher-ellipticity, diskier isophotes, and bulge-to-total ratios of order 0.5 -- all of which suggest that they are `fast-rotators' and rotational motions could have contaminated the velocity dispersion estimate. There are obvious dust features within about 300 pc of the center in about 35% of the sample, predominantly in power-law rather than core galaxies.

We study a sample of 43 early-type galaxies, selected from the SDSS because they appeared to have velocity dispersion > 350 km/s. High-resolution photometry in the SDSS i passband using HRC-ACS on board the HST shows that just less than half of the sample is made up of superpositions of two or three galaxies, so the reported velocity dispersion is incorrect. The other half of the sample is made up of single objects with genuinely large velocity dispersions. None of these objects has sigma larger than 426 +- 30 km/s. These objects define rather different relations than the bulk of the early-type galaxy population: for their luminosities, they are the smallest, most massive and densest galaxies in the Universe. Although the slopes of the scaling relations they define are rather different from those of the bulk of the population, they lie approximately parallel to those of the bulk "at fixed sigma". These objects appear to be of two distinct types: the less luminous (M_r>-23) objects are rather flattened and extremely dense for their luminosities -- their properties suggest some amount of rotational support and merger histories with abnormally large amounts of gaseous dissipation. The more luminous objects (M_r<-23) tend to be round and to lie in or at the centers of clusters. Their properties are consistent with the hypothesis that they are BCGs. Models in which BCGs form from predominantly radial mergers having little angular momentum predict that they should be prolate. If viewed along the major axis, such objects would appear to have abnormally large sigma for their sizes, and to be abnormally round for their luminosities. This is true of the objects in our sample once we account for the fact that the most luminous galaxies (M_r<-23.5), and BCGs, become slightly less round with increasing luminosity.

We use the physics of ellipsoidal collapse to model the probability distribution function of the smoothed dark matter density field in real and redshift space. We provide a simple approximation to the exact collapse model which shows clearly how the evolution can be thought of as a modification of the spherical evolution model as well as of the Zeldovich Approximation. In essence, our model specifies how the initial smoothed overdensity and shear fields can be used to determine the shape and size of the region at later times. We use our parametrization to extend previous work on the real-space PDF so that it predicts the redshift space PDF as well. Our results are in good agreement with measurements of the PDF in simulations of clustering from Gaussian initial conditions down to scales on which the rms fluctuation is slightly greater than unity. We also show how the highly non-Gaussian non-linear redshifted density field in a numerical simulation can be transformed so that it provides an estimate of the shape of the initial real-space PDF. When applied to our simulations, our method recovers the initial Gaussian PDF, provided the variance in the nonlinear smoothed field is less than 4.

We present an analysis of the small-to-intermediate scale clustering of samples of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey and the 2dF-SDSS LRG and QSO (2SLAQ) survey carefully matched to have the same rest-frame colours and luminosity. We study the spatial two-point auto-correlation function in both redshift-space and real-space of a combined sample of over 10,000 LRGs, which represent the most massive galaxies in the universe with stellar masses > 10^11 h^-1 M_sun and space densities 10^-4 h^-3 Mpc^-3. We find no significant evolution in the amplitude r_0 of the correlation function with redshift, but do see a slight decrease in the slope with increasing redshift over 0.19 < z < 0.55 and scales of 0.32 < r < 32 h^-1 Mpc. We compare our measurements with the predicted evolution of dark matter clustering and use the halo model to interpret our results. We find that our clustering measurements are inconsistent (>99.9% significance) with a passive model whereby the LRGs do not merge with one another; a model with a merger rate of 7.5 +/- 2.3% from z = 0.55 to z = 0.19 (i.e. an average rate of 2.4% Gyr^-1) provides a better fit to our observations. Our clustering and number density measurements are consistent with the hypothesis that the merged LRGs were originally central galaxies in different haloes which, following the merger of these haloes, merged to create a single Brightest Cluster Galaxy. In addition, we show that the small-scale clustering signal constrains the scatter in halo merger histories. When combined with measurements of the luminosity function, our results suggest that this scatter is sub-Poisson. While this is a generic prediction of hierarchical models, it has not been tested before.

Halo model interpretations of the luminosity dependence of galaxy clustering assume that there is a central galaxy in every sufficiently massive halo, and that this central galaxy is very different from all the others in the halo. The halo model decomposition makes the remarkable prediction that the mean luminosity of the non-central galaxies in a halo should be almost independent of halo mass: the predicted increase is about 20% while the halo mass increases by a factor of more than 20. In contrast, the luminosity of the central object is predicted to increase approximately linearly with halo mass at low to intermediate masses, and logarithmically at high masses. We show that this weak, almost non-existent mass-dependence of the satellites is in excellent agreement with the satellite population in group catalogs constructed by two different collaborations. This is remarkable, because the halo model prediction was made without ever identifying groups and clusters. The halo model also predicts that the number of satellites in a halo is drawn from a Poisson distribution with mean which depends on halo mass. This, combined with the weak dependence of satellite luminosity on halo mass, suggests that the Scott effect, such that the luminosities of very bright galaxies are merely the statistically extreme values of a general luminosity distribution, may better apply to the most luminous satellite galaxy in a halo than to BCGs. If galaxies are identified with halo substructure at the present time, then central galaxies should be about 4 times more massive than satellite galaxies of the same luminosity, whereas the differences between the stellar M/L ratios should be smaller. Therefore, a comparison of the weak lensing signal from central and satellite galaxies should provide useful constraints. [abridged]

Noisy distance estimates associated with photometric rather than spectroscopic redshifts lead to a mis-estimate of the luminosities, and produce a correlated mis-estimate of the sizes. We consider a sample of early-type galaxies from the SDSS DR6 and apply the generalization of the V_max method to correct for these biases. We show that our technique recovers the true redshift, magnitude and size distributions, as well as the true size-luminosity relation. Regardless the specific application outlined here, our method impacts a broader range of studies, when at least one distance-dependent quantity is involved.

Nonlinear evolution can sometimes be modelled by a deterministic mapping from initial to final of the local smoothed overdensity. Perturbation theory methods base on this deterministic and local mapping and ignore the 'cloud-in-cloud' effect, while the excursion set approach methods take this nonlocality into account. We compared these methods using the spherical collapse mapping and showed that, on scales where the rms fluctuation is small, both models give similar results and they are in good agreement with numerical simulations. If the deterministic mapping depends on quantities other than overdensity, this will also manifest as stochasticity if the other quantities are ignored. We considered the Zeldovich approximation and Ellipsoidal Collapse model, both include the tidal field in the evolution. Our anaylsis shows that the change in cell shape effect should be included on scales where the rms is of order of unity or larger. On scales where the rms is less than 2 methods based on the spherical collapse model allow a rather accurate reconstruction of the shape of the initial distribution from the nonlinear field. This can be used as the basis for constraining the statistical properties of the initial fluctuation field. (Abridge)

We develop an analytic theory for the redshift space bispectrum of dark matter, haloes and galaxies. This is done within the context of the halo model of structure formation, as this allows for the self-consistent inclusion of linear and non-linear redshift space distortions and also for the non-linearity of the halo bias. The model is applicable over a wide range of scales: on the largest scales the predictions reduce to those of the standard perturbation theory (PT); on smaller scales they are determined primarily by the nonlinear virial velocities of galaxies within haloes, and this gives rise to the U-shaped anisotropy in the reduced bispectrum -- a finger print of the Finger-Of-God distortions. We then confront the predictions with bispectrum measurements from an ensemble of numerical simulations. On very large scales, k=0.05 h/Mpc, we find reasonably good agreement between our Halo Model, PT and the data, to within the errors. On smaller scales, k=0.1 h/Mpc, the measured bispectra differ from the PT at the level of 10-20%, especially for colinear triangle configurations. The Halo Model predictions improve over PT, but are accurate to no better than 10%. On smaller scales k=0.5-1.0 h/Mpc, our model provides a significant improvement over PT, which breaks down. This implies that studies which use the lowest order PT to extract galaxy bias information are not robust on scales k>0.1 h/Mpc. The analytic and simulation results also indicate that there is no observable scale for which the configuration dependence of the reduced bispectrum is constant--hierarchical models for the higher order correlation functions in redshift space are unlikely to be useful. It is hoped that our model will facilitate extraction of information from large-scale structure surveys of the Universe.

We apply the Minimum Description Length model selection approach to the detection of extra-solar planets, and use this example to show how specification of the experimental design affects the prior distribution on the model parameter space and hence the posterior likelihood which, in turn, determines which model is regarded as most `correct'. Our analysis shows how conditioning on the experimental design can render a non-compact parameter space effectively compact, so that the MDL model selection problem becomes well-defined.

The baryonic acoustic signature in the large-scale clustering pattern of galaxies has been detected in the two-point correlation function. Its precise spatial scale has been forwarded as a rigid-rod ruler test for the space-time geometry, and hence as a probe for tracking the evolution of Dark Energy. Percent-level shifts in the measured position can bias such a test and erode its power to constrain cosmology. This paper addresses some of the systematic effects that might induce shifts: namely non-linear corrections from matter evolution, redshift space distortions and biasing. We tackle these questions through analytic methods and through a large battery of numerical simulations, with total volume of the order ~100[Gpc\h]^3. A toy-model calculation shows that if the non-linear corrections simply smooth the acoustic peak, then this gives rise to an `apparent' shifting to smaller scales. However if tilts in the broad band power spectrum are induced then this gives rise to more pernicious `physical' shifts. Our numerical simulations show evidence of both: in real space and at z=0, for the dark matter we find percent level shifts; for haloes the shifts depend on halo mass, with larger shifts being found for the most biased samples, up to 3%. From our analysis we find that physical shifts are greater than ~0.4% at z=0. In redshift space these effects are exacerbated, but at higher redshifts are alleviated. We develop an analytical model to understand this, based on solutions to the pair conservation equation using characteristic curves. When combined with modeling of pairwise velocities the model reproduces the main trends found in the data. The model may also help to unbias the acoustic peak.

Many physical properties of galaxies correlate with one another, and these correlations are often used to constrain galaxy formation models. Such correlations include the color-magnitude relation, the luminosity-size relation, the Fundamental Plane, etc. However, the transformation from observable (e.g. angular size, apparent brightness) to physical quantity (physical size, luminosity), is often distance-dependent. Noise in the distance estimate will lead to biased estimates of these correlations, thus compromising the ability of photometric redshift surveys to constrain galaxy formation models. We describe two methods which can remove this bias. One is a generalization of the V_max method, and the other is a maximum likelihood approach. We illustrate their effectiveness by studying the size-luminosity relation in a mock catalog, although both methods can be applied to other scaling relations as well. We show that if one simply uses photometric redshifts one obtains a biased relation; our methods correct for this bias and recover the true relation.

The Sunyaev Zel'dovich (SZ) signal is highly non-Gaussian, so the SZ power spectrum (along with the mean $y$ parameter) does not provide a complete description of the SZ effect. Therefore, SZ-based constraints on cosmological parameters and on cluster gastrophysics which assume Gaussianity will be biased. We derive an analytic expression for the $n$-point joint PDF of the SZ power spectrum. Our derivation, which is based on the halo model, has several advantages: it is expressed in an integral form which allows quick computation; it is applicable to any given survey and any given angular scale; it is straightforward to incorporate many of the complexities which arise when modeling the SZ signal. To illustrate, we use our expression to estimate $p(C_\ell)$, the one-point PDF of the SZ power spectrum. For small sky coverage (applicable to BIMA/CBI and the Sunyaev Zel'dovich Array experiments), our analysis shows that $p(C_\ell)$ on the several arc-minute scale is expected to be strongly skewed, peaking at a value well below the mean and with a long tail which extends to tail high $C_\ell$ values. In the limit of large sky coverage (applicable to the South Pole Telescope and Planck), $p(C_\ell)$ approaches a Gaussian form. However, even in this limit, the variance of the power spectrum is very different from the naive Gaussian-based estimate. This is because different $\ell$ models are strongly correlated, making the cosmic variance of the SZ band-power much larger than the naive estimate. Our analysis should also be useful for modeling the PDF of the power spectrum induced by gravitational lensing at large $\ell$.

We examine the scale dependence of dark matter, halo and galaxy clustering on very large scales (0.01<k[h/Mpc]<0.15), due to non-linear effects from dynamics and halo bias. We pursue a two line offensive: high resolution numerical simulations are used to establish several new results, and an analytic model is developed to understand their origins. Our simulations show: (i) that the z=0 dark matter power spectrum is suppressed relative to linear theory by ~5% on scales (0.05<k[h/Mpc]<0.075); (ii) that, indeed, halo bias is non-linear over the scales we probe and that the scale dependence is a strong function of halo mass. High mass haloes show no suppression of power on scales (k<0.07[h/Mpc]), and only show amplification on smaller scales, whereas low mass haloes show strong, ~5-10%, suppression over the range (0.05 <k[h/Mpc] <0.15). Our results have relevance for studies of the baryon acoustic oscillation features. Non-linear mode-mode coupling: (i) damps these features on progressively larger scales as halo mass increases; (ii) produces small shifts in the positions of the peaks and troughs which depend on halo mass. Our analytic model is described in the language of the `halo-model'. However, for the first time the halo-halo clustering term is propagated into the non-linear regime using `1-loop' perturbation theory and a non-linear halo bias model. We show that, with bias parameters derived from simulations, the model predictions are in agreement with the numerical results. We then use the model to explore the scale dependence of galaxies of different colour and find significant differences between the power spectra of the two populations. Thus understanding the scale dependent bias for a given galaxy sample will be crucial for deriving accurate cosmological constraints. (Abridged)

The luminosity functions of galaxies and quasars provide invaluable information about galaxy and quasar formation. Estimating the luminosity function from magnitude limited samples is relatively straightforward, provided that the distances to the objects in the sample are known accurately; techniques for doing this have been available for about thirty years. However, distances are usually known accurately for only a small subset of the sample. This is true of the objects in the Sloan Digital Sky Survey, and will be increasingly true of the next generation of deep multi-color photometric surveys. Estimating the luminosity function when distances are only known approximately (e.g., photometric redshifts are available, but spectroscopic redshifts are not) is more difficult. I describe two algorithms which can handle this complication: one is a generalization of the V_max algorithm, and the other is a maximum likelihood approach. Because these methods account for uncertainties in the distance estimate, they impact a broader range of studies. For example, they are useful for studying the abundances of galaxies which are sufficiently nearby that the contribution of peculiar velocity to the spectroscopic redshift is not negligible, so only a noisy estimate of the true distance is available. In this respect, peculiar velocities and photometric redshift errors have similar effects. The methods developed here are also useful for estimating the stellar luminosity function in samples where accurate parallax distances are not available.

Black hole masses are tightly correlated with the stellar velocity dispersions of the bulges which surround them, and slightly less-well correlated with the bulge luminosity. It is common to use these correlations to estimate the expected abundance of massive black holes. This is usually done by starting from an observed distribution of velocity dispersions or luminosities and then changing variables. This procedure neglects the fact that there is intrinsic scatter in these black hole mass--observable correlations. Accounting for this scatter results in estimates of black hole abundances which are larger by almost an order of magnitude at masses >10^9 M_sun. Including this scatter is particularly important for models which seek to infer quasar lifetimes and duty cycles from the local black hole mass function. However, even when scatter has been accounted for, the M_bh-sigma relation predicts fewer massive black holes than does the M_bh-L relation. This is because the sigma-L relation in the black hole samples currently available is inconsistent with that in the SDSS sample from which the distributions of L or sigma are based: the black hole samples have smaller L for a given sigma, or larger sigma for a given L. The sigma-L relation in the black hole samples is similarly discrepant with that in other samples of nearby early-type galaxies. This suggests that current black hole samples are biased: if this is a selection rather than physical effect, then the M_bh-sigma and M_bh-L relations currently in the literature are also biased from their true values.

We discuss two slightly counter-intuitive findings about the environmental dependence of clustering in the Sloan Digital Sky Survey. First, we find that the relation between clustering strength and density is not monotonic: galaxies in the densest regions are more strongly clustered than are galaxies in regions of moderate overdensity; galaxies in moderate overdensities are more strongly clustered than are those in moderate underdensities; but galaxies in moderate underdensities are less clustered than galaxies in the least dense regions. We argue that this is natural if clustering evolved gravitationally from a Gaussian field, since the highest peaks and lowest troughs in Gaussian fields are similarly clustered. The precise non-monotonic dependence of galaxy clustering on density is very well reproduced in a mock catalog which is based on a halo-model decomposition of galaxy clustering. Second, the distribution of galaxy counts in our sample is rather well described by a Poisson cluster model. We show that, despite their Poisson nature, correlations with environment are expected in such models. More remarkably, the expected trends are very like those in standard models of halo bias, despite the fact that correlations with environment in these models arise purely from the fact that dense regions are dense because they happen to host more massive halos. This is in contrast to the usual analysis which assumes that it is the large scale environment which determines the halo mass function. [Abridged]

We present measurements of the redshift-space three-point correlation function of 50,967 Luminous Red Galaxies (LRGs) from Data Release 3 (DR3) of the Sloan Digital Sky Survey (SDSS). We have studied the shape dependence of the reduced three-point correlation function (Qz(s,q,theta)) on three different scales, s=4, 7 and 10 h-1 Mpc, and over the range of 1 < q < 3 and 0 < theta < 180. On small scales (s=4 h-1 Mpc), Qz is nearly constant, with little change as a function of q and theta. However, there is evidence for a shallow U-shaped behaviour (with theta) which is expected from theoretical modeling of Qz . On larger scales (s=7 and 10 h-1 Mpc), the U-shaped anisotropy in Qz (with theta) is more clearly detected. We compare this shape-dependence in Qz(s,q,theta) with that seen in mock galaxy catalogues which were generated by populating the dark matter halos in large N-body simulations with mock galaxies using various Halo Occupation Distributions (HOD). We find that the combination of the observed number density of LRGs, the (redshift-space) two-point correlation function and Qz provides a strong constraint on the allowed HOD parameters (M_min, M_1, alpha) and breaks key degeneracies between these parameters. For example, our observed Qz disfavors mock catalogues that overpopulate massive dark matter halos with many LRG satellites. We also estimate the linear bias of LRGs to be b=1.87+/-0.07 in excellent agreement with other measurements.

A dark matter halo is said to have formed when at least half its mass hass been assembled into a single progenitor. With this definition, it is possible to derive a simple but useful analytic estimate of the distribution of halo formation times. The standard estimate of this distribution depends on the shape of the conditional mass function--the distribution of progenitor masses of a halo as a function of time. If the spherical collapse model is used to estimate the progenitor mass function, then the formation times one infers systematically underestimate those seen in numerical simulations of hierarchical gravitational clustering. We provide estimates of halo formation which may be related to an ellipsoidal collapse model. These estimates provide a substantially better description of the simulations. We also provide an alternative derivation of the formation time distribution which is based on the assumption that haloes increase their mass through binary mergers only. Our results are useful for models which relate halo structure to halo formation.

It is common to estimate black hole abundances by using a measured correlation between black hole mass and another more easily measured observable such as the velocity dispersion or luminosity of the surrounding bulge. The correlation is used to transform the distribution of the observable into an estimate of the distribution of black hole masses. However, different observables provide different estimates: the Mbh-sigma relation predicts fewer massive black holes than does the Mbh-L relation. This is because the sigma-L relation in black hole samples currently available is inconsistent with that in the SDSS sample, from which the distributions of L or sigma are based: the black hole samples have smaller L for a given sigma or have larger sigma for a given L. This is true whether L is estimated in the optical or in the NIR. If this is a selection rather than physical effect, then the Mbh-sigma and Mbh-L relations currently in the literature are also biased from their true values. We provide a framework for describing the effect of this bias. We then combine it with a model of the bias to make an estimate of the true intrinsic relations. While we do not claim to have understood the source of the bias, our simple model is able to reproduce the observed trends. If we have correctly modeled the selection effect, then our analysis suggests that the bias in the <Mbh|sigma> relation is likely to be small, whereas the <Mbh|L> relation is biased towards predicting more massive black holes for a given luminosity. In addition, it is likely that the Mbh-L relation is entirely a consequence of more fundamental relations between Mbh and sigma, and between sigma and L. The intrinsic relation we find suggests that at fixed luminosity, older galaxies tend to host more massive black holes.

The size-luminosity relation of early-type Brightest Cluster Galaxies (BCGs), R_e ~ L^0.88, is steeper than that for the bulk of the early-type galaxy population, for which R_e ~ L^0.68. In addition, although BCGs are hardly offset from the Fundamental Plane defined by the bulk of the early-type population, they show considerably smaller scatter. The larger than expected sizes of BCGs, and the increased homogeneity, are qualitatively consistent with models which seek to explain the colors of the most massive galaxies by invoking dry dissipationless mergers, since dissipation tends to reduce the sizes of galaxies, and wet mergers which result in star formation would tend to increase the scatter in luminosity at fixed size and velocity dispersion. Furthermore, BCGs define the same g-r color-magnitude relation as the bulk of the early-type population. If BCGs formed from dry mergers, then BCG progenitors must have been red for their magnitudes, suggesting that they hosted older stellar populations than typical for their luminosities. Our findings have two other consequences. First, the R_e-L relation of the early-type galaxy population as a whole (i.e., normal plus BCG) exhibits some curvature. Some of this curvature must be a consequence of the fact that an increasing fraction of the most luminous galaxies are BCGs. The second consequence is suggested by the fact that, despite following a steeper size-luminosity relation, BCGs tend to define a tight relation between dynamical mass R_e sigma^2/G and luminosity. As consequence, we find that BCGs define a shallower sigma-L relation than the bulk of the early-type galaxy population.

We consider nonlinear biasing models of galaxies with particular attention to a correlation between linear and quadratic biasing coefficients, b_1 and b_2. We first derive perturbative expressions for b_1 and b_2 in halo and peak biasing models. Then we compute power spectra and bispectra of dark matter particles and halos using N-body simulation data and of volume-limited subsamples of Sloan Digital Sky Survey (SDSS) galaxies, and determine their b_1 and b_2. We find that the values of those coefficients at linear regimes (k<0.2h/Mpc) are fairly insensitive to the redshift-space distortion and the survey volume shape. The resulting normalized amplitudes of bispectra, Q, for equilateral triangles, are insensitive to the values of b_1 implying that b_2 indeed correlates with b_1. The present results explain the previous finding of Kayo et al. (2004) for the hierarchical relation of three-point correlation functions of SDSS galaxies. While the relations between b_1 and b_2 are quantitatively different for specific biasing models, their approximately similar correlations indicate a fairly generic outcome of the biasing due to the gravity in primordial Gaussian density fields.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

A generic prediction of hierarchical clustering models is that the mass function of dark haloes in dense regions in the Universe should be top-heavy. We provide a novel test of this prediction using a sample of galaxies drawn from the Sloan Digital Sky Survey. To perform the test, we compare measurements of galaxy clustering in dense and underdense regions. We find that galaxies in dense regions cluster significantly more strongly than those in less dense regions. This is true over the entire 0.1--30 Mpc pair separation range for which we can make accurate measurements. We make similar measurements in realistic mock catalogs in which the only environmental effects are those which arise from the predicted correlation between halo mass and environment. We also provide an analytic halo-model based calculation of the effect. Both the mock catalogs and the analytic calculation provide rather good descriptions of the SDSS measurements. Thus, our results provide strong support for hierarchical models. They suggest that, unless care is taken to study galaxies at fixed mass, correlations between galaxy properties and the surrounding environment are almost entirely due to more fundamental correlations between galaxy properties and host halo mass, and between halo mass and environment.

We present a measurement of the probability distribution function (PDF) of the transmitted flux in the Lya forest from a sample of 3492 quasars included in the SDSS DR3 data release. Our intention is to investigate the sensitivity of the Lya flux PDF as measured from low resolution and low signal-to-noise data to a number of systematic errors such as uncertainties in the mean flux, continuum and noise estimate. The quasar continuum is described by the superposition of a power law and emission lines. We perform a power law continuum fitting on a spectrum-by-spectrum basis, and obtain an average continuum slope of 0.59 +/- 0.36 in the redshift range 2.5<z<3.5. Taking into account the variation in the continuum indices increases the mean flux by 3 and 7 per cent at z=3 and 2.4, respectively, as compared to the values inferred with a single (mean) continuum slope. We compare our measurements to the PDF obtained with mock lognormal spectra, whose statistical properties have been constrained to match the observed Lya flux PDF and power spectrum of high resolution data. Using our power law continuum fitting and the SDSS pipeline noise estimate yields a poor agreement between the observed and mock PDFs. Allowing for a break in the continuum slope and, more importantly, for residual scatter in the continuum level substantially improves the agreement. A decrease of 10-15 per cent in the mean quasar continuum with a typical rms variance at the 20 per cent level can account for the data, provided that the noise excess correction is no larger than 10 per cent.

We present a mark correlation analysis of the galaxies in the Sloan Digital Sky Survey using weights provided by MOPED. The large size of the sample permits statistically significant statements about how galaxies with different metallicities and star formation histories are spatially correlated. Massive objects formed a larger fraction of their stars at higher redshifts and over shorter timescales than did less massive objects (sometimes called down-sizing). We find that those galaxies which dominated the cosmic star formation at z~3 are predominantly in clusters today, whereas galaxies which dominate the star formation at z~0 inhabit substantially lower mass objects in less dense regions today. Hence, our results indicate that star formation and chemical enrichment occured first in the denser regions of the Universe, and moved to less dense regions at later times.

We discuss an analytic approach for modeling structure formation in sheets, filaments and knots. This is accomplished by combining models of triaxial collapse with the excursion set approach: sheets are defined as objects which have collapsed along only one axis, filaments have collapsed along two axes, and halos are objects in which triaxial collapse is complete. In the simplest version of this approach, which we develop here, large scale structure shows a clear hierarchy of morphologies: the mass in large-scale sheets is partitioned up among lower mass filaments, which themselves are made-up of still lower mass halos. Our approach provides analytic estimates of the mass fraction in sheets, filaments and halos, and its evolution, for any background cosmological model and any initial fluctuation spectrum. In the currently popular $\Lambda$CDM model, our analysis suggests that more than 99% of the cosmic mass is in sheets, and 72% in filaments, with mass larger than $10^{10} M_{\odot}$ at the present time. For halos, this number is only 46%. Our approach also provides analytic estimates of how halo abundances at any given time correlate with the morphology of the surrounding large-scale structure, and how halo evolution correlates with the morphology of large scale structure.

We present measurements of the normalised redshift-space three-point correlation function (Q_z) of galaxies from the Sloan Digital Sky Survey (SDSS) main galaxy sample. We have applied our "npt" algorithm to both a volume-limited (36738 galaxies) and magnitude-limited sample (134741 galaxies) of SDSS galaxies, and find consistent results between the two samples, thus confirming the weak luminosity dependence of Q_z recently seen by other authors. We compare our results to other Q_z measurements in the literature and find it to be consistent within the full jack-knife error estimates. However, we find these errors are significantly increased by the presence of the ``Sloan Great Wall'' (at z ~ 0.08) within these two SDSS datasets, which changes the 3-point correlation function (3PCF) by 70% on large scales (s>=10h^-1 Mpc). If we exclude this supercluster, our observed Q_z is in better agreement with that obtained from the 2dFGRS by other authors, thus demonstrating the sensitivity of these higher-order correlation functions to large-scale structures in the Universe. This analysis highlights that the SDSS datasets used here are not ``fair samples'' of the Universe for the estimation of higher-order clustering statistics and larger volumes are required. We study the shape-dependence of Q_z(s,q,theta) as one expects this measurement to depend on scale if the large scale structure in the Universe has grown via gravitational instability from Gaussian initial conditions. On small scales (s <= 6h^-1 Mpc), we see some evidence for shape-dependence in Q_z, but at present our measurements are consistent with a constant within the errors (Q_z ~ 0.75 +/- 0.05). On scales >10h^-1 Mpc, we see considerable shape-dependence in Q_z.

We present measurements of the redshift-space luminosity-weighted or `marked' correlation function in the SDSS. These are compared with a model in which the luminosity function and luminosity dependence of clustering are the same as that observed, and in which the form of the luminosity-weighted correlation function is entirely a consequence of the fact that massive halos populate dense regions. We do this by using mock catalogs which are constrained to reproduce the observed luminosity function and the luminosity dependence of clustering, as well as by using the language of the redshift-space halo-model. These analyses show that marked correlations may show a signal on large scales even if there are no large-scale physical effects--the statistical correlation between halos and their environment will produce a measureable signal. Our model is in good agreement with the measurements, indicating that the halo mass function in dense regions is top-heavy; the correlation between halo mass and large scale environment is the primary driver for correlations between galaxy properties and environment; and the luminosity of the central galaxy in a halo is different from (in general, brighter than) that of the other objects in the halo. Thus our measurement provides strong evidence for the accuracy of these three standard assumptions of galaxy formation models. These assumptions also form the basis of current halo-model based interpretations of galaxy clustering.

We study the photometric and spectral properties of 39320 early-type galaxies within the SDSS, as a function of both local environment and redshift. The distance to the nearest cluster of galaxies and the distance to the 10th nearest luminous neighbor were used to define two extremes in environment. The properties of early-type galaxies are weakly but significantly different in these two extremes: objects in dense environments are less than ~1 Gyr older and alpha-enhanced by ~0.02 relative to their counterparts of the same velocity dispersion in less dense regions, but the metallicities show no dependence on environment. This suggests that, in dense regions, the stars in early-type galaxies formed at slightly earlier times, and on a slightly shorter timescale, than in less dense regions. In addition, objects at lower redshifts are older, but have the same metallicities and alpha-enhancements as their counterparts of the same sigma at higher redshifts, as expected if the low redshift sample is a passively-aged version of the sample at higher redshifts. These significant conclusions are derived from several pieces of evidence discussed in this paper, including the FP of early-types, the observed evolution of absorption linestrengths (measured from a large ensemble of composite spectra which we provide in tabular form) over small lookback times, and comparison with the predictions of single burst (non-solar abundance) stellar population synthesis models (which also indicate that age, metallicity and alpha-enhancement all increase with velocity dispersion).

The two-point correlation function has been the standard statistic for quantifying how galaxies are clustered. The statistic uses the positions of galaxies, but not their properties. Clustering as a function of galaxy property, be it type, luminosity, color, etc., is usually studied by analysing a subset of the full population, the galaxies in the subset chosen because they have a similar range of properties. We explore an alternative technique---marked correlations---in which one weights galaxies by some property or `mark' when measuring clustering statistics. Marked correlations are particularly well-suited to quantifying how the properties of galaxies correlate with their environment. Therefore, measurements of marked statistics, with luminosity, stellar mass, color, star-formation rate, etc. as the mark, permit sensitive tests of galaxy formation models. We make measurements of such marked statistics in semi-analytic galaxy formation models to illustrate their utility. These measurements show that close pairs of galaxies are expected to be red, to have larger stellar masses, and to have smaller star formation rates. We also show that the simplest unbiased estimator of the particular marked statistic we use extensively is very simple to measure---it does not require construction of a random catalog---and provide an estimate of its variance. Large wide-field surveys of the sky are revolutionizing our view of galaxies and how they evolve. Our results indicate that application of marked statistics to this high quantity of high-quality data will provide a wealth of information about galaxy formation.

Marked statistics allow sensitive tests of how galaxy properties correlate with environment, as well as of how correlations between galaxy properties are affected by environment. A halo-model description of marked correlations is developed, which incorporates the effects which arise from the facts that typical galaxy marks (e.g., luminosity, color, star formation rate, stellar mass) depend on the mass of the parent halo, and that massive haloes extend to larger radii and populate denser regions. Comparison with measured marked statistics in semi-analytic galaxy formation models shows good agreement on scales smaller than a Megaparsec, and excellent agreement on larger scales. The halo-model description shows clearly that the behaviour of some low-order marked statistics on these scales encodes information about the mean galaxy mark as a function of halo mass, but is insensitive to mark-gradients within haloes. Higher-order statistics encode information about higher order moments of the distribution of marks within haloes. This information is obtained without ever having to identify haloes or clusters in the galaxy distribution. On scales smaller than a Megaparsec, the halo-model calculation shows that marked statistics allow sensitive tests of whether or not central galaxies in haloes are a special population. A prescription for including more general mark-gradients in the halo-model description is also provided. The formalism developed here is particularly well-suited to interpretation of marked statistics in astrophysical datasets, because it is phrased in the same language that is currently used to interpret more standard measures of galaxy clustering.

We describe the results of a search for galaxies with large (> 350 km/s) velocity dispersions. The largest systems we have found appear to be the extremes of the early-type galaxy population: compared to other galaxies with similar luminosities, they have the largest velocity dispersions and the smallest sizes. However, they are not distant outliers from the Fundamental Plane and mass-to-light scaling relations defined by the bulk of the early-type galaxy population. They may host the most massive black holes in the Universe, and their abundance and properties can be used to constrain galaxy formation models. Clear outliers from the scaling relations tend to be objects in superposition (angular separations smaller than 1 arcsec), evidence for which comes sometimes from the spectra, sometimes from the images, and sometimes from both. The statistical properties of the superposed pairs, e.g., the distribution of pair separations and velocity dispersions, can be used to provide useful information about the expected distribution of image multiplicities, separations and flux ratios due to gravitational lensing by multiple lenses, and may also constrain models of their interaction rates.

(Abridged) We use the triaxial halo model of Smith & Watts to investigate the impact of dark matter halo shapes on the bispectrum. Analytic expressions are derived and evaluated through numerical integration. Two models for the ellipsoidal halo profiles are considered: a toy model designed to isolate the effects of halo shape on the clustering alone; and the realistic model of Jing & Suto. For equilateral k-space triangles, we show that the predictions of the triaxial model are suppressed, relative to the spherical model, by roughly ~7% and ~4% for the two profiles respectively. The reduced bispectrum as a function of triangle configuration is found to be very sensitive to halo shapes on small scales. The generic features are that, relative to the spherical model, the signal is suppressed for k-space triangles that are close to equilateral and boosted for configurations that are colinear. This appears to be a unique signature of halo triaxiality. The galaxy bispectrum is also explored. Two models for the halo occupation distribution (HOD) are considered: the binomial distribution of Scoccimarro et al. and the Poisson satellite model of Kravtsov et al. Our predictions show that the galaxy bispectrum is also sensitive to halo shapes, although relative to the mass, the effects are reduced. The HOD of Kravtsov et al. is found to be more sensitive. This owes to the mean number of galaxies per halo being a steeper function of mass in this model, hence giving stronger weight to the high mass (more triaxial) haloes. Interestingly, the functional form is not strongly sensitive to the details of the HOD, but is mainly determined by the halo shape. However, a combination of measurements made on different scales and for different k-space triangles are sensitive to both halo shape and the HOD.

In hierarchical models, density fluctuations on different scales are correlated. This induces correlations between dark halo masses, their formation histories, and their larger-scale environments. In turn, this produces a correlation between galaxy properties and environment. This correlation is entirely statistical in nature. We show how the observed clustering of galaxies can be used to quantify the importance of this statistical correlation relative to other physical effects which may also give rise to correlations between the properties of galaxies and their surroundings. We also develop a halo model description of this environmental dependence of clustering.

We measure the projected correlation function w_p(r_p) from the Sloan Digital Sky Survey for a flux-limited sample of 118,000 galaxies and for a volume limited subset of 22,000 galaxies with absolute magnitude M_r<-21. Both correlation functions show subtle but systematic departures from the best-fit power law, in particular a change in slope at r_p~1-2 Mpc/h. These departures are stronger for the volume-limited sample, which is restricted to relatively luminous galaxies. We show that the inflection point in w_p(r_p) can be naturally explained by contemporary models of galaxy clustering, according to which it marks the transition from a large scale regime dominated by galaxy pairs in separate dark matter halos to a small scale regime dominated by galaxy pairs in the same dark matter halo. For example, given the dark halo population predicted by an inflationary cold dark matter scenario, the projected correlation function of the volume-limited sample can be well reproduced by a model in which the mean number of M_r<-21 galaxies in a halo of mass M>M_1=4.74 X 10^{13}\msun/h is <N>_M=(M/M_1)^{0.89}, with 75% of the galaxies residing in less massive, single-galaxy halos. This physically motivated model has the same number of free parameters as a power law, and it fits the w_p(r_p) data better, with a \chi^2/d.o.f.=0.93 compared to 6.12 (for 10 degrees of freedom, incorporating the covariance of the correlation function errors). Departures from a power-law correlation function encode information about the relation between galaxies and dark matter halos. Higher precision measurements of these departures for multiple classes of galaxies will constrain galaxy bias and provide new tests of the theory of galaxy formation.

We study the luminosity and color dependence of the galaxy 2-point correlation function in the Sloan Digital Sky Survey, starting from a sample of 200,000 galaxies over 2500 deg^2. We concentrate on the projected correlation function w(r_p), which is directly related to the real space \xi(r). The amplitude of w(r_p) grows continuously with luminosity, rising more steeply above the characteristic luminosity L_*. Redder galaxies exhibit a higher amplitude and steeper correlation function at all luminosities. The correlation amplitude of blue galaxies increases continuously with luminosity, but the luminosity dependence for red galaxies is less regular, with bright red galaxies more strongly clustered at large scales and faint red galaxies more strongly clustered at small scales. We interpret these results using halo occupation distribution (HOD) models assuming concordance cosmological parameters. For most samples, an HOD model with two adjustable parameters fits the w(r_p) data better than a power-law, explaining inflections at r_p ~ 1-3 Mpc/h as the transition between the 1-halo and 2-halo regimes of \xi(r). The implied minimum mass for a halo hosting a central galaxy above a luminosity threshold L grows as M_min ~ L at low luminosities and more steeply above L_*. The mass at which an average halo has one satellite galaxy brighter than L is M_1 ~ 23 M_min(L). These results imply a conditional luminosity function (at fixed halo mass) in which central galaxies lie far above a Schechter function extrapolation of the satellite population. HOD models nicely explain the joint luminosity-color dependence of w(r_p) in terms of the color fractions of central and satellite populations as a function of halo mass. The inferred HOD properties are in good qualitative agreement with theoretical predictions.

We report the first result of the clustering analysis of Sloan Digital Sky Survey (SDSS) quasars. We compute the two-point correlation function (2PCF) of SDSS quasars in redshift space at $8h^{-1}{\rm Mpc} < s < 500h^{-1}{\rm Mpc}$, with particular attention to its baryonic signature. Our sample consists of 19986 quasars extracted from the SDSS Data Release 4 (DR4). The redshift range of the sample is $0.72 \le z \le 2.24$ (the mean redshift is $\bar z = 1.46$) and the reddening-corrected $i$-band apparent magnitude range is $15.0 \le m_{i,{\rm rc}} \le 19.1$. Due to the relatively low number density of the quasar sample, the bump in the power spectrum due to the baryon density, $\Omega_{\rm b}$, is not clearly visible. The effect of the baryon density is, however, to distort the overall shape of the 2PCF.The degree of distortion makes it an interesting alternate measure of the baryonic signature. Assuming a scale-independent linear bias and the spatially flat universe, i.e., $\Omega_{\rm b} + \Omega_{\rm d} + \Omega_\Lambda =1$, where $\Omega_{\rm d}$ and $\Omega_\Lambda$ denote the density parameters of dark matter and the cosmological constant, we combine the observed quasar 2PCF and the predicted matter 2PCF to put constraints on $\Omega_{\rm b}$ and $\Omega_\Lambda$. Our result is fitted as $0.80- 2.8\Omega_{\rm b} < \Omega_\Lambda < 0.90 - 1.4\Omega_{\rm b}$ at the 2$\sigma$ confidence level, which is consistent with results from other cosmological observations such as WMAP. (abridged)

We present an 8 sigma detection of cosmic magnification measured by the variation of quasar density due to gravitational lensing by foreground large scale structure. To make this measurement we used 3800 square degrees of photometric observations from the Sloan Digital Sky Survey (SDSS) containing \~200,000 quasars and 13 million galaxies. Our measurement of the galaxy-quasar cross-correlation function exhibits the amplitude, angular dependence and change in sign as a function of the slope of the observed quasar number counts that is expected from magnification bias due to weak gravitational lensing. We show that observational uncertainties (stellar contamination, Galactic dust extinction, seeing variations and errors in the photometric redshifts) are well controlled and do not significantly affect the lensing signal. By weighting the quasars with the number count slope, we combine the cross-correlation of quasars for our full magnitude range and detect the lensing signal at >4 sigma in all five SDSS filters. Our measurements of cosmic magnification probe scales ranging from 60 kpc/h to 10 Mpc/h and are in good agreement with theoretical predictions based on the WMAP concordance cosmology. As with galaxy-galaxy lensing, future measurements of cosmic magnification will provide useful constraints on the galaxy-mass power spectrum.

We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular $\Lambda$ Cold Dark Matter cosmogony. We identify voids by requiring them to be regions of space with a mean overdensity of -0.8 or less. Each of the simulations contains thousands of voids. The distribution of void sizes in the different simulations shows good agreement. Voids very clearly correspond to minima in the smoothed initial density field. We find a universal void mass profile of the form $\rho(<r)/\rho(r_{\rm eff}) \propto \exp[(r/r_{\rm eff})^\alpha]$ where $r_{\rm eff}$ is the effective radius of a void and $\alpha\sim 2$. The mass function of haloes in voids is steeper than that of haloes that populate denser regions. The abundances of void haloes seem to evolve somewhat more strongly between redshifts $\sim 1$ and 0 than the global abundances of haloes. (abridged)

We present the "C4 Cluster Catalog", a new sample of 748 clusters of galaxies identified in the spectroscopic sample of the Second Data Release (DR2) of the Sloan Digital Sky Survey (SDSS). The C4 cluster--finding algorithm identifies clusters as overdensities in a seven-dimensional position and color space, thus minimizing projection effects which plagued previous optical clusters selection. The present C4 catalog covers ~2600 square degrees of sky with groups containing 10 members to massive clusters having over 200 cluster members with redshifts. We provide cluster properties like sky location, mean redshift, galaxy membership, summed r--band optical luminosity (L_r), velocity dispersion, and measures of substructure. We use new mock galaxy catalogs to investigate the sensitivity to the various algorithm parameters, as well as to quantify purity and completeness. These mock catalogs indicate that the C4 catalog is ~90% complete and 95% pure above M_200 = 1x10^14 solar masses and within 0.03 <=z <= 0.12. The C4 algorithm finds 98% of X-ray identified clusters and 90% of Abell clusters within 0.03 <= z <= 0.12. We show that the L_r of a cluster is a more robust estimator of the halo mass (M_200) than the line-of-sight velocity dispersion or the richness of the cluster. L_r. The final SDSS data will provide ~2500 C4 clusters and will represent one of the largest and most homogeneous samples of local clusters.

We combine the Cosmic Lens All-Sky Survey (CLASS) with new Sloan Digital Sky Survey (SDSS) data on the local velocity dispersion distribution function of E/S0 galaxies, $\phi(\sigma)$, to derive lens statistics constraints on $\Omega_\Lambda$ and $\Omega_m$. Previous studies of this kind relied on a combination of the E/S0 galaxy luminosity function and the Faber-Jackson relation to characterize the lens galaxy population. However, ignoring dispersion in the Faber-Jackson relation leads to a biased estimate of $\phi(\sigma)$ and therefore biased and overconfident constraints on the cosmological parameters. The measured velocity dispersion function from a large sample of E/S0 galaxies provides a more reliable method for probing cosmology with strong lens statistics. Our new constraints are in good agreement with recent results from the redshift-magnitude relation of Type Ia supernovae. Adopting the traditional assumption that the E/S0 velocity function is constant in comoving units, we find a maximum likelihood estimate of $\Omega_\Lambda = 0.74$--0.78 for a spatially flat unvierse (where the range reflects uncertainty in the number of E/S0 lenses in the CLASS sample), and a 95% confidence upper bound of $\Omega_\Lambda<0.86$. If $\phi(\sigma)$ instead evolves in accord with extended Press-Schechter theory, then the maximum likelihood estimate for $\Omega_\Lambda$ becomes 0.72--0.78, with the 95% confidence upper bound $\Omega_\Lambda<0.89$. Even without assuming flatness, lensing provides independent confirmation of the evidence from Type Ia supernovae for a nonzero dark energy component in the universe.

We present an analysis of the color-magnitude-velocity dispersion relation for a sample of 39320 early-type galaxies within the Sloan Digital Sky Survey. We demonstrate that the color-magnitude relation is entirely a consequence of the fact that both the luminosities and colors of these galaxies are correlated with stellar velocity dispersions. Previous studies of the color-magnitude relation over a range of redshifts suggest that the luminosity of an early-type galaxy is an indicator of its metallicity, whereas residuals in color from the relation are indicators of the luminosity-weighted age of its stars. We show that this, when combined with our finding that velocity dispersion plays a crucial role, has a number of interesting implications. First, galaxies with large velocity dispersions tend to be older (i.e., they scatter redward of the color-magnitude relation). Similarly, galaxies with large dynamical mass estimates also tend to be older. In addition, at fixed luminosity, galaxies which are smaller, or have larger velocity dispersions, or are more massive, tend to be older. Second, models in which galaxies with the largest velocity dispersions are also the most metal poor are difficult to reconcile with our data. However, at fixed velocity dispersion, galaxies have a range of ages and metallicities: the older galaxies have smaller metallicities, and vice-versa. Finally, a plot of velocity dispersion versus luminosity can be used as an age indicator: lines of constant age run parallel to the correlation between velocity dispersion and luminosity.

We present measurements of the redshift--space three-point correlation function of galaxies in the Sloan Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy morphology, color and luminosity. In order to control systematics due to selection effects, we used $r$--band, volume-limited samples of galaxies, constructed from the magnitude-limited SDSS data ($14.5<r<17.5$), and further divided the samples into two morphological types (early and late) or two color populations (red and blue). The three-point correlation function of SDSS galaxies follow the hierarchical relation well and the reduced three-point amplitudes in redshift--space are almost scale-independent ($Q_z=0.5\sim1.0$). In addition, their dependence on the morphology, color and luminosity is not statistically significant. Given the robust morphological, color and luminosity dependences of the two-point correlation function, this implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that simple deterministic linear relation with the underlying mass could not explain our measurements on these scales.

According to the current paradigm, galaxies form and reside in extended cold dark matter (CDM) halos and in turn are key tracers of cosmological structure. Understanding how different types of galaxies occupy halos of different masses is one of the major challenges facing extragalactic astrophysics. The observed galaxy properties depend on the environment surrounding the galaxy. Within the framework of most galaxy formation models the environmental dependence of the galaxy population is mainly due to the change of the halo mass function with large-scale environment. Such models make precise predictions for how galaxy clustering should depend on environment. We will illustrate this by presenting analytical models of dark matter and galaxy clustering along with results obtained from numerical simulations. With these results we can hope to obtain a better understanding of the link between galaxies and dark matter and thereby constrain galaxy formation models.

We present a model for the distribution of void sizes and its evolution within the context of hierarchical scenarios of gravitational structure formation. For a proper description of the hierarchical buildup of the system of voids in the matter distribution, not only the "void-in-void" problem should be taken into account, but also that of the "void-in-cloud" issue. Within the context of the excursion set formulation of an evolving void hierarchy is one involving a "two-barrier" excursion problem, unlike the "one-barrier" problem for the dark halo evolution. This leads to voids having a peaked size distribution at any cosmic epoch, centered on a characteristic void size that evolves self-similarly in time, this in distinct contrast to the distribution of virialized halo masses which do not have a small-scale cut-off.

A generic prediction of hierarchical gravitational clustering models is that the distribution of halo formation times should depend relatively strongly on halo mass, massive haloes forming more recently, and depend only weakly, if at all, on the large scale environment of the haloes. We present a novel test of this assumption which uses a statistic that proves to be particularly well-suited to detecting and quantifying weak correlations with environment. We find that close pairs of haloes form at slightly higher redshifts than do more widely separated halo pairs, suggesting that haloes in dense regions form at slightly earlier times than do haloes of the same mass in less dense regions. The environmental trends we find are useful for models which relate the properties of galaxies to the formation histories of the haloes which surround them.

The most commonly used definition of halo formation is the time when a halo's most massive progenitor first contains at least half the final mass of its parent. Reasonably accurate formulae for the distribution of formation times of haloes of fixed mass have been available for some time. We use numerical simulations of hierarchical gravitational clustering to test the accuracy of formulae for the mass at formation. We also derive and test a formula for the joint distribution of formation masses and times. The structure of a halo is expected to be related to its accretion history. Our tests show that our formulae for formation masses and times are reasonably accurate, so we expect that they will aid future analytic studies of halo structure.

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.

We present a model for the distribution of void sizes and its evolution in the context of hierarchical scenarios of gravitational structure formation. We find that at any cosmic epoch the voids have a size distribution which is well-peaked about a characteristic void size which evolves self-similarly in time. This is in distinct contrast to the distribution of virialized halo masses which does not have a small-scale cut-off. In our model, the fate of voids is ruled by two processes. The first process affects those voids which are embedded in larger underdense regions: the evolution is effectively one in which a larger void is made up by the mergers of smaller voids, and is analogous to how massive clusters form from the mergers of less massive progenitors. The second process is unique to voids, and occurs to voids which happen to be embedded within a larger scale overdensity: these voids get squeezed out of existence as the overdensity collapses around them. It is this second process which produces the cut-off at small scales. In the excursion set formulation of cluster abundance and evolution, solution of the cloud-in-cloud problem, i.e., counting as clusters only those objects which are not embedded in larger clusters, requires study of random walks crossing one barrier. We show that a similar formulation of void evolution requires study of a two-barrier problem: one barrier is required to account for voids-in-voids, and the other for voids-in-clouds. Thus, in our model, the void size distribution is a function of two parameters, one of which reflects the dynamics of void formation, and the other the formation of collapsed objects.

Within the context of hierarchical scenarios of gravitational structure formation we describe how an evolving hierarchy of voids evolves on the basis of {\em two} processes, the {\em void-in-void} process and the {\it void-in-cloud} process. The related analytical formulation in terms of a {\it two-barrier} excursion problem leads to a self-similarly evolving peaked void size distribution.

We measure the large-scale real-space power spectrum P(k) using a sample of 205,443 galaxies from the Sloan Digital Sky Survey, covering 2417 square degrees with mean redshift z~0.1. We employ a matrix-based method using pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 22 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.02 h/Mpc < k < 0.3h/Mpc. We pay particular attention to modeling, quantifying and correcting for potential systematic errors, nonlinear redshift distortions and the artificial red-tilt caused by luminosity-dependent bias. Our final result is a measurement of the real-space matter power spectrum P(k) up to an unknown overall multiplicative bias factor. Our calculations suggest that this bias factor is independent of scale to better than a few percent for k<0.1h/Mpc, thereby making our results useful for precision measurements of cosmological parameters in conjunction with data from other experiments such as the WMAP satellite. As a simple characterization of the data, our measurements are well fit by a flat scale-invariant adiabatic cosmological model with h Omega_m =0.201+/- 0.017 and L* galaxy sigma_8=0.89 +/- 0.02 when fixing the baryon fraction Omega_b/Omega_m=0.17 and the Hubble parameter h=0.72; cosmological interpretation is given in a companion paper.

Several measurements of QSO-galaxy correlations have reported signals much larger than predictions of magnification by large-scale structure. We find that the expected signal depends stronly on the properties of the foreground galaxy population. On arcminute scales it can be either larger or smaller by a factor of two for different galaxy types in comparison with a linearly biased version of the mass distribution. Thus the resolution of some of the excess measurements may lie in examining the halo occupation properties of the galaxy population sampled by a given survey; this is also the primary information such measurements will provide. We use the halo model of clustering and simulations to predict the magnification induced cross-correlations and errors for forthcoming surveys. With the full Sloan Digital Sky Survey the statistical errors will be below 1 percent for the galaxy-galaxy correlations and significantly larger for QSO-galaxy correlations. Thus accurate constraints on parameters of the galaxy halo occupation distribution can be obtained from small scale measurements and on the bias parameter from large scales. Since the lensing induced cross-correlation measures the first moment of the halo occupation number of galaxies, these measurements can provide the basis for interpreting galaxy clustering measurements which measure the second and higher order moments.

In previous work on galaxy clusters, several authors reported a discovery of an unusual population of galaxies, which have spiral morphologies, but do not show any star formation activity. These galaxies are called ``passive spirals'', and have been interesting since it has been difficult to understand the existence of such galaxies. Using a volume limited sample (0.05<z<0.1 and Mr<-20.5; 25813 galaxies) of the Sloan Digital Sky Survey data, we have found 73 (0.28$\pm$0.03%) passive spiral galaxies and studied their environments. It is found that passive spiral galaxies live in local galaxy density 1-2 Mpc$^{-2}$ and 1-10 cluster-centric virial radius. Thus the origins of passive spiral galaxies are likely to be cluster related. These characteristic environments coincide with the previously reported environment where galaxy star formation rate suddenly declines and the so-called morphology-density relation turns. It is likely that the same physical mechanism is responsible for all of these observational results. The existence of passive spiral galaxies suggests that a physical mechanism that works calmly is preferred to dynamical origins such as major merger/interaction since such a mechanism can destroy spiral arm structures. Compared with observed cluster galaxy evolution such as the Butcher-Oemler effect and the morphological Butcher-Oemler effect, passive spiral galaxies are likely to be a galaxy population in transition between red, elliptical/S0 galaxies in low redshift clusters and blue, spiral galaxies more numerous in higher redshift clusters.

We present measurements of the angular cross-correlation between luminous red galaxies from the Sloan Digital Sky Survey and the cosmic microwave background temperature maps from the Wilkinson Microwave Anisotropy Probe. We find a statistically significant achromatic positive correlation between these two data sets, which is consistent with the expected signal from the late Integrated Sachs-Wolfe (ISW) effect. We do not detect any anti-correlation on small angular scales as would be produced from a large Sunyaev-Zel'dovich (SZ) effect, although we do see evidence for some SZ effect for our highest redshift samples. Assuming a flat universe, our preliminary detection of the ISW effect provides independent physical evidence for the existence of dark energy.

The distribution of early-type galaxy velocity dispersions, phi(sigma), is measured using a sample drawn from the SDSS database. Its shape differs significantly from that which one obtains by simply using the mean correlation between luminosity, L, and velocity dispersion, sigma, to transform the luminosity function into a velocity function: ignoring the scatter around the mean sigma-L relation is a bad approximation. An estimate of the contribution from late-type galaxies is also made, which suggests that phi(sigma) is dominated by early-type galaxies at velocities larger than ~ 200 km/s.

We investigate the evolution of the fractions of late type cluster galaxies as a function of redshift, using one of the largest, most uniform cluster samples available. The sample consists of 514 clusters of galaxies in the range 0.02<z<0.3 from the Sloan Digital Sky Survey Cut & Enhance galaxy cluster catalog. This catalog was created using a single automated cluster finding algorithm on uniform data from a single telescope, with accurate CCD photometry, thus, minimizing selection biases. We use four independent methods to analyze the evolution of the late type galaxy fraction. Specifically, we select late type galaxies based on: restframe g-r color, u-r color, galaxy profile fitting and concentration index. The first criterion corresponds to the one used in the classical Butcher-Oemler analyses. The last three criteria are more sensitive to the morphological type of the galaxies. In all four cases, we find an increase in the fraction of late type galaxies with increasing redshift, significant at the 99.9% level. The results confirm that cluster galaxies do change colors with redshift (the Butcher-Oemler effect) and, in addition, they change their morphology to later-type toward higher redshift -- indicating a morphological equivalent of the Butcher-Oemler effect. We also find a tendency of richer clusters to have lower fractions of late type galaxies. The trend is consistent with a ram pressure stripping model, where richer clusters have more effective ram pressure due to their higher temperature.

We present here results on the composite luminosity functions of galaxies in the clusters of galaxies selected from the SDSS Cut and Enhance cluster catalog (CE; Goto et al. 2001). We construct the composite luminosity function in the five SDSS bands, u,g,r,i and z, using 204 CE clusters ranging from z=0.02 to z=0.25. We use photometric redshifts to construct composite luminosity functions. The robustness of the weighting scheme was tested using Monte Carlo simulation. We find the slope of composite LFs become flatter toward redder color band. Comparing with field LFs of SDSS (Blanton et al. 2001), cluster LFs have brighter characteristic magnitude and flatter slopes in g, r, i and z bands. These results are consistent with the hypothesis that the cluster LF has two distinct underlying populations i.e. the bright end of the LF is dominated by bright early types that follow a gaussian--like luminosity distribution, while the faint--end of the cluster LF is a steep power--law like function dominated by star--forming (bluer) galaxies. We also study the composite luminosity functions for early--type and late--type galaxies using profile fits, a concentration parameter and u-r color to classify galaxy morphology. The strong dependence of LF on galaxy morphology is found.The faint end slope of the LF is always flatter for early--type galaxies than late--type regardless of passband and methodology. These results are consistent with the hypothesis that the cluster region are dominated by bright elliptical galaxies.

I develop a model for the abundance and spatial distribution of dark matter subclumps. The model shows that subclumps of massive parent halos formed at earlier times than subclumps of the same mass in lower mass parents; equivalently, halos in dense regions at a given time formed earlier than halos of the same mass in less dense regions. This may provide the basis for interpreting recent observations which indicate that the stellar populations of the most massive elliptical galaxies are also the oldest.

A magnitude limited sample of nearly 9000 early-type galaxies, in the redshift range 0.01 < z < 0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The sample was used to study how early-type galaxy observables, including luminosity L, effective radius R_o, surface brightness I_o, color, and velocity dispersion sigma, are correlated with one another. Measurement biases are understood with mock catalogs which reproduce all of the observed scaling relations and their dependences on fitting technique. At any given redshift, the intrinsic distribution of luminosities, sizes and velocity dispersions in our sample are all approximately Gaussian. A maximum likelihood analysis shows that sigma ~ L^{0.25\pm 0.012}, R_o ~ L^{0.63\pm 0.025}, and R_o ~ I^{-0.75\pm 0.02} in the r* band. In addition, the mass-to-light ratio within the effective radius scales as M_o/L ~ L^{0.14\pm 0.02} or M_o/L ~ M_o^{0.22\pm 0.05}, and galaxies with larger effective masses have smaller effective densities: Delta_o ~ M_o^{-0.52\pm 0.03}. These relations are approximately the same in the g*, i* and z* bands. Relative to the population at the median redshift in the sample, galaxies at lower and higher redshifts have evolved only little, with more evolution in the bluer bands. The luminosity function is consistent with weak passive luminosity evolution and a formation time of about 9 Gyrs ago.

A sample of nearly 9000 early-type galaxies, in the redshift range 0.01 < z < 0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. This paper describes how the sample was selected, presents examples of images and seeing corrected fits to the observed surface brightness profiles, describes our method for estimating K-corrections, and shows that the SDSS spectra are of sufficiently high quality to measure velocity dispersions accurately. It also provides catalogs of the measured photometric and spectroscopic parameters. In related papers, these data are used to study how early-type galaxy observables, including luminosity, effective radius, surface brightness, color, and velocity dispersion, are correlated with one another.

The colors and chemical abundances of early-type galaxies at redshifts z<0.3 are studied using a sample of nearly 9000 galaxies, selected from the Sloan Digital Sky Survey using morphological and spectral criteria. In this sample, redder galaxies have larger velocity dispersions. Color also correlates with magnitude and size, but these correlations are entirely due to the L-sigma and R_o-sigma relations: the primary correlation is color-sigma. The red light in early-type galaxies is, on average, slightly more centrally concentrated than the blue. Because of these color gradients, the strength of the color--magnitude relation depends on whether or not the colors are defined using a fixed metric aperture. Chemical evolution and star formation histories of early-type galaxies are investigated using co-added spectra of similar objects. At fixed sigma, the population at z ~ 0.2 had weaker Mg_2 and stronger H_beta absorption compared to the population at z ~ 0. It was also bluer. Comparison of these colors and line-strengths, and their evolution, with single-burst stellar population models suggests a formation time of 9 Gyrs ago, consistent with a Fundamental Plane analysis of this sample.

A magnitude limited sample of nearly 9000 early-type galaxies, in the redshift range 0.01 < z < 0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The Fundamental Plane relation in this sample is R_o ~ sigma^{1.49\pm 0.05} I_o^{-0.75\pm 0.01} in the r* band. It is approximately the same in the g*, i* and z* bands. Relative to the population at the median redshift in the sample, galaxies at lower and higher redshifts have evolved only little. If the Fundamental Plane is used to quantify this evolution then the apparent magnitude limit can masquerade as evolution; once this selection effect has been accounted for, the evolution is consistent with that of a passively evolving population which formed the bulk of its stars about 9 Gyrs ago. One of the principal advangtages of the SDSS sample over previous samples is that the galaxies in it lie in environments ranging from isolation in the field to the dense cores of clusters. The Fundamental Plane shows that galaxies in dense regions are slightly different from galaxies in less dense regions.

(Abridged) We present in this paper a detailed analysis of the effect of environment on the star-formation activity of galaxies within the EDR of the SDSS. We have used the Halpha emission line to derive the star-formation rate (SFR) for each galaxy within a volume-limited sample of 8598 galaxies with 0.05 < z < 0.095 and M(r)<= -20.45. We find that the SFR of galaxies is strongly correlated with the local (projected) galaxy density and thus we present here the density-SFR relation that is analogous to the density-morphology relation. The effect of density on the SFR of galaxies is seen in three ways. First, the overall distribution of SFRs is shifted to lower values in dense environments compared with the field population. Second, the effect is most noticeable for the strongly star-forming galaxies in the 75th percentile of the SFR distribution. Third, there is a ``break'' (or characteristic density) in the density-SFR relation at a local galaxy density of 1h-2 Mpc-2. To understand this break further, we have studied the SFR of galaxies as a function of clustercentric radius from 17 clusters and groups objectively selected from the SDSS EDR data. The distribution of SFRs of cluster galaxies begins to change, compared with the field population, at a clustercentric radius of 3-4 virial radii, which is consistent with the characteristic break in density that we observe in the density-SFR relation. Our tests suggest that the density-morphology relation alone is unlikely to explain the density-SFR relation we observe. Taken all together, these works demonstrate that the decrease in SFR of galaxies in dense environments is a universal phenomenon over a wide range in density (from 0.08 to 10h-2 Mpc-2) and redshift (out to z = 0.5).

The effective optical depth in the Ly-alpha forest region of 1061 low-resolution QSO spectra drawn from the SDSS database decreases with decreasing redshift over the range 2.5 < z < 4. Although the evolution is relatively smooth, tau_eff ~ (1+z)^{3.8 pm 0.2}, at z ~ 3.2 the effective optical depth decreases suddenly, by about ten percent with respect to this smoother evolution. It climbs back to the original smooth scaling again by z ~ 2.9. We describe two techniques, one of which is new, for quantifying this evolution which give consistent results. A variety of tests show that the feature is not likely to be a consequence of how the QSO sample was selected, nor the result of flux calibration or other systematic effects. Other authors have argued that, at this same epoch, the temperature of the IGM also shows a departure from an otherwise smooth decrease with time. These features in the evolution of the temperature and the optical depth are signatures of the reionization of He II.

We develop the formalism to include substructure in the halo model of clustering. Real halos are not likely to be perfectly smooth, but have substructure which has so far been neglected in the halo model -- our formalism allows one to estimate the effects of this substructure on measures of clustering. We derive expressions for the two-point correlation function, the power-spectrum, the cross-correlation between galaxies and mass, as well as higher order correlation functions. Simple forms of the formulae are obtained for the limit in which the size of the subclumps and mass fraction in them is small. Inclusion of substructure allows for a more accurate analysis of the statistical effects of gravitational lensing. It can also bring the halo model predictions into better agreement with the small-scale structure seen in recent high resolution simulations of hierarchical clustering.

The angular distribution of galaxies encodes a wealth of information about large scale structure. Ultimately, the Sloan Digital Sky Survey (SDSS) will record the angular positions of order 10^8 galaxies in five bands, adding significantly to the cosmological constraints. This is the first in a series of papers analyzing a rectangular stripe 2.5x90 degrees from early SDSS data. We present the angular correlation function for galaxies in four separate magnitude bins on angular scales ranging from 0.003 degrees to 15 degrees. Much of the focus of this paper is on potential systematic effects. We show that the final galaxy catalog -- with the mask accounting for regions of poor seeing, reddening, bright stars, etc. -- is free from external and internal systematic effects for galaxies brighter than r* = 22. Our estimator of the angular correlation function includes the effects of the integral constraint and the mask. The full covariance matrix of errors in these estimates is derived using mock catalogs with further estimates using a number of other methods.

We review the formalism and applications of the halo-based description of nonlinear gravitational clustering. In this approach, all mass is associated with virialized dark matter halos; models of the number and spatial distribution of the halos, and the distribution of dark matter within each halo, are used to provide estimates of how the statistical properties of large scale density and velocity fields evolve as a result of nonlinear gravitational clustering. We first describe the model, and demonstrate its accuracy by comparing its predictions with exact results from numerical simulations of nonlinear gravitational clustering. We then present several astrophysical applications of the halo model: these include models of the spatial distribution of galaxies, the nonlinear velocity, momentum and pressure fields, descriptions of weak gravitational lensing, and estimates of secondary contributions to temperature fluctuations in the cosmic microwave background.

The expansion of the Universe leads to a rapid drop in the mean hydrogen Ly-alpha opacity tau of the intergalactic medium (IGM), tau (1+z)^{3.8}, between redshifts 4 and 3. Measurements of the temperature evolution of the IGM and of the He II opacity both suggest that He II reionizes in this redshift range. We use hydrodynamical simulations to show that the temperature increase associated with He II reionization leads to a relatively sudden decrease in tau around the reionization epoch of ~ 10 per cent. We find clear evidence for such a feature in the evolution of tau determined from a sample of ~ 1100 quasars obtained from the SDSS. He II reionization starts at redshift ~ 3.4, and lasts for (Delta z) ~ 0.4. The increase in the IGM temperature also explains the widths of hydrogen absorption lines as measured in high-resolution spectra.

The stochasticity in the distribution of dark haloes in the cosmic density field is reflected in the distribution function $P_V(N_h|\delta_m)$ which gives the probability of finding $N_h$ haloes in a volume $V$ with mass density contrast $\delta_m$. We study the properties of this function using high-resolution $N$-body simulations, and find that $P_V(N_n|\delta_m)$ is significantly non-Poisson. The ratio between the variance and the mean goes from $\sim 1$ (Poisson) at $1+\delta_m\ll 1$ to $<1$ (sub-Poisson) at $1+\delta_m\sim 1$ to $>1$ (super-Poisson) at $1+\delta_m\gg 1$. The mean bias relation is found to be well described by halo bias models based on the Press-Schechter formalism. The sub-Poisson variance can be explained as a result of halo-exclusion while the super-Poisson variance at high $\delta_m$ may be explained as a result of halo clustering. A simple phenomenological model is proposed to describe the behavior of the variance as a function of $\delta_m$. Galaxy distribution in the cosmic density field predicted by semi-analytic models of galaxy formation shows similar stochastic behavior. We discuss the implications of the stochasticity in halo bias to the modelling of higher-order moments of dark haloes and of galaxies.

We compute the angular power spectrum C_l from 1.5 million galaxies in early SDSS data on large angular scales, l<600. The data set covers about 160 square degrees, with a characteristic depth of order 1 Gpc/h in the faintest (21<r<22) of our four magnitude bins. Cosmological interpretations of these results are presented in a companion paper by Dodelson et al (2001). The data in all four magnitude bins are consistent with a simple flat ``concordance'' model with nonlinear evolution and linear bias factors of order unity. Nonlinear evolution is particularly evident for the brightest galaxies. A series of tests suggest that systematic errors related to seeing, reddening, etc., are negligible, which bodes well for the sixtyfold larger sample that the SDSS is currently collecting. Uncorrelated error bars and well-behaved window functions make our measurements a convenient starting point for cosmological model fitting.

A sample of nearly 9000 early-type galaxies, in the redshift range 0.01<z<0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The sample was used to study how early-type galaxy observables, including luminosity L, effective radius R, surface brightness I, color, and velocity dispersion V, are correlated with one another. Measurement biases are understood with mock catalogs which reproduce all of the observed scaling relations. At any given redshift, the intrinsic distribution of luminosities, sizes and velocity dispersions in our sample are all approximately Gaussian. In the r* band L ~ V^3.91, L ~ R^1.58, R ~ I^(-0.75), and the Fundamental Plane relation is R ~ V^(1.49) I^(-0.75). These relations are approximately the same in the g*, i* and z* bands. At fixed luminosity, the mass-to-light ratio scales as M/L ~ L^0.14. The g*-r* color scales as V^0.25. Color also correlates with magnitude and size, but these correlations are entirely due to the L-V and R-V relations. Chemical evolution and star formation histories are investigated using co-added spectra of similar objects in our sample. Chemical abundances correlate primarily with velocity dispersion. At fixed V, the higher redshift population is bluer, is weaker in Mg2, and is stronger in Hbeta than the population nearby. In addition, the population at higher redshifts is slightly more luminous. These differences are consistent with that of a passively evolving population which formed the bulk of its stars about 9 Gyrs ago. The Fundamental Plane suggests that galaxies in dense regions are slightly different from those in less dense regions, but the co-added spectra and color--magnitude relations show no statistically significant dependence on environment.

The decay of very heavy metastable relics of the Early Universe can produce ultra-high energy cosmic rays (UHECRs) in the halo of our own Galaxy. In this model, no Greisen-Zatsepin-Kuzmin cutoff is expected because of the short propagation distances. We show here that, as a consequence of the hierarchical build up of the halo, this scenario predicts the existence of small scale anisotropies in the arrival directions of UHECRs, in addition to a large scale anisotropy, known from previous studies. We also suggest some other observable consequences of this scenario which will be testable with upcoming experiments, as Auger, EUSO and OWL.

We use numerical simulations of a (480 Mpc/h)^3 volume to show that the distribution of peak heights in maps of the temperature fluctuations from the kinematic and thermal Sunyaev-Zeldovich effects will be highly non-Gaussian, and very different from the peak height distribution of a Gaussian random field. We then show that it is a good approximation to assume that each peak in either SZ effect is associated with one and only one dark matter halo. This allows us to use our knowledge of the properties of haloes to estimate the peak height distributions. At fixed optical depth, the distribution of peak heights due to the kinematic effect is Gaussian, with a width which is approximately proportional to optical depth; the non-Gaussianity comes from summing over a range of optical depths. The optical depth is an increasing function of halo mass, and the distribution of halo speeds is Gaussian, with a dispersion which is approximately independent of halo mass. This means that observations of the kinematic effect can be used to put constraints on how the abundance of massive clusters evolves, and on the evolution of cluster velocities. The non-Gaussianity of the thermal effect, on the other hand, comes primarily from the fact that, on average, the effect is larger in more massive haloes, and the distribution of halo masses is highly non-Gaussian. We also show that because haloes of the same mass may have a range of density and velocity dispersion profiles, the relation between halo mass and the amplitude of the thermal effect is not deterministic, but has some scatter.

We present measurements of parameters of the 3-dimensional power spectrum of galaxy clustering from 222 square degrees of early imaging data in the Sloan Digital Sky Survey. The projected galaxy distribution on the sky is expanded over a set of Karhunen-Loeve eigenfunctions, which optimize the signal-to-noise ratio in our analysis. A maximum likelihood analysis is used to estimate parameters that set the shape and amplitude of the 3-dimensional power spectrum. Our best estimates are Gamma=0.188 +/- 0.04 and sigma_8L = 0.915 +/- 0.06 (statistical errors only), for a flat Universe with a cosmological constant. We demonstrate that our measurements contain signal from scales at or beyond the peak of the 3D power spectrum. We discuss how the results scale with systematic uncertainties, like the radial selection function. We find that the central values satisfy the analytically estimated scaling relation. We have also explored the effects of evolutionary corrections, various truncations of the KL basis, seeing, sample size and limiting magnitude. We find that the impact of most of these uncertainties stay within the 2-sigma uncertainties of our fiducial result.

The Sloan Digital Sky Survey is one of the first multicolor photometric and spectroscopic surveys designed to measure the statistical properties of galaxies within the local Universe. In this Letter we present some of the initial results on the angular 2-point correlation function measured from the early SDSS galaxy data. The form of the correlation function, over the magnitude interval 18<r*<22, is shown to be consistent with results from existing wide-field, photographic-based surveys and narrower CCD galaxy surveys. On scales between 1 arcminute and 1 degree the correlation function is well described by a power-law with an exponent of ~ -0.7. The amplitude of the correlation function, within this angular interval, decreases with fainter magnitudes in good agreement with analyses from existing galaxy surveys. There is a characteristic break in the correlation function on scales of approximately 1-2 degrees. On small scales, < 1', the SDSS correlation function does not appear to be consistent with the power-law form fitted to the 1'< theta <0.5 deg data. With a data set that is less than 2% of the full SDSS survey area, we have obtained high precision measurements of the power-law angular correlation function on angular scales 1' < theta < 1 deg, which are robust to systematic uncertainties. Because of the limited area and the highly correlated nature of the error covariance matrix, these initial results do not yet provide a definitive characterization of departures from the power-law form at smaller and larger angles. In the near future, however, the area of the SDSS imaging survey will be sufficient to allow detailed analysis of the small and large scale regimes, measurements of higher-order correlations, and studies of angular clustering as a function of redshift and galaxy type.

Early photometric data from the Sloan Digital Sky Survey (SDSS) contain angular positions for 1.5 million galaxies. In companion papers, the angular correlation function $w(\theta)$ and 2D power spectrum $C_l$ of these galaxies are presented. Here we invert Limber's equation to extract the 3D power spectrum from the angular results. We accomplish this using an estimate of $dn/dz$, the redshift distribution of galaxies in four different magnitude slices in the SDSS photometric catalog. The resulting 3D power spectrum estimates from $w(\theta)$ and $C_l$ agree with each other and with previous estimates over a range in wavenumbers $0.03 < k/{\rm h Mpc}^{-1} < 1$. The galaxies in the faintest magnitude bin ($21 < \rstar < 22$, which have median redshift $z_m=0.43$) are less clustered than the galaxies in the brightest magnitude bin ($18 < \rstar < 19$ with $z_m=0.17$), especially on scales where nonlinearities are important. The derived power spectrum agrees with that of Szalay et al. (2001) who go directly from the raw data to a parametric estimate of the power spectrum. The strongest constraints on the shape parameter $\Gamma$ come from the faintest galaxies (in the magnitude bin $21 < \rstar < 22$), from which we infer $\Gamma = 0.14^{+0.11}_{-0.06}$ (95% C.L.).

Current models of galaxy formation applied to understanding the large-scale structure of the universe have two parts. The first is an accurate solution of the equations of motion for the dark matter due to gravitational clustering. The second consists of making physically reasonable approximations to the behavior of baryons inside dark matter halos. The first uses large, computationally intensive, $n$-body simulations. We argue that because the second step is, at least at present, uncertain, it is possible to obtain similar galaxy distributions without solving the first step exactly. We describe an algorithm which is several orders of magnitude faster than n-body simulations, but which is, nevertheless, rather accurate. The algorithm combines perturbation theory with virialized halo models of the nonlinear density and velocity fields. For two- and three-point statistics the resulting fields are exact on large scales, and rather accurate well into the nonlinear regime, particularly for two-point statistics in real and redshift space. We then show how to use this algorithm to generate mock galaxy distributions from halo occupation numbers. As a first application, we show that it provides a good description of the clustering of galaxies in the PSCz survey. We also discuss applications to the estimation of non-Gaussian contributions to error bars and covariance matrix of the power spectrum, in real and redshift space, for galaxies and dark matter. The results for the latter show good agreement with simulations, supporting the use of our method to constrain cosmological parameters from upcoming galaxy surveys.

The excursion set approach allows one to estimate the abundance and spatial distribution of virialized dark matter haloes efficiently and accurately. The predictions of this approach depend on how the nonlinear processes of collapse and virialization are modelled. We present simple analytic approximations which allow us to compare the excursion set predictions associated with spherical and ellipsoidal collapse. In particular, we present formulae for the universal unconditional mass function of bound objects and the conditional mass function which describes the mass function of the progenitors of haloes in a given mass range today. We show that the ellipsoidal collapse based moving barrier model provides a better description of what we measure in the numerical simulations than the spherical collapse based constant barrier model, although the agreement between model and simulations is better at large lookback times. Our results for the conditional mass function can be used to compute accurate approximations to the local-density mass function which quantifies the tendency for massive haloes to populate denser regions than less massive haloes. This happens because low density regions can be thought of as being collapsed haloes viewed at large lookback times, whereas high density regions are collapsed haloes viewed at small lookback times. Although we have only applied our analytic formulae to two simple barrier shapes, we show that they are, in fact, accurate for a wide variety of moving barriers. We suggest how they can be used to study the case in which the initial dark matter distribution is not completely cold.

A simple model of how objects of different masses stream towards each other as they cluster gravitationally is described. The model shows how the mean streaming velocity of dark matter particles is related to the motions of the parent dark matter haloes. It also provides a reasonably accurate description of how the pairwise velocity dispersion of dark matter particles differs from that of the parent haloes. The analysis is then extended to describe the streaming motions of galaxies. This shows explicitly that the streaming motions measured in a given galaxy sample depend on how the sample was selected, and shows how to account for this dependence on sample selection. In addition,we show that the pairwise dispersion should also depend on sample type. Our model predicts that, on small scales, redshift space distortions should affect red galaxies more strongly than blue.

We present a simple model for the shape of the distribution function of galaxy peculiar velocities. We show how both nonlinear and linear theory terms combine to produce a distribution which has an approximately Gaussian core with exponential wings. The model is easily extended to study how the statistic depends on the type of particle used to trace the velocity field (dark matter particles, dark matter haloes, galaxies), and on the density of the environment in which the test particles are. Comparisons with simulations suggest that our model is accurate. We also show that the evolution of the peculiar velocities depends on the local, rather than the global density. Since clusters populate denser regions on average, using cluster velocities with the linear theory scaling may lead to an overestimate of the global value of Omega. Conversely, using linear theory with the global value of Omega to scale cluster velocities from the initial to the present time results in an underestimate of their true velocities. In general, however, the directions of motions of haloes are rather well described by linear theory. Our results help to simplify models of redshift-space distortions considerably.

We write the correlation function of dark matter particles, xi(r), as the sum of two terms - one which accounts for nonlinear evolution, and dominates on small scales, and another which is essentially the term from linear theory, and dominates on large scales. We use models of the number and spatial distribution of haloes and halo density profiles to describe the nonlinear term and its evolution. The result provides a good description of the evolution of xi(r) in simulations. We then use this decomposition to provide simple and accurate models of how the single particle velocity dispersion evolves with time, and how the first and second moments of the pairwise velocity distribution depend on scale. The key idea is to use the simple physics of linear theory on large scales, the simple physics of the virial theorem on small scales, and our model for the correlation function to tell us how to weight the two types of contributions (linear and nonlinear) to the pairwise velocity statistics. When incorporated into the streaming model, our results will allow a simple accurate description of redshift-space distortions over the entire range of linear to highly nonlinear regimes.

We study galaxy clustering using halo models, where gravitational clustering is described in terms of dark matter halos. At small scales, clustering statistics are dominated by halo density profiles, whereas at large scales, correlations are the result of combining non-linear perturbation theory with halo biasing. Galaxies are assumed to follow the dark matter profiles, and galaxy formation efficiency is given by the number of galaxies as a function of halo mass. This approach leads to generic predictions: the galaxy power spectrum shows a power-law behavior even though the dark matter does not, and the galaxy higher-order correlations show smaller amplitudes at small scales than their dark matter counterparts, as observed in galaxy catalogs. We find that requiring to fit both the second and third order moments of the APM galaxies provides a strong constraint on galaxy formation models. The data at large scales require that galaxy formation be relatively efficient at small masses, m =10^10 Msun/h, whereas data at smaller scales require that the number of galaxies in a halo scale as the mass to the 0.8th power in the high-mass limit. These constraints are independent of those derived from the luminosity function or Tully-Fisher relation. We also predict the power spectrum, bispectrum, and higher-order moments of the mass density field. Although halo models agree well with measurements of the mass power spectrum and the higher order Sp parameters in N-body simulations, the model assumption that halos are spherical leads to disagreement in the configuration dependence of the bispectrum at small scales. We stress the importance of finite volume effects in higher-order statistics and show how they can be estimated in this approach.

The decay of very heavy metastable relics of the Early Universe can produce ultra-high energy cosmic rays (UHECRs) in the halo of our own Galaxy. On distance scales of the order of the halo size, energy losses are negligible---no Greisen-Zatsepin-Kuzmin cutoff is expected. In this letter we show that, as a consequence of the hierarchical build up of the halo, this scenario predicts the existence of small scale anisotropies in the arrival directions of UHECRs. We also suggest some consequences of this scenario which will be testable with upcoming experiments, as Auger.

The Press--Schechter, excursion set approach allows one to make predictions about the shape and evolution of the mass function of bound objects. It combines the assumption that objects collapse spherically with the assumption that the initial density fluctuations were Gaussian and small. While the predicted mass function is reasonably accurate at the high mass end, it has more low mass objects than are seen in simulations of hierarchical clustering. This discrepancy can be reduced substantially if bound structures are assumed to form from an ellipsoidal, rather than a spherical collapse. In the spherical model, a region collapses if the initial density within it exceeds a threshold value, delta_sc. This value is independent of the initial size of the region, and since the mass of the collapsed object is related to its initial size, delta_sc is independent of final mass. In the ellipsoidal model, the collapse of a region depends on the surrounding shear field, as well as on its initial overdensity. Therefore, there is a relation between the density threshold value required for collapse, and the mass of the final object. We provide a fitting function to this delta_ec(m) relation for initially Gaussian fields which simplifies the inclusion of ellipsoidal dynamics in the excursion set approach. We discuss the relation between the excursion set predictions and the halo distribution in high resolution N-body simulations, and show that our simple parametrization of the ellipsoidal collapse model represents a significant improvement on the spherical model on an object-by-object basis. Finally, we show that the associated statistical predictions, the mass function and the large scale halo-to-mass bias relation, are also more accurate than the standard predictions.

We present accurate predictions of the correlation function of hotspots in the microwave background radiation for gaussian theories such as those predicted in most inflation models. The correlation function of peaks above a certain threshold depends only on the threshold and the power spectrum of temperature fluctuations. Since there are both potentially observable quantities in a microwave background map, there are no adjustable parameters in the predictions. These correlations should therefore provide a powerful test of the gaussian hypothesis, and provide a useful discriminant between inflation and topological defect models such as the cosmic string model. The correlations have a number of oscillatory features, which should be detectable at high signal-to-noise with future satellite experiments such as MAP and Planck.

Dark matter haloes are biased tracers of the underlying dark matter distribution. We use a simple model to provide a relation between the abundance of dark matter haloes and their spatial distribution on large scales. Our model shows that knowledge of the unconditional mass function alone is sufficient to provide an accurate estimate of the large scale bias factor. Then we use the mass function measured in numerical simulations of SCDM, OCDM and LCDM to compute this bias. Comparison with these simulations shows that this simple way of estimating the bias relation and its evolution is accurate for less massive haloes as well as massive ones. In particular, we show that haloes which are less/more massive than typical M* haloes at the time they form are more/less strongly clustered than formulae based on the standard Press-Schechter mass function predict.

In hierarchical models of gravitational clustering, virialized haloes are biased tracers of the matter distribution. As discussed by Mo & White (1996), this bias is nonlinear and stochastic. They developed a model which allows one to write down analytic expressions for the mean of the bias relation, in the initial Lagrangian, and the evolved, Eulerian spaces. We provide analytic expressions for the higher order moments as well. We also show how the Mo-White approach can be extended to compute the evolution, not just of the haloes, but of the dark matter distribution itself. The model predictions for the mean and scatter in the Eulerian bias relation, as well as for the Eulerian halo-mass and halo-halo correlation functions, are in reasonable agreement with numerical simulations of hierarchical gravitational clustering for haloes of a wide range of masses, whereas the predictions for the corresponding Lagrangian quantities are accurate only for massive haloes.

We describe a simple efficient algorithm that allows one to construct Monte-Carlo realizations of merger histories of dark matter halos. The algorithm is motivated by the excursion set model (Bond et al. 1991) for the conditional and unconditional halo mass functions. The forest of trees constructed using this algorithm depends on the underlying power spectrum. For Poisson or white-noise initial power-spectra, the forest has exactly the same properties as the ensemble of trees described by Sheth (1996) and Sheth & Pitman (1997). In this case, many ensemble averaged higher order statistics of the tree distribution can be computed analytically. For Gaussian initial conditions with more general power-spectra, mean properties of the ensemble closely resemble the mean properties expected from the excursion set approach. Various statistical quantities associated with the trees constructed using our algorithm are in good agreement with what is measured in numerical simulations of hierarchical gravitational clustering.

A model of the gravitationally evolved dark matter distribution, in the Eulerian space, is developed. It is a simple extension of the excursion set model that is commonly used to estimate the mass function of collapsed dark matter haloes. In addition to describing the evolution of the dark matter itself, the model allows one to describe the evolution of the Eulerian space distribution of the haloes. It can also be used to describe density profiles, on scales larger than the virial radius, of these haloes, and to quantify the way in which matter flows in and out of Eulerian cells. When the initial Lagrangian space distribution is white noise Gaussian, the model suggests that the Inverse Gaussian distribution should provide a reasonably good approximation to the evolved Eulerian density field, in agreement with numerical simulations. Application of this model to clustering from more general Gaussian initial conditions is discussed at the end.

We develop a model for the growth of dark matter halos and use it to study their evolved density profiles. In this model, halos are spherical and form by quiescent accretion of matter in clumps, called satellites. The halo mass as a function of redshift is given by the mass of the most massive progenitor, and is determined from Monte-Carlo realizations of the merger-history tree. Inside the halo, satellites move under the action of the gravitational force of the halo and a dynamical friction drag force. The associated equation of motion is solved numerically. The energy lost to dynamical friction is transferred to the halo in the form of kinetic energy. As they sink into the halo, satellites continually lose matter as a result of tidal stripping. The stripped matter moves inside the halo free of dynamical friction. The evolved density profiles are steeper than those obtained by assuming that, once they have been accreted onto the parent halo, satellites remain at a fixed distance from the halo center. We find that the final density profile depends mainly on the rate of infall of matter onto the halo. This, in turn, depends on the initial fluctuation field as well as on cosmology. For mass scales where the effective spectral index of the initial density field is less than -1, the model predicts a profile which can only approximately be matched by the one parameter family of curves suggested by Navarro, Frenk and White (1997). For scale-free power-spectra with initial slope $n$, the density profile within about 1% of the virial radius is $\rho\propto r^{-\beta}$, with $3(3+n)/(5+n)\le\beta\le 3(3+n)/(4+n)$.

This paper presents a new derivation of the Generalized Poisson distribution. This distribution provides a good fit to the evolved, counts-in-cells distribution measured in numerical simulations of hierarchical clustering from Poisson initial conditions. This new derivation can be used to construct a useful analytic model of the evolution of clustering measured in these simulations. The model is consistent with the assumption that, as the universe expands and the comoving sizes of regions change as a result of gravitational instability, the number of such expanding and contracting regions is conserved. The model neglects the influence of external tides on the evolution of such regions. Indeed, in the context of this model, the GPD can be thought of as arising from a simple variant of the well-studied spherical collapse model, in which tidal effects are also neglected. The model can be understood as a simple generalization of the excursion set model which is used to estimate the number density of collapsed, virialized halos. The generalization allows one to describe the spatial distribution of the halos and the mass simultaneously, so the extent to which these halos are biased tracers of the underlying matter distribution is easy to estimate.

The Press--Schechter description of gravitational clustering from an initially Poisson distribution is shown to be equivalent to the well studied Galton--Watson branching process. This correspondence is used to provide a detailed description of the evolution of hierarchical clustering, including a complete description of the merger history tree. The relation to branching process epidemic models means that the Press--Schechter description can also be understood using the formalism developed in the study of queues. The queueing theory formalism, also, is used to provide a complete description of the merger history of any given Press--Schechter clump. In particular, an analytic expression for the merger history of any given Poisson Press--Schechter clump is obtained. This expression allows one to calculate the partition function of merger history trees. It obeys an interesting scaling relation; the partition function for a given pair of initial and final epochs is the same as that for certain other pairs of initial and final epochs. The distribution function of counts in randomly placed cells, as a function of time, is also obtained using the branching process and queueing theory descriptions. Thus, the Press--Schechter description of the gravitational evolution of clustering from an initially Poisson distribution is now complete. All these interrelations show why the Press--Schechter approach works well in a statistical sense, but cannot provide a detailed description of the dynamics of the clustering particles themselves. One way to extend these results to more general Gaussian initial conditions is discussed.

The correlation function xi(r) of matter in the non-linear regime is assumed to be determined by the density profiles rho(r) and the mass distribution n(M) of virialized halos. The Press--Schechter approach is used to compute n(M), and the stable clustering hypothesis is used to determine the density profiles of these Press--Schechter halos. Thus, the shape and amplitude of xi(r) on small scales is related to the initial power spectrum of density fluctuations. The case of clustering from scale-free initial conditions is treated in detail. If n is the slope of the initial power spectrum of density fluctuations, then stable clustering requires that xi(r)\propto r^{-gamma}, where gamma is a known function of n. If halo--halo correlations can be neglected, then rho(r)\propto r^{-epsilon}, where epsilon = (gamma+3)/2 = 3(4+n)/(5+n). For all values of n of current interest, this slope is steeper than the value 3(3+n)/(4+n) that was obtained by Hoffman & Shaham in their treatment of the shapes of the outer regions of collapsed halos. Our main result is a prediction for the amplitude of the non-linear correlation function. The predicted amplitude and its dependence on n are in good quantitative agreement with N-body simulations of self-similar clustering. If stable clustering is a good approximation only inside the half-mass radii of Press--Schechter halos, then the density contrast required for the onset of stable clustering can be estimated. This density contrast is in the range ~300-600 and increases with the initial slope n, in agreement with estimates from N-body simulations.

With appropriate modifications, the Hoffman--Ribak algorithm that constructs constrained realizations of Gaussian random fields having the correct ensemble properties can also be used to construct constrained realizations of those non-Gaussian random fields that are obtained by transformations of an underlying Gaussian field. For example, constrained realizations of lognormal, generalized Rayleigh, and chi-squared fields having $n$ degrees of freedom constructed this way will have the correct ensemble properties. The lognormal field is considered in detail. For reconstructing Gaussian random fields, constrained realization techniques are similar to reconstructions obtained using minimum variance techniques. A comparison of this constrained realization approach with minimum variance, Wiener filter reconstruction techniques, in the context of lognormal random fields, is also included. The resulting prescriptions for constructing constrained realizations as well as minimum variance reconstructions of lognormal random fields are useful for reconstructing masked regions in galaxy catalogues on smaller scales than previously possible, for assessing the statistical significance of small-scale features in the microwave background radiation, and for generating certain non-Gaussian initial conditions for $N$-body simulations.

The distribution of pairwise, relative peculiar velocities, $f(u;r)$, on small nonlinear scales, $r$, is derived from the Press--Schechter approach. This derivation assumes that Press--Schechter clumps are virialized and isothermal. The virialized assumption requires that the circular velocity, $V_c \propto M^{1/3}$, where $M$ denotes the mass of the clump. The isothermal assumption means that the circular velocity is independent of radius. Further, it is assumed that the velocity distribution within a clump is Maxwellian, that the pairwise relative velocity distribution is isotropic, and that on nonlinear scales clump-clump motions are unimportant when calculating the distribution of velocity differences. Comparison with $N$-body simulations shows that, on small nonlinear scales, all these assumptions are accurate. For most power spectra of interest, the resulting line of sight, pairwise, relative velocity distribution, $f(u_{\rm r})$, is well approximated by an exponential, rather than a Gaussian distribution. This simple Press--Schechter model is also able to provide a natural explanation for the observed, non-Maxwellian shape of $f(v)$, the distribution of peculiar velocities.