The search for primordial gravitational waves in the Cosmic Microwave Background (CMB) will soon be limited by our ability to remove the lensing contamination to $B$-mode polarization. The often-used quadratic estimator for lensing is known to be suboptimal for surveys that are currently operating and will continue to become less and less efficient as instrumental noise decreases. While foregrounds can in principle be mitigated by observing in more frequency bands, progress in delensing hinges entirely on algorithmic advances. We demonstrate here a new inference method that solves this problem by sampling the exact Bayesian posterior of any desired cosmological parameters, of the gravitational lensing potential, and of the delensed CMB maps, given lensed temperature and polarization data. We validate the method using simulated CMB data with non-white noise and masking on up to 650\,deg$^2$ patches of sky. A unique strength of this approach is the ability to jointly estimate cosmological parameters which control both the primordial CMB and the lensing potential, which we demonstrate here for the first time by sampling both the tensor-to-scalar ratio, $r$, and the amplitude of the lensing potential, $A_\phi$. The method allows us to perform the most precise check to-date of several important approximations underlying CMB-S4 $r$ forecasting, and we confirm these yield the correct expected uncertainty on $r$ to better than 10%.

One of the major targets for next-generation cosmic microwave background (CMB) experiments is the detection of the primordial B-mode signal. Planning is under way for Stage-IV experiments that are projected to have instrumental noise small enough to make lensing and foregrounds the dominant source of uncertainty for estimating the tensor-to-scalar ratio $r$ from polarization maps. This makes delensing a crucial part of future CMB polarization science. In this paper we present a likelihood method for estimating the tensor-to-scalar ratio $r$ from CMB polarization observations, which combines the benefits of a full-scale likelihood approach with the tractability of the quadratic delensing technique. This method is a pixel space, all order likelihood analysis of the quadratic delensed B modes, and it essentially builds upon the quadratic delenser by taking into account all order lensing and pixel space anomalies. Its tractability relies on a crucial factorization of the pixel space covariance matrix of the polarization observations which allows one to compute the full Gaussian approximate likelihood profile, as a function of $r$, at the same computational cost of a single likelihood evaluation.

We develop the first algorithm able to jointly compute the maximum {\it a posteriori} estimate of the Cosmic Microwave Background (CMB) temperature and polarization fields, the gravitational potential by which they are lensed, and cosmological parameters such as the tensor-to-scalar ratio, $r$. This is an important step towards sampling from the joint posterior probability function of these quantities, which, assuming Gaussianity of the CMB fields and lensing potential, contains all available cosmological information and would yield theoretically optimal constraints. Attaining such optimal constraints will be crucial for next-generation CMB surveys like CMB-S4, where limits on $r$ could be improved by factors of a few over currently used sub-optimal quadratic estimators. The maximization procedure described here depends on a newly developed lensing algorithm, which we term \textsc{LenseFlow}, and which lenses a map by solving a system of ordinary differential equations. This description has conceptual advantages, such as allowing us to give a simple non-perturbative proof that the lensing determinant is equal to unity in the weak-lensing regime. The algorithm itself maintains this property even on pixelized maps, which is crucial for our purposes and unique to \textsc{LenseFlow} as compared to other lensing algorithms we have tested. It also has other useful properties such as that it can be trivially inverted (i.e. delensing) for the same computational cost as the forward operation, and can be used to compute lensing adjoint, Jacobian, and Hessian operators. We test and validate the maximization procedure on flat-sky simulations covering up to 600\,deg$^2$ with non-uniform noise and masking.

As modern precision cosmological measurements continue to show agreement with the broad features of the standard $\Lambda$-Cold Dark Matter ($\Lambda$CDM) cosmological model, we are increasingly motivated to look for small departures from the standard model's predictions which might not be detected with standard approaches. While searches for extensions and modifications of $\Lambda$CDM have to date turned up no convincing evidence of beyond-the-standard-model cosmology, the list of models compared against $\Lambda$CDM is by no means complete and is often governed by readily-coded modifications to standard Boltzmann codes. Also, standard goodness-of-fit methods such as a naive $\chi^2$ test fail to put strong pressure on the null $\Lambda$CDM hypothesis, since modern datasets have orders of magnitudes more degrees of freedom than $\Lambda$CDM. Here we present a method of tuning goodness-of-fit tests to detect potential sub-dominant extra-$\Lambda$CDM signals present in the data through compressing observations in a way that maximizes extra-$\Lambda$CDM signal variation over noise and $\Lambda$CDM variation. This method, based on a Karhunen-Lo\`{e}ve transformation of the data, is tuned to be maximally sensitive to particular types of variations characteristic of the tuning model; but, unlike direct model comparison, the test is also sensitive to features that only partially mimic the tuning model. As an example of its use, we apply this method in the context of a nonstandard primordial power spectrum compared against the $2015$ $Planck$ CMB temperature and polarization power spectrum. We find weak evidence of extra-$\Lambda$CDM physics, conceivably due to known systematics in the 2015 Planck polarization release.

