The capacity of Stage IV lensing surveys to measure the neutrino mass sum and differentiate between the normal and inverted mass hierarchies depends on the impact of nuisance parameters describing small-scale baryonic astrophysics and intrinsic alignments. For a Euclid-like survey, we perform the first combined weak lensing and galaxy clustering Fisher analysis with baryons, intrinsic alignments, and massive neutrinos for both hierarchies. We use a matter power spectrum generated from a halo model that captures the impact of baryonic feedback and adiabatic contraction. For weak lensing, we find that baryons cause severe degradation to forecasts of the neutrino mass sum, $\Sigma$, approximately doubling $\sigma_{\Sigma}$. We show that including galaxy clustering constraints from Euclid and BOSS, and cosmic microwave background (CMB) Planck priors, can reduce this degradation to $\sigma_{\Sigma}$ to 9% and 16% for the normal and inverted hierarchies respectively. The combined forecasts, $\sigma_{\Sigma_{\rm{NH}}}=0.034\, \rm{eV}$ and $\sigma_{\Sigma_{\rm{IH}}}=0.034\, \rm{eV}$, preclude a meaningful distinction of the hierarchies but could be improved upon with future CMB priors on $n_s$ and information from neutrinoless double beta decay to achieve a 2$\sigma$ distinction. The effect of intrinsic alignments on forecasts is shown to be minimal, with $\sigma_{\Sigma}$ even experiencing mild improvements due to information from the intrinsic alignment signal. We find that while adiabatic contraction and intrinsic alignments will require careful calibration to prevent significant biasing of $\Sigma$, there is less risk presented by feedback from energetic events like AGN and supernovae.

We conduct a thorough investigation into the possibility that residing in an overdense region of the Universe may induce bias in measurements of the large-scale structure. We compute the conditional correlation function and angular power spectrum of density and lensing fluctuations while holding the local spherically averaged density fixed and show that for Gaussian fields this has no effect on the angular power at $l>0$. We identify a range of scales where a perturbative approach allows analytic progress to be made, and we compute leading-order conditional power spectra using an Edgeworth expansion and second-order perturbation theory. We find no evidence for any significant bias to cosmological power spectra from our local density contrast. We show that when smoothed over a large region around the observer, conditioning on the local density typically affects density power spectra by less than a percent at cosmological distances, below cosmic variance. We find that while typical corrections to the lensing angular power spectrum can be at the 10% level on the largest angular scales and for source redshifts $z_s \lesssim 0.1$, for the typical redshifts targeted by upcoming wide imaging surveys the corrections are sub-percent and negligible, in contrast to previous claims in the literature. Using an estimate of the local spherically averaged density from a composite galaxy redshift catalogue we find that the corrections from conditioning on our own local density are below cosmic variance and subdominant to other non-linear effects. We discuss the potential implications of our results for cosmology and point out that a measurement of the local density contrast may be used as a consistency test of cosmological models.

We have calculated the detectable merger rate of primordial black holes, as a function of the redshift, as well as the binary's mass ratio, total mass and chirp mass (observables that have not previously been explored in great detail for PBHs). We consider both the current and design sensitivity of LIGO and five different primordial black hole mass functions, as well as showing a comparison to a predicted astrophysical black hole merger rate. We show that the empirical preference for nearly equal-mass binaries in current LIGO/Virgo data can be consistent with a PBH hypothesis once observational selection effects are taken into account. However, current data do exclude some PBH mass distributions, and future data may be able to rule out the possibility that all observed BH mergers had a primordial origin.

