We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.

Observations show that luminous blue variables (LBVs) are far more dispersed than massive O-type stars, and Smith & Tombleson suggested that these large separations are inconsistent with a single-star evolution model of LBVs. Instead, they suggested that the large distances are most consistent with binary evolution scenarios. To test these suggestions, we modelled young stellar clusters and their passive dissolution, and we find that, indeed, the standard single-star evolution model is mostly inconsistent with the observed LBV environments. If LBVs are single stars, then the lifetimes inferred from their luminosity and mass are far too short to be consistent with their extreme isolation. This implies that there is either an inconsistency in the luminosity- to-mass mapping or the mass-to-age mapping. In this paper, we explore binary solutions that modify the mass-to-age mapping and are consistent with the isolation of LBVs. For the binary scenarios, our crude models suggest that LBVs are rejuvenated stars. They are either the result of mergers or they are mass gainers and received a kick when the primary star exploded. In the merger scenario, if the primary is about 19 solar masses, then the binary has enough time to wander far afield, merge and form a rejuvenated star. In the mass-gainer and kick scenario, we find that LBV isolation is consistent with a wide range of kick velocities, anywhere from 0 to ~ 105 km/s. In either scenario, binarity seems to play a major role in the isolation of LBVs.

We stack maps of the thermal Sunyaev-Zel'dovich effect produced by the Planck Collaboration around the centers of cosmic voids defined by the distribution of galaxies in the CMASS sample of the Baryon Oscillation Spectroscopic Survey, scaled by the void effective radii. We report a first detection of the associated cross-correlation at the $3.4\sigma$ level: voids are under-pressured relative to the cosmic mean. We compare the measured Compton-$y$ profile around voids with a model based solely on the spatial modulation of halo abundance with environmental density. The amplitude of the detected signal is marginally lower than predicted by an overall amplitude $\alpha_v=0.67\pm0.2$. We discuss the possible interpretations of this measurement in terms of modelling uncertainties, excess pressure in low-mass halos, or non-local heating mechanisms.

The cosmic microwave background (CMB) places strong constraints on models of dark matter (DM) that deviate from standard cold DM (CDM), and on initial conditions beyond the scalar adiabatic mode. Here, the full \textit{Planck} data set (including temperature, $E$-mode polarisation, and lensing deflection) is used to test the possibility that some fraction of the DM is composed of ultralight axions (ULAs). This represents the first use of CMB lensing to test the ULA model. We find no evidence for a ULA component in the mass range $10^{-33}\leq m_a\leq 10^{-24}\text{ eV}$. We put percent-level constraints on the ULA contribution to the DM, improving by up to a factor of two compared to the case with temperature anisotropies alone. Axion DM also provides a low-energy window onto the high-energy physics of inflation through the interplay between the vacuum misalignment production of axions and isocurvature perturbations. We perform the first systematic investigation into the parameter space of ULA isocurvature, using an accurate isocurvature transfer function at all $m_{a}$ values. We precisely identify a "window of co-existence" for $10^{-25}\text{ eV}\leq m_a\leq10^{-24}\text{ eV}$ where the data allow, simultaneously, a $\sim10\%$ contribution of ULAs to the DM, and $\sim 1\%$ contributions of isocurvature and tensors to the CMB power. ULAs in this window (and \textit{all} lighter ULAs) are shown to be consistent with a large inflationary Hubble parameter, $H_I\sim 10^{14}\text{ GeV}$. The window of co-existence will be fully probed by proposed CMB-S4 observations with increased accuracy in the high-$\ell$ lensing power and low-$\ell$ $E$ and $B$-mode polarisation. If ULAs in the window exist, this could allow for two independent measurements of $H_I$ in the CMB using the axion DM content and isocurvature, and the tensor contribution to $B$-modes.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

We report a measurement of the power spectrum of cosmic microwave background (CMB) lensing from two seasons of Atacama Cosmology Telescope Polarimeter (ACTPol) CMB data. The CMB lensing power spectrum is extracted from both temperature and polarization data using quadratic estimators. We obtain results that are consistent with the expectation from the best-fit Planck LCDM model over a range of multipoles L=80-2100, with an amplitude of lensing A_lens = 1.06 +/- 0.15 (stat.) +/- 0.06 (sys.) relative to Planck. Our measurement of the CMB lensing power spectrum gives sigma_8 Omega_m^0.25 = 0.643 +/- 0.054; including baryon acoustic oscillation scale data, we constrain the amplitude of density fluctuations to be sigma_8 = 0.831 +/- 0.053. We also update constraints on the neutrino mass sum. We verify our lensing measurement with a number of null tests and systematic checks, finding no evidence of significant systematic errors. This measurement relies on a small fraction of the ACTPol data already taken; more precise lensing results can therefore be expected from the full ACTPol dataset.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

