We quantify the shape and overdensity of the galaxy distribution in the `Local Sheet' within a sphere of $R=8$ Mpc, and compare these properties with the expectations of the $\Lambda$CDM model. We measure ellipsoidal axis ratios of $c/a\approx0.16$ and $b/a\approx0.79$, indicating that the distribution of galaxies in the Local Volume can be approximated by a flattened oblate ellipsoid, consistent with the `sheet'-like configuration noted in previous studies. In contrast with previous estimates that the Local Sheet has a density close to average, we find that the number density of faint and bright galaxies in the Local Volume is $\approx1.7$ and $\approx5.2$ times denser, respectively, than the mean number density of galaxies of the same luminosity. Comparison with simulations shows that the number density contrasts of bright and faint galaxies within $8$ Mpc alone make the Local Volume a $\approx 2.5\sigma$ outlier in the $\Lambda$CDM cosmology. Our results indicate that the cosmic neighbourhood of the Milky Way may be unusual for galaxies of similar luminosity. The impact of the peculiar properties of our neighbourhood on the properties of the Milky Way and other nearby galaxies is not yet understood and warrants further study.

We present a detailed analysis of the physical processes that cause halo assembly bias -- the dependence of halo clustering on proxies of halo formation time. We focus on the origin of assembly bias in the mass range corresponding to the hosts of typical galaxies and use halo concentration as our chief proxy of halo formation time. We also repeat our key analyses across a broad range of halo masses and for alternative formation time definitions. We show that splashback subhaloes are responsible for two thirds of the assembly bias signal, but do not account for the entire effect. After splashback subhaloes have been removed, we find that the remaining assembly bias signal is due to a relatively small fraction ($\lesssim 10\%$) of haloes in dense regions. We test a number of additional physical processes thought to contribute to assembly bias and demonstrate that the two key processes are the slowing of mass growth by large-scale tidal fields and by the high velocities of ambient matter in sheets and filaments. We also rule out several other proposed physical causes of halo assembly bias. Based on our results, we argue that there are three processes that contribute to assembly bias of low-mass halos: large-scale tidal fields, gravitational heating due to the collapse of large-scale structures, and splashback subhaloes located outside the virial radius.

The splashback radius $R_{\rm sp}$, the apocentric radius of particles on their first orbit after falling into a dark matter halo, has recently been suggested as a physically motivated halo boundary that separates accreting from orbiting material. Using the SPARTA code presented in Paper I, we analyze the orbits of billions of particles in cosmological simulations of structure formation and measure $R_{\rm sp}$ for a large sample of halos that span a mass range from dwarf galaxy to massive cluster halos, reach redshift 8, and include WMAP, Planck, and self-similar cosmologies. We analyze the dependence of $R_{\rm sp}/R_{\rm 200m}$ and $M_{\rm sp}/M_{\rm 200m}$ on the mass accretion rate $\Gamma$, halo mass, redshift, and cosmology. The scatter in these relations varies between 0.02 and 0.1 dex. While we confirm the known trend that $R_{\rm sp}/R_{\rm 200m}$ decreases with $\Gamma$, the relationships turn out to be more complex than previously thought, demonstrating that $R_{\rm sp}$ is an independent definition of the halo boundary that cannot trivially be reconstructed from spherical overdensity definitions. We present fitting functions for $R_{\rm sp}/R_{\rm 200m}$ and $M_{\rm sp}/M_{\rm 200m}$ as a function of accretion rate, peak height, and redshift, achieving an accuracy of 5% or better everywhere in the parameter space explored. We discuss the physical meaning of the distribution of particle apocenters and show that the previously proposed definition of $R_{\rm sp}$ as the radius of the steepest logarithmic density slope encloses roughly three-quarters of the apocenters. Finally, we conclude that no analytical model presented thus far can fully explain our results.

The density field in the outskirts of dark matter halos is discontinuous due to a caustic formed by matter at its first apocenter after infall. In this paper, we present an algorithm to identify the "splashback shell" formed by these apocenters in individual simulated halos using only a single snapshot of the density field. We implement this algorithm in the code SHELLFISH (SHELL Finding In Spheroidal Halos) and demonstrate that the code identifies splashback shells correctly and measures their properties with an accuracy of $<5\%$ for halos with more than 50,000 particles and mass accretion rates of $\Gamma_\textrm{DK14}>0.5$. Using SHELLFISH, we present the first estimates for several basic properties of individual splashback shells, such as radius, $R_\textrm{sp}$, mass, and overdensity, and provide fits to the distribution of these quantities as functions of $\Gamma_\textrm{DK14}$, $\nu_\textrm{200m}$, and $z.$ We confirm previous findings that $R_\textrm{sp}$ decreases with increasing $\Gamma_\textrm{DK14}$, but we show that independent of accretion rate, it also decreases with increasing $\nu_\textrm{200m}$. We also study the 3D structures of these shells and find that they generally have non-ellipsoidal oval shapes. We find that splashback radii estimated by SHELLFISH are $20\%-30\%$ larger than those estimated in previous studies from stacked density profiles at high accretion rates. We demonstrate that the latter are biased low due to the contribution of high-mass subhalos to these profiles and show that using the median instead of mean density in each radial bin mitigates the effect of substructure on density profiles and removes the bias.

SCORCH (Simulations and Constructions of the Reionization of Cosmic Hydrogen) is a new project to study the Epoch of Reionization (EoR). In this first paper, we probe the connection between observed high-redshift galaxies and simulated dark matter halos to better understand the abundance and evolution of the primary source of ionizing radiation. High-resolution N-body simulations are run to quantify the abundance of dark matter halos as a function of mass $M$, accretion rate $\dot{M}$, and redshift $z$. A new fit for the halo mass function $dn/dM$ is $\approx 20\%$ more accurate at the high-mass end where bright galaxies are expected to reside. A novel approach is used to fit the halo accretion rate function $dn/d\dot{M}$ in terms of the halo mass function. Abundance matching against the observed galaxy luminosity function is used to estimate the luminosity-mass relation and the luminosity-accretion-rate relation. The inferred star formation efficiency is not monotonic with $M$ nor $\dot{M}$, but reaches a maximum value at a characteristic mass $\sim 2 \times 10^{11}\ M_\odot$ and a characteristic accretion rate $\sim 6 \times 10^2\ M_\odot/{\rm yr}$ at $z \approx 6$. We find a universal EoR luminosity-accretion-rate relation and construct a fiducial model for the galaxy luminosity function. The Schechter parameters evolve such that $\phi_\star$ decreases, $M_\star$ is more positive (fainter), and $\alpha$ is more negative (steeper) at higher redshifts. We forecast for the upcoming James Webb Space Telescope and show that with apparent magnitude limit $m_{\rm AB} \approx 31\ (32)$, it can observe $\gtrsim 11\ (24)$ unlensed galaxies per square degree per unit redshift at least down to $M_\star$ at $z \lesssim 13\ (14)$.