We examine the cosmological implications of measurements of the void-galaxy cross-correlation at redshift $z=0.57$ combined with baryon acoustic oscillation (BAO) data at $0.1<z<2.4$. We find direct evidence of the late-time acceleration due to dark energy at $>10\sigma$ significance from these data alone, independent of the cosmic microwave background and supernovae. Using a nucleosynthesis prior on $\Omega_bh^2$, we measure the Hubble constant to be $H_0=72.3\pm1.9\;{\rm km\,s}^{-1}{\rm Mpc}^{-1}$ from BAO+voids at $z<2$, and $H_0=69.0\pm1.2\;{\rm km\,s}^{-1}{\rm Mpc}^{-1}$ when adding Lyman-$\alpha$ BAO at $z=2.34$, both independent of the CMB. Adding voids to CMB, BAO and supernova data greatly improves measurement of the dark energy equation of state, increasing the figure of merit by >40%, but remaining consistent with flat flat $\Lambda$ cold dark matter.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We develop an estimator to extract the mean radial velocity profile of cosmic voids via the kinematic Sunyaev-Zeldovich effect of pairs of galaxies surrounding them. The estimator is tested with simulated pure kSZ map and void catalogue data from the same simulation. The results show that the recovered signal could be attenuated by low angular resolution of the map or large aperture photometry filter radius size, but the mean radial velocity profile can be fully recovered with our estimator. By applying the estimator to the Planck 2D-ILC CMB map, with galaxy and void catalogues from BOSS DR12, we find that the estimated void velocity profile is $3.31\sigma$ apart from null detection for for voids with continuously rising density profiles asymptoting to the mean density; and $1.75\sigma$ for voids with positive density contrast shell surrounded. By fitting the reconstructed data to the theoretical profile, we find the reduced $\chi^{2}$ to be $1.19$ and $0.62$ for the two types of void, respectively, indicating a good fit of the model to the data. We then forecast the detectability of the radial velocity profile of cosmic voids with future CMB surveys, including SPT-3G, AdvACT, and Simons Observatory. We find that the contamination effect from CMB residuals is negligible with survey area over $2,000~{\rm deg}^2$, especially with aperture photometry size less than $1\,{\rm arcmin}$. But the effect from instrumental noise is non-negligible. For future SPT-3G, AdvACT and Simons Observatory, the detection is potentially achievable from $3\sigma$ to $10\sigma$ C.L., depending on specific instrumental parameters. This opens a new window of probing dynamics of the cosmic structures from the kinematic Sunyaev-Zeldovich effect.

We report a $5.3\sigma$ detection of the gravitational lensing effect of cosmic voids from the Baryon Oscillation Spectroscopic (BOSS) Data Release 12 seen in the $Planck$ 2018 cosmic microwave background (CMB) lensing convergence map. To make this detection, we introduce new optimal techniques for void stacking and filtering of the CMB maps, such as binning voids by a combination of their observed galaxy density and size to separate those with distinctive lensing signatures. We calibrate theoretical expectations for the void-lensing signal using mock catalogs generated in a suite of 108 full-sky lensing simulations from Takahashi et al. (2017). Relative to these templates, we measure the lensing amplitude parameter in the data to be $A_L=1.10 \pm 0.21$ using a matched-filter stacking technique, and confirm it using an alternative Wiener filtering method. We demonstrate that the result is robust against thermal Sunyaev-Zel'dovich contamination and other sources of systematics. We use the lensing measurements to test the relationship between the matter and galaxy distributions within voids, and show that the assumption of linear bias with a value consistent with galaxy clustering results is discrepant with observation at $\sim 3\sigma$; we explain why such a result is consistent with simulations and previous results, and is expected as a consequence of void selection effects. We forecast the potential for void-CMB lensing measurements in future data from the Advanced ACT, Simons Observatory and CMB-S4 experiments, showing that, for the same number of voids, the achievable precision improves by a factor of more than two compared to $Planck$.

