We present a novel analysis technique for liquid xenon time projection chambers that allows for a lower threshold by relying on events with a prompt scintillation signal consisting of single detected photons. The energy threshold of the LUX dark matter experiment is primarily determined by the smallest scintillation response detectable, which previously required a 2-fold coincidence signal in its photomultiplier arrays, enforced in data analysis. The technique presented here exploits the double photoelectron emission effect observed in some photomultiplier models at vacuum ultraviolet wavelengths. We demonstrate this analysis using an electron recoil calibration dataset and place new constraints on the spin-independent scattering cross section of weakly interacting massive particles (WIMPs) down to 2.5 GeV/c$^2$ WIMP mass using the 2013 LUX dataset. This new technique is promising to enhance light WIMP and astrophysical neutrino searches in next-generation liquid xenon experiments.

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.).

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.

The PICASSO dark matter search experiment operated an array of 32 superheated droplet detectors containing 3.0 kg of C$_{4}$F$_{10}$ and collected an exposure of 231.4 kgd at SNOLAB between March 2012 and January 2014. We report on the final results of this experiment which includes for the first time the complete data set and improved analysis techniques including \mbox{acoustic} localization to allow fiducialization and removal of higher activity regions within the detectors. No signal consistent with dark matter was observed. We set limits for spin-dependent interactions on protons of $\sigma_p^{SD}$~=~1.32~$\times$~10$^{-2}$~pb (90\%~C.L.) at a WIMP mass of 20 GeV/c$^{2}$. In the spin-independent sector we exclude cross sections larger than $\sigma_p^{SI}$~=~4.86~$\times$~10$^{-5 }$~pb~(90\% C.L.) in the region around 7 GeV/c$^{2}$. The pioneering efforts of the PICASSO experiment have paved the way forward for a next generation detector incorporating much of this technology and experience into larger mass bubble chambers.

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%.

The DEAP-1 low-background liquid argon detector was used to measure scintillation pulse shapes of electron and nuclear recoil events and to demonstrate the feasibility of pulse-shape discrimination (PSD) down to an electron-equivalent energy of 20 keV. In the surface dataset using a triple-coincidence tag we found the fraction of beta events that are misidentified as nuclear recoils to be $<1.4\times 10^{-7}$ (90% C.L.) for energies between 43-86 keVee and for a nuclear recoil acceptance of at least 90%, with 4% systematic uncertainty on the absolute energy scale. The discrimination measurement on surface was limited by nuclear recoils induced by cosmic-ray generated neutrons. This was improved by moving the detector to the SNOLAB underground laboratory, where the reduced background rate allowed the same measurement with only a double-coincidence tag. The combined data set contains $1.23\times10^8$ events. One of those, in the underground data set, is in the nuclear-recoil region of interest. Taking into account the expected background of 0.48 events coming from random pileup, the resulting upper limit on the electronic recoil contamination is $<2.7\times10^{-8}$ (90% C.L.) between 44-89 keVee and for a nuclear recoil acceptance of at least 90%, with 6% systematic uncertainty on the absolute energy scale. We developed a general mathematical framework to describe PSD parameter distributions and used it to build an analytical model of the distributions observed in DEAP-1. Using this model, we project a misidentification fraction of approx. $10^{-10}$ for an electron-equivalent energy threshold of 15 keV for a detector with 8 PE/keVee light yield. This reduction enables a search for spin-independent scattering of WIMPs from 1000 kg of liquid argon with a WIMP-nucleon cross-section sensitivity of $10^{-46}$ cm$^2$, assuming negligible contribution from nuclear recoil backgrounds.

The Atacama Cosmology Telescope (ACT) is designed to make high angular resolution measurements of anisotropies in the Cosmic Microwave Background (CMB) at millimeter wavelengths. We describe ACTPol, an upgraded receiver for ACT, which uses feedhorn-coupled, polarization-sensitive detector arrays, a 3 degree field of view, 100 mK cryogenics with continuous cooling, and meta material anti-reflection coatings. ACTPol comprises three arrays with separate cryogenic optics: two arrays at a central frequency of 148 GHz and one array operating simultaneously at both 97 GHz and 148 GHz. The combined instrument sensitivity, angular resolution, and sky coverage are optimized for measuring angular power spectra, clusters via the thermal Sunyaev-Zel'dovich and kinetic Sunyaev-Zel'dovich signals, and CMB lensing due to large scale structure. The receiver was commissioned with its first 148 GHz array in 2013, observed with both 148 GHz arrays in 2014, and has recently completed its first full season of operations with the full suite of three arrays. This paper provides an overview of the design and initial performance of the receiver and related systems.

