Normalized to: Sheth, R.

**200** article(s) in total. **440** co-authors, from **1** to **49** common article(s). Median position in authors list is **3,0**.

[1] oai:arXiv.org:1912.07072 [pdf] - 2045381

Non-linear statistics of primordial black holes from gaussian curvature
perturbations

Comments: v2: improved discussions and citations. Constraint on the abundance
related to minimums added. Discussion on the broad power spectrum and related
appendix added. 19 pages, 5 figures

Submitted: **2019-12-15**, last modified: **2020-02-10**

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.

[2] oai:arXiv.org:1907.13311 [pdf] - 2050185

A universal threshold for primordial black hole formation

Comments: 5 pages and 5 figures, v3: clarifications and more numerical
evidences for the lower bound added. Version accepted in PRD

Submitted: **2019-07-31**, last modified: **2020-01-27**

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$.

[3] oai:arXiv.org:1912.06153 [pdf] - 2026544

Probing black hole accretion tracks, scaling relations and radiative
efficiencies from stacked X-ray active galactic nuclei

Comments: 12 pages, 8 Figures. MNRAS, accepted. A discussion around the points
raised by arXiv:1909.10821 is included in the Appendix

Submitted: **2019-12-12**

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.

[4] oai:arXiv.org:1910.10698 [pdf] - 1985419

Linear Point and Sound Horizon as Purely Geometric standard rulers

Comments: 13 pages, 7 figures

Submitted: **2019-10-23**

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.

[5] oai:arXiv.org:1910.10175 [pdf] - 1985368

Constraining black hole-galaxy scaling relations from the large-scale
clustering of Active Galactic Nuclei and implied mean radiative efficiency

Comments: 31 pages, 7 Figures. Accepted

Submitted: **2019-10-22**

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.

[6] oai:arXiv.org:1904.11996 [pdf] - 1975251

Galaxy properties as revealed by MaNGA II. Differences in stellar
populations of slow and fast rotator ellipticals and dependence on
environment

Comments: 22 pages, 28 figures. Accepted for publication in MNRAS

Submitted: **2019-04-26**, last modified: **2019-08-27**

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.

[7] oai:arXiv.org:1904.11992 [pdf] - 1975250

Galaxy properties as revealed by MaNGA. I. Constraints on Initial Mass
Function and M$_{*}$/L gradients in ellipticals

Comments: 22 pages, 29 figures, accepted for publication in MNRAS

Submitted: **2019-04-26**, last modified: **2019-08-27**

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.

[8] oai:arXiv.org:1903.02007 [pdf] - 1948159

Cosmic web anisotropy is the primary indicator of halo assembly bias

Comments: 21 pages, 14 figures; v2: Added discussion, results strengthened.
Conclusions unchanged. Accepted by MNRAS

Submitted: **2019-03-05**, last modified: **2019-08-20**

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.

[9] oai:arXiv.org:1907.09147 [pdf] - 1920217

Excursion set peaks in energy as a model for haloes

Comments: 9 pages, 5 figures

Submitted: **2019-07-22**

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.

[10] oai:arXiv.org:1906.02081 [pdf] - 1896063

Exact enumeration approach to first-passage time distribution of
non-Markov random walks

Comments: 23 pages, 4 figures, 1 table and 5 appendices. Published version

Submitted: **2019-06-05**

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.

[11] oai:arXiv.org:1904.13347 [pdf] - 1966740

Constraints on the redshift evolution of astrophysical feedback with
Sunyaev-Zeldovich effect cross-correlations

Comments: 21 pages, 12 figures, comments welcome

Submitted: **2019-04-30**

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.

[12] oai:arXiv.org:1707.08280 [pdf] - 1866948

Radial Acceleration Relation between Baryons and Dark or Phantom Matter
in the Super-critical Acceleration Regime of Nearly Spherical Galaxies

Comments: ApJ, accepted (22 pages including 19 figures)

Submitted: **2017-07-25**, last modified: **2019-04-12**

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.

[13] oai:arXiv.org:1903.05161 [pdf] - 1848637

Cosmic voids: a novel probe to shed light on our Universe

Comments: Science White Paper for the Astro2020 Decadal Survey, 5 pages text,
14 pages with references and affiliations

Submitted: **2019-03-12**

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.

[14] oai:arXiv.org:1902.09350 [pdf] - 1860044

Modeling Nearly Spherical Pure-Bulge Galaxies with a Stellar
Mass-to-Light Ratio Gradient under the $\Lambda$CDM and MOND Paradigms: II.
The Orbital Anisotropy of Slow Rotators within the Effective Radius

Comments: 11 figures, 3 tables, ApJ, accepted for publication

Submitted: **2019-02-25**

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.

[15] oai:arXiv.org:1901.11036 [pdf] - 1830682

Black hole scaling relations of active and quiescent galaxies:
Addressing selection effects and constraining virial factors

Comments: 15 pages, 5 Figures. MNRAS, accepted

Submitted: **2019-01-30**

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.

[16] oai:arXiv.org:1811.12312 [pdf] - 1901766

Cosmic distance inference from purely geometric BAO methods: Linear
Point standard ruler and Correlation Function Model Fitting

Comments: 12 pages, 2 tables

Submitted: **2018-11-29**

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.

[17] oai:arXiv.org:1811.03132 [pdf] - 1781668

Density and velocity profiles around cosmic voids

Comments: 16 pages, 18 figures

Submitted: **2018-11-07**

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.

[18] oai:arXiv.org:1804.00040 [pdf] - 1705242

Modeling Nearly Spherical Pure-bulge Galaxies with a Stellar
Mass-to-Light Ratio Gradient under the $\Lambda$CDM and MOND Paradigms: I.
Methodology, Dynamical Stellar Mass, and Fundamental Mass Plane

Comments: ApJ, in press (32 pages in preprint style including 22 figures; Fig
21 corrected)

Submitted: **2018-03-30**, last modified: **2018-05-31**

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.

[19] oai:arXiv.org:1802.04207 [pdf] - 1807278

The Excursion set approach: Stratonovich approximation and Cholesky
decomposition

Comments: 6 pages, 2 figures, Accepted for publication in MNRAS

Submitted: **2018-02-12**, last modified: **2018-05-29**

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.

[20] oai:arXiv.org:1712.05414 [pdf] - 1663254

$M_*/L$ gradients driven by IMF variation: Large impact on dynamical
stellar mass estimates

Comments: 13 pages, 7 figures, MNRAS in press

Submitted: **2017-12-14**, last modified: **2018-03-21**

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.

[21] oai:arXiv.org:1712.00324 [pdf] - 1684565

Bimodal Formation Time Distribution for Infall Dark Matter Halos

Comments: 21 pages, 7 figures, accepted for publication by ApJ

Submitted: **2017-12-01**, last modified: **2018-03-14**

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.

[22] oai:arXiv.org:1706.09906 [pdf] - 1648584

Halo assembly bias and the tidal anisotropy of the local halo
environment

Comments: 19 pages, 15 figures; v2: revised in response to referee comments,
added references and discussion, conclusions unchanged. Accepted in MNRAS

Submitted: **2017-06-29**, last modified: **2018-02-25**

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.

[23] oai:arXiv.org:1801.04568 [pdf] - 1656265

The dependence of galaxy clustering on tidal environment in the Sloan
Digital Sky Survey

Comments: 11 pages, 6 figures, submitted to MNRAS

Submitted: **2018-01-14**, last modified: **2018-01-16**

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.

[24] oai:arXiv.org:1710.07296 [pdf] - 1626350

Stellar mass functions and implications for a variable IMF

Comments: 17 pages, 19 figures, 4 tables. Accepted for publication by MNRAS.
Tables 1, C1 and C2 are provided as ancillary files

Submitted: **2017-10-19**, last modified: **2017-12-15**

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.

[25] oai:arXiv.org:1711.09063 [pdf] - 1724849

