We introduce a new method for performing robust Bayesian estimation of the three-dimensional spatial power spectrum at the Epoch of Reionization (EoR), from interferometric observations. The versatility of this technique allows us to present two approaches. First, when the observations span only a small number of independent spatial frequencies ($k$-modes) we sample directly from the spherical power spectrum coefficients that describe the EoR signal realisation. Second, when the number of $k$-modes to be included in the model becomes large, we sample from the joint probability density of the spherical power spectrum and the signal coefficients, using Hamiltonian Monte Carlo methods to explore this high dimensional ($\sim$ 20000) space efficiently. This approach has been successfully applied to simulated observations that include astrophysically realistic foregrounds in a companion publication (Sims et al. 2016). Here we focus on explaining the methodology in detail, and use simple foreground models to both demonstrate its efficacy, and highlight salient features. In particular, we show that including an arbitrary flat spectrum continuum foreground that is $10^8$ times greater in power than the EoR signal has no detectable impact on our parameter estimates of the EoR power spectrum recovered from the data.

A key physical quantity during reionization is the size of HII regions. Previous studies found a characteristic bubble size which increases rapidly during reionization, with apparent agreement between simulations and analytic excursion set theory. Using four different methods, we critically examine this claim. In particular, we introduce the use of the watershed algorithm -- widely used for void finding in galaxy surveys -- which we show to be an unbiased method with the lowest dispersion and best performance on Monte-Carlo realizations of a known bubble size PDF. We find that a friends-of-friends algorithm declares most of the ionized volume to be occupied by a network of volume-filling regions connected by narrow tunnels. For methods tuned to detect the volume-filling regions, previous apparent agreement between simulations and theory is spurious, and due to a failure to correctly account for the window function of measurement schemes. The discrepancy is already obvious from visual inspection. Instead, HII regions in simulations are significantly larger (by factors of 10-1000 in volume) than analytic predictions. The size PDF is narrower, and evolves more slowly with time, than predicted. It becomes more sharply peaked as reionization progresses. These effects are likely caused by bubble mergers, which are inadequately modeled by analytic theory. Our results have important consequences for high-redshift 21cm observations, the mean free path of ionizing photons, and the visibility of Ly-alpha emitters, and point to a fundamental failure in our understanding of the characteristic scales of the reionization process.

The Universe is mostly composed of large and relatively empty domains known as cosmic voids, whereas its matter content is predominantly distributed along their boundaries. The remaining material inside them, either dark or luminous matter, is attracted to these boundaries and causes voids to expand faster and to grow emptier over time. Using the distribution of galaxies centered on voids identified in the Sloan Digital Sky Survey and adopting minimal assumptions on the statistical motion of these galaxies, we constrain the average matter content $\Omega_\mathrm{m}=0.281\pm0.031$ in the Universe today, as well as the linear growth rate of structure $f/b=0.417\pm0.089$ at median redshift $\bar{z}=0.57$, where $b$ is the galaxy bias ($68\%$ C.L.). These values originate from a percent-level measurement of the anisotropic distortion in the void-galaxy cross-correlation function, $\varepsilon = 1.003\pm0.012$, and are robust to consistency tests with bootstraps of the data and simulated mock catalogs within an additional systematic uncertainty of half that size. They surpass (and are complementary to) existing constraints by unlocking cosmological information on smaller scales through an accurate model of nonlinear clustering and dynamics in void environments. As such, our analysis furnishes a powerful probe of deviations from Einstein's general relativity in the low-density regime which has largely remained untested so far. We find no evidence for such deviations in the data at hand.

The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]

We investigate the potential of using cosmic voids as a probe to constrain cosmological parameters through the gravitational lensing effect of the cosmic microwave background (CMB) and make predictions for the next generation surveys. By assuming the detection of a series of $\approx 5 - 10$ voids along a line of sight within a square-degree patch of the sky, we found that they can be used to break the degeneracy direction of some of the cosmological parameter constraints (for example $\omega_b$ and $\Omega_\Lambda$) in comparison with the constraints from random CMB skies with the same size area for a survey with extensive integration time. This analysis is based on our current knowledge of the average void profile and analytical estimates of the void number function. We also provide combined cosmological parameter constraints between a sky patch where series of voids are detected and a patch without voids (a randomly selected patch). The full potential of this technique relies on an accurate determination of the void profile to $\approx 10$% level. For a small-area CMB observation with extensive integration time and a high signal-to-noise ratio, CMB lensing with such series of voids will provide a complementary route to cosmological parameter constraints to the CMB observations. Example of parameter constraints with a series of five voids on a $1.0^{\circ} \times 1.0^{\circ}$ patch of the sky are $100\omega_b = 2.20 \pm 0.27$, $\omega_c = 0.120 \pm 0.022$, $\Omega_\Lambda = 0.682 \pm 0.078$, $\Delta_{\mathcal{R}}^2 = \left(2.22 \pm 7.79\right) \times 10^{-9}$, $n_s = 0.962 \pm 0.097$ and $\tau = 0.925 \pm 1.747$ at 68% C.L.

We present in this paper a new Bayesian semi-blind approach for foreground removal in observations of the 21-cm signal with interferometers. The technique, which we call HIEMICA (HI Expectation-Maximization Independent Component Analysis), is an extension of the Independent Component Analysis (ICA) technique developed for two-dimensional (2D) CMB maps to three-dimensional (3D) 21-cm cosmological signals measured by interferometers. This technique provides a fully Bayesian inference of power spectra and maps and separates the foregrounds from signal based on the diversity of their power spectra. Only relying on the statistical independence of the components, this approach can jointly estimate the 3D power spectrum of the 21-cm signal and, the 2D angular power spectrum and the frequency dependence of each foreground component, without any prior assumptions about foregrounds. This approach has been tested extensively by applying it to mock data from interferometric 21-cm intensity mapping observations under idealized assumptions of instrumental effects. We also discuss the impact when the noise properties are not known completely. As a first step toward solving the 21 cm power spectrum analysis problem we compare the semi-blind HIEMICA technique with the commonly used Principal Component Analysis (PCA). Under the same idealized circumstances the proposed technique provides significantly improved recovery of the power spectrum. This technique can be applied straightforwardly to all 21-cm interferometric observations, including epoch of reionization measurements, and can be extended to single-dish observations as well.

