Scientific analyses often rely on slow, but accurate forward models for observable data conditioned on known model parameters. While various emulation schemes exist to approximate these slow calculations, these approaches are only safe if the approximations are well understood and controlled. This workshop submission reviews and updates a previously published method, which has been used in cosmological simulations, to (1) train an emulator while simultaneously estimating posterior probabilities with MCMC and (2) explicitly propagate the emulation error into errors on the posterior probabilities for model parameters. We demonstrate how these techniques can be applied to quickly estimate posterior distributions for parameters of the $\Lambda$CDM cosmology model, while also gauging the robustness of the emulator approximation.

We present an efficient implementation of Wiener filtering of real-space linear field and optimal quadratic estimator of its power spectrum Band-powers. We first recast the field reconstruction into an optimization problem, which we solve using quasi-Newton optimization. We then recast the power spectrum estimation into the field marginalization problem, from which we obtain an expression that depends on the field reconstruction solution and a determinant term. We develop a novel simulation based method for the latter. We extend the simulations formalism to provide the covariance matrix for the power spectrum. We develop a flexible framework that can be used on a variety of cosmological fields and present results for a variety of test cases, using simulated examples of projected density fields, projected shear maps from galaxy lensing, and observed Cosmic Microwave Background (CMB) temperature anisotropies, with a wide range of map incompleteness and variable noise. For smaller cases where direct numerical inversion is possible, we show that our solution matches that created by direct Wiener Filtering at a fraction of the overall computation cost. Even more significant reduction of computational is achieved by this implementation of optimal quadratic estimator due to the fast evaluation of the Hessian matrix. This technique allows for accurate map and power spectrum reconstruction with complex masks and nontrivial noise properties.

One of the main unsolved problems of cosmology is how to maximize the extraction of information from nonlinear data. If the data are nonlinear the usual approach is to employ a sequence of statistics (N-point statistics, counting statistics of clusters, density peaks or voids etc.), along with the corresponding covariance matrices. However, this approach is computationally prohibitive and has not been shown to be exhaustive in terms of information content. Here we instead develop a Bayesian approach, expanding the likelihood around the maximum posterior of linear modes, which we solve for using optimization methods. By integrating out the modes using perturbative expansion of the likelihood we construct an initial power spectrum estimator, which for a fixed forward model contains all the cosmological information if the initial modes are gaussian distributed. We develop a method to construct the window and covariance matrix such that the estimator is explicitly unbiased and nearly optimal. We then generalize the method to include the forward model parameters, including cosmological and nuisance parameters, and primordial non-gaussianity. We apply the method in the simplified context of nonlinear structure formation, using either simplified 2-LPT dynamics or N-body simulations as the nonlinear mapping between linear and nonlinear density, and 2-LPT dynamics in the optimization steps used to reconstruct the initial density modes. We demonstrate that the method gives an unbiased estimator of the initial power spectrum, providing among other a near optimal reconstruction of linear baryonic acoustic oscillations.

Cosmological inflation generates primordial density perturbations on all scales, including those far too small to contribute to the cosmic microwave background. At these scales, isolated ultracompact minihalos of dark matter can form well before standard structure formation, if the perturbations have sufficient amplitude. Minihalos affect pulsar timing data and are potentially bright sources of gamma rays. The resulting constraints significantly extend the observable window of inflation in the presence of cold dark matter, coupling two of the key problems in modern cosmology.

We test the sensitivity of neutrino parameter constraints from combinations of CMB and LSS data sets to the assumed form of the primordial power spectrum (PPS) using Bayesian model selection. Significantly, none of the tested combinations, including recent high-precision local measurements of $\mathrm{H}_0$ and cluster abundances, indicate a signal for massive neutrinos or extra relativistic degrees of freedom. For PPS models with a large, but fixed number of degrees of freedom, neutrino parameter constraints do not change significantly if the location of any features in the PPS are allowed to vary, although neutrino constraints are more sensitive to PPS features if they are known a priori to exist at fixed intervals in $\log k$. Although there is no support for a non-standard neutrino sector from constraints on both neutrino mass and relativistic energy density, we see surprisingly strong evidence for features in the PPS when it is constrained with data from Planck 2015, SZ cluster counts, and recent high-precision local measurements of $\mathrm{H}_0$. Conversely combining Planck with matter power spectrum and BAO measurements yields a much weaker constraint. Given that this result is sensitive to the choice of data this tension between SZ cluster counts, Planck and $\mathrm{H}_0$ measurements is likely an indication of unmodeled systematic bias that mimics PPS features, rather than new physics in the PPS or neutrino sector.