The Linear Point Standard Ruler for galaxy survey data: validation with
mock catalogues

Comments: The real data distance constraints presented in the companion paper
are improved: arXiv:1703.01275

Submitted: **2017-11-24**

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].

[26] oai:arXiv.org:1703.01275 [pdf] - 1716926

The Linear Point: A cleaner cosmological standard ruler

Comments: The optimal two-point correlation function bin-size is employed.
Results are updated and the distance constraints are improved

Submitted: **2017-03-03**, last modified: **2017-11-24**

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.

[27] oai:arXiv.org:1511.01909 [pdf] - 1597979

Effective Window Function for Lagrangian Halos

Comments: 15 pages, 14 figures, matched the published version, discussion and
presentation significantly improved, conclusions unchanged

Submitted: **2015-11-05**, last modified: **2017-11-22**

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.

[28] oai:arXiv.org:1707.04096 [pdf] - 1585891

Dependence of halo bias on mass and environment

Comments: 9 pages, 6 figures, accepted by MNRAS

Submitted: **2017-07-13**, last modified: **2017-08-31**

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.

[29] oai:arXiv.org:1702.01762 [pdf] - 1575983

Selection bias in dynamically-measured supermassive black hole samples:
Scaling relations and correlations between residuals in semi-analytic galaxy
formation models

Comments: 11 pages, 7 figures. Minor changes to match version to be published
in MNRAS

Submitted: **2017-02-06**, last modified: **2017-04-18**

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.

[30] oai:arXiv.org:1701.01701 [pdf] - 1564005

Consistency relations for the Lagrangian halo bias and their
implications

Comments: 17 pages, 19 figures, match the published version

Submitted: **2017-01-06**, last modified: **2017-04-03**

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.

[31] oai:arXiv.org:1702.08527 [pdf] - 1581662

Comparing PyMorph and SDSS photometry. II. The differences are more than
semantics and are not dominated by intracluster light

Comments: 15 pages, 17 figures, accepted for publication in MNRAS. The PyMorph
luminosities and stellar masses are available at
https://www.physics.upenn.edu/~ameert/SDSS_PhotDec/

Submitted: **2017-02-27**, last modified: **2017-03-20**

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.

[32] oai:arXiv.org:1703.06145 [pdf] - 1574687

Revisiting the bulge-halo conspiracy I: Dependence on galaxy properties
and halo mass

Comments: 23 pages, 10 figures, 3 Appendices (with an extra 7 figures). ApJ,
accepted. Main results in Figures 3, 5, 6, 8

Submitted: **2017-03-17**

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.

[33] oai:arXiv.org:1508.01170 [pdf] - 1548121

Beating non-linearities: improving the Baryon Acoustic Oscillations with
the linear point

Comments: 10 pages, 9 figures. Added discussion on the correlation function
inflection point. This version matches the published MNRAS paper

Submitted: **2015-08-05**, last modified: **2017-03-16**

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.

[34] oai:arXiv.org:1604.01036 [pdf] - 1538894

The high mass end of the stellar mass function: Dependence on stellar
population models and agreement between fits to the light profile

Comments: 18 pages, 17 figures, accepted for publication in MNRAS. The PyMorph
luminosities and stellar masses are available at
https://www.physics.upenn.edu/~ameert/SDSS_PhotDec/

Submitted: **2016-04-04**, last modified: **2017-02-27**

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.

[35] oai:arXiv.org:1702.01722 [pdf] - 1581397

The Halo Boundary of Galaxy Clusters in the SDSS

Comments: 17 pages, 8 figures

Submitted: **2017-02-06**, last modified: **2017-02-07**

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.

[36] oai:arXiv.org:1701.01732 [pdf] - 1534239

Selection bias in dynamically-measured super-massive black hole samples:
dynamical masses and dependence on S\'ersic index

Comments: 12 pages, 8 Figures. MNRAS, in press

Submitted: **2017-01-06**

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.

[37] oai:arXiv.org:1603.01276 [pdf] - 1530513

Selection bias in dynamically-measured super-massive black hole samples:
its consequences and the quest for the most fundamental relation

Comments: 24 pages, 11 Figures, two Appendices. Replaced with accepted version
in MNRAS. Minor revisions. A press release of the Royal Astronomical Society
can be found at
http://www.ras.org.uk/news-and-press/2861-black-holes-and-the-prospects-for-measuring-gravitational-waves

Submitted: **2016-03-03**, last modified: **2017-01-05**

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.

[38] oai:arXiv.org:1611.03619 [pdf] - 1512695

Excursion set peaks: the role of shear

Comments: 17 pages, 13 figures, submitted to MNRAS

Submitted: **2016-11-11**

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.

[39] oai:arXiv.org:1611.03613 [pdf] - 1574568

Constraints on halo formation from cross-correlations with correlated
variables

Comments: 15 pages, 17 figures, submitted to MNRAS

Submitted: **2016-11-11**

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.

[40] oai:arXiv.org:1610.02485 [pdf] - 1531954

An Order Statistics Approach to the Halo Model for Galaxies

Comments: 16 pages, 18 figures, submitted to MNRAS

Submitted: **2016-10-08**

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

[41] oai:arXiv.org:1603.09348 [pdf] - 1475301

Selection bias in dynamically-measured super-massive black hole samples:
consequences for pulsar timing arrays

Comments: 6 pages 4 figures, submitted to MNRAS letters

Submitted: **2016-03-30**

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.

[42] oai:arXiv.org:1507.05627 [pdf] - 1347488

The universality of the virial halo mass function and models for
non-universality of other halo definitions

Comments: 22 pages, 21 figures, accepted for publication in MNRAS - some
changes with respect to the previous version and a new appendix

Submitted: **2015-07-20**, last modified: **2015-12-09**

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.

[43] oai:arXiv.org:1510.07702 [pdf] - 1327522

The massive end of the luminosity and stellar mass functions and
clustering from CMASS to SDSS: Evidence for and against passive evolution

Comments: Accepted for publication in MNRAS, 15 pages, 12 figures

Submitted: **2015-10-26**

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.

[44] oai:arXiv.org:1501.02800 [pdf] - 1449758

Avoiding progenitor bias: The structural and mass evolution of Brightest
Group and Cluster Galaxies in Hierarchical models since z~1

Comments: 9 pages, 5 figures. ApJ, in press. Additional (minor) editing
included in text and figure captions

Submitted: **2015-01-12**, last modified: **2015-02-23**

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.

[45] oai:arXiv.org:1212.1166 [pdf] - 895197

Self-consistency of the Excursion Set Approach

Comments: 5 pages, 2 figures

Submitted: **2012-12-05**, last modified: **2014-11-11**

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.

[46] oai:arXiv.org:1411.2597 [pdf] - 1451815

On the intermediate-redshift central stellar mass-halo mass relation,
and implications for the evolution of the most massive galaxies since z~1

Comments: 6 pages, 4 figures. ApJ Letters, accepted

Submitted: **2014-11-10**

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.

[47] oai:arXiv.org:1401.3185 [pdf] - 848101

Stochasticity in halo formation and the excursion set approach

Comments: 5 pages, 1 figure. Matches published version

Submitted: **2014-01-14**, last modified: **2014-07-08**

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.

[48] oai:arXiv.org:1401.8177 [pdf] - 880175

On the Markovian assumption in the excursion set approach: the
approximation of Markov Velocities

Comments: 14 pages, 7 figures. Matches MNRAS version

Submitted: **2014-01-31**, last modified: **2014-07-08**

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.

[49] oai:arXiv.org:1211.6122 [pdf] - 831837

Systematic effects on the size-luminosity relation: dependence on model
fitting and morphology

Comments: 27 pages, 31 figures, accepted for publication in MNRAS

Submitted: **2012-11-26**, last modified: **2014-06-04**

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.

[50] oai:arXiv.org:1305.0724 [pdf] - 848079