We conduct a thorough study into the feasibility of measuring large-scale correlated proper motions of galaxies with astrometric surveys. We introduce a harmonic formalism for analysing proper motions and their correlation functions on the sphere based on spin-weighted spherical harmonics, and study the statistics of the transverse velocity field induced by large-scale structure. We use a likelihood formalism to derive optimal estimators for the secular parallax due to the Solar System's motion relative to distant objects, and compute the variance and bias due to peculiar velocities and relativistic aberration. We use a simulated catalogue of galaxy proper motions with radial distributions and noise properties similar to those expected from Gaia to forecast the detectability of the proper motion dipole, whose amplitude may be considered a proxy for the Hubble constant. We find cosmic variance to be the limiting source of noise for this measurement, forecasting a detectability of $1$-$2\sigma$ on a single component of the local velocity, increasing to $2$-$4\sigma$ (equivalent to a 25%-50% measurement of the Hubble constant) if the CMB dipole is included as prior information. We conduct a thorough study into the radial dependence of the signal-to-noise, finding that most of the information comes from galaxies closer than a few hundred Mpc. We forecast that the amplitude of peculiar transverse velocities can potentially be measured with 10$\sigma$ significance; such a measurement would offer a unique probe of cosmic flows and a valuable test of the cosmological model.

Accurate and precise covariance matrices will be important in enabling planned cosmological surveys to detect new physics. Standard methods imply either the need for many N-body simulations in order to obtain an accurate estimate, or a precise theoretical model. We combine these approaches by constructing a likelihood function conditioned on simulated and theoretical covariances, consistently propagating noise from the finite number of simulations and uncertainty in the theoretical model itself using an informative Inverse-Wishart prior. Unlike standard methods, our approach allows the required number of simulations to be less than the number of summary statistics. We recover the linear 'shrinkage' covariance estimator in the context of a Bayesian data model, and test our marginal likelihood on simulated mock power spectrum estimates. We conduct a thorough investigation into the impact of prior confidence in different choices of covariance models on the quality of model fits and parameter variances. In a simplified setting we find that the number of simulations required can be reduced if one is willing to accept a mild degradation in the quality of model fits, finding that even weakly informative priors can help to reduce the simulation requirements. We identify the correlation matrix of the summary statistics as a key quantity requiring careful modelling. Our approach can be easily generalized to any covariance model or set of summary statistics, and elucidates the role of hybrid estimators in cosmological inference.

This paper reports the first results of a direct dark matter search with the DEAP-3600 single-phase liquid argon (LAr) detector. The experiment was performed 2 km underground at SNOLAB (Sudbury, Canada) utilizing a large target mass, with the LAr target contained in a spherical acrylic vessel of 3600 kg capacity. The LAr is viewed by an array of PMTs, which would register scintillation light produced by rare nuclear recoil signals induced by dark matter particle scattering. An analysis of 4.44 live days (fiducial exposure of 9.87 tonne-days) of data taken with the nearly full detector during the initial filling phase demonstrates the detector performance and the best electronic recoil rejection using pulse-shape discrimination in argon, with leakage $<1.2\times 10^{-7}$ (90% C.L.) between 16 and 33 keV$_{ee}$. No candidate signal events are observed, which results in the leading limit on WIMP-nucleon spin-independent cross section on argon, $<1.2\times 10^{-44}$ cm$^2$ for a 100 GeV/c$^2$ WIMP mass (90% C.L.).

Baryonic effects on large-scale structure, if not accounted for, can significantly bias dark energy constraints. As the detailed physics of the baryons is not yet well-understood, correcting for baryon effects introduces additional parameters which must be marginalized over, increasing the uncertainties on the inferred cosmological parameters. Forthcoming weak lensing surveys are aiming for percent-level precision on the dark energy equation of state, so the problem must be thoroughly examined. We use a halo model with analytic modifications which capture the impact of adiabatic contraction of baryons and feedback on the matter power spectrum, and generalize the Navarro-Frenk-White profile to account for a possible inner core. A Fisher analysis predicts degradations of 40% in the $w_0$-$w_a$ Figure of Merit for a Euclid-like survey, and up to 80% for other cosmological parameters. We forecast potential inner core constraints of a few $\mathrm{kpc}$, while for a fixed inner core, adiabatic concentration and feedback parameters are constrained to a few percent. We explore the scales where baryons and dark energy contribute most to the Fisher information, finding that probing to increasingly non-linear scales does little to reduce degradation. Including external baryon information improves our forecasts, but limiting degradation to 1% requires strong priors. Adding Planck cosmic microwave background priors improves the Figure of Merit by a factor of 2.7 and almost completely recovers the individual marginalized errors on $w_0$ and $w_a$. We also quantify the calibration of baryon modelling required to reduce biases of dark energy forecasts to acceptable levels for forthcoming lensing surveys.