We present the temperature and polarization angular power spectra measured by the Atacama Cosmology Telescope Polarimeter (ACTPol). We analyze night-time data collected during 2013-14 using two detector arrays at 149 GHz, from 548 deg$^2$ of sky on the celestial equator. We use these spectra, and the spectra measured with the MBAC camera on ACT from 2008-10, in combination with Planck and WMAP data to estimate cosmological parameters from the temperature, polarization, and temperature-polarization cross-correlations. We find the new ACTPol data to be consistent with the LCDM model. The ACTPol temperature-polarization cross-spectrum now provides stronger constraints on multiple parameters than the ACTPol temperature spectrum, including the baryon density, the acoustic peak angular scale, and the derived Hubble constant. Adding the new data to planck temperature data tightens the limits on damping tail parameters, for example reducing the joint uncertainty on the number of neutrino species and the primordial helium fraction by 20%.

Measurements of cosmic microwave background (CMB) anisotropies provide strong evidence for the existence of dark matter and dark energy. They can also test its composition, probing the energy density and particle mass of different dark-matter and dark-energy components. CMB data have already shown that ultra-light axions (ULAs) with mass in the range $10^{-32}~{\rm eV} \to 10^{-26}~{\rm eV}$ compose a fraction $< 0.01$ of the cosmological critical density. Here, the sensitivity of a proposed CMB-Stage IV (CMB-S4) experiment (assuming a 1 arcmin beam and $< 1~\mu K{\rm-arcmin}$ noise levels over a sky fraction of 0.4) to the density of ULAs and other dark-sector components is assessed. CMB-S4 data should be $\sim 10$ times more sensitive to the ULA energy-density than Planck data alone, across a wide range of ULA masses $10^{-32}< m_{a}< 10^{-23}~{\rm eV}$, and will probe axion decay constants of $f_{a}\approx 10^{16}~{\rm GeV}$, at the grand unified scale. CMB-S4 could improve the CMB lower bound on the ULA mass from $\sim 10^{-25}~{\rm eV}$ to $10^{-23}~{\rm eV}$, nearing the mass range probed by dwarf galaxy abundances and dark-matter halo density profiles. These improvements will allow for a multi-$\sigma$ detection of percent-level departures from CDM over a wide range of masses. Much of this improvement is driven by the effects of weak gravitational lensing on the CMB, which breaks degeneracies between ULAs and neutrinos. We also find that the addition of ULA parameters does not significantly degrade the sensitivity of the CMB to neutrino masses. These results were obtained using the axionCAMB code (a modification to the CAMB Boltzmann code), presented here for public use.

We present the detection of a giant radio halo (GRH) in the Sunyaev-Zel'dovich (SZ)-selected merging galaxy cluster ACT-CL J0256.5+0006 ($z = 0.363$), observed with the Giant Metrewave Radio Telescope at 325 MHz and 610 MHz. We find this cluster to host a faint ($S_{610} = 5.6 \pm 1.4$ mJy) radio halo with an angular extent of 2.6 arcmin, corresponding to 0.8 Mpc at the cluster redshift, qualifying it as a GRH. J0256 is one of the lowest-mass systems, $M_{\rm 500,SZ} = (5.0 \pm 1.2) \times 10^{14} M_\odot$, found to host a GRH. We measure the GRH at lower significance at 325 MHz ($S_{325} = 10.3 \pm 5.3$ mJy), obtaining a spectral index measurement of $\alpha^{610}_{325} = 1.0^{+0.7}_{-0.9}$. This result is consistent with the mean spectral index of the population of typical radio halos, $\alpha = 1.2 \pm 0.2$. Adopting the latter value, we determine a 1.4 GHz radio power of $P_{1.4\text{GHz}} = (1.0 \pm 0.3) \times 10^{24}$ W Hz$^{-1}$, placing this cluster within the scatter of known scaling relations. Various lines of evidence, including the ICM morphology, suggest that ACT-CL J0256.5+0006 is composed of two subclusters. We determine a merger mass ratio of 7:4, and a line-of-sight velocity difference of $v_\perp = 1880 \pm 280$ km s$^{-1}$. We construct a simple merger model to infer relevant time-scales in the merger. From its location on the $P_{\rm 1.4GHz}{-}L_{\rm X}$ scaling relation, we infer that we observe ACT-CL J0256.5+0006 approximately 500 Myr before first core crossing.