The excursion set approach in non-Gaussian random fields

Comments: 14 pages, new sections and figures describing new results, discussion
and references added

Submitted: **2013-05-03**, last modified: **2014-01-20**

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.

[51] oai:arXiv.org:1401.2460 [pdf] - 791876

Environmental dependence of bulge-dominated galaxy sizes in hierarchical
models of galaxy formation. Comparison with the local Universe

Comments: 27 pages, 14 figures, 2 tables. MNRAS, in press

Submitted: **2014-01-10**

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.

[52] oai:arXiv.org:1311.1212 [pdf] - 791768

Cosmology with massive neutrinos II: on the universality of the halo
mass function and bias

Comments: 21 pages, 9 figures

Submitted: **2013-11-05**

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.

[53] oai:arXiv.org:1305.5830 [pdf] - 1171576

Bias deconstructed: Unravelling the scale dependence of halo bias using
real space measurements

Comments: 13 pages, 9 figures; v3 -- Matches published version

Submitted: **2013-05-24**, last modified: **2013-10-09**

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.

[54] oai:arXiv.org:1304.7778 [pdf] - 1166300

The massive end of the luminosity and stellar mass functions: Dependence
on the fit to the light profile

Comments: 8 pages, 6 figures; accepted for publication in MNRAS. Electronic
versions of Tables 2 and 3, which list the binned counts, are provided as
"ancillary" files

Submitted: **2013-04-29**, last modified: **2013-09-09**

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.

[55] oai:arXiv.org:1212.0964 [pdf] - 706877

Peak-peak correlations in the cosmic background radiation from cosmic
strings

Comments: 11 pages, 7 figures and 1 table; V2: Added one figure, comments,
references and moderate corrections. Accepted for publication in MNRAS

Submitted: **2012-12-05**, last modified: **2013-08-14**

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.

[56] oai:arXiv.org:1306.0551 [pdf] - 848081

The importance of stepping up in the excursion set approach

Comments: 10 pages, 5 figures

Submitted: **2013-06-03**

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.

[57] oai:arXiv.org:1305.1505 [pdf] - 1166457

An accurate tool for the fast generation of dark matter halo catalogs

Comments: 14 pages, 9 PDF figures, accepted by Monthly Notices of the Royal
Astronomical Society. The PINOCCHIO code can be downloaded here:
http://adlibitum.oats.inaf.it/monaco/Homepage/Pinocchio/index.html

Submitted: **2013-05-07**

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.

[58] oai:arXiv.org:1207.7117 [pdf] - 648333

Non-local Lagrangian bias

Comments: 13 pages, 6 figures; new section shows agreement with bispectrum
measurements in simulations

Submitted: **2012-07-30**, last modified: **2013-04-05**

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.

[59] oai:arXiv.org:1212.4157 [pdf] - 640461

Ellipsoidal halo finders and implications for models of triaxial halo
formation

Comments: 18 pages, 17 figures, accepted for publication in MNRAS

Submitted: **2012-12-17**, last modified: **2013-03-17**

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]

[60] oai:arXiv.org:1208.0229 [pdf] - 637162

Lagrangian bias in the local bias model

Comments: 14 pages, published on JCAP

Submitted: **2012-08-01**, last modified: **2013-03-11**

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.

[61] oai:arXiv.org:1210.1483 [pdf] - 1151880

Excursion set peaks: a self-consistent model of dark halo abundances and
clustering

Comments: 12 pages, 8 figures; v2 -- revised version accepted in MNRAS. Added
clarifications, new figures and discussion of further tests of the framework,
and of consequences for non-universality of the mass function. Original
conclusions unchanged

Submitted: **2012-10-04**, last modified: **2013-02-13**

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.

[62] oai:arXiv.org:1211.0287 [pdf] - 619960

Measures of galaxy environment -- II. Rank-ordered mark correlations

Comments: 12 pages, 15 figures; published in MNRAS. The mock galaxy catalogues
and environment measures used in this paper and Paper I (arXiv:1109.6328) are
available here:
http://tao.it.swin.edu.au/mock-galaxy-factory/precomputed-mock-catalogues/

Submitted: **2012-11-01**, last modified: **2013-01-29**

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.

[63] oai:arXiv.org:1301.5128 [pdf] - 1159196

Stochastic bias in multi-dimensional excursion set approaches

Comments: 9 pages, 5 figures, submitted to MNRAS

Submitted: **2013-01-22**

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.

[64] oai:arXiv.org:1207.1853 [pdf] - 602618

On the anisotropic density distribution on large scales

Comments: 6 pages, 3 figures; With a moderate revision of the original version,
as published in MNRAS

Submitted: **2012-07-08**, last modified: **2012-12-13**

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.

[65] oai:arXiv.org:1206.3506 [pdf] - 1124143

Peaks theory and the excursion set approach

Comments: 10 pages, 4 figures; v2 -- minor changes, added discussion in sec2.2,
fixed a typo. Accepted in MNRAS

Submitted: **2012-06-15**, last modified: **2012-08-10**

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.

[66] oai:arXiv.org:1205.3401 [pdf] - 1123408

Scale dependent halo bias in the excursion set approach

Comments: 16 pages, 10 figures; v2 -- minor changes, typos fixed, references
added, accepted in MNRAS

Submitted: **2012-05-15**, last modified: **2012-08-09**

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.

[67] oai:arXiv.org:1203.6609 [pdf] - 1117645

The clustering of galaxies in the SDSS-III Baryon Oscillation
Spectroscopic Survey: a large sample of mock galaxy catalogues

Comments: 21 pages, 18 figures, 2 tables. See companion papers that share the
"The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic
Survey:" title. Mocks and covariance matrices will be available at
http://www.marcmanera.net/mocks when the paper is accepted and DR9 released.
Minor edits. Amended author's names. Submitted to MNRAS

Submitted: **2012-03-29**, last modified: **2012-04-01**

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.

[68] oai:arXiv.org:1107.3652 [pdf] - 525786

The luminosities of the brightest cluster galaxies and brightest
satellites in SDSS groups

Comments: 12 pages, 9 figures; v2 -- Revised to match version accepted in
MNRAS. Includes a new section on using mark correlations to test extreme
value statistics

Submitted: **2011-07-19**, last modified: **2012-03-30**

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.

[69] oai:arXiv.org:1202.5804 [pdf] - 481272

Cosmic Tides

Comments: 4 pages, 4 figures

Submitted: **2012-02-26**

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.

[70] oai:arXiv.org:1201.3614 [pdf] - 496728

Gravity and Large-Scale Non-local Bias

Comments: 26 pages, 9 figures. v2: improved appendix A

Submitted: **2012-01-17**, last modified: **2012-02-21**

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)

[71] oai:arXiv.org:1202.3438 [pdf] - 1116675

Plain fundamentals of Fundamental Planes: Analytics and algorithms

Comments: 11 pages, 1 figure, associated IDL code, MNRAS accepted

Submitted: **2012-02-15**

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.

[72] oai:arXiv.org:1111.6977 [pdf] - 1092011

Formation times, mass growth histories and concentrations of dark matter
haloes

Comments: 15 pages, 15 figures, MNRAS accepted (minor changes + 1 figure and
some references)

Submitted: **2011-11-29**, last modified: **2012-01-19**

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.

[73] oai:arXiv.org:1201.3876 [pdf] - 618541

One step beyond: The excursion set approach with correlated steps

Comments: 5 pages, 1 figure. Uses mn2e class style

Submitted: **2012-01-18**

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.

[74] oai:arXiv.org:1106.2041 [pdf] - 480017

A hierarchy of voids: More ado about nothing

Comments: 9 pages, 3 figures; v2 -- improved Eulerian void-assignment
algorithm, new figures (including LCDM walks) and clarified discussion.
Conclusions regarding walks with correlated steps unchanged. Accepted in
MNRAS