The Dark matter Experiment using Argon Pulse-shape discrimination (DEAP) has been designed for a direct detection search for particle dark matter using a single-phase liquid argon target. The projected cross section sensitivity for DEAP-3600 to the spin-independent scattering of Weakly Interacting Massive Particles (WIMPs) on nucleons is $10^{-46}~\rm{cm}^{2}$ for a 100 GeV/$c^2$ WIMP mass with a fiducial exposure of 3 tonne-years. This paper describes the physical properties and construction of the DEAP-3600 detector.

We investigate the feasibility of measuring the effects of peculiar velocities in large-scale structure using the dipole of the redshift-space cross-correlation function. We combine number counts of galaxies with brightness-temperature fluctuations from 21cm intensity mapping, demonstrating that the dipole may be measured at modest significance ($\lesssim 2\sigma$) by combining the upcoming radio survey CHIME with the future redshift surveys of DESI and Euclid. More significant measurements ($\lesssim~10\sigma$) will be possible by combining intensity maps from the SKA with these of DESI or Euclid, and an even higher significance measurement ($\lesssim 100\sigma$) may be made by combining observables completely internally to the SKA. We account for effects such as contamination by wide-angle terms, interferometer noise and beams in the intensity maps, non-linear enhancements to the power spectrum, stacking multiple populations, sensitivity to the magnification slope, and the possibility that number counts and intensity maps probe the same tracers. We also derive a new expression for the covariance matrix of multi-tracer redshift-space correlation function estimators with arbitrary orientation weights, which may be useful for upcoming surveys aiming at measuring redshift-space clustering with multiple tracers.

Constraints on cosmological parameters from large-scale structure have traditionally been obtained from two-point statistics. However, non-linear structure formation renders these statistics insufficient in capturing the full information content available, necessitating the measurement of higher-order moments to recover information which would otherwise be lost. We construct quantities based on non-linear and non-local transformations of weakly non-Gaussian fields that Gaussianize the full multivariate distribution at a given order in perturbation theory. Our approach does not require a model of the fields themselves and takes as input only the first few polyspectra, which could be modelled or measured from simulations or data, making our method particularly suited to observables lacking a robust perturbative description such as the weak-lensing shear. We apply our method to simulated density fields, finding a significantly reduced bispectrum and an enhanced correlation with the initial field. We demonstrate that our method reconstructs a large proportion of the linear baryon acoustic oscillations, improving the information content over the raw field by 35%. We apply the transform to toy 21cm intensity maps, showing that our method still performs well in the presence of complications such as redshift-space distortions, beam smoothing, pixel noise, and foreground subtraction. We discuss how this method might provide a route to constructing a perturbative model of the fully non-Gaussian multivariate likelihood function.

Lensing of the CMB is an important effect, and is usually modelled by remapping the unlensed CMB fields by a lensing deflection. However the lensing deflections also change the photon path so that the emission angle is no longer orthogonal to the background last-scattering surface. We give the first calculation of the emission-angle corrections to the standard lensing approximation from dipole (Doppler) sources for temperature and quadrupole sources for temperature and polarization. We show that while the corrections are negligible for the temperature and E-mode polarization, additional large-scale B-modes are produced with a white spectrum that dominates those from post-Born field rotation (curl lensing). On large scales about one percent of the total lensing-induced B-mode amplitude is expected to be due to this effect. However, the photon emission angle does remain orthogonal to the perturbed last-scattering surface due to time delay, and half of the large-scale emission-angle B modes cancel with B modes from time delay to give a total contribution of about half a percent. While not important for planned observations, the signal could ultimately limit the ability of delensing to reveal low amplitudes of primordial gravitational waves.We also derive the rotation of polarization due to multiple deflections between emission and observation. The rotation angle is of quadratic order in the deflection angle, and hence negligibly small: polarization typically rotates by less than an arcsecond, orders of magnitude less than a small-scale image rotates due to post-Born field rotation (which is quadratic in the shear). The field-rotation B modes dominate the other effects on small scales.