Cosmic voids gravitationally lens the cosmic microwave background (CMB) radiation, resulting in a distinct imprint on degree scales. We use the simulated CMB lensing convergence map from the MICE N-body simulation to calibrate our detection strategy for a given void definition and galaxy tracer density. We then identify cosmic voids in DES Year 1 data and stack the Planck 2015 lensing convergence map on their locations, probing the consistency of simulated and observed void lensing signals. When fixing the shape of the stacked convergence profile to that calibrated from simulations, we find imprints at the $3{\sigma}$ significance level for various analysis choices. The best measurement strategies based on the MICE calibration process yield $S/N \sim 4$ for DES Y1, and the best-fit amplitude recovered from the data is consistent with expectations from MICE ($A \sim 1$). Given these results as well as the agreement between them and N-body simulations, we conclude that the previously reported excess integrated Sachs-Wolfe (ISW) signal associated with cosmic voids in DES Y1 has no counterpart in the Planck CMB lensing map.

We present catalogues of cosmic voids identified in the distribution of Luminous Red Galaxies (LRGs) and Quasi Stellar Objects (QSOs) in the fourteenth data release (DR14) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We perform a multivariate analysis to assess the level of contamination in these catalogues by spurious Poisson underdensities. We find that the LRG void catalogue is largely free from contamination but that the QSO catalogue may be heavily contaminated. We analyse the multipoles of the void-galaxy cross-correlation function in these catalogues to obtain constraints on the growth rate of structure around voids. We find a value of $\beta(z=0.703)=0.58^{+0.33}_{-0.28}$ for the LRG voids and $\beta(z=1.53)=0.15^{+0.13}_{-0.12}$ for the QSO voids.

We revisit a cosmological constraint on dark matter decaying into dark radiation at late times. In Enqvist et al. (2015), we mainly focused on the effects of decaying dark matter (DDM) on the cosmic microwave background (CMB) and nonlinear matter power spectrum. Extending our previous analysis, here we use N-body simulation to investigate how DDM affects the halo mass function. This allows us to incorporate the cluster counts observed by the Sunyaev-Zel'dovich effect to study a bound on the lifetime of DDM. We also update the data of CMB and cosmic shear power spectrum with the Planck 2015 results and KiDS450 observations, respectively. From these cosmological observations, we obtain an lower bound on the lifetime $\Gamma^{-1}\ge 175\,$Gyr from the Planck2015 results (CMB+SZ cluster count) combined with the KiDS450 and the recent measurements of the baryon acoustic scale.

We present a measurement of the anisotropic void-galaxy cross-correlation function in the CMASS galaxy sample of the BOSS DR12 data release. We perform a joint fit to the data for redshift space distortions (RSD) due to galaxy peculiar velocities and anisotropies due to the Alcock-Paczynski (AP) effect, for the first time using a velocity field reconstruction technique to remove the complicating effects of RSD in the void centre positions themselves. Fits to the void-galaxy function give a 1% measurement of the AP parameter combination $D_A(z)H(z)/c = 0.4367\pm 0.0045$ at redshift $z=0.57$, where $D_A$ is the angular diameter distance and $H$ the Hubble parameter, exceeding the precision obtainable from baryon acoustic oscillations (BAO) by a factor of ~3.5 and free of systematic errors. From voids alone we also obtain a 10% measure of the growth rate, $f\sigma_8(z=0.57)=0.501\pm0.051$. The parameter degeneracies are orthogonal to those obtained from galaxy clustering. Combining void information with that from BAO and galaxy RSD in the same CMASS sample, we measure $D_A(0.57)/r_s=9.383\pm 0.077$ (at 0.8% precision), $H(0.57)r_s=(14.05\pm 0.14)\;10^3$ kms$^{-1}$Mpc$^{-1}$ (1%) and $f\sigma_8=0.453\pm0.022$ (4.9%), consistent with cosmic microwave background (CMB) measurements from Planck. These represent a factor \sim2 improvement in precision over previous results through the inclusion of void information. Fitting a flat cosmological constant $\Lambda$CDM model to these results in combination with Planck CMB data, we find up to an 11% reduction in uncertainties on $H_0$ and $\Omega_m$ compared to use of the corresponding BOSS consensus values. Constraints on extended models with non-flat geometry and a dark energy of state that differs from $w=-1$ show an even greater improvement.