Submitted: **2011-06-10**, last modified: **2011-11-08**

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.

[75] oai:arXiv.org:1105.2261 [pdf] - 480013

Halo bias in the excursion set approach with correlated steps

Comments: 7 pages, 3 figures; v2 -- minor clarifications, accepted in MNRAS

Submitted: **2011-05-11**, last modified: **2011-11-03**

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.

[76] oai:arXiv.org:1105.1990 [pdf] - 480012

Halo abundances and counts-in-cells: The excursion set approach with
correlated steps

Comments: 14 pages, 9 figures; v2 -- revised version with explicit
demonstration that the original conclusions hold for LCDM, expanded
discussion on stochasticity of barrier. Accepted in MNRAS

Submitted: **2011-05-10**, last modified: **2011-11-03**

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)

[77] oai:arXiv.org:1006.2243 [pdf] - 426054

The cosmological free-free signal from galaxy groups and clusters

Comments: 11 pages, 7 figures

Submitted: **2010-06-11**, last modified: **2011-10-17**

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.

[78] oai:arXiv.org:1105.3378 [pdf] - 1076691

How unusual are the Shapley Supercluster and the Sloan Great Wall?

Comments: 13 pages, 8 figures

Submitted: **2011-05-17**

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.

[79] oai:arXiv.org:1010.2839 [pdf] - 1041291

Modelling the shapes of the largest gravitationally bound objects

Comments: 15 pages, 8 figures. New appendix added, extended discussion. Matches
version accepted by MNRAS

Submitted: **2010-10-14**, last modified: **2011-05-16**

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.

[80] oai:arXiv.org:1011.1501 [pdf] - 1041785

Evidence of major dry mergers at M* > 2 x 10^11 Msun from curvature in
early-type galaxy scaling relations?

Comments: 5 pages, 3 figures. Accepted for publication in MNRAS

Submitted: **2010-11-05**, last modified: **2011-01-19**

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.

[81] oai:arXiv.org:1012.1335 [pdf] - 333497

Scale-dependent halo bias from scale-dependent growth

Comments: 14 pages, submitted to Phys. Rev. D

Submitted: **2010-12-06**, last modified: **2011-01-12**

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.

[82] oai:arXiv.org:1007.3500 [pdf] - 1033801

Optimal linear reconstruction of dark matter from halo catalogs

Comments: Added Figure 3 to show the scatter between the weighted halo field vs
the mass field, Accepted for publication in MNRAS

Submitted: **2010-07-20**, last modified: **2010-11-08**

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)

[83] oai:arXiv.org:1005.3770 [pdf] - 1032637

Curvature in the color-magnitude relation but not in color-sigma: Major
dry mergers at M* > 2 x 10^11 Msun?

Comments: 22 pages, 22 figures and 3 tables. Accepted for publication in MNRAS

Submitted: **2010-05-20**, last modified: **2010-11-05**

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.

[84] oai:arXiv.org:1009.3449 [pdf] - 282653

Modeling scale-dependent bias on the baryonic acoustic scale with the
statistics of peaks of Gaussian random fields

Comments: (v1): 23 pages text, 8 figures + appendix (v2): typos fixed,
references added, accepted for publication in PRD

Submitted: **2010-09-17**, last modified: **2010-10-19**

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.

[85] oai:arXiv.org:0908.2285 [pdf] - 1003482

The nonlinear redshift space probability distribution function in models
with local primordial non-Gaussianity

Comments: Matched version accepted by MNRAS

Submitted: **2009-08-17**, last modified: **2010-08-02**

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}.

[86] oai:arXiv.org:1003.4740 [pdf] - 1025923

Halo model description of the non-linear dark matter power spectrum at
$k \gg 1$ Mpc$^{-1}$

Comments: Accepted for publication in MNRAS - minor changes

Submitted: **2010-03-24**, last modified: **2010-06-30**

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.

[87] oai:arXiv.org:0910.1093 [pdf] - 1017853

Galaxy luminosities, stellar masses, sizes, velocity dispersions as a
function of morphological type

Comments: 42 pages, 34 figures, 9 tables. Accepted for publication in MNRAS

Submitted: **2009-10-07**, last modified: **2010-01-29**

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.

[88] oai:arXiv.org:0909.4544 [pdf] - 142229

Redshift space correlations and scale-dependent stochastic biasing of
density peaks

Comments: 24 pages, 7 figures (v3): references added (v4): added
figure+appendix. In press in PRD

Submitted: **2009-09-24**, last modified: **2010-01-18**

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).

[89] oai:arXiv.org:0911.0436 [pdf] - 1018189

The Substructure Hierarchy in Dark Matter Haloes

Comments: 18 pages, 24 figures, accepted for publication on MNRAS

Submitted: **2009-11-02**, last modified: **2010-01-07**

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.

[90] oai:arXiv.org:0910.1214 [pdf] - 1017860

Convolution and deconvolution based estimates of galaxy scaling
relations from photometric redshift surveys

Comments: 7 pages, 4 figures, accepted for publication in MNRAS; v2 includes a
new section and other minor changes

Submitted: **2009-10-07**, last modified: **2010-01-06**

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.

[91] oai:arXiv.org:0912.0012 [pdf] - 1018565

Sizes and ages of SDSS ellipticals: Comparison with hierarchical galaxy
formation models

Comments: 13 pages, 9 Figures, 1 Table. Accepted by MNRAS

Submitted: **2009-11-30**

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.

[92] oai:arXiv.org:0911.1829 [pdf] - 30482

Density profiles and voids in modified gravity models

Comments: 11 pages, 10 figures

Submitted: **2009-11-10**

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.

[93] oai:arXiv.org:0906.1314 [pdf] - 1002479

Large scale bias and the inaccuracy of the peak-background split

Comments: MNRAS accepted. New section with Spherical Overdensity identified
halos included. Appendix enlarged

Submitted: **2009-06-06**, last modified: **2009-11-09**

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.

[94] oai:arXiv.org:0906.0424 [pdf] - 1002427

Reconstructing galaxy fundamental distributions and scaling relations
from photometric redshift surveys. Applications to the SDSS early-type sample

Comments: 14 pages, 12 figures, 1 table, new section 3.1 and appendix added,
MNRAS in press

Submitted: **2009-06-02**, last modified: **2009-10-07**

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.

[95] oai:arXiv.org:0905.1702 [pdf] - 1002184

Halo abundances in the f_{nl} model

Comments: 10 pages, 3 figures, match version accepted by MNRAS

Submitted: **2009-05-11**, last modified: **2009-08-01**

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.

[96] oai:arXiv.org:0905.1706 [pdf] - 1002185

The initial shear field in models with primordial local non-Gaussianity
and implications for halo and void abundances

Comments: 16 pages, 8 figures, match version accepted by MNRAS

Submitted: **2009-05-11**, last modified: **2009-08-01**

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.

[97] oai:arXiv.org:0712.4100 [pdf] - 8468

Dark matter halo creation in moving barrier models

Comments: MNRAS accepted, 17 pages, 10 figures

Submitted: **2007-12-25**, last modified: **2009-07-10**

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.

[98] oai:arXiv.org:0906.2190 [pdf] - 25172

Non-Gaussian Distribution and Clustering of Hot and Cold Pixels in the
WMAP Five-Year Sky

Comments: 16 pages, 14 figures, submitted to MNRAS

Submitted: **2009-06-11**

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.

[99] oai:arXiv.org:0812.0200 [pdf] - 19027

Spherical Collapse and Cluster Counts in Modified Gravity Models

Comments: 9 pages, 8 figures Accepted by Phys. Rev. D

Submitted: **2008-11-30**, last modified: **2009-03-30**

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.

[100] oai:arXiv.org:0901.0757 [pdf] - 20042

On the equivalence between the effective cosmology and excursion set
treatments of environment