We investigate the problem of noise bias in maximum likelihood and maximum a posteriori estimators for cosmic shear. We derive the leading and next-to-leading order biases and compute them in the context of galaxy ellipticity measurements, extending previous work on maximum likelihood inference for weak lensing. We show that a large part of the bias on these point estimators can be removed using information already contained in the likelihood when a galaxy model is specified, without the need for external calibration. We test these bias-corrected estimators on simulated galaxy images similar to those expected from planned space-based weak lensing surveys, with promising results. We find that the introduction of an intrinsic shape prior can help with mitigation of noise bias, such that the maximum a posteriori estimate can be made less biased than the maximum likelihood estimate. Second-order terms offer a check on the convergence of the estimators, but are largely sub-dominant. We show how biases propagate to shear estimates, demonstrating in our simple setup that shear biases can be reduced by orders of magnitude and potentially to within the requirements of planned space-based surveys at mild signal-to-noise. We find that second-order terms can exhibit significant cancellations at low signal-to-noise when Gaussian noise is assumed, which has implications for inferring the performance of shear-measurement algorithms from simplified simulations. We discuss the viability of our point estimators as tools for lensing inference, arguing that they allow for the robust measurement of ellipticity and shear.

We analyse the feasibility of detecting the polarization of the CMB caused by scattering of the remote temperature quadrupole by galaxy clusters with forthcoming CMB polarization surveys. For low-redshift clusters, the signal is strongly correlated with the local large-scale temperature and polarization anisotropies, and the best prospect for detecting the cluster signal is via cross-correlation. For high-redshift clusters, the correlation with the local temperature is weaker and the power in the uncorrelated component of the cluster polarization can be used to enhance detection. We derive linear and quadratic maximum-likelihood estimators for these cases, and forecast signal-to-noise values for the SZ surveys of a Planck-like mission and SPTPol. Our estimators represent an optimal `stacking' analysis of the polarization from clusters. We find that the detectability of the effect is sensitive to the cluster gas density distribution, as well as the telescope resolution, cluster redshift distribution, and sky coverage. We find that the effect is too small to be detected in current and near-future SZ surveys without dedicated polarization follow-up, and that an r.m.s. noise on the Stokes parameters of roughly 1 uK-arcmin for each cluster field is required for a 2-sigma detection, assuming roughly 550 clusters are observed. We discuss and quantify potential biases from the kinetic part of the signal caused by the relative motion of the cluster with respect to the CMB, and from the background CMB polarization behind the cluster, discussing ways in which these biases might be mitigated. Our formalism should be important for next-generation CMB polarization missions, which we argue will be able to measure this effect with high signal-to-noise. This will allow for an important consistency test of the LCDM model on scales that are inaccessible to other probes.

We demonstrate that the intrinsic alignment of galaxies with large-scale tidal fields sources an extra contribution to the recently-detected cross-correlation of galaxy shear and weak lensing of the microwave background. The extra term is the analogy of the 'GI' term in standard cosmic shear studies, and results in a reduction in the amplitude of the cross-correlation. We compute the intrinsic alignment contribution in linear and non-linear theory, and show that it can be at roughly the 15% level for the CFHT Stripe 82 redshift distribution, if the canonical amplitude of intrinsic alignments is assumed. The new term can therefore potentially reconcile the apparently low value of the measured cross-correlation with standard LCDM. We discuss various small-scale effects in the signal and the dependence on the source redshift distribution. We discuss the exciting possibility of self-calibrating intrinsic alignments with a joint analysis of cosmic shear and weak lensing of the microwave background

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in response to the Call for White Papers for the definition of the L2 and L3 Missions in the ESA Science Programme. PRISM would have two instruments: (1) an imager with a 3.5m mirror (cooled to 4K for high performance in the far-infrared---that is, in the Wien part of the CMB blackbody spectrum), and (2) an Fourier Transform Spectrometer (FTS) somewhat like the COBE FIRAS instrument but over three orders of magnitude more sensitive. Highlights of the new science (beyond the obvious target of B-modes from gravity waves generated during inflation) made possible by these two instruments working in tandem include: (1) the ultimate galaxy cluster survey gathering 10e6 clusters extending to large redshift and measuring their peculiar velocities and temperatures (through the kSZ effect and relativistic corrections to the classic y-distortion spectrum, respectively) (2) a detailed investigation into the nature of the cosmic infrared background (CIB) consisting of at present unresolved dusty high-z galaxies, where most of the star formation in the universe took place, (3) searching for distortions from the perfect CMB blackbody spectrum, which will probe a large number of otherwise inaccessible effects (e.g., energy release through decaying dark matter, the primordial power spectrum on very small scales where measurements today are impossible due to erasure from Silk damping and contamination from non-linear cascading of power from larger length scales). These are but a few of the highlights of the new science that will be made possible with PRISM.

We investigate the prospects for constraining alternative theories of gravity with a typical near-term low-budget 21 cm intensity mapping experiment. We derive the 21 cm brightness temperature perturbation consistently in linear theory including all line-of-sight and relativistic effects. We uncover new terms that are a small correction on large scales, analogous to those recently found in the context of galaxy surveys. We then perform a Fisher matrix analysis of the B_0 parametrization of f(R) gravity, where B_0 is proportional to the square of Compton wavelength of the scalaron. We find that our 21 cm survey, in combination with CMB information from Planck, will be able to place a 95% upper limit of 7 x 10^{-5} on B_0 in flat models with a LCDM expansion history, improving on current cosmological constraints by several orders of magnitude. We argue that this constraint is limited by our ability to model the mildly non-linear regime of structure formation in modified gravity. We also perform a model-independent principal component analysis on the free functions introduced into the field equations by modified gravity, mu and Sigma. We find that 20--30 modes of the free functions will be `well-constrained' by our combination of observables, the lower and upper limits dependent on the criteria used to define the `goodness' of the constraint. These constraints are found to be robust to uncertainties in the time-dependence of the bias. Our analysis reveals that our observables are sensitive primarily to temporal variations in Sigma and scale variations in mu. We argue that the inclusion of 21 cm intensity maps will significantly improve constraints on any cosmological deviations from General Relativity in large-scale structure in a very cost-effective manner.