The largest structures in the cosmic web probe the dynamical nature of dark energy through their integrated Sachs-Wolfe imprints. In the strength of the signal, typical cosmic voids have shown good consistency with expectation $A_{\rm ISW}=\Delta T^{\rm data} / \Delta T^{\rm theory}=1$, given the substantial cosmic variance. Discordantly, large-scale hills in the gravitational potential, or supervoids, have shown excess signals. In this study, we mapped out 87 new supervoids in the total 5000 deg$^2$ footprint of the Dark Energy Survey at $0.2<z<0.9$ to probe these anomalous claims. We found an excess imprinted profile with $ A_{\rm ISW}\approx4.1\pm2.0$ amplitude. The combination with independent BOSS data reveals an ISW imprint of supervoids at the $3.3\sigma$ significance level with an enhanced $A_{\rm ISW}\approx5.2\pm1.6$ amplitude. The tension with $\Lambda$CDM predictions is equivalent to $2.6\sigma$ and remains unexplained.

Redshift space distortions within voids provide a unique method to test for environmental dependence of the growth rate of structures in low density regions, where effects of modified gravity theories might be important. We derive a linear theory model for the redshift space void-galaxy correlation that is valid at all pair separations, including deep within the void, and use this to obtain expressions for the monopole $\xi^s_0$ and quadrupole $\xi^s_2$ contributions. Our derivation highlights terms that have previously been neglected but are important within the void interior. As a result our model differs from previous works and predicts new physical effects, including a change in the sign of the quadrupole term within the void radius. We show how the model can be generalised to include a velocity dispersion. We compare our model predictions to measurements of the correlation function using mock void and galaxy catalogues modelled after the BOSS CMASS galaxy sample using the Big MultiDark $N$-body simulation, and show that the linear model with dispersion provides an excellent fit to the data at all scales, $0\leq s\leq120\;h^{-1}$Mpc. While the RSD model matches simulations, the linear bias approximation does not hold within voids, and care is needed in fitting for the growth rate. We show that fits to the redshift space correlation recover the growth rate $f(z=0.52)$ to a precision of $2.7\%$ using the simulation volume of $(2.5\;h^{-1}\mathrm{Gpc})^3$.

Redshift space distortions (RSD) in the void-galaxy correlation $\xi^s$ provide information on the linear growth rate of structure in low density environments. Accurate modelling of these RSD effects can also allow the use of voids in competitive Alcock-Paczynski measurements. Linear theory models of $\xi^s$ are able to provide extremely good descriptions of simulation data on all scales provided the real space void positions are known. However, by reference to simulation data we demonstrate the failure of the assumptions implicit in current models of $\xi^s$ for voids identified directly in redshift space, as would be simplest using real observational data. To overcome this problem we instead propose using a density-field reconstruction method based on the Zeldovich approximation to recover the real space void positions from redshift space data. We show that this recovers the excellent agreement between theory and data for $\xi^s$. Performing the reconstruction requires an input cosmological model so, to be self-consistent, we have to perform reconstruction for every model to be tested. We apply this method to mock galaxy and void catalogues in the Big MultiDark $N$-body simulation and consistently recover the fiducial growth rate to a precision of $3.4\%$ using the simulation volume of $(2.5\;h^{-1}\mathrm{Gpc})^3$.

The properties of large underdensities in the distribution of galaxies in the Universe, known as cosmic voids, are potentially sensitive probes of fundamental physics. We use data from the MultiDark suite of N-body simulations and multiple halo occupation distribution mocks to study the relationship between galaxy voids, identified using a watershed void-finding algorithm, and the gravitational potential $\Phi$. We find that the majority of galaxy voids correspond to local density minima in larger-scale overdensities, and thus lie in potential wells. However, a subset of voids can be identified that closely trace maxima of the gravitational potential and thus stationary points of the velocity field. We identify a new void observable, $\lambda_v$, which depends on a combination of the void size and the average galaxy density contrast within the void, and show that it provides a good proxy indicator of the potential at the void location. A simple linear scaling of $\Phi$ as a function of $\lambda_v$ is found to hold, independent of the redshift and properties of the galaxies used as tracers of voids. We provide an accurate fitting formula to describe the spherically averaged potential profile $\Phi(r)$ about void centre locations. We discuss the importance of these results for the understanding of the evolution history of voids, and for their use in precision measurements of the integrated Sachs-Wolfe effect, gravitational lensing and peculiar velocity distortions in redshift space.