We present galaxy velocity dispersions and dynamical mass estimates for 44 galaxy clusters selected via the Sunyaev-Zel'dovich (SZ) effect by the Atacama Cosmology Telescope. Dynamical masses for 18 clusters are reported here for the first time. Using \Nbody\ simulations, we model the different observing strategies used to measure the velocity dispersions and account for systematic effects resulting from these strategies. We find that the galaxy velocity distributions may be treated as isotropic, and that an aperture correction of up to 7 per cent in the velocity dispersion is required if the spectroscopic galaxy sample is sufficiently concentrated towards the cluster centre. Accounting for the radial profile of the velocity dispersion in simulations enables consistent dynamical mass estimates regardless of the observing strategy. Cluster masses $M_{200}$ are in the range $(1-15)\times10^{14}M_\odot$. Comparing with masses estimated from the SZ distortion assuming a gas pressure profile derived from X-ray observations gives a mean SZ-to-dynamical mass ratio of $1.10\pm0.13$, but there is an additional 0.14 systematic uncertainty due to the unknown velocity bias; the statistical uncertainty is dominated by the scatter in the mass-velocity dispersion scaling relation. This ratio is consistent with previous determinations at these mass scales.

We use microwave temperature maps from two seasons of data from the Atacama Cosmology Telescope (ACTPol) at 146 GHz, together with the Constant Mass CMASS galaxy sample from the Baryon Oscillation Spectroscopic Survey to measure the kinematic Sunyaev-Ze\v{l}dovich (kSZ) effect over the redshift range z = 0.4 - 0.7. We use galaxy positions and the continuity equation to obtain a reconstruction of the line-of-sight velocity field. We stack the cosmic microwave background temperature at the location of each halo, weighted by the corresponding reconstructed velocity. The resulting best fit kSZ model is preferred over the no-kSZ hypothesis at 3.3sigma and 2.9sigma for two independent velocity reconstruction methods, using 25,537 galaxies over 660 square degrees. The effect of foregrounds that are uncorrelated with the galaxy velocities is expected to be well below our signal, and residual thermal Sunyaev-Ze\v{l}dovich contamination is controlled by masking the most massive clusters. Finally, we discuss the systematics involved in converting our measurement of the kSZ amplitude into the mean free electron fraction of the halos in our sample.

There are two redshifts in cosmology: $z_{obs}$, the observed redshift computed via spectral lines, and the model redshift, $z$, defined by the effective FLRW scale factor. In general these do not coincide. We place observational constraints on the allowed distortions of $z$ away from $z_{obs}$ - a possibility we dub redshift remapping. Remapping is degenerate with cosmic dynamics for either $d_L(z)$ or $H(z)$ observations alone: for example, the simple remapping $z = \alpha_1 z_{obs} +\alpha_2 z_{obs}^2$ allows a decelerating Einstein de Sitter universe to fit the observed supernova Hubble diagram as successfully as $\Lambda$CDM, highlighting that supernova data alone cannot prove that the universe is accelerating. We show however, that redshift remapping leads to apparent violations of cosmic distance duality that can be used to detect its presence even when neither a specific theory of gravity nor the Copernican Principle are assumed. Combining current data sets favours acceleration but does not yet rule out redshift remapping as an alternative to dark energy. Future surveys, however, will provide exquisite constraints on remapping and any models -- such as backreaction -- that predict it.