Comments: 4 pages, 0 figures, accepted to MNRAS

Submitted: **2009-01-07**

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.

[101] oai:arXiv.org:0712.3800 [pdf] - 8427

Merger history trees of dark matter haloes in moving barrier models

Comments: MNRAS accepted, 15 pages, 12 figures

Submitted: **2007-12-21**, last modified: **2009-01-02**

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.

[102] oai:arXiv.org:0812.2253 [pdf] - 19423

The nonlinear probability distribution function in models with local
primordial non-Gaussianity

Comments: 7 pages, 2 figures, MNRAS submitted

Submitted: **2008-12-11**

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.

[103] oai:arXiv.org:0812.1037 [pdf] - 19186

Linear theory and velocity correlations of clusters

Comments: 5 pages, 2 figures, MNRAS accepted

Submitted: **2008-12-04**

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.

[104] oai:arXiv.org:0811.2805 [pdf] - 18632

MgII absorption systems and their neighbouring galaxies from a
background subtraction technique

Comments: Submitted to MNRAS

Submitted: **2008-11-17**

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.

[105] oai:arXiv.org:0805.0310 [pdf] - 12284

A Halo Model of Galaxy Colors and Clustering in the SDSS

Comments: 13 pages, 7 figures, accepted to MNRAS. Minor revisions were made,
and a plot was added

Submitted: **2008-05-02**, last modified: **2008-09-30**

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.

[106] oai:arXiv.org:0809.2609 [pdf] - 16357

A Search for the Most Massive Galaxies. III. Global and Central
Structure

Comments: 27 Pages, 22 Figures, 2 Tables, Accepted for Publication in MNRAS

Submitted: **2008-09-15**

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.

[107] oai:arXiv.org:0809.2602 [pdf] - 16353

A Search for the Most Massive Galaxies. II. Structure, Environment and
Formation

Comments: 21 pages, 19 figures, accepted for publication in MNRAS

Submitted: **2008-09-15**

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.

[108] oai:arXiv.org:0805.1238 [pdf] - 12460

Ellipsoidal collapse and the redshift space probability distribution
function of dark matter

Comments: 10 pages, 7 figures. Accepted for publication in MNRAS

Submitted: **2008-05-08**, last modified: **2008-08-05**

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.

[109] oai:arXiv.org:0802.4288 [pdf] - 10535

The 2dF-SDSS LRG and QSO survey: evolution of the clustering of luminous
red galaxies since z = 0.6

Comments: 20 pages, replaced with version accepted for publication in MNRAS

Submitted: **2008-02-28**, last modified: **2008-05-19**

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.

[110] oai:arXiv.org:0707.3218 [pdf] - 3311

Satellite Luminosities in Galaxy Groups

Comments: 8 pages, 3 figures. Matches version accepted by MNRAS

Submitted: **2007-07-21**, last modified: **2008-05-05**

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]

[111] oai:arXiv.org:0804.0818 [pdf] - 11504

Effect of Distance Errors: Applications to SDSS Early-Type Galaxies

Comments: 5 pages, 2 figures, to appear in the proceedings of the XLIIIst
Rencontres de Moriond, Cosmology, La Thuile, Italy, 15-22 March 2008

Submitted: **2008-04-04**

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.

[112] oai:arXiv.org:0711.5029 [pdf] - 7606

Perturbation theory and excursion set estimates of the probability
distribution function of dark matter, and a method for reconstructing the
initial distribution function

Comments: 12 pages, 6 figures; Figures and texts modified; accepted for
publication in MNRAS

Submitted: **2007-11-30**, last modified: **2008-03-18**

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)

[113] oai:arXiv.org:0712.0017 [pdf] - 7611

An analytic model for the bispectrum of galaxies in redshift space

Comments: Submitted to the PRD. 28 pages, 7 Figs. Comparison with numerical
simulations improved and expanded. Analytic model predictions improved on
large scales. Errors corrected in bispectrum estimator for redshift space Q.
Other changes made in line with referees report

Submitted: **2007-12-03**, last modified: **2008-03-18**

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.

[114] oai:arXiv.org:0802.0498 [pdf] - 9749

Experimental design and model selection: The example of exoplanet
detection

Comments: 12 pages, 2 figures; To appear in the Festschrift for Jorma Rissanen;
UPR-1187

Submitted: **2008-02-04**

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.

[115] oai:arXiv.org:astro-ph/0703620 [pdf] - 90429

`Eppur Si Muove': On The Motion of the Acoustic Peak in the Correlation
Function

Comments: Revisions made in line with referee reports. 22 Pages, 9 Figures. PRD
submitted

Submitted: **2007-03-26**, last modified: **2007-12-02**

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.

[116] oai:arXiv.org:0710.1165 [pdf] - 5717

Unbiased estimates of galaxy scaling relations from photometric redshift
surveys

Comments:

Submitted: **2007-10-05**

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.

[117] oai:arXiv.org:astro-ph/0701879 [pdf] - 89008

The probability distribution function of the SZ power spectrum: an
analytical approach

Comments: v2: added one figure and discussions. Will appear in ApJ December
issue. 12 pages, 8 figures

Submitted: **2007-01-30**, last modified: **2007-09-10**

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$.

[118] oai:arXiv.org:astro-ph/0609547 [pdf] - 85139

The Scale Dependence of Halo and Galaxy Bias: Effects in Real Space

Comments: Halo power spectra at z=1 and z=2 added. Matches published version:
PRD. 75, 063512 (2007)

Submitted: **2006-09-20**, last modified: **2007-03-27**

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)

[119] oai:arXiv.org:astro-ph/0703537 [pdf] - 90346

On estimating redshift and luminosity distributions in photometric
redshift surveys

Comments: 9 pages, 6 figures, submitted to MNRAS

Submitted: **2007-03-20**

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.

[120] oai:arXiv.org:astro-ph/0609297 [pdf] - 84889

On the inconsistency between the black hole mass function inferred from
M_bh-sigma and M_bh-L correlation

Comments: 9 pages, 7 figures. ApJ in press. Expanded comparison with
literature. Reference removed as required by Journal

Submitted: **2006-09-12**, last modified: **2007-03-20**

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.

[121] oai:arXiv.org:astro-ph/0703340 [pdf] - 90149

The Three-Point Correlation Function of Luminous Red Galaxies in the
Sloan Digital Sky Survey

Comments: 14 pages. Accepted for publication to the MNRAS. Data accompanying
paper can be found at http://www.dsg.port.ac.uk/~nicholb/3pt/kulkarni/

Submitted: **2007-03-14**, last modified: **2007-03-15**

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.

[122] oai:arXiv.org:astro-ph/0703391 [pdf] - 90200

Strong clustering of underdense regions and the environmental dependence
of clustering from Gaussian initial conditions

Comments: 10 pages, 5 figures, accepted for publication in MNRAS

Submitted: **2007-03-15**

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]

[123] oai:arXiv.org:astro-ph/0611221 [pdf] - 86607

An improved model for the formation times of dark matter haloes

Comments: 8 pages, 5 figures. Replaced to match version accepted for
publication in MNRAS

Submitted: **2006-11-07**, last modified: **2007-01-18**

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.

[124] oai:arXiv.org:astro-ph/0609300 [pdf] - 316521

Selection bias in the M_BH-sigma and M_BH-L correlations and its
consequences

Comments: 12 pages, 7 figures. Accepted by ApJ. We have added a figure showing
that a similar bias is also seen in the K-band. A new appendix describes the
BH samples as well as the fits used in the main text

Submitted: **2006-09-12**, last modified: **2007-01-03**

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.

[125] oai:arXiv.org:astro-ph/0607117 [pdf] - 83320

The luminosities, sizes and velocity dispersions of Brightest Cluster
Galaxies: Implications for formation history

Comments: 16 pages, 16 figures, AJ in press

Submitted: **2006-07-07**, last modified: **2006-12-15**

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.

[126] oai:arXiv.org:astro-ph/0609740 [pdf] - 1233893

Bispectrum and Nonlinear Biasing of Galaxies: Perturbation Analysis,
Numerical Simulation and SDSS Galaxy Clustering

Comments: 14 pages, 8 figures, accepted for publication in PASJ

Submitted: **2006-09-27**, last modified: **2006-12-12**

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.

[127] oai:arXiv.org:astro-ph/0608632 [pdf] - 84509

Cosmological Constraints from the SDSS Luminous Red Galaxies

Comments: Matches accepted PRD version. SDSS data, likelihood code, Markov
chains and ppt figures available at
http://space.mit.edu/home/tegmark/sdss.html 36 journal pages, 25 figs.
CosmoMC plugin at http://cosmologist.info/cosmomc/

Submitted: **2006-08-30**, last modified: **2006-10-30**

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.

[128] oai:arXiv.org:astro-ph/0601407 [pdf] - 79265

The environmental dependence of galaxy clustering in the Sloan Digital
Sky Survey

Comments: 6 pages, 4 figures, paper replaced with accepted MNRAS version

Submitted: **2006-01-18**, last modified: **2006-08-21**

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.

[129] oai:arXiv.org:astro-ph/0608167 [pdf] - 84044

The Probability Distribution of the Lya transmitted flux from a sample
of SDSS quasars

Comments: 16 pages. submitted to MNRAS

Submitted: **2006-08-08**

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.

[130] oai:arXiv.org:astro-ph/0604581 [pdf] - 81650

Environment and the cosmic evolution of star formation

Comments: 4 pages, 4 figures, submitted to ApJL

Submitted: **2006-04-27**

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.

[131] oai:arXiv.org:astro-ph/0511365 [pdf] - 77723

An excursion set model of the cosmic web: The abundance of sheets,
filaments and halos

Comments: 22 pages, 7 figures, Accepted for publication in ApJ

Submitted: **2005-11-11**, last modified: **2006-03-27**

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.

[132] oai:arXiv.org:astro-ph/0602548 [pdf] - 80124

The Effect of Large-Scale Structure on the SDSS Galaxy Three-Point
Correlation Function

Comments: Accepted for publication in MNRAS

Submitted: **2006-02-24**

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.

[133] oai:arXiv.org:astro-ph/0512463 [pdf] - 78663

The Luminosity-Weighted or `Marked' Correlation Function

