We present a public code to generate random fields with an arbitrary probability distribution function (PDF) and an arbitrary correlation function. The algorithm is cosmology-independent, applicable to any stationary stochastic process over a three dimensional grid. We implement it in the case of the matter density field, showing its benefits over the lognormal approximation, which is often used in cosmology for generation of mock catalogues. We find that the covariance of the power spectrum from the new fast realizations is more accurate than that from a lognormal model. As a proof of concept, we also apply the new simulation scheme to the divergence of the Lagrangian displacement field. We find that information from the correlation function and the PDF of the displacement-divergence provides modest improvement over other standard analytical techniques to describe the particle field in the simulation. This suggests that further progress in this direction should come from multi-scale or non-local properties of the initial matter distribution.

Matter in the Universe is distributed in a cosmic web, with a filament of matter typically connecting neighboring galaxy pairs, separated by tens of millions of light-years. Filaments are known to influence the spins of galaxies around them, but it remains unknown whether filaments themselves spin. Here, we measure the angular momentum around intergalactic filaments in a dark-matter cosmological simulation, finding that matter generally rotates substantially around them. Our results suggest that spin is a fundamental property of filaments, part of a full picture of structure formation. Filaments are likely to be the longest coherently rotating objects in the Universe. This finding helps to explain how galaxies spin in the cosmic web, and could provide a mechanism to generate intergalactic magnetic fields.

Dark matter as a Bose-Einstein condensate, such as the axionic scalar field particles of String Theory, can explain the coldness of dark matter on large scales. Pioneering simulations in this context predict a rich wave-like structure, with a ground state soliton core in every galaxy surrounded by a halo of excited states that interfere on the de Broglie scale. This de Broglie scale is largest for low mass galaxies as momentum is lower, providing a simple explanation for the wide cores of dwarf spheroidal galaxies. Here we extend these "wave dark matter" ($\psi$DM) predictions to the newly discovered class of "Ultra Diffuse Galaxies" (UDG) that resemble dwarf spheroidal galaxies but with more extended stellar profiles. Currently the best studied example, DF44, has a uniform velocity dispersion of $\simeq 33$km/s, extending to a least 3 kpc, that we show is reproduced by our $\psi$DM simulations with a soliton radius of $\simeq 0.5$ kpc. In the $\psi$DM context, we show the relatively flat dispersion profile of DF44 lies between massive galaxies with compact dense solitons, as may be present in the Milky Way on a scale of 100pc and lower mass galaxies where the velocity dispersion declines centrally within a wide, low density soliton, like Antlia II, of radius 3Kpc. In contrast, standard CDM requires excessive tangential stellar motions predominate to counter the inherent central cusp which would otherwise enhance the central velocity dispersion.

The standard explanation for galaxy spin starts with the tidal-torque theory (TTT), in which an ellipsoidal dark-matter protohalo, which comes to host the galaxy, is torqued up by the tidal gravitational field around it. We discuss a complementary picture, using the relatively familiar velocity field, instead of the tidal field, whose intuitive connection to the surrounding, possibly faraway matter arrangement is more obscure. In this 'spin from primordial inner motions' (SPIM) concept, implicit in TTT derivations but not previously emphasized, the angular momentum from the gravity-sourced velocity field inside a protohalo largely cancels out, but has some excess from the aspherical outskirts. At first, the net spin scales according to linear theory, a sort of comoving conservation of familiar angular momentum. Then, at collapse, it is conserved in physical coordinates. Small haloes are then typically subject to secondary exchanges of angular momentum. The TTT is useful for analytic estimates. But a literal interpretation of the TTT is inaccurate in detail, without some implicit concepts about smoothing of the velocity and tidal fields. This could lead to misconceptions, for those first learning about how galaxies come to spin. Protohaloes are not perfectly ellipsoidal and do not uniformly torque up, as in a naive interpretation of the TTT; their inner velocity fields retain substantial dispersion. Furthermore, quantitatively, given initial conditions and protohalo boundaries, SPIM is more direct and accurate than the TTT to predict halo spins. We also discuss how SPIM applies to rotating filaments, and the relation between halo mass and spin, in which the total spin of a halo can be thought of as a sum of random contributions.

Here we propose a new scheme to reconstruct the baryon acoustic oscillations (BAO) signal, with key cosmological information, based on deep convolutional neural networks. After training the network with almost no fine-tuning, in the test set, the network recovers large-scale modes accurately: the correlation coefficient between the ground truth and recovered initial conditions still reach $90\%$ at $k \leq 0.2~ h\mathrm{Mpc}^{-1}$, which significantly improves the BAO signal-to-noise ratio until the scale $k=0.4~ h\mathrm{Mpc}^{-1}$. Furthermore, our scheme is independent of the survey boundary since it reconstructs initial condition based on local density distribution in configuration space, which means that we can gain more information from the whole survey space. Finally, we found our trained network is not sensitive to the cosmological parameters and works very well in those cosmologies close to that of our training set. This new scheme will possibly help us dig out more information from the current, on-going and future galaxy surveys.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

We propose the Cosmic Web Detachment (CWD) model, a framework to interpret the star-formation history of galaxies in a cosmological context. The CWD model unifies several starvation mechanisms known to disrupt or stop star formation into one single physical framework. Galaxies begin accreting star-forming gas at early times via a network of primordial filaments, simply related to the pattern of density fluctuations in the initial conditions. But when shell-crossing occurs on intergalactic scales, this pattern is disrupted, and the galaxy detaches from its primordial filaments, ending the accretion of cold gas. We argue that CWD encompasses known external processes halting star formation, such as harassment, strangulation and starvation. On top of these external processes, internal feedback processes such as AGN contribute to stop in star formation as well. By explicitly pointing out the non-linear nature of CWD events we introduce a simple formalism to identify CWD events in N-body simulations. With it we reproduce and explain, in the context of CWD, several observations including downsizing, the cosmic star formation rate history, the galaxy mass-color diagram and the dependence of the fraction of red galaxies with mass and local density.

The cosmic web (the arrangement of matter in the universe), spider's webs, and origami tessellations are linked by their geometry (specifically, of sectional-Voronoi tessellations). This motivates origami and textile artistic representations of the cosmic web. It also relates to the scientific insights origami can bring to the cosmic web; we show results of some cosmological computer simulations, with some origami-tessellation properties. We also adapt software developed for cosmic-web research to provide an interactive tool for general origami-tessellation design.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

The most common statistic used to analyze large-scale structure surveys is the correlation function, or power spectrum. Here, we show how `slicing' the correlation function on local density brings sensitivity to interesting non-Gaussian features in the large-scale structure, such as the expansion or contraction of baryon acoustic oscillations (BAO) according to the local density. The sliced correlation function measures the large-scale flows that smear out the BAO, instead of just correcting them as reconstruction algorithms do. Thus, we expect the sliced correlation function to be useful in constraining the growth factor, and modified gravity theories that involve the local density. Out of the studied cases, we find that the run of the BAO peak location with density is best revealed when slicing on a $\sim 40$ Mpc/$h$ filtered density. But slicing on a $\sim100$ Mpc/$h$ filtered density may be most useful in distinguishing between underdense and overdense regions, whose BAO peaks are separated by a substantial $\sim 5$ Mpc/$h$ at $z=0$. We also introduce `curtain plots' showing how local densities drive particle motions toward or away from each other over the course of an $N$-body simulation.