Cosmic voids in the large-scale structure of the Universe affect the peculiar motions of objects in their vicinity. Although these motions are difficult to observe directly, the clustering pattern of their surrounding tracers in redshift space is influenced in a unique way. This allows to investigate the interplay between densities and velocities around voids, which is solely dictated by the laws of gravity. With the help of $N$-body simulations and derived mock-galaxy catalogs we calculate the average density fluctuations around voids identified with a watershed algorithm in redshift space and compare the results with the expectation from general relativity and the $\Lambda$CDM model. We find linear theory to work remarkably well in describing the dynamics of voids. Adopting a Bayesian inference framework, we explore the full posterior of our model parameters and forecast the achievable accuracy on measurements of the growth rate of structure and the geometric distortion through the Alcock-Paczynski effect. Systematic errors in the latter are reduced from $\sim15\%$ to $\sim5\%$ when peculiar velocities are taken into account. The relative parameter uncertainties in galaxy surveys with number densities comparable to the SDSS MAIN (CMASS) sample probing a volume of $1h^{-3}{\rm Gpc}^3$ yield $\sigma_{f/b}\left/(f/b)\right.\sim2\%$ ($20\%$) and $\sigma_{D_AH}/D_AH\sim0.2\%$ ($2\%$), respectively. At this level of precision the linear-theory model becomes systematics dominated, with parameter biases that fall beyond these values. Nevertheless, the presented method is highly model independent; its viability lies in the underlying assumption of statistical isotropy of the Universe.

We present a non-parametric, model-independent method to reconstruct the spherical density profiles of void stacks in real space, without redshift-space distortions. Our method uses the expected spherical symmetry of stacked voids to build the shape of the spherical density profile of cosmic voids in real space without any assumption about the cosmological model. We test the reconstruction algorithm with both a toy model and a full dark matter simulation. We present the result for the simulation: the reconstruction of the spherical density profile for a simulated stacked void in real space. We also present a first application of the algorithm to reconstruct real cosmic void stacks density profiles in real space from the Sloan Digital Sky Survey (Sutter et al. 2012b). We discuss capabilities of the algorithm and possible future improvements. Reconstructed density profiles from real voids open the way to the study of the spherically averaged dynamical structure of voids.

We show that the number of observed voids in galaxy redshift surveys is a sensitive function of the equation of state of dark energy. Using the Fisher matrix formalism we find the error ellipses in the $w_0-w_a$ plane when the equation of state of dark energy is assumed to be of the form $w_{CPL}(z)=w_0 +w_a z/(1+z)$. We forecast the number of voids to be observed with the ESA Euclid satellite and the NASA WFIRST mission, taking into account updated details of the surveys to reach accurate estimates of their power. The theoretical model for the forecast of the number of voids is based on matches between abundances in simulations and the analytical prediction. To take into account the uncertainties within the model, we marginalize over its free parameters when calculating the Fisher matrices. The addition of the void abundance constraints to the data from Planck, HST and supernova survey data noticeably tighten the $w_0-w_a$ parameter space. We thus quantify the improvement in the constraints due to the use of voids and demonstrate that the void abundance is a sensitive new probe for the dark energy equation of state.

Cosmic voids are a promising environment to characterize neutrino-induced effects on the large-scale distribution of matter in the universe. We perform a comprehensive numerical study of the statistical properties of voids, identified both in the matter and galaxy distributions, in massive and massless neutrino cosmologies. The matter density field is obtained by running several independent $N$-body simulations with cold dark matter and neutrino particles, while the galaxy catalogs are modeled by populating the dark matter halos in simulations via a halo occupation distribution (HOD) model to reproduce the clustering properties observed by the Sloan Digital Sky Survey (SDSS) II Data Release 7. We focus on the impact of massive neutrinos on the following void statistical properties: number density, ellipticities, two-point statistics, density and velocity profiles. Considering the matter density field, we find that voids in massive neutrino cosmologies are less evolved than those in the corresponding massless neutrinos case: there is a larger number of small voids and a smaller number of large ones, their profiles are less evacuated, and they present a lower wall at the edge. Moreover, the degeneracy between $\sigma_8$ and $\Omega_{\nu}$ is broken when looking at void properties. In terms of the galaxy density field, we find that differences among cosmologies are difficult to detect because of the small number of galaxy voids in the simulations. Differences are instead present when looking at the matter density and velocity profiles around these voids.

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

How do peculiar velocities affect observed voids? To answer this question we use the VIDE toolkit to identify voids in mock galaxy populations embedded within an N-body simulation both with and without peculiar velocities included. We compare the resulting void populations to assess the impact on void properties. We find that void abundances and spherically-averaged radial density profiles are mildly affected by peculiar velocities. However, peculiar velocities can distort by up to 10% the shapes for a particular subset of voids depending on the void size and density contrast, which can lead to increased variance in Alcock-Paczy\'nski test. We offer guidelines for performing optimal cuts on the void catalogue to reduce this variance by removing the most severely affected voids while preserving the unaffected ones. In addition, since this shape distortion is largely limited to the line of sight, we show that the void radii are only affected at the $\sim$ 10% level and the macrocenter positions at the $\sim$ 20% (even before performing cuts), meaning that cosmological probes based on the Integrated Sachs-Wolfe and gravitational lensing are not severely impacted by peculiar velocities.

What do we know about voids in the dark matter distribution given the Sloan Digital Sky Survey (SDSS) and assuming the $\Lambda\mathrm{CDM}$ model? Recent application of the Bayesian inference algorithm BORG to the SDSS Data Release 7 main galaxy sample has generated detailed Eulerian and Lagrangian representations of the large-scale structure as well as the possibility to accurately quantify corresponding uncertainties. Building upon these results, we present constrained catalogs of voids in the Sloan volume, aiming at a physical representation of dark matter underdensities and at the alleviation of the problems due to sparsity and biasing on galaxy void catalogs. To do so, we generate data-constrained reconstructions of the presently observed large-scale structure using a fully non-linear gravitational model. We then find and analyze void candidates using the VIDE toolkit. Our methodology therefore predicts the properties of voids based on fusing prior information from simulations and data constraints. For usual void statistics (number function, ellipticity distribution and radial density profile), all the results obtained are in agreement with dark matter simulations. Our dark matter void candidates probe a deeper void hierarchy than voids directly based on the observed galaxies alone. The use of our catalogs therefore opens the way to high-precision void cosmology at the level of the dark matter field. We will make the void catalogs used in this work available at http://www.cosmicvoids.net.