Small temperature anisotropies in the Cosmic Microwave Background can be sourced by density perturbations via the late-time integrated Sachs-Wolfe effect. Large voids and superclusters are excellent environments to make a localized measurement of this tiny imprint. In some cases excess signals have been reported. We probed these claims with an independent data set, using the first year data of the Dark Energy Survey in a different footprint, and using a different super-structure finding strategy. We identified 52 large voids and 102 superclusters at redshifts $0.2 < z < 0.65$. We used the Jubilee simulation to a priori evaluate the optimal ISW measurement configuration for our compensated top-hat filtering technique, and then performed a stacking measurement of the CMB temperature field based on the DES data. For optimal configurations, we detected a cumulative cold imprint of voids with $\Delta T_{f} \approx -5.0\pm3.7~\mu K$ and a hot imprint of superclusters $\Delta T_{f} \approx 5.1\pm3.2~\mu K$ ; this is $\sim1.2\sigma$ higher than the expected $|\Delta T_{f}| \approx 0.6~\mu K$ imprint of such super-structures in $\Lambda$CDM. If we instead use an a posteriori selected filter size ($R/R_{v}=0.6$), we can find a temperature decrement as large as $\Delta T_{f} \approx -9.8\pm4.7~\mu K$ for voids, which is $\sim2\sigma$ above $\Lambda$CDM expectations and is comparable to previous measurements made using SDSS super-structure data.

We present a new method for detection of the integrated Sachs-Wolfe (ISW) imprints of cosmic superstructures on the cosmic microwave background, based on a matched filtering approach. The expected signal-to-noise ratio for this method is comparable to that obtained from the full cross-correlation, and unlike other stacked filtering techniques it is not subject to an a posteriori bias. We apply this method to Planck CMB data using voids and superclusters identified in the CMASS galaxy data from the Sloan Digital Sky Survey Data Release 12, and measure the ISW amplitude to be $A_\mathrm{ISW}=1.64\pm0.53$ relative to the $\Lambda$CDM expectation, corresponding to a $3.1\sigma$ detection. In contrast to some previous measurements of the ISW effect of superstructures, our result is in agreement with the $\Lambda$CDM model.

Galaxies and their dark matter halos populate a complicated filamentary network around large, nearly empty regions known as cosmic voids. Cosmic voids are usually identified in spectroscopic galaxy surveys, where 3D information about the large-scale structure of the Universe is available. Although an increasing amount of photometric data is being produced, its potential for void studies is limited since photometric redshifts induce line-of-sight position errors of $\sim50$ Mpc/$h$ or more that can render many voids undetectable. In this paper we present a new void finder designed for photometric surveys, validate it using simulations, and apply it to the high-quality photo-$z$ redMaGiC galaxy sample of the Dark Energy Survey Science Verification (DES-SV) data. The algorithm works by projecting galaxies into 2D slices and finding voids in the smoothed 2D galaxy density field of the slice. Fixing the line-of-sight size of the slices to be at least twice the photo-$z$ scatter, the number of voids found in these projected slices of simulated spectroscopic and photometric galaxy catalogs is within 20% for all transverse void sizes, and indistinguishable for the largest voids of radius $\sim 70$ Mpc/$h$ and larger. The positions, radii, and projected galaxy profiles of photometric voids also accurately match the spectroscopic void sample. Applying the algorithm to the DES-SV data in the redshift range $0.2<z<0.8$, we identify 87 voids with comoving radii spanning the range 18-120 Mpc/$h$, and carry out a stacked weak lensing measurement. With a significance of $4.4\sigma$, the lensing measurement confirms the voids are truly underdense in the matter field and hence not a product of Poisson noise, tracer density effects or systematics in the data. It also demonstrates, for the first time in real data, the viability of void lensing studies in photometric surveys.