Ultra-light axions (ULAs) with masses in the range 10^{-33} eV <m <10^{-20} eV are motivated by string theory and might contribute to either the dark-matter or dark-energy density of the Universe. ULAs could suppress the growth of structure on small scales, or lead to an enhanced integrated Sachs-Wolfe effect on large-scale cosmic microwave-background (CMB) anisotropies. In this work, cosmological observables over the full ULA mass range are computed, and then used to search for evidence of ULAs using CMB data from the Wilkinson Microwave Anisotropy Probe (WMAP), Planck satellite, Atacama Cosmology Telescope, and South Pole Telescope, as well as galaxy clustering data from the WiggleZ galaxy-redshift survey. In the mass range 10^{-32} eV < m <10^{-25.5} eV, the axion relic-density \Omega_{a} (relative to the total dark-matter relic density \Omega_{d}) must obey the constraints \Omega_{a}/\Omega_{d} < 0.05 and \Omega_{a}h^{2} < 0.006 at 95%-confidence. For m> 10^{-24} eV, ULAs are indistinguishable from standard cold dark matter on the length scales probed, and are thus allowed by these data. For m < 10^{-32} eV, ULAs are allowed to compose a significant fraction of the dark energy.

We present the first attempt to use a combination of CMB, LIGO, and PPTA data to constrain both the tilt and the running of primordial tensor power spectrum through constraints on the gravitational wave energy density generated in the early universe. Combining measurements at different cosmological scales highlights how complementary data can be used to test the predictions of early universe models including the inflationary consistency relation. Current data prefers a slightly positive tilt ($n_t = 0.06^{+0.63}_{-0.89}$) and a negative running ($n_{t, {\rm run}} < -0.22$) for the tensor power spectrum spectrum. Interestingly, the addition of direct gravitational wave detector data alone puts strong bounds on the tensor-to-scalar ratio $r < 0.2$ since the large positive tensor tilt preferred by the \textit{Planck} temperature power spectrum is no longer allowed. Adding the recently released BICEP2/KECK and Planck 353 GHz Polarization cross correlation data puts an even stronger bound $r<0.1$ . We comment on possible effects of a large positive tilt on the background expansion and show that depending on the assumptions regarding the UV cutoff ($k_{\rm UV}/k_* = 10^{24}$) of the primordial spectrum of gravitational waves, the strongest bounds on $n_t = 0.07^{+0.52}_{-0.80}$ are derived from this effect.

We present evidence of the gravitational lensing of the cosmic microwave background by $10^{13}$ solar mass dark matter halos. Lensing convergence maps from the Atacama Cosmology Telescope Polarimeter (ACTPol) are stacked at the positions of around 12,000 optically-selected CMASS galaxies from the SDSS-III/BOSS survey. The mean lensing signal is consistent with simulated dark matter halo profiles, and is favored over a null signal at 3.2 sigma significance. This result demonstrates the potential of microwave background lensing to probe the dark matter distribution in galaxy group and galaxy cluster halos.

We correlate the positions of radio galaxies in the FIRST survey with the CMB lensing convergence estimated from the Atacama Cosmology Telescope over 470 square degrees to determine the bias of these galaxies. We remove optically cross-matched sources below redshift $z=0.2$ to preferentially select Active Galactic Nuclei (AGN). We measure the angular cross-power spectrum $C_l^{\kappa g}$ at $4.4\sigma$ significance in the multipole range $100<l<3000$, corresponding to physical scales between $\approx$ 2--60 Mpc at an effective redshift $z_{\rm eff}= 1.5$. Modelling the AGN population with a redshift-dependent bias, the cross-spectrum is well fit by the Planck best-fit $\Lambda$CDM cosmological model. Fixing the cosmology we fit for the overall bias model normalization, finding $b(z_{\rm eff}) = 3.5 \pm 0.8$ for the full galaxy sample, and $b(z_{\rm eff})=4.0\pm1.1 (3.0\pm1.1)$ for sources brighter (fainter) than 2.5 mJy. This measurement characterizes the typical halo mass of radio-loud AGN: we find $\log(M_{\rm halo} / M_\odot) = 13.6^{+0.3}_{-0.4}$.