We use properties of void populations identified in $N$-body simulations to forecast the ability of upcoming galaxy surveys to differentiate models of f(R) gravity from \lcdm~cosmology. We analyze multiple simulation realizations, which were designed to mimic the expected number densities, volumes, and redshifts of the upcoming Euclid satellite and a lower-redshift ground-based counterpart survey, using the public {\tt VIDE} toolkit. We examine void abundances, ellipicities, radial density profiles, and radial velocity profiles at redshifts 1.0 and 0.43. We find that stronger f(R) coupling strengths eliminates small voids and produces voids up to $\sim 20\%$ larger in radius, leading to a significant tilt in the void number function. Additionally, under the influence of modified gravity, voids at all scales tend to be measurably emptier with correspondingly higher compensation walls. The velocity profiles reflect this, showing increased outflows inside voids and increased inflows outside voids. Using the void number function as an example, we forecast that future surveys can constrain the modified gravity coupling strength to $\sim 3 \times 10^{-5}$ using voids.

We present a general method for accelerating by more than an order of magnitude the convolution of pixelated functions on the sphere with a radially-symmetric kernel. Our method splits the kernel into a compact real-space component and a compact spherical harmonic space component. These components can then be convolved in parallel using an inexpensive commodity GPU and a CPU. We provide models for the computational cost of both real-space and Fourier space convolutions and an estimate for the approximation error. Using these models we can determine the optimum split that minimizes the wall clock time for the convolution while satisfying the desired error bounds. We apply this technique to the problem of simulating a cosmic microwave background (CMB) anisotropy sky map at the resolution typical of the high resolution maps produced by the Planck mission. For the main Planck CMB science channels we achieve a speedup of over a factor of ten, assuming an acceptable fractional rms error of order 1.e-5 in the power spectrum of the output map.

Infrared and radio observations of the Epoch of Reionization promise to revolutionize our understanding of the cosmic dawn, and major efforts with the JWST, MWA and HERA are underway. While measurements of the ionizing sources with infrared telescopes and the effect of these sources on the intergalactic medium with radio telescopes \emph{should} be complementary, to date the wildly disparate angular resolutions and survey speeds have made connecting proposed observations difficult. In this paper we develop a method to bridge the gap between radio and infrared studies. While the radio images may not have the sensitivity and resolution to identify individual bubbles with high fidelity, by leveraging knowledge of the measured power spectrum we are able to separate regions that are likely ionized from largely neutral, providing context for the JWST observations of galaxy counts and properties in each. By providing the ionization context for infrared galaxy observations, this method can significantly enhance the science returns of JWST and other infrared observations.

We propose a novel technique to probe the expansion history of the Universe based on the clustering statistics of cosmic voids. In particular, we compute their two-point statistics in redshift space on the basis of realistic mock galaxy catalogs and apply the Alcock-Paczynski test. In contrast to galaxies, we find void auto-correlations to be marginally affected by peculiar motions, providing a model-independent measure of cosmological parameters without systematics from redshift-space distortions. Because only galaxy-galaxy and void-galaxy correlations have been considered in these types of studies before, the presented method improves both statistical and systematic uncertainties on the product of angular diameter distance and Hubble rate, furnishing the potentially cleanest probe of cosmic geometry available to date.

We estimate the galaxy density along lines of sight to hard extragalactic gamma-ray sources by correlating source positions on the sky with a void catalog based on the Sloan Digital Sky Survey (SDSS). Extragalactic gamma-ray sources that are detected at very high energy (VHE; E>100 GeV) or have been highlighted as VHE-emitting candidates in the Fermi Large Area Telescope hard source catalog (together referred to as "VHE-like" sources) are distributed along underdense lines of sight at the 2.4 sigma level. There is also a less suggestive correlation for the Fermi hard source population (1.7 sigma). A correlation between 10-500 GeV flux and underdense fraction along the line of sight for VHE-like and Fermi hard sources is found at 2.4 sigma and 2.6 sigma, respectively. The preference for underdense sight lines is not displayed by gamma-ray emitting galaxies within the second Fermi catalog, containing sources detected above 100 MeV, or the SDSS DR7 quasar catalog. We investigate whether this marginal correlation might be a result of lower extragalactic background light (EBL) photon density within the underdense regions and find that, even in the most extreme case of a entirely underdense sight line, the EBL photon density is only 2% less than the nominal EBL density. Translating this into gamma-ray attenuation along the line of sight for a highly attenuated source with opacity tau(E,z) ~5, we estimate that the attentuation of gamma-rays decreases no more than 10%. This decrease, although non-neglible, is unable to account for the apparent hard source correlation with underdense lines of sight.

We present VIDE, the Void IDentification and Examination toolkit, an open-source Python/C++ code for finding cosmic voids in galaxy redshift surveys and N-body simulations, characterizing their properties, and providing a platform for more detailed analysis. At its core, VIDE uses a substantially enhanced version of ZOBOV (Neyinck 2008) to calculate a Voronoi tessellation for estimating the density field and a performing a watershed transform to construct voids. Additionally, VIDE provides significant functionality for both pre- and post-processing: for example, vide can work with volume- or magnitude-limited galaxy samples with arbitrary survey geometries, or dark matter particles or halo catalogs in a variety of common formats. It can also randomly subsample inputs and includes a Halo Occupation Distribution model for constructing mock galaxy populations. VIDE uses the watershed levels to place voids in a hierarchical tree, outputs a summary of void properties in plain ASCII, and provides a Python API to perform many analysis tasks, such as loading and manipulating void catalogs and particle members, filtering, plotting, computing clustering statistics, stacking, comparing catalogs, and fitting density profiles. While centered around ZOBOV, the toolkit is designed to be as modular as possible and accommodate other void finders. VIDE has been in development for several years and has already been used to produce a wealth of results, which we summarize in this work to highlight the capabilities of the toolkit. VIDE is publicly available at http://bitbucket.org/cosmicvoids/vide public and http://www.cosmicvoids.net.

We use properties of void populations identified in N-body simulations to forecast the ability of upcoming galaxy surveys to differentiate models of f(R) gravity from \Lambda CDM cosmology. We analyze simulations designed to mimic the densities, volumes, and clustering statistics of upcoming surveys, using the public VIDE toolkit. We examine void abundances as a basic probe at redshifts 1.0 and 0.4. We find that stronger f(R) coupling strengths produce voids up to ~20% larger in radius, leading to a significant shift in the void number function. As an initial estimate of the constraining power of voids, we use this change in the number function to forecast a constraint on the coupling strength of $\Delta f_{R_{0}} = 10^{-5}$.