We present a public catalogue of voids in the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 LOWZ and CMASS galaxy surveys. This catalogue contains information on the location, sizes, densities, shapes and bounding surfaces of 8956 independent, disjoint voids, making it the largest public void catalogue to date. Voids are identified using a version of the ZOBOV algorithm, the operation of which has been calibrated through tests on mock galaxy populations in N-body simulations, as well as on a suite of 4096 mock catalogues which fully reproduce the galaxy clustering, survey masks and selection functions. Based on this, we estimate a false positive detection rate of 3%. Comparison with mock catalogues limits deviations of the void size distribution from that predicted in the $\Lambda$CDM model to be less than 6% for voids with effective radius $8<R_v<60\,h^{-1}$Mpc and in the redshift range $0.15<z<0.7$. This could tightly constrain modified gravity scenarios and models with a varying equation of state, but we identify systematic biases which must be accounted for to reduce the theoretical uncertainty in the predictions for these models to the current level of precision attained from the data. We also examine the distribution of void densities and identify a deficit of the deepest voids relative to $\Lambda$CDM expectations, which is significant at more than the $3\sigma$ equivalent level. We discuss possible explanations for this discrepancy but at present its cause remains unknown.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

The statistical study of voids in the matter distribution promises to be an important tool for precision cosmology, but there are known discrepancies between theoretical models of voids and the voids actually found in large simulations or galaxy surveys. The empirical properties of observed voids are also not well understood. In this paper, we study voids in an N-body simulation, using the ZOBOV watershed algorithm. As in other studies, we use sets of subsampled dark matter particles as tracers to identify voids, but we use the full-resolution simulation output to measure dark matter densities at the identified locations. Voids span a wide range of sizes and densities, but there is a clear trend towards larger voids containing deeper density minima, a trend which is expected for all watershed void finders. We also find that the tracer density at void locations is usually smaller than the true density, and that this relationship depends on the sampling density of tracers. We show that fits given in the literature fail to match the observed density profiles of voids. The average enclosed density contrast within watershed voids varies widely with both the size of the void and the minimum density within it, but is always far from the shell crossing threshold expected from theoretical models. Voids with deeper density minima also show much broader density profiles. We discuss the implications of these results for the excursion set approach to modelling such voids.

We consider decaying dark matter (DDM) as a resolution to the possible tension between cosmic microwave background (CMB) and weak lensing (WL) based determinations of the amplitude of matter fluctuations, $\sigma_8$. We perform N-body simulations in a model where dark matter decays into dark radiation and develop an accurate fitting formula for the non-linear matter power spectrum, which enables us to test the DDM model by the combined measurements of CMB, WL and the baryon acoustic oscillation (BAO). We employ a Markov chain Monte Carlo analysis to examine the overlap of posterior distributions of the cosmological parameters, comparing CMB alone with WL+BAO. We find an overlap that is significantly larger in the DDM model than in the standard CDM model. This may be hinting at DDM, although current data is not constraining enough to unambiguously favour a non-zero dark matter decay rate $\Gamma$. From the combined CMB+WL data, we obtain a lower bound $\Gamma^{-1}\ge 97$ Gyr at 95 % C.L, which is less tight than the constraint from CMB alone.

We study how the properties of cosmic voids depend on those of the tracer galaxy populations in which they are identified. We use a suite of halo occupation distribution (HOD) mocks in a simulation, identify voids in these populations using the ZOBOV void finder and measure their abundances, sizes, tracer densities, and dark matter content. To separate the effects of bias from those of sampling density, we do the same for voids traced by randomly down-sampled subsets of the simulation dark matter particles. At the same sampling density, galaxy bias reduces the total number of voids by ~50% and can dramatically change their size distribution. The matter content of voids in biased and unbiased tracers also differs. Deducing void properties from simulation therefore requires the use of realistic galaxy mocks. We discuss how the void observables can be related to their matter content. In particular we consider the compensation of the total mass deficit in voids and find that the distinction between over- and under-compensated voids is not a function of void size alone, as has previously been suggested. However, we find a simple linear relationship between the average density of tracers in the void and the total mass compensation on much larger scales. The existence of this linear relationship holds independent of the bias and sampling density of the tracers. This provides a universal tool to classify void environments and will be important for the use of voids in observational cosmology.

