Landscape cosmology posits the existence of a convoluted, multidimensional, scalar potential -- the "landscape" -- with vast numbers of metastable minima. Random matrices and random functions in many dimensions provide toy models of the landscape, allowing the exploration of conceptual issues associated with these scenarios. We compute the relative number and slopes of minima as a function of the vacuum energy $\Lambda$ in an $N$-dimensional Gaussian random potential, quantifying the associated probability density, $p(\Lambda)$. After normalisations $p(\Lambda)$ depends only on the dimensionality $N$ and a single free parameter $\gamma$, which is related to the power spectrum of the random function. For a Gaussian landscape with a Gaussian power spectrum, the fraction of positive minima shrinks super-exponentially with $N$; at $N=100$, $p(\Lambda>0) \approx 10^{-780}$. Likewise, typical eigenvalues of the Hessian matrices reveal that the flattest approaches to typical minima grow flatter with $N$, while the ratio of the slopes of the two flattest directions grows with $N$. We discuss the implications of these results for both swampland and conventional anthropic constraints on landscape cosmologies. In particular, for parameter values when positive minima are extremely rare, the flattest approaches to minima where $\Lambda \approx 0$ are much flatter than for typical minima, increasingly the viability of quintessence solutions.

In simple inflationary cosmological scenarios the near-exponential growth can be followed by a long period in which the Universe is dominated by the oscillating inflaton condensate. The condensate is initially almost homogeneous, but perturbations grow gravitationally, eventually fragmenting the condensate if it is not disrupted more quickly by resonance or prompt reheating. We show that the gravitational fragmentation of the condensate is well-described by the Schr\"odinger-Poisson equations and use numerical solutions to show that large overdensities form quickly after the onset of nonlinearity. This is the first exploration of this phase of nonlinear dynamics in the very early universe, which can affect the detailed form of the inflationary power spectrum and the dark matter fraction when the dark sector is directly coupled to the inflaton.

We recently showed that postulated ultracompact minihalos with a steep density profile do not form in realistic simulations with enhanced initial perturbations. In this paper we assume that a small fraction of the dark matter consists of primordial black holes (PBHs) and simulate the formation of structures around them. We find that in this scenario halos with steep density profiles do form, consistent with theoretical predictions. If the rest of the dark matter consists of weakly interacting massive particles (WIMPs), we also show that WIMPs in the dense innermost part of these halos would annihilate and produce a detectable gamma-ray signal. The non-detection of this signal implies that PBHs make up at most one billionth of the dark matter, provided that their mass is greater than one millionth of the mass of the Sun. Similarly, a detection of PBHs would imply that the remaining dark matter could not be WIMPs.

PyUltraLight simulates the dynamics of ultralight dark matter in a non-expanding background. PyUltraLight can describe the evolution of several interacting ultralight dark matter halos or one or more halos orbiting a central, fixed Newtonian potential, the latter scenario corresponding to dwarf galaxies orbiting a massive central galaxy. We verify PyUltraLight by showing that it reproduces qualitative dynamical features of previously published simulations and demonstrate that it has excellent energy-conservation properties. PyUltraLight is implemented in a Python-based Jupyter notebook, solving the Schrodinger-Poisson equation governing ultralight scalar field dark matter dynamics in the non-relativistic regime using a symmetrised split-step pseudospectral algorithm. The notebook interface makes it simple to specify simulation parameters and visualise the resulting output but performance-critical routines are managed via calls to computationally efficient compiled libraries. PyUltraLight runs on standard desktop hardware with support for shared memory mutlithreading and is available on GitHub.

We perform three-dimensional simulations of structure formation in the early Universe, when boosting the primordial power spectrum on approximately kpc scales. We demonstrate that our simulations are capable of producing power-law profiles close to the steep $\rho\propto r^{-9/4}$ halo profiles that are commonly assumed to be a good approximation to ultracompact minihalos (UCMHs). However, we show that for more realistic initial conditions in which halos are neither perfectly symmetric nor isolated, the steep power-law profile is disrupted and we find that the Navarro-Frenk-White profile is a better fit to most halos. In the presence of background fluctuations even extreme, nearly spherical initial conditions do not remain exceptional. Nonetheless, boosting the amplitude of initial fluctuations causes all structures to form earlier and thus at larger densities. With sufficiently large amplitude of fluctuations we find that values for the concentration of typical halos in our simulations can become very large. However, despite the signal coming from dark matter annihilation inside the cores of these halos being enhanced, it is still orders-of-magnitude smaller compared to the usually assumed UCMH profile. The upper bound on the primordial power spectrum from the non-observation of UCMHs should therefore be re-evaluated.