Taking N-body simulations with volumes and particle densities tuned to match the SDSS DR7 spectroscopic main sample, we assess the ability of current void catalogs (e.g., Sutter et al. 2012b) to distinguish a model of coupled dark matter-dark energy from {\Lambda}CDM cosmology using properties of cosmic voids. Identifying voids with the VIDE toolkit, we find no statistically significant differences in the ellipticities, but find that coupling produces a population of significantly larger voids, possibly explaining the recent result of Tavasoli et al. (2013). In addition, we use the universal density profile of Hamaus et al. (2014) to quantify the relationship between coupling and density profile shape, finding that the coupling produces broader, shallower, undercompensated profiles for large voids by thinning the walls between adjacent medium-scale voids. We find that these differences are potentially measurable with existing void catalogs once effects from survey geometries and peculiar velocities are taken into account.

Understanding the internal structure and spatial distribution of cosmic voids is crucial when considering them as probes of cosmology. We present recent advances in modeling void density- and velocity-profiles in real space, as well as void two-point statistics in redshift space, by examining voids identified via the watershed transform in state-of-the-art $\Lambda$CDM n-body simulations and mock galaxy catalogs. The simple and universal characteristics that emerge from these statistics indicate the self-similarity of large-scale structure and suggest cosmic voids to be among the most pristine objects to consider for future studies on the nature of dark energy, dark matter and modified gravity.

We discuss various applications of VIDE, the Void IDentification and Examination toolkit, an open-source Python/C++ code for finding cosmic voids in galaxy redshift surveys and N-body simulations. Based on a substantially enhanced version of ZOBOV, VIDE not only finds voids, but also summarizes their properties, extracts statistical information, and provides a Python-based platform for more detailed analysis, such as manipulating void catalogs and particle members, filtering, plotting, computing clustering statistics, stacking, comparing catalogs, and fitting density profiles. VIDE also provides significant additional functionality for pre-processing inputs: for example, VIDE can work with volume- or magnitude-limited galaxy samples with arbitrary survey geometries, or dark matter particles or halo catalogs in a variety of common formats. It can also randomly subsample inputs and includes a Halo Occupation Distribution model for constructing mock galaxy populations. VIDE has been used for a wide variety of applications, from discovering a universal density profile to estimating primordial magnetic fields, and is publicly available at http://bitbucket.org/cosmicvoids/vide_public and http://www.cosmicvoids.net.

We investigate the formation, growth, merger history, movement, and destruction of cosmic voids detected via the watershed transform code VIDE in a cosmological N-body dark matter {\Lambda}CDM simulation. By adapting a method used to construct halo merger trees, we are able to trace individual voids back to their initial appearance and record the merging and evolution of their progenitors at high redshift. For the scales of void sizes captured in our simulation, we find that the void formation rate peaks at scale factor 0.3, which coincides with a growth in the void hierarchy and the emergence of dark energy. Voids of all sizes appear at all scale factors, though the median initial void size decreases with time. When voids become detectable they have nearly their present-day volumes. Almost all voids have relatively stable growth rates and suffer only infrequent minor mergers. Dissolution of a void via merging is very rare. Instead, most voids maintain their distinct identity as annexed subvoids of a larger parent. The smallest voids are collapsing at the present epoch, but void destruction ceases after scale factor 0.3. In addition, voids centers tend to move very little, less than 0.01 of their effective radii per ln a, over their lifetimes. Overall, most voids exhibit little radical dynamical evolution; their quiet lives make them pristine probes of cosmological initial conditions and the imprint of dark energy.

We perform an Alcock-Paczynski test using stacked cosmic voids identified in the SDSS Data Release 7 main sample and Data Release 10 LOWZ and CMASS samples. We find ~1,500 voids out to redshift $0.6$ using a heavily modified and extended version of the watershed algorithm ZOBOV, which we call VIDE (Void IDentification and Examination). To assess the impact of peculiar velocities we use the mock void catalogs presented in Sutter et al. (2013). We find a constant uniform flattening of 14% along the line of sight when peculiar velocities are included. This flattening appears universal for all void sizes at all redshifts and for all tracer densities. We also use these mocks to identify an optimal stacking strategy. After correcting for systematic effects we find that our Alcock-Paczynski measurement leads to a preference of our best-fit value of $\Omega_{\rm M}\sim 0.15$ over $\Omega_{\rm M} = 1.0$ by a likelihood ratio of 10. Likewise, we find a factor of $4.5$ preference of the likelihood ratio for a $\Lambda$CDM $\Omega_{\rm M} = 0.3$ model and a null measurement. Taken together, we find substantial evidence for the Alcock-Paczynski signal in our sample of cosmic voids. Our assessment using realistic mocks suggests that measurements with future SDSS releases and other surveys will provide tighter cosmological parameter constraints. The void-finding algorithm and catalogs used in this work will be made publicly available at http://www.cosmicvoids.net.

We present results from an end-to-end simulation pipeline interferometric observations of cosmic microwave background polarization. We use both maximum-likelihood and Gibbs sampling techniques to estimate the power spectrum. In addition, we use Gibbs sampling for image reconstruction from interferometric visibilities. The results indicate the level to which various systematic errors (e.g., pointing errors, gain errors, beam shape errors, cross- polarization) must be controlled in order to successfully detect and characterize primordial B modes as well as other scientific goals. In addition, we show that Gibbs sampling is an effective method of image reconstruction for interferometric data in other astrophysical contexts.

The Dark Sky Simulations are an ongoing series of cosmological N-body simulations designed to provide a quantitative and accessible model of the evolution of the large-scale Universe. Such models are essential for many aspects of the study of dark matter and dark energy, since we lack a sufficiently accurate analytic model of non-linear gravitational clustering. In July 2014, we made available to the general community our early data release, consisting of over 55 Terabytes of simulation data products, including our largest simulation to date, which used $1.07 \times 10^{12}~(10240^3)$ particles in a volume $8h^{-1}\mathrm{Gpc}$ across. Our simulations were performed with 2HOT, a purely tree-based adaptive N-body method, running on 200,000 processors of the Titan supercomputer, with data analysis enabled by yt. We provide an overview of the derived halo catalogs, mass function, power spectra and light cone data. We show self-consistency in the mass function and mass power spectrum at the 1% level over a range of more than 1000 in particle mass. We also present a novel method to distribute and access very large datasets, based on an abstraction of the World Wide Web (WWW) as a file system, remote memory-mapped file access semantics, and a space-filling curve index. This method has been implemented for our data release, and provides a means to not only query stored results such as halo catalogs, but also to design and deploy new analysis techniques on large distributed datasets.