We examine the integrated Sachs-Wolfe (ISW) imprint of voids and superclusters on the cosmic microwave background. We first study results from the Jubilee $N$-body simulation. From Jubilee, we obtain the full-sky ISW signal from structures out to redshift $z=1.4$ and a mock luminous red galaxy (LRG) catalogue. We confirm that the expected signal in the concordance \Lambda CDM model is very small and likely to always be much smaller than the anisotropies arising at the last scattering surface. Any current detections of such an imprint must, therefore, predominantly arise from something other than an ISW effect in a \Lambda CDM universe. Using the simulation as a guide, we then look for the signal using a catalogue of voids and superclusters from the Sloan Digital Sky Survey. We find a result that is consistent with the \Lambda CDM model, i.e. a signal consistent with zero.

We discuss the universality and self-similarity of void density profiles, for voids in realistic mock luminous red galaxy (LRG) catalogues from the Jubilee simulation, as well as in void catalogues constructed from the SDSS LRG and Main Galaxy samples. Voids are identified using a modified version of the ZOBOV watershed transform algorithm, with additional selection cuts. We find that voids in simulation are self-similar, meaning that their average rescaled profile does not depend on the void size, or -- within the range of the simulated catalogue -- on the redshift. Comparison of the profiles obtained from simulated and real voids shows an excellent match. The profiles of real voids also show a universal behaviour over a wide range of galaxy luminosities, number densities and redshifts. This points to a fundamental property of the voids found by the watershed algorithm, which can be exploited in future studies of voids.

The discovery of a void of size $\sim200\;h^{-1}$Mpc and average density contrast of $\sim-0.1$ aligned with the Cold Spot direction has been recently reported. It has been argued that, although the first-order integrated Sachs-Wolfe (ISW) effect of such a void on the CMB is small, the second-order Rees-Sciama (RS) contribution exceeds this by an order of magnitude and can entirely explain the observed Cold Spot temperature profile. In this paper we examine this surprising claim using both an exact calculation with the spherically symmetric Lema\^itre-Tolman-Bondi metric, and perturbation theory about a background Friedmann-Robertson-Walker (FRW) metric. We show that both approaches agree well with each other, and both show that the dominant temperature contribution of the postulated void is an unobservable dipole anisotropy. If this dipole is subtracted, we find that the remaining temperature anisotropy is dominated by the linear ISW signal, which is orders of magnitude larger than the second-order RS effect, and that the total magnitude is too small to explain the observed Cold Spot profile. We calculate the density and size of a void that would be required to explain the Cold Spot, and show that the probability of existence of such a void is essentially zero in $\Lambda$CDM. We identify the importance of \emph{a posteriori} selection effects in the identification of the Cold Spot, but argue that even after accounting for them, a supervoid explanation of the Cold Spot is always disfavoured relative to a random statistical fluctuation on the last scattering surface.

The stacked density profile of cosmic voids in the galaxy distribution provides an important tool for the use of voids for precision cosmology. We study the density profiles of voids identified using the ZOBOV watershed transform algorithm in realistic mock luminous red galaxy (LRG) catalogues from the Jubilee simulation, as well as in void catalogues constructed from the SDSS LRG and Main Galaxy samples. We compare different methods for reconstructing density profiles scaled by the void radius and show that the most commonly used method based on counts in shells and simple averaging is statistically flawed as it underestimates the density in void interiors. We provide two alternative methods that do not suffer from this effect; one based on Voronoi tessellations is also easily able to account from artefacts due to finite survey boundaries and so is more suitable when comparing simulation data to observation. Using this method we show that voids in simulation are exactly self-similar, meaning that their average rescaled profile does not depend on the void size. Within the range of our simulation we also find no redshift dependence of the profile. Comparison of the profiles obtained from simulated and real voids shows an excellent match. The profiles of real voids also show a universal behaviour over a wide range of galaxy luminosities, number densities and redshifts. This points to a fundamental property of the voids found by the watershed algorithm, which can be exploited in future studies of voids.