We forecast constraints on cosmological parameters with primary CMB anisotropy information and weak lensing reconstruction with a future post-Planck CMB experiment, the Cosmic Origins Explorer (COrE), using oscillation data on the neutrino mass splittings as prior information. Our MCMC simulations in flat models with a non-evolving equation-of-state of dark energy w give typical 68% upper bounds on the total neutrino mass of 0.136 eV and 0.098 eV for the inverted and normal hierarchies respectively, assuming the total summed mass is close to the minimum allowed by the oscillation data for the respective hierarchies (0.10 eV and 0.06 eV). Including information from future baryon acoustic oscillation measurements with the complete BOSS, Type 1a supernovae distance moduli from WFIRST, and a realistic prior on the Hubble constant, these upper limits shrink to 0.118 eV and 0.080 eV for the inverted and normal hierarchies, respectively. Addition of these distance priors also yields percent-level constraints on w. We find tension between our MCMC results and the results of a Fisher matrix analysis, most likely due to a strong geometric degeneracy between the total neutrino mass, the Hubble constant, and w in the unlensed CMB power spectra. If the minimal-mass, normal hierarchy were realised in nature, the inverted hierarchy should be disfavoured by the full data combination at typically greater than the 2-sigma level. For the minimal-mass inverted hierarchy, we compute the Bayes' factor between the two hierarchies for various combinations of our forecast datasets, and find that the future probes considered here should be able to provide `strong' evidence (odds ratio 12:1) for the inverted hierarchy. Finally, we consider potential biases of the other cosmological parameters from assuming the wrong hierarchy and find that all biases on the parameters are below their 1-sigma marginalised errors.