We present a simple empirical function for the average density profile of cosmic voids, identified via the watershed technique in $\Lambda$CDM N-body simulations. This function is universal across void size and redshift, accurately describing a large radial range of scales around void centers with only two free parameters. In analogy to halo density profiles, these parameters describe the scale radius and the central density of voids. While we initially start with a more general four-parameter model, we find two of its parameters to be redundant, as they follow linear trends with the scale radius in two distinct regimes of the void sample, separated by its compensation scale. Assuming linear theory, we derive an analytic formula for the velocity profile of voids and find an excellent agreement with the numerical data as well. In our companion paper [Sutter et al., Mon. Not. R. Astron. Soc. 442, 462 (2014)] the presented density profile is shown to be universal even across tracer type, properly describing voids defined in halo and galaxy distributions of varying sparsity, allowing us to relate various void populations by simple rescalings. This provides a powerful framework to match theory and simulations with observational data, opening up promising perspectives to constrain competing models of cosmology and gravity.

We present and study cosmic voids identified using the watershed void finder VIDE in the Sloan Digital Sky Survey Data Release 9, compare these voids to ones identified in mock catalogs, and assess the impact of the survey mask on void statistics such as number functions, ellipticity distributions, and radial density profiles. The nearly 1,000 identified voids span three nearly volume-limited samples from redshift z = 0.43 to 0.7. For comparison we use 98 of the publicly available 2LPT-based mock galaxy catalogs of Manera et al., and also generate our own mock catalogs by applying a Halo Occupation Distribution model to an N-body simulation. We find that the mask reduces the number density of voids at all scales by a factor of three and slightly skews the relative size distributions. This engenders an increase in the mean ellipticity by roughly 30%. However, we find that radial density profiles are largely robust to the effects of the mask. We see excellent agreement between the data and both mock catalogs, and find no tension between the observed void properties and the properties derived from {\Lambda}CDM simulations. We have added the void catalogs from both data and mock galaxy populations discussed in this work to the Public Cosmic Void Catalog at http://www.cosmicvoids.net.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

To study the impact of sparsity and galaxy bias on void statistics, we use a single large-volume, high-resolution N-body simulation to compare voids in multiple levels of subsampled dark matter, halo populations, and mock galaxies from a Halo Occupation Distribution model tuned to different galaxy survey densities. We focus our comparison on three key observational statistics: number functions, ellipticity distributions, and radial density profiles. We use the hierarchical tree structure of voids to interpret the impacts of sampling density and galaxy bias, and theoretical and empirical functions to describe the statistics in all our sample populations. We are able to make simple adjustments to theoretical expectations to offer prescriptions for translating from analytics to the void properties measured in realistic observations. We find that sampling density has a much larger effect on void sizes than galaxy bias. At lower tracer density, small voids disappear and the remaining voids are larger, more spherical, and have slightly steeper profiles. When a proper lower mass threshold is chosen, voids in halo distributions largely mimic those found in galaxy populations, except for ellipticities, where galaxy bias leads to higher values. We use the void density profile of Hamaus et al. (2014) to show that voids follow a self-similar and universal trend, allowing simple translations between voids studied in dark matter and voids identified in galaxy surveys. We have added the mock void catalogs used in this work to the Public Cosmic Void Catalog at http://www.cosmicvoids.net.

We report the first measurement of the diminutive lensing signal arising from matter underdensities associated with cosmic voids. While undetectable individually, by stacking the weak gravitational shear estimates around 901 voids detected in SDSS DR7 by Sutter et al. (2012a), we find substantial evidence for a depression of the lensing signal compared to the cosmic mean. This depression is most pronounced at the void radius, in agreement with analytical models of void matter profiles. Even with the largest void sample and imaging survey available today, we cannot put useful constraints on the radial dark-matter void profile. We invite independent investigations of our findings by releasing data and analysis code to the public at https://github.com/pmelchior/void-lensing

Galaxy bias, the unknown relationship between the clustering of galaxies and the underlying dark matter density field is a major hurdle for cosmological inference from large-scale structure. While traditional analyses focus on the absolute clustering amplitude of high-density regions mapped out by galaxy surveys, we propose a relative measurement that compares those to the underdense regions, cosmic voids. On the basis of realistic mock catalogs we demonstrate that cross correlating galaxies and voids opens up the possibility to calibrate galaxy bias and to define a static ruler thanks to the observable geometric nature of voids. We illustrate how the clustering of voids is related to mass compensation and show that volume-exclusion significantly reduces the degree of stochasticity in their spatial distribution. Extracting the spherically averaged distribution of galaxies inside voids from their cross correlations reveals a remarkable concordance with the mass-density profile of voids.

The two fundamental assumptions of the standard cosmological model - that the initial fluctuations are statistically isotropic and Gaussian - are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice and frequency dependence. Many of these anomalies were previously observed in the WMAP data, and are now confirmed at similar levels of significance (about 3 sigma). However, we find little evidence for non-Gaussianity, with the exception of a few statistical signatures that seem to be associated with specific anomalies. In particular, we find that the quadrupole-octopole alignment is also connected to a low observed variance of the CMB signal. A power asymmetry is now found to persist to scales corresponding to about l=600, and can be described in the low-l regime by a phenomenological dipole modulation model. However, any primordial power asymmetry is strongly scale-dependent and does not extend to arbitrarily small angular scales. Finally, it is plausible that some of these features may be reflected in the angular power spectrum of the data, which shows a deficit of power on similar scales. Indeed, when the power spectra of two hemispheres defined by a preferred direction are considered separately, one shows evidence for a deficit in power, while its opposite contains oscillations between odd and even modes that may be related to the parity violation and phase correlations also detected in the data. Although these analyses represent a step forward in building an understanding of the anomalies, a satisfactory explanation based on physically motivated models is still lacking.

How do observed voids relate to the underlying dark matter distribution? To examine the spatial distribution of dark matter contained within voids identified in galaxy surveys, we apply Halo Occupation Distribution models representing sparsely and densely sampled galaxy surveys to a high-resolution N-body simulation. We compare these galaxy voids to voids found in the halo distribution, low-resolution dark matter, and high-resolution dark matter. We find that voids at all scales in densely sampled surveys - and medium- to large-scale voids in sparse surveys - trace the same underdensities as dark matter, but they are larger in radius by ~20%, they have somewhat shallower density profiles, and they have centers offset by ~0.4Rv rms. However, in void-to-void comparison we find that shape estimators are less robust to sampling, and the largest voids in sparsely sampled surveys suffer fragmentation at their edges. We find that voids in galaxy surveys always correspond to underdensities in the dark matter, though the centers may be offset. When this offset is taken into account, we recover almost identical radial density profiles between galaxies and dark matter. All mock catalogs used in this work are available at http://www.cosmicvoids.net.