The study of the interesting cosmological properties of voids in the Universe depends on the efficient and robust identification of such voids in galaxy redshift surveys. Recently, Sutter et al. (2012) have published a public catalogue of voids in the Sloan Digital Sky Survey Data Release 7 main galaxy and luminous red galaxy samples, using the void-finding algorithm ZOBOV, which is based on the watershed transform. We examine the properties of this catalogue and show that it suffers from several problems and inconsistencies, including the identification of some extremely overdense regions as voids. As a result, cosmological results obtained using this catalogue need to be reconsidered. We provide instead an alternative, self-consistent, public catalogue of voids in the same galaxy data, obtained from using an improved version of the same watershed transform algorithm. We provide a more robust method of dealing with survey boundaries and masks, as well as with a radially varying selection function, which means that our method can be applied to any other survey. We discuss some basic properties of the voids thus discovered, and describe how further information may be obtained from the catalogue. In addition, we apply an inversion of the algorithm to the same data to obtain a corresponding catalogue of large-scale overdense structures, or "superclusters". Our catalogues are available for public download at www.hip.fi/nadathur/download/dr7catalogue.

Sutter et al. have responded to the criticisms we made of their cosmic void catalogue in our recent paper presenting an alternative catalogue (Nadathur & Hotchkiss 2013). Unfortunately, their response contains several statements which are incorrect, as we point out in this note.

We present initial results from the Jubilee ISW project, which models the expected \LambdaCDM Integrated Sachs-Wolfe (ISW) effect in the Jubilee simulation. The simulation volume is (6 Gpc/h)^3, allowing power on very large-scales to be incorporated into the calculation. Haloes are resolved down to a mass of 1.5x10^12 M_sun/h, which allows us to derive a catalogue of mock Luminous Red Galaxies (LRGs) for cross-correlation analysis with the ISW signal. We find the ISW effect observed on a projected sky to grow stronger at late times with the evolution of the ISW power spectrum matching expectations from linear theory. Maps of the gravitational lensing effect, including the convergence and deflection fields, are calculated using the same potential as for the ISW. We calculate the redshift dependence of the ISW-LRG cross-correlation signal for a full sky survey with no noise considerations. For l < 30, the signal is strongest for lower redshift bins (z ~ 0.2 to 0.5), whereas for l > 30 the signal is best observed with surveys covering z ~ 0.6-1.0.