The use of Eulerian 'standard perturbation theory' to describe mass assembly in the early universe has traditionally been limited to modes with k $\lesssim$ 0.1 h/Mpc at z = 0. At larger k the SPT power spectrum deviates from measurements made using N-body simulations. Recently, there has been progress in extending the reach of perturbation theory to larger k using ideas borrowed from effective field theory. We revisit the computation of the redshift-space matter power spectrum within this framework, including for the first time the full one-loop time dependence. We use a resummation scheme proposed by Vlah et al. to account for damping of baryonic acoustic oscillations due to large-scale random motions and show that this has a significant effect on the multipole power spectra. We renormalize by comparison to a suite of custom N-body simulations matching the MultiDark MDR1 cosmology. At z = 0 and for scales k $\lesssim$ 0.4 h/Mpc we find that the EFT furnishes a description of the real-space power spectrum up to $\sim$2%, for the $\ell$ = 0 mode up to $\sim$5%, and for the $\ell$ = 2, 4 modes up to $\sim$25%. We argue that, in the MDR1 cosmology, positivity of the $\ell$ = 0 mode gives a firm upper limit of k $\approx$ 0.75 h/Mpc for the validity of the one-loop EFT prediction in redshift space using only the lowest-order counterterm. We show that replacing the one-loop growth factors by their Einstein-de Sitter counterparts is a good approximation for the $\ell$ = 0 mode, but can induce deviations as large as 2% for the $\ell$ = 2, 4 modes. An accompanying software bundle, distributed under open source licenses, includes Mathematica notebooks describing the calculation, together with parallel pipelines capable of computing both the necessary one-loop SPT integrals and the effective field theory counterterms.

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.

We consider methods with which to answer the question "is any observed galaxy cluster too unusual for Lambda-CDM?" After emphasising that many previous attempts to answer this question will overestimate the confidence level at which Lambda-CDM can be ruled out, we outline a consistent approach to these rare clusters, which allows the question to be answered. We define three statistical measures, each of which are sensitive to changes in cluster populations arising from different modifications to the cosmological model. We also use these properties to define the 'equivalent mass at redshift zero' for a cluster --- the mass of an equally unusual cluster today. This quantity is independent of the observational survey in which the cluster was found, which makes it an ideal proxy for ranking the relative unusualness of clusters detected by different surveys. These methods are then used on a comprehensive sample of observed galaxy clusters and we confirm that all are less than 2-sigma deviations from the Lambda-CDM expectation. Whereas we have only applied our method to galaxy clusters, it is applicable to any isolated, collapsed, halo. As motivation for future surveys, we also calculate where in the mass redshift plane the rarest halo is most likely to be found, giving information as to which objects might be the most fruitful in the search for new physics.

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.

Within a cosmological context, we study the behaviour of collisionless particles in the weak field approximation to General Relativity, allowing for large gradients of the fields and relativistic velocities for the particles. We consider a spherically symmetric setup such that high resolution simulations are possible with minimal computational resources. We test our formalism by comparing it to two exact solutions: the Schwarzschild solution and the Lema\^itre-Tolman-Bondi model. In order to make the comparison we consider redshifts and lensing angles of photons passing through the simulation. These are both observable quantities and hence are gauge independent. We demonstrate that our scheme is more accurate than a Newtonian scheme, correctly reproducing the leading-order post-Newtonian correction. In addition, our setup is able to handle shell-crossings, which is not possible within a fluid model. Furthermore, by introducing angular momentum, we find configurations corresponding to bound objects which may prove useful for numerical studies of the effects of modified gravity, dynamical dark energy models or even compact bound objects within General Relativity.

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 investigate the hemispherical power asymmetry in the cosmic microwave background on small angular scales. We find an anomalously high asymmetry in the multipole range l=601-2048, with a naive statistical significance of 6.5 sigma. However, we show that this extreme anomaly is simply a coincidence of three other effects, relativistic power modulation, edge effects from the mask applied, and inter-scale correlations. After correcting for all of these effects, the significance level drops to ~1 sigma, i.e., there is no anomalous intrinsic asymmetry in the small angular scales. Using this null result, we derive a constraint on a potential dipolar modulation amplitude, A(k)<0.0045 on the ~10 Mpc-scale, at 95% C.L. This new constraint must be satisfied by any theoretical model attempting to explain the hemispherical asymmetry at large angular scales.

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.

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.

We address the issue of cosmological backreaction from non-linear structure formation by constructing an approximation for the time evolved metric of a dust dominated universe based on a gradient expansion. Our metric begins as a perturbation of a flat Friedmann-Robertson-Walker state described by a nearly scale invariant, Gaussian, power-law distribution, and evolves in time until non-linear structures have formed. After describing and attempting to control for certain complications in the implementation of this approach, this metric then forms a working model of the universe. We numerically calculate the evolution of the average scale factor in this model and hence the backreaction. We argue that, despite its limitations, this model is more realistic than previous models that have confronted the issue of backreaction. We find that the \emph{instantaneous} effects of backreaction in this model could be as large as $\sim10%$ of the background. This suggests that a proper understanding of the \emph{cumulative} effects of backreaction could be crucial for precision cosmology and any future exploration of the dark sector.