The Planck nominal mission cosmic microwave background (CMB) maps yield unprecedented constraints on primordial non-Gaussianity (NG). Using three optimal bispectrum estimators, separable template-fitting (KSW), binned, and modal, we obtain consistent values for the primordial local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result fNL^local= 2.7+/-5.8, fNL^equil= -42+/-75, and fNL^ortho= -25+\-39 (68% CL statistical). NG is detected in the data; using skew-C_l statistics we find a nonzero bispectrum from residual point sources, and the ISW-lensing bispectrum at a level expected in the LambdaCDM scenario. The results are based on comprehensive cross-validation of these estimators on Gaussian and non-Gaussian simulations, are stable across component separation techniques, pass an extensive suite of tests, and are confirmed by skew-C_l, wavelet bispectrum and Minkowski functional estimators. Beyond estimates of individual shape amplitudes, we present model-independent, 3-dimensional reconstructions of the Planck CMB bispectrum and thus derive constraints on early-Universe scenarios that generate primordial NG, including general single-field models of inflation, excited initial states (non-Bunch-Davies vacua), and directionally-dependent vector models. We provide an initial survey of scale-dependent feature and resonance models. These results bound both general single-field and multi-field model parameter ranges, such as the speed of sound, c_s \geq 0.02 (95% CL), in an effective field theory parametrization, and the curvaton decay fraction r_D \geq 0.15 (95% CL). The Planck data significantly limit the viable parameter space of the ekpyrotic/cyclic scenarios. The amplitude of the 4-point function in the local model tauNL < 2800 (95% CL). These constraints represent the highest precision tests to date of physical mechanisms for the origin of cosmic structure.

We present a novel, general-purpose method for deconvolving and denoising images from gridded radio interferometric visibilities using Bayesian inference based on a Gaussian process model. The method automatically takes into account incomplete coverage of the uv-plane, signal mode coupling due to the primary beam, and noise mode coupling due to uv sampling. Our method uses Gibbs sampling to efficiently explore the full posterior distribution of the underlying signal image given the data. We use a set of widely diverse mock images with a realistic interferometer setup and level of noise to assess the method. Compared to results from a proxy for point source- based CLEAN method we find that in terms of RMS error and signal-to-noise ratio our approach performs better than traditional deconvolution techniques, regardless of the structure of the source image in our test suite. Our implementation scales as O(np log np), provides full statistical and uncertainty information of the reconstructed image, requires no supervision, and provides a robust, consistent framework for incorporating noise and parameter marginalizations and foreground removal.

Based on CMB maps from the 2013 Planck Mission data release, this paper presents the detection of the ISW effect, i.e., the correlation between the CMB and large-scale evolving gravitational potentials. The significance of detection ranges from 2 to 4 sigma, depending on which method is used. We investigate three separate approaches, which cover essentially all previous studies, as well as breaking new ground. (i) Correlation of the CMB with the Planck reconstructed gravitational lensing potential (for the first time). This detection is made using the lensing-induced bispectrum; the correlation between lensing and the ISW effect has a significance close to 2.5 sigma. (ii) Cross-correlation with tracers of LSS, yielding around 3 sigma significance, based on a combination of radio (NVSS) and optical (SDSS) data. (iii) Aperture photometry on stacked CMB fields at the locations of known large-scale structures, which yields a 4 sigma signal when using a previously explored catalogue, but shows strong discrepancies in amplitude and scale compared to expectations. More recent catalogues give more moderate results, ranging from negligible to 2.5 sigma at most, but with a more consistent scale and amplitude, the latter being still slightly above what is expected from numerical simulations within LCMD. Where they can be compared, these measurements are compatible with previous work using data from WMAP, which had already mapped these scales to the limits of cosmic variance. Planck's broader frequency coverage confirms that the signal is achromatic, bolstering the case for ISW detection. As a final step we use tracers of large-scale structure to filter the CMB data, presenting maps of the ISW temperature perturbation. These results provide complementary and independent evidence for the existence of a dark energy component that governs the current accelerated expansion of the Universe.

Recently, Nadathur & Hotchkiss (2013) submitted a paper discussing a new cosmic void catalog. This paper includes claims about the void catalog described in Sutter et al. (2012). In this note, we respond to those claims, clarify some discrepancies between the text of Sutter et al. (2012) and the most recent version of the catalog, and provide some comments on the differences between our catalog and that of Nadathur & Hotchkiss (2013). All updates and documentation for our catalog are available at http://www.cosmicvoids.net.

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 investigate the impact of instrumental systematic errors in interferometric measurements of the cosmic microwave background (CMB) temperature and polarization power spectra. We simulate interferometric CMB observations to generate mock visibilities and estimate power spectra using the statistically optimal maximum likelihood technique. We define a quadratic error measure to determine allowable levels of systematic error that do not induce power spectrum errors beyond a given tolerance. As an example, in this study we focus on differential pointing errors. The effects of other systematics can be simulated by this pipeline in a straightforward manner. We find that, in order to accurately recover the underlying B-modes for r=0.01 at 28<l<384, Gaussian-distributed pointing errors must be controlled to 0.7^\circ rms for an interferometer with an antenna configuration similar to QUBIC, in agreement with analytical estimates. Only the statistical uncertainty for 28<l<88 would be changed at ~10% level. With the same instrumental configuration, we find the pointing errors would slightly bias the 2-\sigma upper limit of the tensor-to-scalar ratio r by ~10%. We also show that the impact of pointing errors on the TB and EB measurements is negligibly small.

The detection of the primordial $B$-mode spectrum of the polarized cosmic microwave background (CMB) signal may provide a probe of inflation. However, observation of such a faint signal requires excellent control of systematic errors. Interferometry proves to be a promising approach for overcoming such a challenge. In this paper we present a complete simulation pipeline of interferometric observations of CMB polarization, including systematic errors. We employ two different methods for obtaining the power spectra from mock data produced by simulated observations: the maximum likelihood method and the method of Gibbs sampling. We show that the results from both methods are consistent with each other, as well as, within a factor of 6, with analytical estimates. Several categories of systematic errors are considered: instrumental errors, consisting of antenna gain and antenna coupling errors, and beam errors, consisting of antenna pointing errors, beam cross-polarization and beam shape (and size) errors. In order to recover the tensor-to-scalar ratio, $r$, within a 10% tolerance level, which ensures the experiment is sensitive enough to detect the $B$-signal at $r=0.01$ in the multipole range $28 < \ell < 384$, we find that, for a QUBIC-like experiment, Gaussian-distributed systematic errors must be controlled with precisions of $|g_{rms}| = 0.1$ for antenna gain, $|\epsilon_{rms}| = 5 \times 10^{-4}$ for antenna coupling, $\delta_{rms} \approx 0.7^\circ$ for pointing, $\zeta_{rms} \approx 0.7^\circ$ for beam shape, and $\mu_{rms} = 5 \times 10^{-4}$ for beam cross-polarization.