Comments: 10 pages, 9 figures, submitted to MNRAS

Submitted: **2005-12-17**

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.

[134] oai:arXiv.org:astro-ph/0509360 [pdf] - 75942

Evolution and environment of early-type galaxies

Comments: 58 pages, 22 figures. Accepted by AJ. The full version of Tables 1,
2, 3, 4 and 7, 8 are available at
http://www.physics.upenn.edu/~bernardm/PAPERS/EnvEvol/TABLES/ A README file
(bernardi_TABLES.README) explains the content of the directory. The 925
composite spectra can be downloaded from
http://www.physics.upenn.edu/~bernardm/PAPERS/EnvEvol/SPECTRA/ A README file
(bernardi_COMPOSITES.README) explains the content of the directory. Table 9:
the ZP of the Age-sigma relation has been corrected

Submitted: **2005-09-14**, last modified: **2005-11-30**

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).

[135] oai:arXiv.org:astro-ph/0511773 [pdf] - 78131

Marked correlations in galaxy formation models

Comments: 13 pages, 10 figures, submitted to MNRAS in July 2005

Submitted: **2005-11-28**

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.

[136] oai:arXiv.org:astro-ph/0511772 [pdf] - 78130

The halo-model description of marked statistics

Comments: 12 pages, 8 figures, MNRAS, in press

Submitted: **2005-11-28**

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.

[137] oai:arXiv.org:astro-ph/0510696 [pdf] - 77192

A search for the most massive galaxies: Double Trouble?

Comments: 20 pages, 8 figures. Accepted by AJ. The full set of figures in
Appendix B is available at
http://www.physics.upenn.edu/~bernardm/PAPERS/BIGEtypes/bernardi.FIG-B.ps.gz
Figure 8 did not show the set of galaxies described in the text of the
appendix. This has now been corrected

Submitted: **2005-10-24**, last modified: **2005-10-31**

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.

[138] oai:arXiv.org:astro-ph/0508382 [pdf] - 75271

The impact of halo shapes on the bispectrum in cosmology

Comments: Accepted for publication in MNRAS; minor changes; 19 pages, 10
figures

Submitted: **2005-08-17**, last modified: **2005-10-05**

(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.

[139] oai:arXiv.org:astro-ph/0509593 [pdf] - 76174

The environmental dependence of clustering in hierarchical models

Comments: 11 pages, 6 figures, MNRAS in press

Submitted: **2005-09-20**

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.

[140] oai:arXiv.org:astro-ph/0408569 [pdf] - 67060

The Luminosity and Color Dependence of the Galaxy Correlation Function

Comments: 64 pages, 24 figures. Minor changes to match accepted ApJ version

Submitted: **2004-08-30**, last modified: **2005-08-31**

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.

[141] oai:arXiv.org:astro-ph/0301280 [pdf] - 54251

On Departures From a Power Law in the Galaxy Correlation Function

Comments: 23 pages, 5 figures. Minor changes to match accepted ApJ version

Submitted: **2003-01-14**, last modified: **2005-08-31**

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.

[142] oai:arXiv.org:astro-ph/0412631 [pdf] - 1233470

Large Scale Clustering of Sloan Digital Sky Survey Quasars: Impact of
the Baryon Density and the Cosmological Constant

Comments: 26 pages, 12 figures, Accepted for publication in the PASJ

Submitted: **2004-12-26**, last modified: **2005-05-13**

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)

[143] oai:arXiv.org:astro-ph/0504510 [pdf] - 72607

Detection of Cosmic Magnification with the Sloan Digital Sky Survey

Comments: 12 pages, 8 figures, 2 tables; accepted for publication in ApJ

Submitted: **2005-04-22**

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.

[144] oai:arXiv.org:astro-ph/0409162 [pdf] - 67243

Voids in a $\Lambda$CDM Universe

Comments: 12 pages, 16 figures, submitted to MNRAS; 13 pages, 17 figures,
revised version, accepted for publication

Submitted: **2004-09-07**, last modified: **2005-04-05**

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)

[145] oai:arXiv.org:astro-ph/0503713 [pdf] - 142625

The C4 Clustering Algorithm: Clusters of Galaxies in the Sloan Digital
Sky Survey

Comments: 32 pages of figures and text accepted in AJ. Electronic version with
additional tables, links, and figures is available at
http://www.ctio.noao.edu/~chrism/c4

Submitted: **2005-03-31**

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.

[146] oai:arXiv.org:astro-ph/0401138 [pdf] - 358245

Improved Cosmological Constraints from Gravitational Lens Statistics

Comments: 35 pages, 15 figures, to be published in ApJ

Submitted: **2004-01-08**, last modified: **2004-12-02**

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.

[147] oai:arXiv.org:astro-ph/0409571 [pdf] - 67652

Colors, magnitudes and velocity dispersions in early-type galaxies:
Implications for galaxy ages and metallicities

Comments: 12 pages, 9 figures. Accepted by AJ

Submitted: **2004-09-23**

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.

[148] oai:arXiv.org:astro-ph/0403638 [pdf] - 1233347

Three-point Correlation Functions of SDSS Galaxies in Redshift Space:
Morphology, Color, and Luminosity Dependence

Comments: order of figures are changed. 9 pages, 15 figures

Submitted: **2004-03-27**, last modified: **2004-07-16**

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.

[149] oai:arXiv.org:astro-ph/0405261 [pdf] - 64781

How does clustering depend on environment?