We investigate the possibility of constraining the dark energy equation of state by measuring the ratio of Einstein radii in a strong gravitational lens system with two source planes. This quantity is independent of the Hubble parameter and directly measures the growth of angular diameter distances as a function of redshift. We investigate the prospects for a single double source plane system and for a forecast population of systems discovered by re-observing a population of single source lenses already known from a photometrically selected catalogue such as CASSOWARY or from a spectroscopically selected catalogue such as SLACS. We find that constraints comparable to current data-sets (15% uncertainty on the dark equation of state at 68%CL) are possible with a handful of double source plane systems. We also find that the method's degeneracy between Omega_M and w is almost orthogonal to that of CMB and BAO measurements, making this method highly complimentary to current probes.

Cosmological parameter constraints from the CMB power spectra alone suffer several well-known degeneracies. These degeneracies can be broken by numerical artefacts and also a variety of physical effects that become quantitatively important with high-accuracy data e.g. from the Planck satellite. We study degeneracies in models with flat and non-flat spatial sections, non-trivial dark energy and massive neutrinos, and investigate the importance of various physical degeneracy-breaking effects. We test the CAMB power spectrum code for numerical accuracy, and demonstrate that the numerical calculations are accurate enough for degeneracies to be broken mainly by true physical effects (the integrated Sachs-Wolfe effect, CMB lensing and geometrical and other effects through recombination) rather than numerical artefacts. We quantify the impact of CMB lensing on the power spectra, which inevitably provides degeneracy-breaking information even without using information in the non-Gaussianity. Finally we check the numerical accuracy of sample-based parameter constraints using CAMB and CosmoMC. In an appendix we document recent changes to CAMB's numerical treatment of massive neutrino perturbations, which are tested along with other recent improvements by our degeneracy exploration results.

We have undertaken SED modeling of discs around low mass T Tauri stars in order to delineate regions of the infrared two colour plane (K - [8] versus K - [24]) that correspond to discs in different evolutionary stages. This provides a ready tool for classifying the nature of star-disc systems based on infrared photometry. In particular we demonstrate the distinct loci followed by discs that undergo `uniform draining' (reduction in surface density by a spatially uniform factor) from those that clear from the inside out. We draw attention to the absence of objects on this `draining locus' in those star forming regions where the 24um sensitivity would permit their detection, as compared with the ~20 objects in these regions with colours suggestive of inner holes. We thus conclude that discs predominantly clear from the inside out. We also apply our classification of the infrared two colour plane to stars of spectral type M3-M5 in the IC 348 cluster and show that some of its members (dubbed `weak excess' sources by Muzerolle et al 2010) that had previously been assumed to be in a state of clearing are instead likely to be optically thick discs in which the dust is well settled towards the mid-plane. Nevertheless, there are many discs in a state of partial clearing in IC 348, with their abundance (relative to the total population of disc bearing stars) being around four times higher than for comparable stars in Taurus. However, the number of partially cleared discs relative to the total number of late type stars is similarly low in both regions (10 and 20 % respectively). We therefore conclude that IC 348 represents a more evolved version of the Taurus population (with more of its discs being highly settled or partially cleared) but that the timescale for clearing is similarly short (a few times 10^5 years) in both cases.

Markarian 501 is only the second extragalactic source to be detected with high statistical certainty at TeV energies; it is similar in many ways to Markarian 421. The Whipple Observatory gamma-ray telescope has been used to observe the AGN Markarian 501 in 1996 and 1997, the years subsequent to its initial detection. The apparent variability on the one-day time-scale observed in TeV gamma rays in 1995 is confirmed and compared with the variability in Markarian 421. Observations at X-ray and optical wavelengths from 1997 are also presented.

Markarian 421 was the first extragalactic source to be detected with high statistical certainty at TeV energies. The Whipple Observatory gamma-ray telescope has been used to observe the Active Galactic Nucleus, Markarian 421 in 1996 and 1997. The rapid variability observed in TeV gamma rays in previous years is confirmed. Doubling times as short as 15 minutes are reported with flux levels reaching 15 photons per minute. The TeV energy spectrum is derived using two independent methods. The implications for the intergalactic infra-red medium of an observed unbroken power law spectrum up to energies of 5 TeV is discussed.