Detection of B-mode polarization of the cosmic microwave background (CMB) radiation is one of the frontiers of observational cosmology. Because they are an order of magnitude fainter than E-modes, it is quite a challenge to detect B-modes. Having more manageable systematics, interferometers prove to have a substantial advantage over imagers in detecting such faint signals. Here, we present a method for Bayesian inference of power spectra and signal reconstruction from interferometric data of the CMB polarization signal by using the technique of Gibbs sampling. We demonstrate the validity of the method in the flat-sky approximation for a simulation of an interferometric observation on a finite patch with incomplete uv-plane coverage, a finite beam size and a realistic noise model. With a computational complexity of O(n^{3/2}), n being the data size, Gibbs sampling provides an efficient method for analyzing upcoming cosmology observations.

We report on the first application of the Alcock-Paczynski test to stacked voids in spectroscopic galaxy redshift surveys.We use voids from the Sutter et al. (2012) void catalog, which was derived from the Sloan Digital Sky Survey Data Release 7 main sample and luminous red galaxy catalogs. The construction of that void catalog removes potential shape measurement bias by using a modified version of the ZOBOV algorithm and by removing voids near survey boundaries and masks. We apply the shape-fitting procedure presented in Lavaux & Wandelt (2012) to ten void stacks out to redshift z=0.36. Combining these measurements, we determine the mean cosmologically induced "stretch" of voids in three redshift bins, with 1-sigma errors of 5-15%. The mean stretch is consistent with unity, providing no indication of a distortion induced by peculiar velocities. While the statistical errors are too large to detect the Alcock-Paczynski effect over our limited redshift range, this proof-of-concept analysis defines procedures that can be applied to larger spectroscopic galaxy surveys at higher redshifts to constrain dark energy using the expected statistical isotropy of structures that are minimally affected by uncertainties in galaxy velocity bias.

We present a general method for accelerating by more than an order of magnitude the convolution of pixelated function on the sphere with a radially-symmetric kernel. Our method splits the kernel into a compact real-space, and a compact spherical harmonic space component that can then be convolved in parallel using an inexpensive commodity GPU and a CPU, respectively. We provide models for the computational cost of both real-space and Fourier space convolutions and an estimate for the approximation error. Using these models we can determine the optimum split that minimizes the wall clock time for the convolution while satisfying the desired error bounds. We apply this technique to the problem of simulating a cosmic microwave background sky map at the resolution typical of the high resolution maps of the cosmic microwave background anisotropies produced by the Planck space craft. For the main Planck CMB science channels we achieve a speedup of over a factor of ten, assuming an acceptable fractional rms error of order 10^-5 in the (power spectrum of the) output map.

We produce the most comprehensive public void catalog to date using the Sloan Digital Sky Survey Data Release 7 main sample out to redshift z=0.2 and the luminous red galaxy sample out to z=0.44. Using a modified version of the parameter-free void finder ZOBOV, we fully take into account the presence of the survey boundary and masks. Our strategy for finding voids is thus appropriate for any survey configuration. We produce two distinct catalogs: a complete catalog including voids near any masks, which would be appropriate for void galaxy surveys, and a bias-free catalog of voids away from any masks, which is necessary for analyses that require a fair sampling of void shapes and alignments. Our discovered voids have effective radii from 5 to 135 h^-1 Mpc. We discuss basic catalog statistics such as number counts and redshift distributions and describe some additional data products derived from our catalog, such as radial density profiles and projected density maps. We find that radial profiles of stacked voids show a qualitatively similar behavior across nearly two decades of void radii and throughout the full redshift range.

We present a first estimate based on a cosmological gasdynamics simulation of galaxy cluster radio halo counts to be expected in forthcoming low-frequency radio surveys. Our estimate is based on a FLASH simulation of the LCDM model for which we have assigned radio power to clusters via a model that relates radio emissivity to cluster magnetic field strength, intracluster turbulence, and density. We vary several free parameters of this model and find that radio halo number counts vary by up to a factor of two for average magnetic fields ranging from 0.2 to 3.1 uG. However, we predict significantly fewer low-frequency radio halos than expected from previous semi-analytic estimates, although this discrepancy could be explained by frequency-dependent radio halo probabilities as predicted in reacceleration models. We find that upcoming surveys will have difficulty in distinguishing models because of large uncertainties and low number counts. Additionally, according to our modeling we find that expected number counts can be degenerate with both reacceleration and hadronic secondary models of cosmic ray generation. We find that relations between radio power and mass and X-ray luminosity may be used to distinguish models, and by building mock radio sky maps we demonstrate that surveys such as LOFAR may have sufficient resolution and sensitivity to break this model degeneracy by imaging many individual clusters.

We present a Bayesian angular power spectrum and signal map inference engine which can be adapted to interferometric observations of anisotropies inthe cosmic microwave background, 21 cm emission line mapping of galactic brightness fluctuations, or 21 cm absorption line mapping of neutral hydrogen in the dark ages. The method uses Gibbs sampling to generate a sampled representation of the angular power spectrum posterior and the posterior of signal maps given a set of measured visibilities in the uv-plane. We use a mock interferometric CMB observation to demonstrate the validity of this method in the flat-sky approximation when adapted to take into account arbitrary coverage of the uv-plane, mode-mode correlations due to observations on a finite patch, and heteroschedastic visibility errors. The computational requirements scale as O(n_p log n_p) where n_p measures the ratio of the size of the detector array to the inter-detector spacing, meaning that Gibbs sampling is a promising technique for meeting the data analysis requirements of future cosmology missions.