Comments: 3 pages, 1 figure, uses iaus.cls. To appear in the proceedings of IAU
Colloquium 195: "Outskirts of galaxy clusters: intense life in the suburbs",
Torino, Italy, March 12-16, 2004

Submitted: **2004-05-13**

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.

[150] oai:arXiv.org:astro-ph/0404397 [pdf] - 64305

A Hierarchy of Voids

Comments: 6 pages, 4 figures, to appear in the Proceedings of IAU Colloquium
195: "Outskirts of Galaxy Clusters: intense life in the suburbs", Torino
Italy, March 2004

Submitted: **2004-04-20**

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.

[151] oai:arXiv.org:astro-ph/0402237 [pdf] - 62751

On the environmental dependence of halo formation

Comments: 7 pages, 10 figures, submitted to MNRAS

Submitted: **2004-02-11**

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.

[152] oai:arXiv.org:astro-ph/0402055 [pdf] - 62569

Formation times and masses of dark matter haloes

Comments: 5 pages, 5 figures, MNRAS, in press

Submitted: **2004-02-02**

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.

[153] oai:arXiv.org:astro-ph/0310723 [pdf] - 1356594

Cosmological parameters from SDSS and WMAP

Comments: Minor revisions to match accepted PRD version. SDSS data and ppt
figures available at http://www.hep.upenn.edu/~max/sdsspars.html

Submitted: **2003-10-27**, last modified: **2004-01-15**

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.

[154] oai:arXiv.org:astro-ph/0311260 [pdf] - 60840

A hierarchy of voids: Much ado about nothing

Comments: 23 pages, 9 figures, submitted to MNRAS

Submitted: **2003-11-11**

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.

[155] oai:arXiv.org:astro-ph/0310914 [pdf] - 60577

Void Hierarchy and Cosmic Structure

Comments: Contribution to ``Multi-Wavelength Cosmology'', conference Mykonos,
Greece, June 2003, eds. M. Plionis et al. 4 pages, 3 figs

Submitted: **2003-10-31**

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.

[156] oai:arXiv.org:astro-ph/0310725 [pdf] - 1938866

The 3D power spectrum of galaxies from the SDSS

Comments: 41 ApJ pages, 40 figs. Power spectrum data, likelihood software and
ppt pics available at http://www.hep.upenn.edu/~max/sdsspower.html

Submitted: **2003-10-27**, last modified: **2003-10-29**

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.

[157] oai:arXiv.org:astro-ph/0304203 [pdf] - 56111

Quasar-Galaxy and Galaxy-Galaxy Cross-Correlations: Model Predictions
with Realistic Galaxies

Comments: 11 pages, 4 figures, matches version in MNRAS

Submitted: **2003-04-10**, last modified: **2003-10-16**

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.

[158] oai:arXiv.org:astro-ph/0301303 [pdf] - 1233050

The Environment of Passive Spiral Galaxies in the SDSS

Comments: 31 pages, 13 figures, PASJ in press

Submitted: **2003-01-15**, last modified: **2003-07-22**

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.

[159] oai:arXiv.org:astro-ph/0307335 [pdf] - 723278

Physical Evidence for Dark Energy

Comments: Submitted to Phys. Rev. Letters (6 pages)

Submitted: **2003-07-18**, last modified: **2003-07-20**

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.

[160] oai:arXiv.org:astro-ph/0303092 [pdf] - 55331

The velocity dispersion function of early-type galaxies

Comments: Minor changes, matches version to appear in ApJ, 1 September 2003

Submitted: **2003-03-04**, last modified: **2003-07-15**

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.

[161] oai:arXiv.org:astro-ph/0205413 [pdf] - 1232895

Composite Luminosity Functions Based on the Sloan Digital Sky Survey
"Cut & Enhance" Galaxy Cluster Catalog

Comments: 21 pages, 5 figures, PASJ, 54, 515, miner typos fixed

Submitted: **2002-05-23**, last modified: **2003-06-11**

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.

[162] oai:arXiv.org:astro-ph/0301302 [pdf] - 1233049

The Morphological Butcher-Oemler effect in the SDSS Cut&Enhance Galaxy
Cluster Catalog

Comments: 44 pages, 15 figures, accepted for PASJ

Submitted: **2003-01-15**, last modified: **2003-06-11**

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.

[163] oai:arXiv.org:astro-ph/0304355 [pdf] - 56263

Substructure in dark matter halos: Towards a model of the abundance and
spatial distribution of subclumps

Comments: 5 pages, 2 figures, submitted to MNRAS

Submitted: **2003-04-20**

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.

[164] oai:arXiv.org:astro-ph/0301626 [pdf] - 54597

Early-type galaxies in the SDSS. III. The Fundamental Plane

Comments: 27 pages, 10 figures. Accepted by AJ (scheduled for April 2003). This
paper is part III of a revised version of astro-ph/0110344

Submitted: **2003-01-30**

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.

[165] oai:arXiv.org:astro-ph/0301624 [pdf] - 54595

Early-type galaxies in the SDSS. II. Correlations between observables

Comments: 29 pages, 11 figures. Accepted by AJ (scheduled for April 2003). This
paper is part II of a revised version of astro-ph/0110344

Submitted: **2003-01-30**

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.

[166] oai:arXiv.org:astro-ph/0301629 [pdf] - 354103

Early-type galaxies in the SDSS. IV. Colors and chemical evolution

Comments: 38 pages, 16 figures. Accepted by AJ (scheduled for April 2003). This
paper is part IV of a revised version of astro-ph/0110344. The full version
of Table 3 and the 182 composite spectra can be downloaded from
http://astrophysics.phys.cmu.edu/~bernardi/SDSS/Etypes/SPECFITS A README file
(bernardi4.README) explains the content of the directory

Submitted: **2003-01-30**

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.

[167] oai:arXiv.org:astro-ph/0301631 [pdf] - 142559

Early-type galaxies in the SDSS. I. The sample

Comments: 63 pages, 21 figures. Accepted by AJ (scheduled for April 2003). This
paper is part I of a revised version of astro-ph/0110344. The full version of
Tables 2 and 3, i.e. the tables listing the photometric and spectroscopic
parameters of ~ 9000 galaxies, are available at
http://astrophysics.phys.cmu.edu/~bernardi/SDSS/Etypes/TABLES

Submitted: **2003-01-30**

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.

[168] oai:arXiv.org:astro-ph/0210193 [pdf] - 52231

Galaxy Star-Formation as a Function of Environment in the Early Data
Release of the Sloan Digital Sky Survey

Comments: Accepted for publication in ApJ, 24 pages in emulateapj.sty

Submitted: **2002-10-08**

(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).

[169] oai:arXiv.org:astro-ph/0206293 [pdf] - 49924

A feature at z ~ 3.2 in the evolution of the Ly-alpha forest optical
depth

Comments: 45 pages, 25 figures, final version. AJ, accepted

Submitted: **2002-06-17**, last modified: **2002-09-25**

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.

[170] oai:arXiv.org:astro-ph/0208353 [pdf] - 51169

Substructure and the halo model of large-scale structure

Comments: 9 pages, 1 figure. Submitted to MNRAS

Submitted: **2002-08-19**

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.

[171] oai:arXiv.org:astro-ph/0107416 [pdf] - 43787

Analysis of Systematic Effects and Statistical Uncertainties in Angular
Clustering of Galaxies from Early SDSS Data

Comments: 64 pages, 31 figures, new version to match that accepted by ApJ

Submitted: **2001-07-20**, last modified: **2002-07-03**

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.

[172] oai:arXiv.org:astro-ph/0206508 [pdf] - 387059

Halo Models of Large Scale Structure

Comments: 166 pages, 60 figures (low resolution in most color figures); Invited
review article for Physics Reports (final version in press). Compact versions
with two pages on one and full resolution figures are available from
http://www.its.caltech.edu/~asante

Submitted: **2002-06-28**

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.