Cosmological $N$-body simulations play a vital role in studying models for the evolution of the Universe. To compare to observations and make a scientific inference, statistic analysis on large simulation datasets, e.g., finding halos, obtaining multi-point correlation functions, is crucial. However, traditional in-memory methods for these tasks do not scale to the datasets that are forbiddingly large in modern simulations. Our prior paper proposes memory-efficient streaming algorithms that can find the largest halos in a simulation with up to $10^9$ particles on a small server or desktop. However, this approach fails when directly scaling to larger datasets. This paper presents a robust streaming tool that leverages state-of-the-art techniques on GPU boosting, sampling, and parallel I/O, to significantly improve performance and scalability. Our rigorous analysis of the sketch parameters improves the previous results from finding the centers of the $10^3$ largest halos to $\sim 10^4-10^5$, and reveals the trade-offs between memory, running time and number of halos. Our experiments show that our tool can scale to datasets with up to $\sim 10^{12}$ particles while using less than an hour of running time on a single GPU Nvidia GTX 1080.

For over twenty years, the term 'cosmic web' has guided our understanding of the large-scale arrangement of matter in the cosmos, accurately evoking the concept of a network of galaxies linked by filaments. But the physical correspondence between the cosmic web and structural-engineering or textile 'spiderwebs' is even deeper than previously known, and extends to origami tessellations as well. Here we explain that in a good structure-formation approximation known as the adhesion model, threads of the cosmic web form a spiderweb, i.e. can be strung up to be entirely in tension. The correspondence is exact if nodes sampling voids are included, and if structure is excluded within collapsed regions (walls, filaments and haloes), where dark-matter multistreaming and baryonic physics affect the structure. We also suggest how concepts arising from this link might be used to test cosmological models: for example, to test for large-scale anisotropy and rotational flows in the cosmos.

The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public.

In view of future high-precision large-scale structure surveys, it is important to quantify the percent and subpercent level effects in cosmological $N$-body simulations from which theoretical predictions are drawn. One such effect involves deciding whether to zero all modes above the one-dimensional Nyquist frequency, the so-called "corner" modes, in the initial conditions. We investigate this effect by comparing power spectra, density distribution functions, halo mass functions, and halo profiles in simulations with and without these modes. For a simulation with a mass resolution of $m_p \sim 10^{11}\,h^{-1}\,M_{\odot}$, we find that at $z>6$, the difference in the matter power spectrum is large at wavenumbers above $\sim 80$\% of $k_{\rm{Ny}}$, reducing to below 2\% at all scales by $z\sim 3$. Including corner modes results in a better match between low- and high-resolution simulations at wavenumbers around the Nyquist frequency of the low-resolution simulation, but the effect of the corner modes is smaller than the effect of particle discreteness. The differences in mass functions are 3\% for the smallest halos at $z=6$ for the $m_p \sim 10^{11}\,h^{-1}\,M_{\odot}$ simulation, but we find no significant difference in the stacked profiles of well-resolved halos at $z \leq 6$. Thus removing power at $|\mathbf{k}|>k_{\rm{Ny}}$ in the initial conditions of cosmological simulations has a small effect on small scales and high redshifts, typically below a few percent.

We present a Bayesian phase-space reconstruction of the cosmic large-scale matter density and velocity fields from the SDSS-III Baryon Oscillations Spectroscopic Survey Data Release 12 (BOSS DR12) CMASS galaxy clustering catalogue. We rely on a given $\Lambda$CDM cosmology, a mesh resolution in the range of 6-10 $h^{-1}$ Mpc, and a lognormal-Poisson model with a redshift dependent nonlinear bias. The bias parameters are derived from the data and a general renormalised perturbation theory approach. We use combined Gibbs and Hamiltonian sampling, implemented in the \textsc{argo} code, to iteratively reconstruct the dark matter density field and the coherent peculiar velocities of individual galaxies, correcting hereby for coherent redshift space distortions (RSD). Our tests relying on accurate $N$-body based mock galaxy catalogues, show unbiased real space power spectra of the nonlinear density field up to $k\sim0.2\, h$ Mpc$^{-1}$, and vanishing quadrupoles down to $r\sim20\,h^{-1}$ Mpc. We also demonstrate that the nonlinear cosmic web can be obtained from the tidal field tensor based on the Gaussian component of the reconstructed density field. We find that the reconstructed velocities have a statistical correlation coefficient compared to the true velocities of each individual lightcone mock galaxy of $r\sim0.68$ including about 10% of satellite galaxies with virial motions (about $r=0.75$ without satellites). The power spectra of the velocity divergence agree well with theoretical predictions up to $k\sim0.2\,h\,{\rm Mpc}^{-1}$. This work will be especially useful to improve, e.g. BAO reconstructions, kinematic Sunyaev-Zeldovich (kSZ), integrated Sachs-Wolfe (ISW) measurements, or environmental studies.

We have searched for the signature of cosmic voids in the CMB, in both the Planck temperature and lensing-convergence maps; voids should give decrements in both. We use zobov voids from the DR12 SDSS CMASS galaxy sample. We base our analysis on N-body simulations, to avoid a posteriori bias. For the first time, we detect the signature of voids in CMB lensing: the significance is $3.2\sigma$, close to $\Lambda$CDM in both amplitude and projected density-profile shape. A temperature dip is also seen, at modest significance ($2.3\sigma$), with amplitude about 6 times the prediction. This temperature signal is induced mostly by voids with radius between 100 and 150 Mpc/h, while the lensing signal is mostly contributed by smaller voids -- as expected; lensing relates directly to density, while ISW depends on gravitational potential. The void abundance in observations and simulations agree, as well. We also repeated the analysis excluding lower-significance voids: no lensing signal is detected, with an upper limit of about twice the $\Lambda$CDM prediction. But the mean temperature decrement now becomes non-zero at the $3.7\sigma$ level (similar to that found by Granett et al.), with amplitude about 20 times the prediction. However, the observed dependence of temperature on void size is in poor agreement with simulations, whereas the lensing results are consistent with $\Lambda$CDM theory. Thus, the overall tension between theory and observations does not favour non-standard theories of gravity, despite the hints of an enhanced amplitude for the ISW effect from voids.

We apply the Alcock-Paczynski (AP) test to the stacked voids identified using the large-scale structure galaxy catalog from the Baryon Oscillation Spectroscopic Survey (BOSS). This galaxy catalog is part of the Sloan Digital Sky Survey (SDSS) Data Release 12 and is the final catalog of SDSS-III. We also use 1000 mock galaxy catalogs that match the geometry, density, and clustering properties of the BOSS sample in order to characterize the statistical uncertainties of our measurements and take into account systematic errors such as redshift space distortions. For both BOSS data and mock catalogs, we use the ZOBOV algorithm to identify voids, we stack together all voids with effective radii of 30-100Mpc/h in the redshift range 0.43-0.7, and we accurately measure the shape of the stacked voids. Our tests with the mock catalogs show that we measure the stacked void ellipticity with a statistical precision of 2.6%. We find that the stacked voids in redshift space are slightly squashed along the line of sight, which is consistent with previous studies. We repeat this measurement of stacked void shape in the BOSS data assuming several values of Omega_m within the flat LCDM model, and we compare to the mock catalogs in redshift space in order to perform the AP test. We obtain a constraint of $\Omega_m = 0.38^{+0.18}_{-0.15}$ at the 68% confidence level from the AP test. We discuss the various sources of statistical and systematic noise that affect the constraining power of this method. In particular, we find that the measured ellipticity of stacked voids scales more weakly with cosmology than the standard AP prediction, leading to significantly weaker constraints. We discuss how AP constraints will improve in future surveys with larger volumes and densities.