Clowes et al. (2013) have recently reported the discovery of a Large Quasar Group (LQG), dubbed the Huge-LQG, at redshift z~1.3 in the DR7 quasar catalogue of the Sloan Digital Sky Survey. On the basis of its characteristic size ~500 Mpc and longest dimension >1 Gpc, it is claimed that this structure is incompatible with large-scale homogeneity and the cosmological principle. If true, this would represent a serious challenge to the standard cosmological model. However, the homogeneity scale is an average property which is not necessarily affected by the discovery of a single large structure. I clarify this point and provide the first fractal dimension analysis of the DR7 quasar catalogue to demonstrate that it is in fact homogeneous above scales of at most 130 Mpc/h, which is much less than the upper limit for \Lambda CDM. In addition, I show that the algorithm used to identify the Huge-LQG regularly finds even larger clusters of points, extending over Gpc scales, in explicitly homogeneous simulations of a Poisson point process with the same density as the quasar catalogue. This provides a simple null test to be applied to any cluster thus found in a real catalogue, and suggests that the interpretation of LQGs as `structures' is misleading.

A detection of the stacked integrated Sachs-Wolfe (ISW) signal in the CMB of rare superstructures identified in the SDSS Luminous Red Galaxy catalogue has been reported at very high statistical significance. The magnitude of the observed signal has previously been argued to be more than 3 standard deviations larger than the theoretical \Lambda CDM expectation. However, this calculation was made in the linear approximation, and relied on assumptions that may potentially have caused the \Lambda CDM expectation to be underestimated. Here we update the theoretical model calculation and compare it with an analysis of ISW maps obtained from N-body simulations of a \Lambda CDM universe. The differences between model predictions and the map analyses are found to be small and cannot explain the discrepancy with observation, which remains at >3 s.d. significance. We discuss the cosmological significance of this anomaly and speculate on the potential of alternative models to explain it.

A crucial diagnostic of the \Lambda CDM cosmological model is the integrated Sachs-Wolfe (ISW) effect of large-scale structure on the cosmic microwave background (CMB). The ISW imprint of superstructures of size \sim100\;h^{-1} Mpc at redshift $z\sim0.5$ has been detected with $>4\sigma$ significance, however it has been noted that the signal is much larger than expected. We revisit the calculation using linear theory predictions in \Lambda CDM cosmology for the number density of superstructures and their radial density profile, and take possible selection effects into account. While our expected signal is larger than previous estimates, it is still inconsistent by $>3\sigma$ with the observation. If the observed signal is indeed due to the ISW effect then huge, extremely underdense voids are far more common in the observed universe than predicted by \Lambda CDM.

It is well known that large Hubble-induced supergravity corrections to the inflaton field can ruin the flatness of the potential, thus creating a tension between slow-roll inflation and supergravity. In this paper we show that it is possible to obtain a {\it cosmologically flat} direction, embedded within the minimal supersymmetric standard model, despite very large super-Hubble corrections. As an illustration, we show that a flat direction which is lifted by an $n=6$ operator matches the current cosmic microwave background data for a wide range of the Hubble parameter, $10^{5} {\rm GeV}\lesssim H_{inf}\lesssim 10^{8.5}$ GeV. Our approach can be applied to any $F$-term inflationary model.

We provide a model in which both the inflaton and the curvaton are obtained from within the minimal supersymmetric Standard Model, with known gauge and Yukawa interactions. Since now both the inflaton and curvaton fields are successfully embedded within the same sector, their decay products thermalize very quickly before the electroweak scale. This results in two important features of the model: firstly, there will be no residual isocurvature perturbations, and secondly, observable non-Gaussianities can be generated with the non-Gaussianity parameter $f_\mathrm{NL}\sim {\cal O}(5-1000)$ being determined solely by the combination of weak-scale physics and the Standard Model Yukawas.

The detection of primordial gravitational waves, or tensor perturbations, would be regarded as compelling evidence for inflation. The canonical measure of this is the ratio of tensor to scalar perturbations, r. For single-field slow-roll models of inflation with small field excursions, the Lyth bound dictates that if the evolution of the slow-roll parameter epsilon is monotonic, the tensor-to-scalar ratio must be below observationally detectable levels. We describe how non-monotonic evolution of epsilon can evade the Lyth bound and generate observationally large r, even with small field excursions. This has consequences for the scalar power spectrum as it necessarily predicts an enhancement in the spectrum at very small scales and significant scale-dependent running at CMB scales. This effect has not been appropriately accounted for in previous analyses. We describe a mechanism that will generically produce the required behaviour in epsilon and give an example of this mechanism arising in a well-motivated small-field model. This model can produce r\geq0.05 while satisfying all current observational constraints.

We consider observational constraints and fine-tuning issues in a renormalizable model of inflection point inflation, with two independent parameters. We derive constraints on the parameter space of this model arising from the WMAP 7-year power spectrum. It has previously been shown that it is possible to successfully embed this potential in the MSSM. Unfortunately, to do this requires severe fine-tuning. We address this issue by introducing a hybrid field to dynamically uplift the potential with a subsequent smooth phase transition to end inflation at the necessary point. Large parameter regions exist where this drastically reduces the fine-tuning required without ruining the viability of the model. A side effect of this mechanism is that it increases the width of the slow-roll region of the potential, thus also alleviating the problem of the fine-tuning of initial conditions. The MSSM embedding we study has been previously shown to be able to explain the smallness of the neutrino masses. The hybrid transition does not spoil this feature as there exist parameter regions where the fine-tuning parameter is as large as $10^{-1}$ and the neutrino masses remain small.

In the standard cosmological model, the dimming of distant Type Ia supernovae is explained by invoking the existence of repulsive `dark energy' which is causing the Hubble expansion to accelerate. However this may be an artifact of interpreting the data in an (oversimplified) homogeneous model universe. In the simplest inhomogeneous model which fits the SNe Ia Hubble diagram without dark energy, we are located close to the centre of a void modelled by a Lema\'itre-Tolman-Bondi metric. It has been claimed that such models cannot fit the CMB and other cosmological data. This is however based on the assumption of a scale-free spectrum for the primordial density perturbation. An alternative physically motivated form for the spectrum enables a good fit to both SNe Ia (Constitution/Union2) and CMB (WMAP 7-yr) data, and to the locally measured Hubble parameter. Constraints from baryon acoustic oscillations and primordial nucleosynthesis are also satisfied.