[173] oai:arXiv.org:astro-ph/0206319 [pdf] - 49950

Detection of He II reionization in the SDSS quasar sample

Comments: 4 pages, 3 figures. Minor changes, to be published in ApJL

Submitted: **2002-06-18**, last modified: **2002-06-25**

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.

[174] oai:arXiv.org:astro-ph/0105008 [pdf] - 42237

On the Distribution of Haloes, Galaxies and Mass

Comments: 10 pages, 6 figures, Latex using MN2e style. Minor changes. Accepted
for publication in MNRAS

Submitted: **2001-05-01**, last modified: **2002-02-06**

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.

[175] oai:arXiv.org:astro-ph/0107418 [pdf] - 43789

The Angular Power Spectrum of Galaxies from Early SDSS Data

Comments: Replaced to match accepted ApJ version (14 pages). Data, window
functions etc available at http://www.hep.upenn.edu/~max/sdss.html or from
max@physics.upenn.edu

Submitted: **2001-07-20**, last modified: **2002-01-19**

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.

[176] oai:arXiv.org:astro-ph/0110344 [pdf] - 45401

Early-type galaxies in the SDSS

Comments: 108 pages, 50 figures. Submitted to AJ. Higher resolution figures
available at http://astro.uchicago.edu/~bernardi/SDSS_Early_types/

Submitted: **2001-10-15**

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.

[177] oai:arXiv.org:astro-ph/0108288 [pdf] - 44264

Small Scale Anisotropies of UHECRs from Super-Heavy Halo Dark Matter

Comments: Contribution given at ICRC 2001 - August 7-15, 2001 - Hamburg

Submitted: **2001-08-17**

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.

[178] oai:arXiv.org:astro-ph/0104332 [pdf] - 42071

Non-gaussian CMB temperature fluctuations from peculiar velocities of
clusters

Comments: Revised, citation added. To appear in MNRAS

Submitted: **2001-04-20**, last modified: **2001-08-15**

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.

[179] oai:arXiv.org:astro-ph/0107419 [pdf] - 43790

KL Estimation of the Power Spectrum Parameters from the Angular
Distribution of Galaxies in Early SDSS Data

Comments: Fig 1 postscript problem corrected

Submitted: **2001-07-20**, last modified: **2001-07-26**

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.

[180] oai:arXiv.org:astro-ph/0107417 [pdf] - 43788

The Angular Correlation Function of Galaxies from Early SDSS Data

Comments:

Submitted: **2001-07-20**, last modified: **2001-07-22**

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.

[181] oai:arXiv.org:astro-ph/0107421 [pdf] - 1938807

The 3D Power Spectrum from Angular Clustering of Galaxies in Early SDSS
Data

Comments: 25 pages, 19 figures

Submitted: **2001-07-20**

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.).

[182] oai:arXiv.org:astro-ph/0106120 [pdf] - 42916

PTHalos: A fast method for generating mock galaxy distributions

Comments: 13 pages, 10 figures. (references added)

Submitted: **2001-06-07**, last modified: **2001-06-14**

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.

[183] oai:arXiv.org:astro-ph/0105113 [pdf] - 42342

An excursion set model of hierarchical clustering: Ellipsoidal collapse
and the moving barrier

Comments: 18 pages, submitted to MNRAS

Submitted: **2001-05-07**

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.

[184] oai:arXiv.org:astro-ph/0010137 [pdf] - 38484

On the streaming motions of haloes and galaxies

Comments: 10 pages, submitted to MNRAS

Submitted: **2000-10-06**, last modified: **2000-10-08**

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.

[185] oai:arXiv.org:astro-ph/0009167 [pdf] - 38009

Linear and nonlinear contributions to pairwise peculiar velocities

Comments: 16 pages, 8 figures, submitted to MNRAS 2000 August 15

Submitted: **2000-09-11**

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.

[186] oai:arXiv.org:astro-ph/0009166 [pdf] - 38008

Peculiar velocities of galaxies and clusters

Comments: 18 pages, 18 figures, submitted to MNRAS 2000 July 15

Submitted: **2000-09-11**

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.

[187] oai:arXiv.org:astro-ph/0006319 [pdf] - 36699

How Many Galaxies Fit in a Halo? Constraints on Galaxy Formation
Efficiency from Spatial Clustering

Comments: 32 pages, 10 figures, submitted to ApJ. References added

Submitted: **2000-06-22**, last modified: **2000-06-26**

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.

[188] oai:arXiv.org:astro-ph/0006316 [pdf] - 401120

Halo Dark Matter and Ultra-High Energy Cosmic Rays

Comments: 9 pages, 2 figures. Accepted for Publication in Phys. Lett. B

Submitted: **2000-06-21**

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.

[189] oai:arXiv.org:astro-ph/9907024 [pdf] - 107249

Ellipsoidal collapse and an improved model for the number and spatial
distribution of dark matter haloes

Comments: 12 pages, 6 figures, submitted to MNRAS

Submitted: **1999-07-02**

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.

[190] oai:arXiv.org:astro-ph/9904307 [pdf] - 106204

The correlation of peaks in the microwave background

Comments: 8 pages, 5 postscript figures, submitted to MNRAS

Submitted: **1999-04-22**

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.

[191] oai:arXiv.org:astro-ph/9901122 [pdf] - 104716

Large scale bias and the peak background split

Comments: 8 pages, 6 figures, submitted to MNRAS corrected y-label for fig.4
(newlabel = 1 + oldlabel)

Submitted: **1999-01-11**, last modified: **1999-01-13**

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.

[192] oai:arXiv.org:astro-ph/9808138 [pdf] - 102482

Biasing and the distribution of dark matter haloes

Comments: revised to match MNRAS accepted version (28 pages, 23 figures)

Submitted: **1998-08-13**, last modified: **1999-01-11**

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.

[193] oai:arXiv.org:astro-ph/9805319 [pdf] - 101540

An excursion set model for the distribution of dark matter and dark
matter haloes

Comments: 15 pages, 5 figures, submitted to MNRAS Sept. 1997

Submitted: **1998-05-26**

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.

[194] oai:arXiv.org:astro-ph/9805322 [pdf] - 101543

The forest of merger history trees associated with the formation of dark
matter halos

Comments: 13 pages, 10 figures, submitted to MNRAS

Submitted: **1998-05-26**

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.

[195] oai:arXiv.org:astro-ph/9803281 [pdf] - 100816

Mass growth and density profiles of dark matter halos in hierarchical
clustering

Comments: submitted to MNRAS, 19pages, 8 figures

Submitted: **1998-03-24**

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)$.

[196] oai:arXiv.org:astro-ph/9803045 [pdf] - 100580

The Generalized Poisson distribution and a model of clustering from
Poisson initial conditions

Comments: 9 pages, 2 figures, accepted for publication in MNRAS

Submitted: **1998-03-05**

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.

[197] oai:arXiv.org:astro-ph/9602112 [pdf] - 1234593

Galton--Watson branching processes and the growth of gravitational
clustering

Comments: 13 pages, uuencoded, gzipped, postscript, submitted to MN

Submitted: **1996-02-21**

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.

[198] oai:arXiv.org:astro-ph/9602103 [pdf] - 1234591

The non-linear correlation function and the shapes of virialized halos

Comments: 8 pages, uuencoded, gzipped, postscript, submitted to MN

Submitted: **1996-02-20**

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.

[199] oai:arXiv.org:astro-ph/9511096 [pdf] - 1234517

Constrained realizations and minimum variance reconstruction of
non-Gaussian random fields

Comments: 12 pages, gzipped postscript, MNRAS, in press

Submitted: **1995-11-20**

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.

[200] oai:arXiv.org:astro-ph/9511068 [pdf] - 1234514

The distribution of pairwise peculiar velocities in the nonlinear regime

Comments: (MNRAS, in press) 16 pages, uuencoded

Submitted: **1995-11-15**

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.