We discuss an idealized model of halo formation, in which a collapsing halo node is tetrahedral, with a filament extruding from each of its four faces, and with a wall connecting each pair of filaments. In the model, filaments generally spin when they form, and the halo spins if and only if there is some rotation in filaments. This is the simplest-possible fully three-dimensional halo collapse in the 'origami approximation,' in which voids are irrotational, and the dark-matter sheet out of which dark-matter structures form is allowed to fold in position-velocity phase space, but not stretch (i.e., it cannot vary in density along a stream). Up to an overall scaling, the four filament directions, and only three other quantities, such as filament spins, suffice to determine all of the collapse's properties: the shape, mass, and spin of the halo; the densities per unit length and spins of all filaments; and masses per unit area of the walls. If the filaments are arranged regular-tetrahedrally, filament properties obey simple laws, reminiscent of angular-momentum conservation. The model may be most useful in understanding spin correlations between neighboring galaxies joined by filaments; these correlations would give intrinsic alignments between galaxies, essential to understand for accurate cosmological weak-lensing measurements.

We present a cosmic void catalog using the large-scale structure galaxy catalog from the Baryon Oscillation Spectroscopic Survey (BOSS). This galaxy catalog is part of the Sloan Digital Sky Survey (SDSS) Data Release 12 and is the final catalog of SDSS-III. We take into account the survey boundaries, masks, and angular and radial selection functions, and apply the ZOBOV void finding algorithm to the galaxy catalog. We identify a total of 10,643 voids. After making quality cuts to ensure that the voids represent real underdense regions, we obtain 1,228 voids with effective radii spanning the range 20-100Mpc/h and with central densities that are, on average, 30% of the mean sample density. We release versions of the catalogs both with and without quality cuts. We discuss the basic statistics of voids, such as their size and redshift distributions, and measure the radial density profile of the voids via a stacking technique. In addition, we construct mock void catalogs from 1000 mock galaxy catalogs, and find that the properties of BOSS voids are in good agreement with those in the mock catalogs. We compare the stellar mass distribution of galaxies living inside and outside of the voids, and find no significant difference. These BOSS and mock void catalogs are useful for a number of cosmological and galaxy environment studies.

We study the parity-odd part (that we shall call Doppler term) of the linear galaxy two-point correlation function that arises from wide-angle, velocity, Doppler lensing and cosmic acceleration effects. As it is important at low redshift and at large angular separations, the Doppler term is usually neglected in the current generation of galaxy surveys. For future wide-angle galaxy surveys such as Euclid, SPHEREx and SKA, however, we show that the Doppler term must be included. The effect of these terms is dominated by the magnification due to relativistic aberration effects and the slope of the galaxy redshift distribution and it generally mimics the effect of the local type primordial non-Gaussianity with the effective nonlinearity parameter $f_{\rm NL}^{\rm eff}$ of a few, we show that this would affect forecasts on measurements of $f_{\rm NL}$ at low-redshift. Our results show that a survey at low redshift with large number density over a wide area of the sky could detect the Doppler term with a signal-to-noise ratio of $\sim 1-20$, depending on survey specifications.

We reproduce the galaxy clustering catalogue from the SDSS-III Baryon Oscillation Spectroscopic Survey Final Data Release (BOSS DR11&DR12) with high fidelity on all relevant scales in order to allow a robust analysis of baryon acoustic oscillations and redshift space distortions. We have generated (6,000) 12,288 MultiDark PATCHY BOSS (DR11) DR12 light-cones corresponding to an effective volume of $\sim192,000\,[h^{-1}\,{\rm Gpc}]^3$ (the largest ever simulated volume), including cosmic evolution in the redshift range from 0.15 to 0.75. The mocks have been calibrated using a reference galaxy catalogue based on the halo abundance matching modelling of the BOSS DR11&DR12 galaxy clustering data and on the data themselves. The production follows three steps. First, we apply the PATCHY code to generate a dark matter field and an object distribution including nonlinear stochastic galaxy bias. Secondly, we run the halo/stellar distribution reconstruction HADRON code to assign masses to the various objects. This step uses the mass distribution as a function of local density and non-local indicators (i.e., tidal field tensor eigenvalues and relative halo exclusion separation for massive objects) from the reference simulation applied to the corresponding PATCHY dark matter and galaxy distribution. Finally, we apply the SUGAR code to build the light cones. The resulting MultiDark PATCHY mock light cones reproduce the number density, selection function, survey geometry, and in general within 1 $\sigma$, for arbitrary stellar mass bins, the power spectrum up to $k=0.3\,h\,{\rm Mpc}^{-1}$, the two-point correlation functions down to a few Mpc scales, and the three-point statistics of the BOSS DR11&DR12 galaxy samples.

The Gaussianization transform has been proposed as a method to remove the issues of scale-dependent galaxy bias and nonlinearity from galaxy clustering statistics, but these benefits have yet to be thoroughly tested for realistic galaxy samples. In this paper, we test the effectiveness of the Gaussianization transform for different galaxy types by applying it to realistic simulated blue and red galaxy samples. We show that in real space, the shapes of the Gaussianized power spectra of both red and blue galaxies agree with that of the underlying dark matter, with the initial power spectrum, and with each other to smaller scales than do the statistics of the usual (untransformed) density field. However, we find that the agreement in the Gaussianized statistics breaks down in redshift space. We attribute this to the fact that red and blue galaxies exhibit very different fingers of god in redshift space. After applying a finger-of-god compression, the agreement on small scales between the Gaussianized power spectra is restored. We also compare the Gaussianization transform to the clipped galaxy density field and find that while both methods are effective in real space, they have more complicated behaviour in redshift space. Overall, we find that Gaussianization can be useful in recovering the shape of the underlying dark matter power spectrum to k ~ 0.5 h/Mpc and of the initial power spectrum to k ~ 0.4 h/Mpc in certain cases at z = 0.