The Planck satellite, along with several ground based telescopes, have mapped the cosmic microwave background (CMB) at sufficient resolution and signal-to-noise so as to allow a detection of the subtle distortions due to the gravitational influence of the intervening matter distribution. A natural modeling approach is to write a Bayesian hierarchical model for the lensed CMB in terms of the unlensed CMB and the lensing potential. So far there has been no feasible algorithm for inferring the posterior distribution of the lensing potential from the lensed CMB map. We propose a solution that allows efficient Markov Chain Monte Carlo sampling from the joint posterior of the lensing potential and the unlensed CMB map using the Hamiltonian Monte Carlo technique. The main conceptual step in the solution is a re-parameterization of CMB lensing in terms of the lensed CMB and the "inverse lensing" potential. We demonstrate a fast implementation on simulated data including noise and a sky cut, that uses a further acceleration based on a very mild approximation of the inverse lensing potential. We find that the resulting Markov Chain has short correlation lengths and excellent convergence properties, making it promising for application to high resolution CMB data sets of the future.

The reconstruction of the CMB lensing potential is based on a Taylor expansion of lensing effects which is known to have poor convergence properties. For lensing of temperature fluctuations, an understanding of the higher order terms in this expansion which is accurate enough for current experimental sensitivity levels has been developed in Hanson et. al. (2010), as well as a slightly modified Okamoto and Hu quadratic estimator which incorporates lensed rather than unlensed spectra into the estimator weights to mitigate the effect of higher order terms. We extend these results in several ways: (1) We generalize this analysis to the full set of quadratic temperature/polarization lensing estimators, (2) We study the effect of higher order terms for more futuristic experimental noise levels, (3) We show that the ability of the modified quadratic estimator to mitigate the effect of higher order terms relies on a delicate cancellation which occurs only when the true lensed spectra are known. We investigate the sensitivity of this cancellation to uncertainties in or knowledge of these spectra. We find that higher order terms in the Taylor expansion can impact projected error bars at experimental sensitivities similar to those found in future ACTpol/SPTpol experiments.

We study a regression characterization for the quadratic estimator of weak lensing, developed by Hu and Okamoto (2001,2002), for cosmic microwave background observations. This characterization motivates a modification of the quadratic estimator by an adaptive Wiener filter which uses the robust Bayesian techniques described in Strawderman (1971) and Berger (1980). This technique requires the user to propose a fiducial model for the spectral density of the unknown lensing potential but the resulting estimator is developed to be robust to misspecification of this model. The role of the fiducial spectral density is to give the estimator superior statistical performance in a "neighborhood of the fiducial model" while controlling the statistical errors when the fiducial spectral density is drastically wrong. Our estimate also highlights some advantages provided by a Bayesian analysis of the quadratic estimator.

We present a new estimation method for mapping the gravitational lensing potential from observed CMB intensity and polarization fields. Our method uses Bayesian techniques to estimate the average curvature of the potential over small local regions. These local curvatures are then used to construct an estimate of a low pass filter of the gravitational potential. By utilizing Bayesian/likelihood methods one can easily overcome problems with missing and/or non-uniform pixels and problems with partial sky observations (E and B mode mixing, for example). Moreover, our methods are local in nature which allow us to easily model spatially varying beams and are highly parallelizable. We note that our estimates do not rely on the typical Taylor approximation which is used to construct estimates of the gravitational potential by Fourier coupling. We present our methodology with a flat sky simulation under nearly ideal experimental conditions with a noise level of 1 $\mu K$-arcmin for the temperature field, $\sqrt{2}$ $\mu K$-arcmin for the polarization fields, with an instrumental beam full width at half maximum (FWHM) of 0.25 arcmin.