Cosmological constraints derived from galaxy clusters rely on accurate predictions of cluster observable properties, in which feedback from active galactic nuclei (AGN) is a critical component. In order to model the physical effects due to supermassive black holes (SMBH) on cosmological scales, subgrid modeling is required, and a variety of implementations have been developed in the literature. However, theoretical uncertainties due to model and parameter variations are not yet well understood, limiting the predictive power of simulations including AGN feedback. By performing a detailed parameter sensitivity study in a single cluster using several commonly-adopted AGN accretion and feedback models with FLASH, we quantify the model uncertainties in predictions of cluster integrated properties. We find that quantities that are more sensitive to gas density have larger uncertainties (~20% for Mgas and a factor of ~2 for Lx at R500), whereas Tx, Ysz, and Yx are more robust (~10-20% at R500). To make predictions beyond this level of accuracy would require more constraints on the most relevant parameters: the accretion model, mechanical heating efficiency, and size of feedback region. By studying the impact of AGN feedback on the scaling relations, we find that an anti-correlation exists between Mgas and Tx, which is another reason why Ysz and Yx are excellent mass proxies. This anti-correlation also implies that AGN feedback is likely to be an important source of intrinsic scatter in the Mgas-Tx and Lx-Tx relations.

We present 3D adaptive mesh refinement MHD simulations of an isolated galaxy cluster that include injection of kinetic, thermal, and magnetic energy via a central active galactic nucleus (AGN) in order to study and evaluate the role that AGN may play in producing the observed cluster-wide magnetic fields. Using the MHD solver in FLASH 3.3, we compare several sub-resolution approaches to the evolution of AGN, specifically focusing on large-scale jet and bubble models. We examine the effects of magnetized outflows on the accretion history of the black hole and cluster thermodynamic properties, discuss the ability of various models to magnetize the cluster medium, and assess the sensitivity of these models to their underlying subgrid parameters. We find that magnetized jet-based models suffer a severe reduction in accretion rate compared to hydrodynamic jets; however, bubble models remain largely unaffected. While both jets and sporadically-placed bubbles have difficulty reproducing the observed strength and topology of cluster magnetic fields, models based on centrally-located bubbles come closest to observations. Finally, whereas jet models are relatively insensitive to changes in their subgrid parameters, the accretion rate and average magnetic field produced by the bubbles vary by as much as an order of magnitude depending on the grid resolution and accretion strength.

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

While supermassive black holes (SMBHs) play an important role in galaxy and cluster evolution, at present they can only be included in large-scale cosmological simulation via subgrid techniques. However, these subgrid models have not been studied in a systematic fashion. Using a newly-developed fast, parallel spherical overdensity halo finder built into the simulation code FLASH, we perform a suite of dark matter-only cosmological simulations to study the effects of subgrid model choice on relations between SMBH mass and dark matter halo mass and velocity dispersion. We examine three aspects of SMBH subgrid models: the choice of initial black hole seed mass, the test for merging two black holes, and the frequency of applying the subgrid model. We also examine the role that merging can play in determining the relations, ignoring the complicating effects of SMBH-driven accretion and feedback. We find that the choice of subgrid model can dramatically affect the black hole merger rate, the cosmic SMBH mass density, and the low-redshift relations to halo properties. We also find that it is possible to reproduce observations of the low-redshift relations without accretion and feedback, depending on the choice of subgrid model.

Active galactic nuclei (AGN) found at the centers of clusters of galaxies are a possible source for weak cluster-wide magnetic fields. To evaluate this scenario, we present 3D adaptive mesh refinement MHD simulations of a cool-core cluster that include injection of kinetic, thermal, and magnetic energy via an AGN-powered jet. Using the MHD solver in FLASH 2, we compare several sub-resolution approaches that link the estimated accretion rate as measured on the simulation mesh to the accretion rate onto the central black hole and the resulting feedback. We examine the effects of magnetized outflows on the accretion history of the black hole and discuss the ability of these models to magnetize the cluster medium.

Galaxy clusters are valuable cosmological probes. However, cluster mass estimates rely on observable quantities that are affected by complicated baryonic physics in the intracluster medium (ICM), including feedback from active galactic nuclei (AGN). Cosmological simulations have started to include AGN feedback using subgrid models. In order to make robust predictions, the systematics of different implementations and parametrizations need to be understood. We have developed an AGN subgrid model in FLASH that supports a few different black hole accretion models and feedback models. We use this model to study the effect of AGN on X-ray cluster observables and its dependence on model variations.

Using a hydrodynamics plus N-body simulation of galaxy cluster formation within a large volume and mock Chandra X-ray observations, we study the form and evolution of the intrinsic scatter about the best-fit X-ray temperature-mass relation for clusters. We investigate the physical origin of the scatter by correlating it with quantities that are closely related to clusters' formation and merging histories. We also examine the distribution of the scatter for merging and nonmerging populations, identified using halo merger trees derived from the simulation as well as X-ray substructure measures. We find a strong correlation between the scatter in the M-T_X relation and the halo concentration, in the sense that more concentrated clusters tend to be cooler than clusters with similar masses. No bias is found between the merging and relaxed clusters, but merging clusters generally have greater scatter, which is related to the properties of the distribution of halo concentrations. We also detect a signature of non-lognormality in the distribution of scatter for our simulated clusters both at z=0 and at z=1. A detailed comparison of merging clusters identified by substructure measures and by halo merger trees is given in the discussion. We conclude that, when cooling-related effects are neglected, the variation in halo concentrations is a more important factor for driving the intrinsic scatter in the M-T_X relation, while departures from hydrostatic equilibrium due to cluster mergers have a minor effect.

We use high-resolution simulations of large-scale structure formation to analyze the effects of interacting dark matter and dark energy on the evolution of the halo mass function. Using a chi-square likelihood analysis, we find significant differences in the mass function between models of coupled dark matter-dark energy and standard concordance cosmology Lambda-CDM out to redshift z=1.5. We also find a preliminary indication that the Dark Energy Survey should be able to distinguish these models from Lambda-CDM within its mass and redshift contraints. While we can distinguish the effects of these models from Lambda-CDM cosmologies with different fundamental parameters, DES will require independent measurements of sigma-8 to confirm these effects.

We examine how coupled dark matter and dark energy modify the development of Zel'dovich pancakes. We study how the various effects of these theories, such as a fifth force in the dark sector and a modified particle Hubble drag, produce variations in the redshifts of caustic formation and the present-day density profiles of pancakes. We compare our results in direct simulation to a perturbation theory approach for the dark energy scalar field. We determine the range of initial scalar field amplitudes for which perturbation theory is accurate in describing the development of the pancakes. Notably, we find that perturbative methods which neglect kinetic terms in the scalar field equation of motion are not valid for arbitrarily small perturbations. We also examine whether models that have been tuned to match the constraints of current observations can produce new observable effects in the nonlinear structure of pancakes. Our results suggest that a fully realistic three-dimensional simulation will produce significant new observable features, such as modifications to the mass function and halo radial density profile shapes, that can be used to distinguish these models from standard concordance cosmology and from each other.