Here we present a simple, parameter-free, non-perturbative algorithm that gives low-redshift cosmological particle realizations accurate to few-Megaparsec scales, called MUSCLE (MUltiscale Spherical ColLapse Evolution). It has virtually the same cost as producing N-body-simulation initial conditions, since it works with the 'stretch' parameter {\psi}, the Lagrangian divergence of the displacement field. It promises to be useful in quickly producing mock catalogs, and to simplify computationally intensive reconstructions of galaxy surveys. MUSCLE applies a spherical-collapse prescription on multiple Gaussian-smoothed scales. It achieves higher accuracy than perturbative schemes (Zel'dovich and 2LPT), and, by including the void-in-cloud process (voids in large-scale collapsing regions), solves problems with a single-scale spherical-collapse scheme. Slight further improvement is possible by mixing in the 2LPT estimate on large scales. Additionally, we show the behavior of {\psi} for different morphologies (voids, walls, filaments, and haloes). A Python code to produce these realizations is available at http://skysrv.pha.jhu.edu/~neyrinck/muscle.html.

We propose a common terminology for use in describing both temporal merger trees and spatial structure trees for dark-matter halos. We specify a unified data format in HDF5 and provide example I/O routines in C, FORTRAN and PYTHON.

We compare analytical predictions of void volume functions to those measured from N-body simulations, detecting voids with the zobov void finder. We push to very small, nonlinear voids, below few Mpc radius, by considering the unsampled DM density field. We also study the case where voids are identified using halos. We develop analytical formula for the void abundance of both the excursion set approach and the peaks formalism. These formula are valid for random walks smoothed with a top-hat filter in real space, with a large class of realistic barrier models. We test the extent to which the spherical evolution approximation, which forms the basis of the analytical predictions, models the highly aspherical voids that occur in the cosmic web, and are found by a watershed-based algorithm such as zobov. We show that the volume function returned by zobov is quite sensitive to the choice of treatment of sub-voids, a fact that has not been appreciated previously. For reasonable choices of sub-void exclusion, we find that the Lagrangian density delta_v of the zobov voids -- which is predicted to be a constant delta_v = -2.7 in the spherical evolution model -- is different from the predicted value, showing substantial scatter and scale dependence. This result applies to voids identified at z=0 with effective radius between 1 and 10 Mpc/h. Our analytical approximations are flexible enough to give a good description of the resulting volume function; however, this happens for choices of parameter values that are different from those suggested by the spherical evolution assumption. We conclude that analytical models for voids must move away from the spherical approximation in order to be applied successfully to observations, and we discuss some possible ways forward.

In the current cosmological paradigm, an initially flat three-dimensional manifold that pervades space (the `dark-matter sheet') folds up to build concentrations of mass (galaxies), and a cosmic web between them. Galaxies are nodes, connected by a network of filaments and walls. The folding is in six-dimensional (3D position, plus 3D velocity) phase space. The positions of creases, or caustics, mark the edges of structures. Here, I introduce an origami approximation to cosmological structure formation, in which the dark-matter sheet is not allowed to stretch. But it still produces an idealized cosmic web, with nodes, filaments, walls and voids. In 2D, nodes form in `polygonal collapse' (a twist-fold in origami), necessarily generating filaments simultaneously. In 3D, nodes form in `polyhedral collapse,' simultaneously generating filaments and walls. The masses, spatial arrangement, and angular momenta of nodes and filaments are related in the model. I describe some `tetrahedral collapse', or tetrahedral twist-fold, models.

We discuss the Kolmogorov complexity of primordial patches that collapse to form galaxies like the Milky Way; this complexity quantifies the amount of initial data available to form the structure. We also speculate on how the quantity changes with time. Because of dark-matter and baryonic collapse processes, it likely decreases with time, i.e.\ information sinks dominate sources. But sources of new random information do exist; e.g., a central black hole with an accretion disk and jets could in principle broadcast small-scale quantum fluctuations over a substantial portion of a galaxy. A speculative example of how this concept might be useful is in differentiating between warm (WDM) and cold (CDM) dark matter. With WDM, the initial patch that formed the Milky Way would have had few features, making the present high degree of structure a curiosity. The primordial patch would have had only several billion independent information-carrying `pixels' if the WDM particle had a mass of 1 keV. This number of `pixels' is much less than even the number of stars in the Milky Way. If the dark matter is proven to be warm, the high degree of structure in the Milky Way could have arisen in two ways: (1) from a high sensitivity to initial conditions, like an intricate fractal arising from a relatively simple computer code; or (2) from random information generated after the galaxy formed, i.e.\ not entirely deterministically from the initial conditions.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

Warm dark matter (WDM) has been proposed as an alternative to cold dark matter (CDM), to resolve issues such as the apparent lack of satellites around the Milky Way. Even if WDM is not the answer to observational issues, it is essential to constrain the nature of the dark matter. The effect of WDM on haloes has been extensively studied, but the small-scale initial smoothing in WDM also affects the present-day cosmic web and voids. It suppresses the cosmic "sub-web" inside voids, and the formation of both void haloes and subvoids. In $N$-body simulations run with different assumed WDM masses, we identify voids with the ZOBOV algorithm, and cosmic-web components with the ORIGAMI algorithm. As dark-matter warmth increases (i.e., particle mass decreases), void density minima grow shallower, while void edges change little. Also, the number of subvoids decreases. The density field in voids is particularly insensitive to baryonic physics, so if void density profiles and minima could be measured observationally, they would offer a valuable probe of the nature of dark matter. Furthermore, filaments and walls become cleaner, as the substructures in between have been smoothed out; this leads to a clear, mid-range peak in the density PDF.

We study the nature of voids defined as single-stream regions that have not undergone shell-crossing. We use ORIGAMI to determine the cosmic web morphology of each dark matter particle in a suite of cosmological $N$-body simulations, which explicitly calculates whether a particle has crossed paths with others along multiple sets of axes and does not depend on a parameter or smoothing scale. The theoretical picture of voids is that of expanding underdensities with borders defined by shell-crossing. We find instead that locally underdense single-stream regions are not bounded on all sides by multi-stream regions, thus they percolate, filling the simulation volume; we show that the set of multi-stream particles also percolates. This percolation persists to high resolution, where the mass fraction of single-stream voids is low, because the volume fraction remains high; we speculate on the fraction of collapsed mass in the continuum limit of infinite resolution. By introducing a volume threshold parameter to define underdense void "cores", we create a catalog of ORIGAMI voids which consist entirely of single-stream particles and measure their percolation properties, volume functions, and average densities.

The powerful Lagrangian view of structure formation was essentially introduced to cosmology by Zel'dovich. In the current cosmological paradigm, a dark-matter-sheet 3D manifold, inhabiting 6D position-velocity phase space, was flat (with vanishing velocity) at the big bang. Afterward, gravity stretched and bunched the sheet together in different places, forming a cosmic web when projected to the position coordinates. Here, I explain some properties of an origami approximation, in which the sheet does not stretch or contract (an assumption that is false in general), but is allowed to fold. Even without stretching, the sheet can form an idealized cosmic web, with convex polyhedral voids separated by straight walls and filaments, joined by convex polyhedral nodes. The nodes form in 'polygonal' or 'polyhedral' collapse, somewhat like spherical/ellipsoidal collapse, except incorporating simultaneous filament and wall formation. The origami approximation allows phase-space geometries of nodes, filaments, and walls to be more easily understood, and may aid in understanding spin correlations between nearby galaxies. This contribution explores kinematic origami-approximation models giving velocity fields for the first time.

The star formation history of galaxies is a complex process usually considered to be stochastic in nature, for which we can only give average descriptions such as the color-density relation. In this work we follow star-forming gas particles in a hydrodynamical N-body simulation back in time in order to study their initial spatial configuration. By keeping record of the time when a gas particle started forming stars we can produce gas-star isochrone surfaces delineating the surfaces of accreting gas that begin producing stars at different times. These accretion surfaces are closely packed inside dense regions, intersecting each other, and as a result galaxies inside proto-clusters stop accreting gas early, naturally explaining the color dependence on density. The process described here has a purely gravitational / geometrical origin, arguably operating at a more fundamental level than complex processes such as AGN and supernovae, and providing a conceptual origin for the color-density relation.

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.

As an alternative way of describing the cosmological velocity field, we discuss the evolution of rotational invariants constructed from the velocity gradient tensor. Compared with the traditional divergence-vorticity decomposition, these invariants, defined as coefficients of characteristic equation of the velocity gradient tensor, enable a complete classification of all possible flow patterns in the dark-matter comoving frame, including both potential and vortical flows. We show that this tool, first introduced in turbulence two decades ago, proves to be very useful in understanding the evolution of the cosmic web structure, and in classifying its morphology. Before shell-crossing, different categories of potential flow are highly associated with cosmic web structure, because of the coherent evolution of density and velocity. This correspondence is even preserved at some level when vorticity is generated after shell-crossing. The evolution from the potential to vortical flow can be traced continuously by these invariants. With the help of this tool, we show that the vorticity is generated in a particular way that is highly correlated with the large-scale structure. This includes a distinct spatial distribution and different types of alignment between cosmic web and vorticity direction for various vortical flows. Incorporating shell-crossing into closed dynamical systems is highly non-trivial, but we propose a possible statistical explanation for some of these phenomena relating to the internal structure of the three-dimensional invariants space.

Gaussianizing transformations are used statistically in many non-cosmological fields, but in cosmology, we are only starting to apply them. Here I explain a strategy of analyzing the 1-point function (PDF) of a spatial field, together with the 'essential' clustering statistics of the Gaussianized field, which are invariant to a local transformation. In cosmology, if the tracer sampling is sufficient, this achieves two important goals. First, it can greatly multiply the Fisher information, which is negligible on nonlinear scales in the usual $\delta$ statistics. Second, it decouples clustering statistics from a local bias description for tracers such as galaxies.

We study the relationship between dark-matter haloes and matter in the MIP $N$-body simulation ensemble, which allows precision measurements of this relationship, even deeply into voids. What enables this is a lack of discreteness, stochasticity, and exclusion, achieved by averaging over hundreds of possible sets of initial small-scale modes, while holding fixed large-scale modes that give the cosmic web. We find (i) that dark-matter-halo formation is greatly suppressed in voids; there is an exponential downturn at low densities in the otherwise power-law matter-to-halo density bias function. Thus, the rarity of haloes in voids is akin to the rarity of the largest clusters, and their abundance is quite sensitive to cosmological parameters. The exponential downturn appears both in an excursion-set model, and in a model in which fluctuations evolve in voids as in an open universe with an effective $\Omega_m$ proportional to a large-scale density. We also find that (ii) haloes typically populate the average halo-density field in a super-Poisson way, i.e. with a variance exceeding the mean; and (iii) the rank-order-Gaussianized halo and dark-matter fields are impressively similar in Fourier space. We compare both their power spectra and cross-correlation, supporting the conclusion that one is roughly a strictly-increasing mapping of the other. The MIP ensemble especially reveals how halo abundance varies with `environmental' quantities beyond the local matter density; (iv) we find a visual suggestion that at fixed matter density, filaments are more populated by haloes than clusters.

We study the late-time Integrated Sachs-Wolfe (ISW) effect in $f(R)$ gravity using N-body simulations. In the $f(R)$ model under study, the linear growth rate is larger than that in general relativity (GR). This slows down the decay of the cosmic potential and induces a smaller ISW effect on large scales. Therefore, the $\dot\Phi$ (time derivative of the potential) power spectrum at $k<0.1h$/Mpc is suppressed relative to that in GR. In the non-linear regime, relatively rapid structure formation in $f(R)$ gravity boosts the non-linear ISW effect relative to GR, and the $\dot\Phi$ power spectrum at $k>0.1h$/Mpc is increased (100$\%$ greater on small scales at $z=0$). We explore the detectability of the ISW signal via stacking supercluster and supervoids. The differences in the corresponding ISW cold or hot spots are $\sim 20\%$ for structures of $\sim 100$Mpc/$h$. Such differences are greater for smaller structures, but the amplitude of the signal is lower. The high amplitude of ISW signal detected by Granett et al. can not explained in the $f(R)$ model. We find relatively big differences between $f(R)$ and GR in the transverse bulk motion of matter, and discuss its detectability via the relative frequency shifts of photons from multiple lensed images.

We measure the average temperature decrement on the cosmic microwave background (CMB) produced by voids selected in the SDSS DR7 spectroscopic redshift galaxy catalog, spanning redshifts $0<z<0.44$. We find an imprint of amplitude between 2.6 and 2.9$\mu K$ as viewed through a compensated top-hat filter scaled to the radius of each void; we assess the statistical significance of the imprint at ~2$\sigma$. We make crucial use of $N$-body simulations to calibrate our analysis. As expected, we find that large voids produce cold spots on the CMB through the Integrated Sachs-Wolfe (ISW) effect. However, we also find that small voids in the halo density field produce hot spots, because they reside in contracting, larger-scale overdense regions. This is an important effect to consider when stacking CMB imprints from voids of different radius. We have found that the same filter radius that gives the largest ISW signal in simulations also yields close to the largest detected signal in the observations. However, although it is low in significance, our measured signal is much higher-amplitude than expected from ISW in the concordance $\Lambda$CDM universe. The discrepancy is also at the ~2$\sigma$ level. We have demonstrated that our result is robust against the varying of thresholds over a wide range.

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

In the context of the study of the Integrated Sachs Wolfe effect (ISW), we construct a template of the projected density distribution up to $z\simeq 0.7$ by using the Luminous Galaxies (LGs) from the Sloan Digital Sky Survey DR8. We use a photo-z catalogue trained with more than a hundred thousand galaxies from BOSS in the SDSS DR8 imaging area. We consider two different LG samples whose selection matches that of SDSS-III/BOSS: the LOWZ sample ($z\in [0.15,0.5]$) and the CMASS sample ($z\in[0.4,0.7]$). When building the LG density maps we use the information from star density, survey footprint, seeing conditions, sky emission, dust extinction and airmass to explore the impact of these artifacts on the two LG samples. In agreement with previous studies, we find that the CMASS sample is particularly sensitive to Galactic stars, which dominate the contribution to the auto-angular power spectrum below $\ell=7$. Other potential systematics affect mostly the low multipole range ($\ell\in[2,7]$), but leave fluctuations on smaller scales practically unchanged. The resulting power spectra in the multipole range $\ell\in[2,100]$ for the LOWZ, CMASS and LOWZ+CMASS samples are compatible with linear $\Lambda$CDM expectations and constant bias values of $b=1.98 \pm 0.11$, $2.08\pm0.14$ and $1.88\pm 0.11$, respectively, with no traces of non-Gaussianity: $f_{\rm NL}^{\rm local}=59\pm 75$ at 95% confidence level for the full LOWZ+CMASS sample in the range $\ell\in[4,100]$. After cross-correlating WMAP-9yr data with the LOWZ+CMASS LG density field, the ISW signal is detected at the level of 1.62--1.69$\,\sigma$. While this result is in close agreement with predictions from Monte Carlo simulations in the concordance $\Lambda$CDM model, it cannot rule out by itself an Einstein-de Sitter scenario, and has a moderately low signal compared to previous studies conducted on subsets of this LG sample.

Structures like galaxies and filaments of galaxies in the Universe come about from the origami-like folding of an initially flat three-dimensional manifold in 6D phase space. The ORIGAMI method identifies these structures in a cosmological simulation, delineating the structures according to their outer folds. Structure identification is a crucial step in comparing cosmological simulations to observed maps of the Universe. The ORIGAMI definition is objective, dynamical and geometric: filament, wall and void particles are classified according to the number of orthogonal axes along which dark-matter streams have crossed. Here, we briefly review these ideas, and speculate on how ORIGAMI might be useful to find cosmic voids.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

We present an experiment in which we 'ring' a set of cosmological N-body-simulation initial conditions, placing spikes in the initial power spectrum at different wavenumber bins. We then measure where these spikes end up in the final conditions. In the usual overdensity power spectrum, most sensitive to contracting and collapsing dense regions, initial power on slightly non-linear scales (k ~ 0.3 h/Mpc) smears to smaller scales, coming to dominate the initial power once there. Log-density and Gaussianized-density power spectra, sensitive to low-density (expanding) and high-density regions, respond differently: initial spikes spread symmetrically in scale, both upward and downward. In fact, in the power spectrum of 1/(1 + {\delta}), spikes migrate to larger scales, showing the magnifying effect of voids on small-scale modes. These power spectra show much greater sensitivity to small-scale initial features. We also test the difference between an approximation of the Ly-{\alpha} flux field, and its Gaussianized form, and give a toy model that qualitatively explains the symmetric power spreading in Gaussianized-density power spectra. Also, we discuss how to use this framework to estimate power-spectrum covariance matrices. This can be used to track the fate of information in the Universe, that takes the form of initial degrees of freedom, one random spike per initial mode.

Large-scale bulk flows in the Universe distort the initial density field, broadening the baryon-acoustic-oscillation (BAO) feature that was imprinted when baryons were strongly coupled to photons. Additionally, there is a small shift inward in the peak of the conventional overdensity correlation function, a mass-weighted statistic. This shift occurs when high density peaks move toward each other. We explore whether this shift can be removed by applying to the density field a transform (such as a logarithm) that gives fairer statistical weight to fluctuations in underdense regions. Using configuration-space perturbation theory in the Zel'dovich approximation, we find that the log-density correlation function shows a much smaller inward shift in the position of the BAO peak at low redshift than is seen in the overdensity correlation function. We also show that if the initial, Lagrangian density of matter parcels could be estimated at their Eulerian positions, giving a displaced-initial-density field, its peak shift would be even smaller. In fact, a transformed field that accentuates underdensities, such as the reciprocal of the density, pushes the peak the other way, outward. In our model, these shifts in the peak position can be attributed to shift terms, involving the derivative of the linear correlation function, that entirely vanish in this displaced-initial-density field.

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles.We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

We examine the Lagrangian divergence of the displacement field, arguably a more natural object than the density in a Lagrangian description of cosmological large-scale structure. This quantity, which we denote \psi, quantifies the stretching and distortion of the initially homogeneous lattice of dark-matter particles in the universe. \psi\ encodes similar information as the density, but the correspondence has subtleties. It corresponds better to the log-density A than the overdensity \delta. A Gaussian distribution in \psi\ produces a distribution in A with slight skewness; in \delta, we find that in many cases the skewness is further increased by 3. A local spherical-collapse-based (SC) fit found by Bernardeau gives a formula for \psi's particle-by-particle behavior that works quite well, better than applying Lagrangian perturbation theory (LPT) at first or second (2LPT) order. In 2LPT, there is a roughly parabolic relation between initial and final \psi\ that can give overdensities in deep voids, so low-redshift, high-resolution 2LPT realizations should be used with caution. The SC fit excels at predicting \psi\ until streams cross; then, for particles forming haloes, \psi\ plummets as in a waterfall to -3. This gives a new method for producing N-particle realizations. Compared to LPT realizations, such SC realizations give reduced stream-crossing, and better visual and 1-point-PDF correspondence to the results of full gravity. LPT, on the other hand, predicts large-scale flows and the large-scale power-spectrum amplitude better, unless an empirical correction is added to the SC formula.

In a cold-dark-matter universe, cosmological structure formation proceeds in rough analogy to origami folding. Dark matter occupies a three-dimensional 'sheet' of free- fall observers, non-intersecting in six-dimensional velocity-position phase space. At early times, the sheet was flat like an origami sheet, i.e. velocities were essentially zero, but as time passes, the sheet folds up to form cosmic structure. The present paper further illustrates this analogy, and clarifies a Lagrangian definition of caustics and streams: caustics are two-dimensional surfaces in this initial sheet along which it folds, tessellating Lagrangian space into a set of three-dimensional regions, i.e. streams. The main scientific result of the paper is that streams may be colored by only two colors, with no two neighbouring streams (i.e. streams on either side of a caustic surface) colored the same. The two colors correspond to positive and negative parities of local Lagrangian volumes. This is a severe restriction on the connectivity and therefore arrangement of streams in Lagrangian space, since arbitrarily many colors can be necessary to color a general arrangement of three-dimensional regions. This stream two-colorability has consequences from graph theory, which we explain. Then, using N-body simulations, we test how these caustics correspond in Lagrangian space to the boundaries of haloes, filaments and walls. We also test how well outer caustics correspond to a Zel'dovich-approximation prediction.

Tessellations are valuable both conceptually and for analysis in the study of the large-scale structure of the universe. They provide a conceptual model for the 'cosmic web,' and are of great use to analyze cosmological data. Here we describe tessellations in another set of coordinates, of the initially flat sheet of dark matter that gravity folds up in rough analogy to origami. The folds that develop are called caustics, and they tessellate space into stream regions. Tessellations of the dark-matter sheet are also useful in simulation analysis, for instance for density measurement, and to identify structures where streams overlap.

We present the ORIGAMI method of identifying structures, particularly halos, in cosmological N-body simulations. Structure formation can be thought of as the folding of an initially flat three-dimensional manifold in six-dimensional phase space. ORIGAMI finds the outer folds that delineate these structures. Halo particles are identified as those that have undergone shell-crossing along 3 orthogonal axes, providing a dynamical definition of halo regions that is independent of density. ORIGAMI also identifies other morphological structures: particles that have undergone shell-crossing along 2, 1, or 0 orthogonal axes correspond to filaments, walls, and voids respectively. We compare this method to a standard Friends-of-Friends halo-finding algorithm and find that ORIGAMI halos are somewhat larger, more diffuse, and less spherical, though the global properties of ORIGAMI halos are in good agreement with other modern halo-finding algorithms.

We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub-)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non-homogeneous background density (e.g. a host halo). This includes real-space, phase-space, velocity-space and time- space based finders, as well as finders employing a Voronoi tessellation, friends-of-friends techniques, or refined meshes as the starting point for locating substructure.A common post-processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well on the presence and location of substructure and even for properties relating to the inner part part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that basic properties (mass, maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub-)subhalo hierarchy become progressively less important.

Motivated by recent results on lognormal statistics showing that the moment hierarchy of a lognormal variable completely fails at capturing its information content in the large variance regime, we discuss in this work the inadequacy of the hierarchy of correlation functions to describe a correlated lognormal field, which provides a roughly accurate description of the non-linear cosmological matter density field. We present families of fields having the same hierarchy of correlation functions than the lognormal field at all orders. This explicitly demonstrates the little studied though known fact that the correlation function hierarchy never provides a complete description of a lognormal field, and that it fails to capture information in the non-linear regime, where other simple observables are left totally unconstrained. We discuss why perturbative, Edgeworth-like approaches to statistics in the non-linear regime, common in cosmology, can never reproduce or predict that effect, and why it is however generic for tailed fields, hinting at a breakdown of the perturbation theory based on the field fluctuations. We make a rough but successful quantitative connection to N-body simulations results, that showed that the spectrum of the log-density field carries more information than the spectrum of the field entering the non-linear regime.

We investigate the use of a logarithmic density variable in estimating the Lagrangian displacement field, motivated by the success of a logarithmic transformation in restoring information to the matter power spectrum. The logarithmic relation is an extension of the linear relation, motivated by the continuity equation, in which the density field is assumed to be proportional to the divergence of the displacement field; we compare the linear and logarithmic relations by measuring both of these fields directly in a cosmological N-body simulation. The relative success of the logarithmic and linear relations depends on the scale at which the density field is smoothed. Thus we explore several ways of measuring the density field, including Cloud-In-Cell smoothing, adaptive smoothing, and the (scale-independent) Delaunay tessellation, and we use both a Fourier space and a geometrical tessellation approach to measuring the divergence. We find that the relation between the divergence of the displacement field and the density is significantly tighter with a logarithmic density variable, especially at low redshifts and for very small (~2 Mpc/h) smoothing scales. We find that the grid-based methods are more reliable than the tessellation-based method of calculating both the density and the divergence fields, though in both cases the logarithmic relation works better in the appropriate regime, which corresponds to nonlinear scales for the grid-based methods and low densities for the tessellation-based method.

It was recently shown that applying a Gaussianizing transform, such as a logarithm, to the nonlinear matter density field extends the range of useful applicability of the power spectrum by a factor of a few smaller. Such a transform dramatically reduces nonlinearities in both the covariance and the shape of the power spectrum. Here, analyzing Coyote Universe real-space dark matter density fields, we investigate the consequences of these transforms for cosmological parameter estimation. The power spectrum of the log-density provides the tightest cosmological parameter error bars (marginalized or not), giving a factor of 2-3 improvement over the conventional power spectrum in all five parameters tested. For the tilt, n_s, the improvement reaches a factor of 5. Similar constraints are achieved if the log-density power spectrum and conventional power spectrum are analyzed together. Rank-order Gaussianization seems just as useful as a log transform to constrain n_s, but not other parameters. Dividing the overdensity by its dispersion in few-Mpc cells, while it diagonalizes the covariance matrix, does not seem to help with parameter constraints. We also provide a code that emulates these power spectra over a range of concordance cosmological models.

The power spectrum is widely used in astronomy, to analyze temporal or spatial structure. In cosmology, it is used to quantify large-scale structure (LSS) and the cosmic microwave background (CMB). This is because the power spectrum completely quantifies Gaussian random fields, which the CMB and LSS fields seem to be at early epochs. However, at late epochs and small scales, cosmological density fields become highly non-Gaussian. The power spectrum loses power to describe LSS and CMB fields on small scales, most obviously through high covariance in the power spectrum as a function of scale. Practically, this significantly degrades constraints that observations can place on cosmological parameters. However, if a nonlinear transformation that produces a (more) Gaussian 1-point distribution is applied to a field, the extra covariance in the field's power spectrum can be dramatically reduced. In the case of the roughly lognormal low-redshift matter density field, a log transform accomplishes this. Applying a log transform to the density field before measuring the power spectrum also tightens cosmological parameter constraints by a factor of several.

We find a simple, accurate model for the covariance matrix of the real-space cosmological matter power spectrum on slightly nonlinear scales (k~0.1-0.8 h/Mpc at z=0), where off-diagonal matrix elements become substantial. The model includes a multiplicative, scale-independent modulation of the power spectrum. It has only one parameter, the variance (among realizations) of the variance of the nonlinear density field in cells, with little dependence on the cell size between 2-8 Mpc/h. Furthermore, we find that this extra covariance can be modeled out by instead measuring the power spectrum of (delta/sigma_cell), i.e. the ratio of the overdensity to its dispersion in cells a few Mpc in size. Dividing delta by sigma_cell essentially removes the non-Gaussian part of the covariance matrix, nearly diagonalizing it.

The matter density field exhibits a nearly lognormal probability density distribution (PDF) after entering into the nonlinear regime. Recently, it has been shown that the shape of the power spectrum of a logarithmically transformed density field is very close to the linear density power spectrum, motivating an analytic study of it. In this paper, we develop cosmological perturbation theory for the power spectrum of this field. Our formalism is developed in the context of renormalized perturbation theory, which helps to regulate the convergence behavior of the perturbation series, and of the Taylor- series expansion we use of the logarithmic mapping. This approach allows us to handle the critical issue of density smoothing in a straightforward way. We also compare our perturbative results with simulation measurements.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

We find that, even in the presence of discreteness noise, a Gaussianizing transform (producing a more-Gaussian one-point distribution) reduces nonlinearities in the power spectra of cosmological matter and galaxy density fields, in many cases drastically. Although Gaussianization does increase the effective shot noise, it also increases the power spectrum's fidelity to the linear power spectrum on scales where the shot noise is negligible. Gaussianizing also increases the Fisher information in the power spectrum in all cases and resolutions, although the gains are smaller in redshift space than in real space. We also find that the gain in cumulative Fisher information from Gaussianizing peaks at a particular grid resolution that depends on the sampling level.

We show that redshift-space distortions of galaxy correlations have a strong effect on correlation functions with distinct, localized features, like the signature of the baryon acoustic oscillations (BAO). Near the line of sight, the features become sharper as a result of redshift-space distortions. We demonstrate this effect by measuring the correlation function in Gaussian simulations and the Millennium Simulation. We also analyze the SDSS DR7 main-galaxy sample (MGS), splitting the sample into slices 2.5 degrees on the sky in various rotations. Measuring 2D correlation functions in each slice, we do see a sharp bump along the line of sight. Using Mexican-hat wavelets, we localize it to (110 +/- 10) Mpc/h. Averaging only along the line of sight, we estimate its significance at a particular wavelet scale and location at 2.2 sigma. In a flat angular weighting in the (pi,r_p) coordinate system, the noise level is suppressed, pushing the bump's significance to 4 sigma. We estimate that there is about a 0.2% chance of getting such a signal anywhere in the vicinity of the BAO scale from a power spectrum lacking a BAO feature. However, these estimates of the significances make some use of idealized Gaussian simulations, and thus are likely a bit optimistic.

The Cold Spot on the Cosmic Microwave Background could arise due to a supervoid at low redshift through the integrated Sachs-Wolfe effect. We imaged the region with MegaCam on the Canada-France-Hawai'i Telescope and present galaxy counts in photometric redshift bins. We rule out the existence of a 100Mpc radius spherical supervoid with underdensity delta=-0.3 at 0.5<z<0.9 at high significance. The data are consistent with an underdensity at low redshift, but the fluctuations are within the range of cosmic variance and the low density areas are not contiguous on the sky. Thus, we find no strong evidence for a supervoid. We cannot resolve voids smaller than 50Mpc radius; however, these can only make a minor contribution to the CMB temperature decrement.

We construct a map of the time derivative of the gravitational potential traced by SDSS Luminous Red Galaxies. The potential decays on large scales due to cosmic acceleration, leaving an imprint on cosmic microwave background (CMB) radiation through the integrated Sachs-Wolfe (ISW) effect. With a template fit, we directly measure this signature on the CMB at a 2-sigma confidence level. The measurement is consistent with the cross-correlation statistic, strengthening the claim that dark energy is indeed the cause of the correlation. This new approach potentially simplifies the cosmological interpretation. Our constructed linear ISW map shows no evidence for degree-scale cold and hot spots associated with supervoid and supercluster structures. This suggests that the linear ISW effect in a concordance Lambda-CDM cosmology is insufficient to explain the strong CMB imprints from these structures that we previously reported.

We find that nonlinearities in the dark-matter power spectrum are dramatically smaller if the density field first undergoes a logarithmic mapping. In the Millennium simulation, this procedure produces a power spectrum with a shape hardly departing from the linear power spectrum for k <~ 1 h/Mpc at all redshifts. Also, this procedure unveils pristine Fisher information on a range of scales reaching a factor of 2-3 smaller than in the standard power spectrum, yielding 10 times more cumulative signal-to-noise at z=0.

We measure hot and cold spots on the microwave background associated with supercluster and supervoid structures identified in the Sloan Digital Sky Survey Luminous Red Galaxy catalog. The structures give a compelling visual imprint, with a mean temperature deviation of 9.6 +/- 2.2 microK, i.e. above 4 sigma. We interpret this as a detection of the late-time Integrated Sachs-Wolfe (ISW) effect, in which cosmic acceleration from dark energy causes gravitational potentials to decay, heating or cooling photons passing through density crests or troughs. In a flat universe, the linear ISW effect is a direct signal of dark energy.

The observed apparent acceleration of the universe is usually attributed to negative pressure from a mysterious dark energy. This acceleration causes the gravitational potential to decay, heating or cooling photons travelling through crests or troughs of large-scale matter density fluctuations. This phenomenon, the late-time integrated Sachs-Wolfe (ISW) effect, has been detected, albeit at low significance, by cross-correlating various galaxy surveys with the Cosmic Microwave Background (CMB). Recently, the best evidence has come from the statistical combination of results from multiple correlated galaxy data sets. Here we show that vast structures identified in a galaxy survey project an image onto the CMB; stacking regions aligned with superclusters produces a hot spot, and supervoids, a cold spot. At over 4 sigma, this is the clearest evidence of the ISW effect to date. For the first time, our findings pin the effect to discrete structures. The ISW signal from supervoids and superclusters can be combined with other cosmological probes to constrain dark energy and cosmological parameters. In addition, our findings make it more plausible that the extreme Cold Spot and other anomalies in the CMB are caused by supervoids.

Despite a history that dates back at least a quarter of a century studies of voids in the large--scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void--finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recent increased interest in voids, both in very large galaxy surveys and in detailed simulations of cosmic structure formation, this situation is very unfortunate. We here present the first systematic comparison study of thirteen different void finders constructed using particles, haloes, and semi--analytical model galaxies extracted from a subvolume of the Millennium simulation. The study includes many groups that have studied voids over the past decade. We show their results and discuss their differences and agreements. As it turns out, the basic results of the various methods agree very well with each other in that they all locate a major void near the centre of our volume. Voids have very underdense centres, reaching below 10 percent of the mean cosmic density. In addition, those void finders that allow for void galaxies show that those galaxies follow similar trends. For example, the overdensity of void galaxies brighter than $m_B = -20 $ is found to be smaller than about -0.8 by all our void finding algorithms.

ZOBOV (ZOnes Bordering On Voidness) is an algorithm that finds density depressions in a set of points, without any free parameters, or assumptions about shape. It uses the Voronoi tessellation to estimate densities, which it uses to find both voids and subvoids. It also measures probabilities that each void or subvoid arises from Poisson fluctuations. This paper describes the ZOBOV algorithm, and the results from its application to the dark-matter particles in a region of the Millennium Simulation. Additionally, the paper points out an interesting high-density peak in the probability distribution of dark-matter particle densities.

We investigate large-amplitude baryon acoustic oscillations (BAO's) in off-diagonal entries of cosmological power-spectrum covariance matrices. These covariance-matrix BAO's describe the increased attenuation of power-spectrum BAO's caused by upward fluctuations in large-scale power. We derive an analytic approximation to covariance-matrix entries in the BAO regime, and check the analytical predictions using N-body simulations. These BAO's look much stronger than the BAO's in the power spectrum, but seem detectable only at about a one-sigma level in gigaparsec-scale galaxy surveys. In estimating cosmological parameters using matter or galaxy power spectra, including the covariance-matrix BAO's can have a several-percent effect on error-bar widths for some parameters directly related to the BAO's, such as the baryon fraction. Also, we find that including the numerous galaxies in small haloes in a survey can reduce error bars in these cosmological parameters more than the simple reduction in shot noise might suggest.

We investigate the cosmological Fisher information in the non-linear dark-matter power spectrum in the context of the halo model. We find that there is a plateau in information content on translinear scales which is generic to all cosmological parameters we tried. There is a rise in information on smaller scales, but we find that it is quite degenerate among different cosmological parameters (except, perhaps, the tilt). This suggests that it could be difficult to constrain cosmological parameters using the non-linear regime of the dark-matter power spectrum. We suggest ways to get around this problem, such as removing the largest haloes from consideration in survey analysis.

Using the halo model, we investigate the cosmological Fisher information in the non-linear dark-matter power spectrum about the initial amplitude of linear power. We find that there is little information on `translinear' scales (where the one- and two-halo terms are both significant) beyond what is on linear scales, but that additional information is present on small scales, where the one-halo term dominates. This behavior agrees with the surprising results that Rimes & Hamilton (2005, 2006) found using N-body simulations. We argue that the translinear plateau in cumulative information arises largely from fluctuations in the numbers of large haloes in a finite volume. This implies that more information could be extracted on non-linear scales if the masses of the largest haloes in a survey are known.

We present a new, galaxy-halo model of large-scale structure, in which the galaxies entering a given sample are the fundamental objects. Haloes attach to galaxies, in contrast to the standard halo model, in which galaxies attach to haloes. The galaxy-halo model pertains mainly to the relationships between the power spectra of galaxies and mass, and their cross-power spectrum. With surprisingly little input, an intuition-aiding approximation to the galaxy-matter cross-correlation coefficient R(k) emerges, in terms of the halo mass dispersion. This approximation seems valid to mildly non-linear scales (k < ~3 h/Mpc), allowing measurement of the bias and the matter power spectrum from measurements of the galaxy and galaxy-matter power spectra (or correlation functions). This is especially relevant given the recent advances in precision in measurements of the galaxy-matter correlation function from weak gravitational lensing. The galaxy-halo model also addresses the issue of interpreting the galaxy-matter correlation function as an average halo density profile, and provides a simple description of galaxy bias as a function of scale.

We have developed an algorithm to find haloes in an N-body dark matter simulation, called VOBOZ (VOronoi BOund Zones), which has as little dependence on free parameters as we can manage. By using the Voronoi diagram, we achieve nonparametric, `natural' measurements of each particle's density and set of neighbors. We then eliminate much of the ambiguity in merging sets of particles together by identifying every possible density peak, and measuring the probability that each does not arise from Poisson noise. This probability could be interpreted as a measure of a halo's sub-haloness, and might be of use to researchers of the halo model. The first parameter in VOBOZ controls the subtlety of particle unbinding, and may be eliminated if one is cavalier with processor time; even if one is not, the results saturate to the parameter-free answer when the parameter is sufficiently small. The only parameter which remains, an outer density cut-off, does not influence whether or not haloes are identified, nor does it have any effect on subhaloes; it only affects the masses returned for supercluster haloes.

By comparing the PSCz galaxy power spectrum with the results of nested pure dark matter N-body simulations, we try to understand how infrared-selected galaxies populate dark-matter haloes, paying special attention to the method of halo identification in the simulations. We thus test the hypothesis that baryonic physics negligibly affects the distribution of galaxies down to the smallest scales yet observed. We are successful in reproducing the PSCz power spectrum on scales < ~40 h/Mpc, near our resolution limit, by imposing a central density cut-off on simulated haloes, which gives a rough minimum mass and circular velocity of haloes in which PSCz galaxies formed.

By comparing the PSCz galaxy power spectrum with haloes from nested and phased N-body simulations, we try to understand how IRAS infrared-selected galaxies populate dark-matter haloes. We pay special attention to the way we identify haloes in the simulations.

We investigate the relationship between the dark matter and baryons in the linear regime. This relation is quantified by the so-called ``filtering scale''. We show that a simple gaussian ansatz which uses the filtering scale provides a good approximation to the exact solution.