Cosmological analyses can be accelerated by approximating slow calculations using a training set, which is either precomputed or generated dynamically. However, this approach is only safe if the approximations are well understood and controlled. This paper surveys issues associated with the use of machine-learning based emulation strategies for accelerating cosmological parameter estimation. We describe a learn-as-you-go algorithm that is implemented in the Cosmo++ code and (1) trains the emulator while simultaneously estimating posterior probabilities; (2) identifies unreliable estimates, computing the exact numerical likelihoods if necessary; and (3) progressively learns and updates the error model as the calculation progresses. We explicitly describe and model the emulation error and show how this can be propagated into the posterior probabilities. We apply these techniques to the Planck likelihood and the calculation of $\Lambda$CDM posterior probabilities. The computation is significantly accelerated without a pre-defined training set and uncertainties in the posterior probabilities are subdominant to statistical fluctuations. We have obtained a speedup factor of $6.5$ for Metropolis-Hastings and $3.5$ for nested sampling. Finally, we discuss the general requirements for a credible error model and show how to update them on-the-fly.

A short inflationary phase may not erase all traces of the primordial universe. Associated observables include both spatial curvature and "anomalies" in the microwave background or large scale structure. The present curvature $\Omega_{K,0}$ reflects the initial curvature, $\Omega_{K,\mathrm{start}}$, and the angular size of anomalies depends on $k_\mathrm{start}$, the comoving horizon size at the onset of inflation. We estimate posteriors for $\Omega_{K,\mathrm{start}}$ and $k_\mathrm{start}$ using current data and simulations, and show that if either quantity is measured to have a non-zero value, both are likely to be observable. Mappings from $\Omega_{K,\mathrm{start}}$ and $k_\mathrm{start}$ to present-day observables depend strongly on the primordial equation of state; $\Omega_{K,0}$ spans ten orders of magnitude for a given $\Omega_{K,\mathrm{start}}$ while a simple and general relationship connects $\Omega_{K,0}$ and $k_\mathrm{start}$. We show that current bounds on $\Omega_{K,0}$ imply that if $k_\mathrm{start}$ is measurable, the curvature was already small when inflation began. Finally, since the energy density changes slowly during inflation, primordial gravitational wave constraints require that a short inflationary phase is preceded by a nontrivial pre-inflationary phase with critical implications for the expected value of $\Omega_{K,\mathrm{start}}$.

Using the temperature data from Planck we search for departures from a power-law primordial power spectrum, employing Bayesian model-selection and posterior probabilities. We parametrize the spectrum with $n$ knots located at arbitrary values of $\log{k}$, with both linear and cubic splines. This formulation recovers both slow modulations and sharp transitions in the primordial spectrum. The power spectrum is well-fit by a featureless, power-law at wavenumbers $k>10^{-3} \, \mathrm{Mpc}^{-1}$. A modulated primordial spectrum yields a better fit relative to $\Lambda$CDM at large scales, but there is no strong evidence for a departure from a power-law spectrum. Moreover, using simulated maps we show that a local feature at $k \sim 10^{-3} \, \mathrm{Mpc}^{-1}$ can mimic the suppression of large-scale power. With multi-knot spectra we see only small changes in the posterior distributions for the other free parameters in the standard $\Lambda$CDM universe. Lastly, we investigate whether the hemispherical power asymmetry is explained by independent features in the primordial power spectrum in each ecliptic hemisphere, but find no significant differences between them.

An inflationary gravitational wave background consistent with BICEP2 is difficult to reconcile with a simple power-law spectrum of primordial scalar perturbations. Tensor modes contribute to the temperature anisotropies at multipoles with $l\lesssim 100$, and this effect --- together with a prior on the form of the scalar perturbations --- was the source of previous bounds on the tensor-to-scalar ratio. We compute Bayesian evidence for combined fits to BICEP2 and Planck for three nontrivial primordial spectra: a) a running spectral index, b) a cutoff at fixed wavenumber, and c) a spectrum described by a linear spline with a single internal knot. We find no evidence for a cutoff, weak evidence for a running index, and significant evidence for a "broken" spectrum. Taken at face-value, the BICEP2 results require two new inflationary parameters in order to describe both the broken scale invariance in the perturbation spectrum and the observed tensor-to-scalar ratio. Alternatively, this tension may be resolved by additional data and more detailed analyses.