We present the design and the preliminary on sky performance with respect to beams and pass-bands of a multichroic polarimeter array covering the 90 and 146 GHz Cosmic Microwave Background (CMB) bands and its enabling broadband optical system recently deployed on the Atacama Cosmology Telescope (ACT). The constituent pixels are feedhorn-coupled multichroic polarimeters fabricated at NIST. This array is coupled to the ACT telescope via a set of three silicon lenses incorporating novel broad-band metamaterial anti-reflection coatings. This receiver represents the first multichroic detector array deployed for a CMB experiment and paves the way for the extensive use of multichroic detectors and broadband optical systems in the next generation of CMB experiments.

Advanced ACTPol is a polarization-sensitive upgrade for the 6 m aperture Atacama Cosmology Telescope (ACT), adding new frequencies and increasing sensitivity over the previous ACTPol receiver. In 2016, Advanced ACTPol will begin to map approximately half the sky in five frequency bands (28-230 GHz). Its maps of primary and secondary cosmic microwave background (CMB) anisotropies -- imaged in intensity and polarization at few arcminute-scale resolution -- will enable precision cosmological constraints and also a wide array of cross-correlation science that probes the expansion history of the universe and the growth of structure via gravitational collapse. To accomplish these scientific goals, the Advanced ACTPol receiver will be a significant upgrade to the ACTPol receiver, including four new multichroic arrays of cryogenic, feedhorn-coupled AlMn transition edge sensor (TES) polarimeters (fabricated on 150 mm diameter wafers); a system of continuously rotating meta-material silicon half-wave plates; and a new multiplexing readout architecture which uses superconducting quantum interference devices (SQUIDs) and time division to achieve a 64-row multiplexing factor. Here we present the status and scientific goals of the Advanced ACTPol instrument, emphasizing the design and implementation of the Advanced ACTPol cryogenic detector arrays.

We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons/second, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to: 1) construct kilometer scale phase screens necessary for wide-field telescopes, 2) reproduce atmospheric point-spread-function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, 3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and 4) recover system effective area predicted from analytic photometry integrals. This new code, the photon simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope (LSST) design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

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.

Many current and future dark matter and neutrino detectors are designed to measure scintillation light with a large array of photomultiplier tubes (PMTs). The energy resolution and particle identification capabilities of these detectors depend in part on the ability to accurately identify individual photoelectrons in PMT waveforms despite large variability in pulse amplitudes and pulse pileup. We describe a Bayesian technique that can identify the times of individual photoelectrons in a sampled PMT waveform without deconvolution, even when pileup is present. To demonstrate the technique, we apply it to the general problem of particle identification in single-phase liquid argon dark matter detectors. Using the output of the Bayesian photoelectron counting algorithm described in this paper, we construct several test statistics for rejection of backgrounds for dark matter searches in argon. Compared to simpler methods based on either observed charge or peak finding, the photoelectron counting technique improves both energy resolution and particle identification of low energy events in calibration data from the DEAP-1 detector and simulation of the larger MiniCLEAN dark matter detector.

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 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.

The direct search for dark matter is entering a period of increased sensitivity to the hypothetical Weakly Interacting Massive Particle (WIMP). One such technology that is being examined is a scintillation only noble liquid experiment, MiniCLEAN. MiniCLEAN utilizes over 500 kg of liquid cryogen to detect nuclear recoils from WIMP dark matter and serves as a demonstration for a future detector of order 50 to 100 tonnes. The liquid cryogen is interchangeable between argon and neon to study the A$^{2}$ dependence of the potential signal and examine backgrounds. MiniCLEAN utilizes a unique modular design with spherical geometry to maximize the light yield using cold photomultiplier tubes in a single-phase detector. Pulse shape discrimination techniques are used to separate nuclear recoil signals from electron recoil backgrounds. MiniCLEAN will be spiked with additional $^{39}$Ar to demonstrate the effective reach of the pulse shape discrimination capability. Assembly of the experiment is underway at SNOLAB and an update on the project is given.