The tension between the best fit parameters derived by the Planck team and a number of other astronomical measurements suggests either systematics in the astronomical measurements, systematics in the Planck data, the need for new physics, or a combination thereof. We re-analyze the Planck data and find that the $217\,\text{GHz}\times 217\,\text{GHz}$ detector set spectrum used in the Planck analysis is responsible for some of this tension. We use a map-based foreground cleaning procedure, relying on a combination of 353 GHz and 545 GHz maps to reduce residual foregrounds in the intermediate frequency maps used for cosmological inference. For our baseline data analysis, which uses 47% of the sky and makes use of both 353 and 545 GHz data for foreground cleaning, we find the $\Lambda$CDM cosmological parameters $\Omega_c h^2 = 0.1170 \pm 0.0025$, $n_s = 0.9686 \pm 0.0069$, $H_0 = 68.0 \pm 1.1\,\mathrm{km} \mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, $\Omega_b h^2 = 0.02197 \pm 0.00026$, $\ln 10^{10}A_s = 3.082 \pm 0.025$, and $\tau = 0.090 \pm 0.013 $. While in broad agreement with the results reported by the Planck team, these revised parameters imply a universe with a lower matter density of $\Omega_m=0.302\pm0.015$, and parameter values generally more consistent with pre-Planck CMB analyses and astronomical observations. We compare our cleaning procedure with the foreground modeling used by the Planck team and find good agreement. The difference in parameters between our analysis and that of the Planck team is mostly due to our use of cross-spectra from the publicly available survey maps instead of their use of the detector set cross-spectra which include pixels only observed in one of the surveys. We show evidence suggesting residual systematics in the detector set spectra used in the Planck likelihood code, which is substantially reduced for our spectra.

We present a measurement of the gravitational lensing of the Cosmic Microwave Background (CMB) temperature and polarization fields obtained by cross-correlating the reconstructed convergence signal from the first season of ACTPol data at 146 GHz with Cosmic Infrared Background (CIB) fluctuations measured using the Planck satellite. Using an overlap area of 206 square degrees, we detect gravitational lensing of the CMB polarization by large-scale structure at a statistical significance of 4.5 sigma. Combining both CMB temperature and polarization data gives a lensing detection at 9.1 sigma significance. A B-mode polarization lensing signal is present with a significance of 3.2 sigma. We also present the first measurement of CMB lensing--CIB correlation at small scales corresponding to l > 2000. Null tests and systematic checks show that our results are not significantly biased by astrophysical or instrumental systematic effects, including Galactic dust. Fitting our measurements to the best-fit lensing-CIB cross power spectrum measured in Planck data, scaled by an amplitude A, gives A=1.02 +0.12/-0.18 (stat.) +/-0.06(syst.), consistent with the Planck results.

We present a measurement of the one-point probability distribution function (PDF) of the thermal Sunyaev-Zel'dovich (tSZ) decrement in the pixel temperature histogram of filtered 148 GHz sky maps from the Atacama Cosmology Telescope (ACT). The PDF includes the signal from all galaxy clusters in the map, including objects below the signal-to-noise threshold for individual detection, making it a particularly sensitive probe of the amplitude of matter density perturbations, $\sigma_8$. We use a combination of analytic halo model calculations and numerical simulations to compute the theoretical tSZ PDF and its covariance matrix, accounting for all noise sources and including relativistic corrections. From the measured ACT 148 GHz PDF alone, we find $\sigma_8 = 0.793 \pm 0.018$, with additional systematic errors of $\pm 0.017$ due to uncertainty in intracluster medium gas physics and $\pm 0.006$ due to uncertainty in infrared point source contamination. Using effectively the same data set, the statistical error here is a factor of two lower than that found in ACT's previous $\sigma_8$ determination based solely on the skewness of the tSZ signal. In future temperature maps with higher sensitivity, the tSZ PDF will break the degeneracy between intracluster medium gas physics and cosmological parameters.

For future surveys, spectroscopic follow-up for all supernovae will be extremely difficult. However, one can use light curve fitters, to obtain the probability that an object is a Type Ia. One may consider applying a probability cut to the data, but we show that the resulting non-Ia contamination can lead to biases in the estimation of cosmological parameters. A different method, which allows the use of the full dataset and results in unbiased cosmological parameter estimation, is Bayesian Estimation Applied to Multiple Species (BEAMS). BEAMS is a Bayesian approach to the problem which includes the uncertainty in the types in the evaluation of the posterior. Here we outline the theory of BEAMS and demonstrate its effectiveness using both simulated datasets and SDSS-II data. We also show that it is possible to use BEAMS if the data are correlated, by introducing a numerical marginalisation over the types of the objects. This is largely a pedagogical introduction to BEAMS with references to the main BEAMS papers.