This paper introduces a new publicly available numerical library for cosmology, Cosmo++. The library has been designed using object-oriented programming techniques, and fully implemented in C++. Cosmo++ introduces a unified interface for using most of the frequently used numerical methods in cosmology. Most of the features are implemented in Cosmo++ itself, while a part of the functionality is implemented by linking to other publicly available libraries. The most important features of the library are Cosmic Microwave Background anisotropies power spectrum and transfer function calculations, likelihood calculations, parameter space sampling tools, sky map simulations, and mask apodization. Cosmo++ also includes a few mathematical tools that are frequently used in numerical research in cosmology and beyond. A few simple examples are included in Cosmo++ to help the user understand the key features. The library has been fully tested, and we describe some of the important tests in this paper. Cosmo++ is publicly available at http://cosmo.grigoraslanyan.com

Most treatments of large scale anomalies in the microwave sky are a posteriori, with unquantified look-elsewhere effects. We contrast these with physical models of specific inhomogeneities in the early universe which then generate apparent anomalies. Physical models predict correlations between candidate anomalies, as well as the corresponding signals in polarization and large scale structure, reducing the impact of cosmic variance. We compute the apparent spatial curvature associated with large-scale inhomogeneities and show that it is typically small, allowing for a self-consistent analysis. Focussing on a single large plane wave inhomogeneity, we show this can contribute to low-l mode alignment and odd-even asymmetry in the power spectra and the best fit inhomogeneity accounts for a significant part of the claimed odd-even asymmetry. We argue that this approach can be generalized to provide a more quantitative assessment of potential large scale anomalies in the universe.

We analyze seven year and nine year WMAP temperature maps for signatures of three finite flat topologies M_0=T^3, M_1=T^2 x R^1, and M_2=S^1 x R^2. We use Monte-Carlo simulations with the Feldman-Cousins method to obtain confidence intervals for the size of the topologies considered. We analyze the V, W, and Q frequency bands along with the ILC map and find no significant difference in the results. The 95.5% confidence level lower bound on the size of the topology is 1.5L_0 for M_0, 1.4L_0 for M_1, and 1.1L_0 for M_2, where L_0 is the radius of the last scattering surface. Our results agree very well with the recently released results from the Planck temperature data. We show that the likelihood function is not Gaussian in the size, and therefore simulations are important for obtaining accurate bounds on the size. We then introduce the formalism for including polarization data in the analysis. The improvement that we find from WMAP polarization maps is small because of the high level of instrumental noise, but our forecast for Planck maps shows a much better improvement on the lower bound for L. For the M_0 topology we expect an improvement on the lower bound of L from 1.7L_0 to 1.9L_0 at 95.5% confidence level. Using both polarization and temperature data is important because it tests the hypothesis that deviations in the TT spectrum at small l originate in the primordial perturbation spectrum.

We place limits on semiclassical fluctuations that might be present in the primordial perturbation spectrum. These can arise if some signatures of pre-inflationary features survive the expansion, or could be created by whatever mechanism ends inflation. We study two possible models for such remnant fluctuations, both of which break the isotropy of CMB on large scales. We first consider a semiclassical fluctuation in one Fourier mode of primordial perturbations. The second scenario we analyze is a semiclassical Gaussian bump somewhere in space. These models are tested with the seven-year WMAP data using a Markov Chain Monte Carlo Bayesian analysis, and we place limits on these fluctuations. The upper bound for the amplitude of a fluctuation in a single Fourier mode is a<=10^(-4), while for the Gaussian bump a<=10^(-3).

We cross correlate the gravitational lensing map extracted from cosmic microwave background measurements by the Wilkinson Microwave Anisotropy Probe (WMAP) with the radio galaxy distribution from the NRAO VLA Sky Survey (NVSS) by using a quadratic estimator technique. We use the full covariance matrix to filter the data, and calculate the cross-power spectra for the lensing-galaxy correlation. We explore the impact of changing the values of cosmological parameters on the lensing reconstruction, and obtain statistical detection significances at $>3\sigma$. The results of all cross correlations pass the curl null test as well as a complementary diagnostic test using the NVSS data in equatorial coordinates. We forecast the potential for Planck and NVSS to constrain the lensing-galaxy cross correlation as well as the galaxy bias. The lensing-galaxy cross-power spectra are found to be Gaussian distributed.

We study the possibility that the universe has compact topologies T^3, T^2 x R^1, or S^1 x R^2 using the seven-year WMAP data. The maximum likelihood 95% confidence intervals for the size L of the compact direction are 1.7 < L/L_0 < 2.1, 1.8 < L/L_0 < 2.0, 1.2 < L/L_0 < 2.1 for the three cases, respectively, where L_0=14.4 Gpc is the distance to the last scattering surface. An infinite universe is compatible with the data at 4.3 sigma. We find using a Bayesian analysis that the most probable universe has topology T^2 x R^1, with L/L_0=1.9.