
We present weaklensing measurements using the firstyear data from the Hyper
SuprimeCam Strategic Survey Program on the Subaru telescope for eight galaxy
clusters selected through their thermal SunyaevZel'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 signaltonoise ratio of individual cluster lensing
measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked
cluster weaklensing signal. We fit for an average weaklensing mass
distribution using three different profiles, a NavarroFrenkWhite profile, a
darkmatteronly 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 SZestimated mass to the
lensingestimated mass (the socalled hydrostatic mass bias $1b$) of
$0.74^{+0.13}_{0.12}$, which is comparable to previous SZselected 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 Otype stars, and Smith & Tombleson suggested that these large
separations are inconsistent with a singlestar 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 singlestar 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 tomass mapping or the masstoage mapping. In this paper, we
explore binary solutions that modify the masstoage 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
massgainer 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 SunyaevZel'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 crosscorrelation at the $3.4\sigma$ level: voids are
underpressured 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 lowmass
halos, or nonlocal 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 percentlevel 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
lowenergy window onto the highenergy 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 coexistence" 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
coexistence will be fully probed by proposed CMBS4 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
dynamicallyevolving 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, sciencebased assessment of any
given cadence simulation, that will enable the scheduling of the LSST to be as
wellinformed 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 bestfit Planck LCDM model
over a range of multipoles L=802100, 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,
selforganized, 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
nextgeneration groundbased cosmic microwave background experiment, CMBS4,
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. CMBS4 will dramatically advance cosmological
studies by crossing critical thresholds in the search for the Bmode
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 nighttime
data collected during 201314 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 200810, in combination with Planck
and WMAP data to estimate cosmological parameters from the temperature,
polarization, and temperaturepolarization crosscorrelations. We find the new
ACTPol data to be consistent with the LCDM model. The ACTPol
temperaturepolarization crossspectrum 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
darkmatter and darkenergy components. CMB data have already shown that
ultralight 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 CMBStage IV (CMBS4) experiment
(assuming a 1 arcmin beam and $< 1~\mu K{\rmarcmin}$ noise levels over a sky
fraction of 0.4) to the density of ULAs and other darksector components is
assessed. CMBS4 data should be $\sim 10$ times more sensitive to the ULA
energydensity 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. CMBS4 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 darkmatter halo density profiles. These improvements will allow for a
multi$\sigma$ detection of percentlevel 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
SunyaevZel'dovich (SZ)selected merging galaxy cluster ACTCL 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 lowestmass
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 ACTCL J0256.5+0006 is composed of two subclusters. We
determine a merger mass ratio of 7:4, and a lineofsight velocity difference
of $v_\perp = 1880 \pm 280$ km s$^{1}$. We construct a simple merger model to
infer relevant timescales in the merger. From its location on the $P_{\rm
1.4GHz}{}L_{\rm X}$ scaling relation, we infer that we observe ACTCL
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 SunyaevZel'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
$(115)\times10^{14}M_\odot$. Comparing with masses estimated from the SZ
distortion assuming a gas pressure profile derived from Xray observations
gives a mean SZtodynamical 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 massvelocity
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 SunyaevZe\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 lineofsight 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 nokSZ 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
SunyaevZe\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.

Ultralight 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
darkmatter or darkenergy density of the Universe. ULAs could suppress the
growth of structure on small scales, or lead to an enhanced integrated
SachsWolfe effect on largescale cosmic microwavebackground (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 galaxyredshift survey. In the mass range 10^{32} eV < m
<10^{25.5} eV, the axion relicdensity \Omega_{a} (relative to the total
darkmatter 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 tensortoscalar 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 opticallyselected CMASS galaxies from the
SDSSIII/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
crossmatched sources below redshift $z=0.2$ to preferentially select Active
Galactic Nuclei (AGN). We measure the angular crosspower spectrum $C_l^{\kappa
g}$ at $4.4\sigma$ significance in the multipole range $100<l<3000$,
corresponding to physical scales between $\approx$ 260 Mpc at an effective
redshift $z_{\rm eff}= 1.5$. Modelling the AGN population with a
redshiftdependent bias, the crossspectrum is well fit by the Planck bestfit
$\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
radioloud 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 reanalyze 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 mapbased
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
prePlanck 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 crossspectra from the publicly
available survey maps instead of their use of the detector set crossspectra
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
crosscorrelating 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 largescale
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 Bmode polarization lensing signal is present with a
significance of 3.2 sigma. We also present the first measurement of CMB
lensingCIB 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 bestfit lensingCIB 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 onepoint probability distribution function
(PDF) of the thermal SunyaevZel'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 signaltonoise 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 followup 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 nonIa
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 SDSSII 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, PanSTARRS 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 nonIa 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$Karcmin. We present TT, TE, EE, TB, EB, and BB power spectra from
three of these regions. The observed Emode polarization power spectrum,
displaying six acoustic peaks in the range $200<\ell<3000$, is an excellent fit
to the prediction of the bestfit 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 groundbased 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
angularscale temperature measurements. Our projections are based on combining
expected smallscale Emode 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
multifrequency channels of most foregrounds in the temperature data, we
investigate the sensitivity to the only residual secondary component, the
kinematic SunyaevZel'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 arcminutescale resolution for future CMB surveys.