New supernova surveys such as the Dark Energy Survey, Pan-STARRS and the LSST will produce an unprecedented number of photometric supernova candidates, most with no spectroscopic data. Avoiding biases in cosmological parameters due to the resulting inevitable contamination from non-Ia supernovae can be achieved with the BEAMS formalism, allowing for fully photometric supernova cosmology studies. Here we extend BEAMS to deal with the case in which the supernovae are correlated by systematic uncertainties. The analytical form of the full BEAMS posterior requires evaluating 2^N terms, where N is the number of supernova candidates. This `exponential catastrophe' is computationally unfeasible even for N of order 100. We circumvent the exponential catastrophe by marginalising numerically instead of analytically over the possible supernova types: we augment the cosmological parameters with nuisance parameters describing the covariance matrix and the types of all the supernovae, \tau_i, that we include in our MCMC analysis. We show that this method deals well even with large, unknown systematic uncertainties without a major increase in computational time, whereas ignoring the correlations can lead to significant biases and incorrect credible contours. We then compare the numerical marginalisation technique with a perturbative expansion of the posterior based on the insight that future surveys will have exquisite light curves and hence the probability that a given candidate is a Type Ia will be close to unity or zero, for most objects. Although this perturbative approach changes computation of the posterior from a 2^N problem into an N^2 or N^3 one, we show that it leads to biases in general through a small number of misclassifications, implying that numerical marginalisation is superior.

We report on measurements of the cosmic microwave background (CMB) and celestial polarization at 146 GHz made with the Atacama Cosmology Telescope Polarimeter (ACTPol) in its first three months of observing. Four regions of sky covering a total of 270 square degrees were mapped with an angular resolution of $1.3'$. The map noise levels in the four regions are between 11 and 17 $\mu$K-arcmin. We present TT, TE, EE, TB, EB, and BB power spectra from three of these regions. The observed E-mode polarization power spectrum, displaying six acoustic peaks in the range $200<\ell<3000$, is an excellent fit to the prediction of the best-fit cosmological models from WMAP9+ACT and Planck data. The polarization power spectrum, which mainly reflects primordial plasma velocity perturbations, provides an independent determination of cosmological parameters consistent with those based on the temperature power spectrum, which results mostly from primordial density perturbations. We find that without masking any point sources in the EE data at $\ell<9000$, the Poisson tail of the EE power spectrum due to polarized point sources has an amplitude less than $2.4$ $\mu$K$^2$ at $\ell = 3000$ at 95\% confidence. Finally, we report that the Crab Nebula, an important polarization calibration source at microwave frequencies, has 8.7\% polarization with an angle of $150.7^\circ \pm 0.6^\circ$ when smoothed with a $5'$ Gaussian beam.

Future arcminute resolution polarization data from ground-based Cosmic Microwave Background (CMB) observations can be used to estimate the contribution to the temperature power spectrum from the primary anisotropies and to uncover the signature of reionization near $\ell=1500$ in the small angular-scale temperature measurements. Our projections are based on combining expected small-scale E-mode polarization measurements from Advanced ACTPol in the range $300<\ell<3000$ with simulated temperature data from the full Planck mission in the low and intermediate $\ell$ region, $2<\ell<2000$. We show that the six basic cosmological parameters determined from this combination of data will predict the underlying primordial temperature spectrum at high multipoles to better than $1\%$ accuracy. Assuming an efficient cleaning from multi-frequency channels of most foregrounds in the temperature data, we investigate the sensitivity to the only residual secondary component, the kinematic Sunyaev-Zel'dovich (kSZ) term. The CMB polarization is used to break degeneracies between primordial and secondary terms present in temperature and, in effect, to remove from the temperature data all but the residual kSZ term. We estimate a $15 \sigma$ detection of the diffuse homogeneous kSZ signal from expected AdvACT temperature data at $\ell>1500$, leading to a measurement of the amplitude of matter density fluctuations, $\sigma_8$, at $1\%$ precision. Alternatively, by exploring the reionization signal encoded in the patchy kSZ measurements, we bound the time and duration of the reionization with $\sigma(z_{\rm re})=1.1$ and $\sigma(\Delta z_{\rm re})=0.2$. We find that these constraints degrade rapidly with large beam sizes, which highlights the importance of arcminute-scale resolution for future CMB surveys.