The statistics of peak counts in reconstructed shear maps contain information beyond the power spectrum, and can improve cosmological constraints from measurements of the power spectrum alone if systematic errors can be controlled. We study the effect of galaxy shape measurement errors on predicted cosmological constraints from the statistics of shear peak counts with the Large Synoptic Survey Telescope (LSST). We use the LSST image simulator in combination with cosmological N-body simulations to model realistic shear maps for different cosmological models. We include both galaxy shape noise and, for the first time, measurement errors on galaxy shapes. We find that the measurement errors considered have relatively little impact on the constraining power of shear peak counts for LSST.

A main science goal for the Large Synoptic Survey Telescope (LSST) is to measure the cosmic shear signal from weak lensing to extreme accuracy. One difficulty, however, is that with the short exposure time ($\simeq$15 seconds) proposed, the spatial variation of the Point Spread Function (PSF) shapes may be dominated by the atmosphere, in addition to optics errors. While optics errors mainly cause the PSF to vary on angular scales similar or larger than a single CCD sensor, the atmosphere generates stochastic structures on a wide range of angular scales. It thus becomes a challenge to infer the multi-scale, complex atmospheric PSF patterns by interpolating the sparsely sampled stars in the field. In this paper we present a new method, PSFent, for interpolating the PSF shape parameters, based on reconstructing underlying shape parameter maps with a multi-scale maximum entropy algorithm. We demonstrate, using images from the LSST Photon Simulator, the performance of our approach relative to a 5th-order polynomial fit (representing the current standard) and a simple boxcar smoothing technique. Quantitatively, PSFent predicts more accurate PSF models in all scenarios and the residual PSF errors are spatially less correlated. This improvement in PSF interpolation leads to a factor of 3.5 lower systematic errors in the shear power spectrum on scales smaller than $\sim13'$, compared to polynomial fitting. We estimate that with PSFent and for stellar densities greater than $\simeq1/{\rm arcmin}^{2}$, the spurious shear correlation from PSF interpolation, after combining a complete 10-year dataset from LSST, is lower than the corresponding statistical uncertainties on the cosmic shear power spectrum, even under a conservative scenario.

The complete 10-year survey from the Large Synoptic Survey Telescope (LSST) will image $\sim$ 20,000 square degrees of sky in six filter bands every few nights, bringing the final survey depth to $r\sim27.5$, with over 4 billion well measured galaxies. To take full advantage of this unprecedented statistical power, the systematic errors associated with weak lensing measurements need to be controlled to a level similar to the statistical errors. This work is the first attempt to quantitatively estimate the absolute level and statistical properties of the systematic errors on weak lensing shear measurements due to the most important physical effects in the LSST system via high fidelity ray-tracing simulations. We identify and isolate the different sources of algorithm-independent, \textit{additive} systematic errors on shear measurements for LSST and predict their impact on the final cosmic shear measurements using conventional weak lensing analysis techniques. We find that the main source of the errors comes from an inability to adequately characterise the atmospheric point spread function (PSF) due to its high frequency spatial variation on angular scales smaller than $\sim10'$ in the single short exposures, which propagates into a spurious shear correlation function at the $10^{-4}$--$10^{-3}$ level on these scales. With the large multi-epoch dataset that will be acquired by LSST, the stochastic errors average out, bringing the final spurious shear correlation function to a level very close to the statistical errors. Our results imply that the cosmological constraints from LSST will not be severely limited by these algorithm-independent, additive systematic effects.

New data are reported from the operation of a 4.0 kg CF$_{3}$I bubble chamber in the 6800-foot-deep SNOLAB underground laboratory. The effectiveness of ultrasound analysis in discriminating alpha-decay background events from single nuclear recoils has been confirmed, with a lower bound of $>$99.3% rejection of alpha-decay events. Twenty single nuclear recoil event candidates and three multiple bubble events were observed during a total exposure of 553 kg-days distributed over three different bubble nucleation thresholds. The effective exposure for single bubble recoil-like events was 437.4 kg-days. A neutron background internal to the apparatus, of known origin, is estimated to account for five single nuclear recoil events and is consistent with the observed rate of multiple bubble events. This observation provides world best direct detection constraints on WIMP-proton spin-dependent scattering for WIMP masses $>$20 GeV/c$^{2}$ and demonstrates significant sensitivity for spin-independent interactions.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.