We critically investigate current statistical tests applied to high redshift clusters of galaxies in order to test the standard cosmological model and describe their range of validity. We carefully compare a sample of high-redshift, massive, galaxy clusters with realistic Poisson sample simulations of the theoretical mass function, which include the effect of Eddington bias. We compare the observations and simulations using the following statistical tests: the distributions of ensemble and individual existence probabilities (in the >M,>z sense), the redshift distributions, and the 2d Kolmogorov-Smirnov test. Using seemingly rare clusters from Hoyle et al. (2011), and Jee et al. (2011) and assuming the same survey geometry as in Jee et al. (2011, which is less conservative than Hoyle et al. 2011), we find that the (>M,>z) existence probabilities of all clusters are fully consistent with LCDM. However assuming the same survey geometry, we use the 2d K-S test probability to show that the observed clusters are not consistent with being the least probable clusters from simulations at >95% confidence, and are also not consistent with being a random selection of clusters, which may be caused by the non-trivial selection function and survey geometry. Tension can be removed if we examine only a X-ray selected sub sample, with simulations performed assuming a modified survey geometry.

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.

There are observations of 15 high-redshift massive galaxy clusters, which have an extremely small probability with a purely Gaussian initial curvature perturbation. Here we revisit the estimation of the contribution of non-Gaussianities to the cluster mass function and point out serious problems that have resulted in the application of the mass function out of the range of its validity. We remedy the situation and show that the values of f_NL previously claimed to completely reconcile (i.e. at ~100% confidence) the existence of the clusters with LambdaCDM are unphysically small. However, for WMAP cosmology and at 95% confidence, we arrive at the limit f_NL>411, which is similar to previous estimates. We also explore the possibility of a large g_NL as the reason for the observed excess of the massive galaxy clusters. This scenario, g_NL>2*10^6, appears to be in more agreement with CMB and LSS limits for the non-Gaussianity parameters and could also provide an explanation for the overabundance of large voids in the early universe.

Number counts of massive high-redshift clusters provide a window to study primordial non-Gaussianity. The current quality of data, however, forces the statistical analysis to probe a region of parameter space -- the extreme tail of the mass function -- which is neither accessible in any of the currently available theoretical prescriptions for calculating the mass function, nor calibrated in N-body simulations. In this work we present a new analytical prescription for calculating a "resummed" non-Gaussian halo mass function, which is constructed to remain stable in the extreme tail. We show that the prescription works well in the parameter regime that has been currently explored in simulations. We then use Fisher matrix techniques to compare our prescription with an extrapolated fit to N-body simulations, which has recently been used to obtain constraints from data collected by the South Pole Telecope. We show that for the current data, both prescriptions would lead to statistically consistent constraints. As the data improve, however, there is a possibility of introducing a statistically significant bias in the constraints due to the choice of prescription, especially if non-Gaussianity is scale dependent and becomes relatively large on cluster scales. It would then be necessary to test the accuracy of the prescriptions in N-body simulations that can probe clusters with high masses and redshifts in the presence of large non-Gaussianity.

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.

I show that the most common method of quantifying the likelihood that an extreme galaxy cluster could exist is biased and can result in false claims of tension with LambdaCDM. This common method uses the probability that at least one cluster could exist above the mass and redshift of an observed cluster. I demonstrate the existence of the bias using sample cluster populations, describe its origin and explain how to remove it. I then suggest potentially more suitable and unbiased measures of the rareness of individual clusters. Each different measure will be most sensitive to different possible types of new physics. I show how to generalise these measures to quantify the total `rareness' of a set of clusters. It is seen that, when mass uncertainties are marginalised over, there is no tension between the standard LambdaCDM cosmological model and the existence of any observed set of clusters. As a case study, I apply these rareness measures to sample cluster populations generated using primordial density perturbations with a non-Gaussian spectrum.

Multiple inflation is a model based on N=1 supergravity wherein there are sudden changes in the mass of the inflaton because it couples to 'flat direction' scalar fields which undergo symmetry breaking phase transitions as the universe cools. The resulting brief violations of slow-roll evolution generate a non-gaussian signal which we find to be oscillatory and yielding f_NL ~ 5-20. This is potentially detectable by e.g. Planck but would require new bispectrum estimators to do so. We also derive a model-independent result relating the period of oscillations of a phase transition during inflation to the period of oscillations in the primordial curvature perturbation generated by the inflaton.

The form of the inflationary potential is severely restricted if one requires that it be natural in the technical sense, i.e. terms of unrelated origin are not required to be correlated. We determine the constraints on observables that are implied in such natural inflationary models, in particular on $r$, the ratio of tensor to scalar perturbations. We find that the naturalness constraint does not require $r$ to be lare enough to be detectable by the forthcoming searches for B-mode polarisation in CMB maps. We show also that the value of $r$ is a sensitive discriminator between inflationary models.