We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations of the first $10^4$ orbits, thus achieving speed-ups of up to $10^5$ over full simulations. This computationally opens up the stability constrained characterization of multi-planet systems. Our model, trained on $\approx 100,000$ three-planet systems sampled at discrete resonances, generalizes both to a sample spanning a continuous period-ratio range, as well as to a large five-planet sample with qualitatively different configurations to our training dataset. Our approach significantly outperforms previous methods based on systems' angular momentum deficit, chaos indicators, and parametrized fits to numerical integrations. We use SPOCK to constrain the free eccentricities between the inner and outer pairs of planets in the Kepler-431 system of three approximately Earth-sized planets to both be below 0.05. Our stability analysis provides significantly stronger eccentricity constraints than currently achievable through either radial velocity or transit duration measurements for small planets, and within a factor of a few of systems that exhibit transit timing variations (TTVs). Given that current exoplanet detection strategies now rarely allow for strong TTV constraints (Hadden et al., 2019), SPOCK enables a powerful complementary method for precisely characterizing compact multi-planet systems. We publicly release SPOCK for community use.

Any successful alternative gravity theory that obviates the need for dark matter must fit our cosmological observations. Measurements of microwave background polarization trace the large-scale baryon velocity field at recombination and show very strong, $O(1)$, baryon acoustic oscillations. Measurements of the large-scale structure of galaxies at low redshift show much weaker features in the spectrum. If the alternative gravity theory's dynamical equations for the growth rate of structure are linear, then the density field growth can be described by a Green's function: $\delta(\vec x,t) = \delta(\vec x,t')G(x,t,t')$. We show that the Green function, $G(x,t,t')$, must have dramatic features that erase the initial baryon oscillations. This implies an acceleration law that changes sign on the $\sim 150$ Mpc scale. On the other hand, if the alternative gravity theory has a large nonlinear term that couples modes on different scales, then the theory would predict large-scale non-Gaussian features in large-scale structure. These are not seen in the distribution of galaxies nor in the distribution of quasars. No proposed alternative gravity theory for dark matter seems to satisfy these constraints.

We develop a general approach to distill symbolic representations of a learned deep model by introducing strong inductive biases. We focus on Graph Neural Networks (GNNs). The technique works as follows: we first encourage sparse latent representations when we train a GNN in a supervised setting, then we apply symbolic regression to components of the learned model to extract explicit physical relations. We find the correct known equations, including force laws and Hamiltonians, can be extracted from the neural network. We then apply our method to a non-trivial cosmology example-a detailed dark matter simulation-and discover a new analytic formula which can predict the concentration of dark matter from the mass distribution of nearby cosmic structures. The symbolic expressions extracted from the GNN using our technique also generalized to out-of-distribution data better than the GNN itself. Our approach offers alternative directions for interpreting neural networks and discovering novel physical principles from the representations they learn.

The marked power spectrum is capable of placing far tighter constraints on cosmological parameters (particularly the neutrino mass) than the conventional power spectrum. What new information does it contain beyond conventional statistics? Through the development of a perturbative model, we find that the mark induces a significant coupling between small-scale non-Gaussianities and large scales, leading to the additional information content. The model is derived in the context of one-loop perturbation theory and validated by comparison to $N$-body simulations across a variety of mark parameters. At moderate redshifts, including for massive neutrino cosmologies, the theory is in good agreement with the simulations. The mixing of large and small scales complicates the modeling as there is no well-defined convergence radius of the theory at low $z$. Extension to higher perturbative order and biased tracers is possible via a similar approach, and a simple model of the latter is shown to yield promising results. The theory becomes non-perturbative at redshift zero for small smoothing scales. This suggests that marked spectrum may have important contributions from baryonic effects even at low $k$: these will need to be studied before the full power of this tool can be realized.

We introduce a physically-motivated model of the matter power spectrum, based on the halo model and perturbation theory. This model achieves 1\% accuracy on all $k-$scales between $k=0.02h\,\mathrm{Mpc}^{-1}$ to $k=1h\,\mathrm{Mpc}^{-1}$. Our key ansatz is that the number density of halos depends on the non-linear density contrast filtered on some unknown scale $R$. Using the Effective Field Theory of Large Scale Structure to evaluate the two-halo term, we obtain a model for the power spectrum with only two fitting parameters: $R$ and the effective `sound speed', which encapsulates small-scale physics. This is tested with two suites of cosmological simulations across a broad range of cosmologies and found to be highly accurate. Due to its physical motivation, the statistics can be easily extended beyond the power spectrum; we additionally derive the one-loop covariance matrices of cluster counts and their combination with the matter power spectrum. This yields a significantly better fit to simulations than previous models, and includes a new model for super-sample effects, which is rigorously tested with separate universe simulations. At low redshift, we find a significant ($\sim 10\%$) exclusion covariance from accounting for the finite size of halos which has not previously been modeled. Such power spectrum and covariance models will enable joint analysis of upcoming large-scale structure surveys, gravitational lensing surveys and cosmic microwave background maps on scales down to the non-linear scale. We provide a publicly released Python code.

We present a new method for joint cosmological parameter inference and cluster mass calibration from a combination of weak lensing measurements and the abundance of thermal Sunyaev-Zel'dovich (tSZ) selected galaxy clusters. We combine cluster counts with the spherical harmonic cosmic shear power spectrum and the cross-correlation between cluster overdensity and cosmic shear. These correlations constrain the cluster mass-observable relation. We model the observables using a halo model framework, including their full non-Gaussian covariance. Forecasting constraints on cosmological and mass calibration parameters for a combination of LSST cosmic shear and Simons Observatory tSZ cluster counts, we find competitive constraints for cluster cosmology, with a factor two improvement in the dark energy figure of merit compared to LSST cosmic shear alone. We find most of the mass calibration information will be in the large and intermediate scales of the cross-correlation between cluster overdensity and cosmic shear. Finally, we find broadly comparable constraints to traditional analyses based on calibrating masses using stacked cluster lensing measurements, with the benefit of consistently accounting for the correlations with cosmic shear.

We present new constraints on anisotropic birefringence of the cosmic microwave background polarization using two seasons of data from the Atacama Cosmology Telescope covering $456$ square degrees of sky. The birefringence power spectrum, measured using a curved-sky quadratic estimator, is consistent with zero. Our results provide the tightest current constraint on birefringence over a range of angular scales between $5$ arcminutes and $9$ degrees. We improve previous upper limits on the amplitude of a scale-invariant birefringence power spectrum by a factor of between $2$ and $3$. Assuming a nearly-massless axion field during inflation, our result is equivalent to a $2\,\sigma$ upper limit on the Chern-Simons coupling constant between axions and photons of $g_{\alpha\gamma}<4.0\times 10^{-2}/H_I$ where $H_I$ is the inflationary Hubble scale.

We explore synergies between the space-based Wide-Field Infrared Survey Telescope (WFIRST) and the ground-based Rubin Observatory Legacy Survey of Space and Time (LSST). In particular, we consider a scenario where the currently envisioned survey strategy for WFIRST's High Latitude Survey (HLS), i.e., 2000 square degrees in four narrow photometric bands is altered in favor of a strategy that combines rapid coverage of the LSST area (to full LSST depth) in one band. We find that a 5-month WFIRST survey in the W-band can cover the full LSST survey area providing high-resolution imaging for >95% of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where WFIRST spends 1.5 years covering the LSST area. For this second scenario we quantify the constraining power on dark energy equation of state parameters from a joint weak lensing and galaxy clustering analysis, and compare it to an LSST-only survey and to the Reference WFIRST HLS survey. Our survey simulations are based on the WFIRST exposure time calculator and redshift distributions from the CANDELS catalog. Our statistical uncertainties account for higher-order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modeling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. Assuming the 5-month WFIRST wide scenario, we find a significant increase in constraining power for the joint LSST+WFIRST wide survey compared to LSST Y10 (FoM(Wwide)= 2.4 FoM(LSST)) and compared to LSST+WFIRST HLS (FoM(Wwide)= 5.5 FoM(HLS)).

We simulate the scientific performance of the Wide-Field Infrared Survey Telescope (WFIRST) High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6 year HLS Reference survey is currently envisioned to image 2000 deg$^2$ in multiple bands to a depth of $\sim$26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z=3. The combination of deep, multi-band photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g., weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper we explore multi-probe strategies that can be implemented given WFIRST's instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the WFIRST survey (photo-z and shear calibration) and ultimately run a joint likelihood analysis in N-dim parameter space. We find that the WFIRST reference survey alone (no external data sets) can achieve a standard dark energy FoM of >300 when including all probes. This assumes no information from external data sets and realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community driven effort to simulate and optimize the science return of WFIRST.

We construct cosmic microwave background lensing mass maps using data from the 2014 and 2015 seasons of observations with the Atacama Cosmology Telescope (ACT). These maps cover 2100 square degrees of sky and overlap with a wide variety of optical surveys. The maps are signal dominated on large scales and have fidelity such that their correlation with the cosmic infrared background is clearly visible by eye. We also create lensing maps with thermal Sunyaev-Zel'dovich contamination removed using a novel cleaning procedure that only slightly degrades the lensing signal-to-noise ratio. The cross-spectrum between the cleaned lensing map and the BOSS CMASS galaxy sample is detected at $10$-$\sigma$ significance, with an amplitude of $A=1.02 \pm 0.10$ relative to the Planck best-fit LCDM cosmological model with fiducial linear galaxy bias. Our measurement lays the foundation for lensing cross-correlation science with current ACT data and beyond.

The cross-correlation between fluctuations in the electron scattering optical depth $\tau_{\rm es}$ as probed by future Cosmic Microwave Background (CMB) experiments, and fluctuations in the 21cm differential brightness temperature $\Delta T_{\rm 21cm}$ as probed by ground-based radio interferometers, will trace the reionization history of the Universe. In particular, the $\tau_{\rm es}-$21cm cross-correlation should yield a determination of the characteristic bubble size distribution and ionization fraction as a function of redshift. When assuming that the cross-correlation signal is limited by instrumental noise rather than by foregrounds, we estimate its potential detectability by upcoming experiments. Specifically, the combination of HERA and Simons Observatory, CMB-S4 and PICO should yield a signal-to-noise ratio around 3 - 6, while and the exploitation of the SKA should increase it to 10-20. Finally, we have discussed how such levels of detectability can be affected when (simply modeled) 21cm foregrounds are present. For the most promising PICO$\times$SKA configuration, an efficiency of foreground removal to a level of $7\times 10^{-4}$ is needed to achieve a $5\sigma$ detection of the cross-correlation signal; in addition, safe avoidance of foreground contamination in the line-of-sight Fourier modes above $0.03 \,h\,\rm Mpc^{-1}$ would guarantee a detection significance around $3\sigma$.

The density profiles of dwarf galaxies are a highly varied set. If the dark matter is an ultra-light particle such as axions, then simulations predict a distinctive and unique profile. If the axion mass is large enough to fit the ultra-faint dwarf (UFD) satellites($m\gtrapprox 10^{-21}$ eV), then the models do not fit the density profile of Fornax and Sculptor and are ruled out by more than $3-\sigma$ confidence. If the axion mass is in the mass range that can fit mass profiles of Fornax and Sculptor dwarf spheroidals, then its extended profile implies enormous masses ($\approx10^{11}-10^{12}M_{\odot}$) for the UFDs. These large masses for the UFDS are ruled out by more than $3-\sigma$ confidence by dynamical friction arguments. The tension would increase further considering star formation histories and stellar masses of the UFDs. Unless future ultra-light dark matter (ULDM) simulations with baryonic feedback show a significant change in the density structure of the halos, the current data is incompatible with the ULDM scenario. Relaxing the slope constraint from classical dwarf galaxies would lead to excluding ULDM with mass less than $6\times10^{-22}$ eV.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

We propose using the upcoming Cosmic Microwave Background (CMB) ground based experiments to detect the signal of ALPs (Axion like particles) interacting with magnetic fields in galaxy clusters. The conversion between CMB photons and ALPs in the presence of the cluster magnetic field can cause a polarized spectral distortion in the CMB around a galaxy cluster. The strength of the signal depends upon the redshift of the galaxy cluster and will exhibit a distinctive spatial profile around it depending upon the structure of electron density and magnetic field. This distortion produces a different shape from the other known spectral distortions like $y$-type and $\mu$-type and hence are separable from the multi-frequency CMB observation. The spectrum is close to kinematic Sunyaev-Zeldovich (kSZ) signal but can be separated from it using the polarization information. For the future ground-based CMB experiments such as Simons Observatory and CMB-S4, we estimate the measurability of this signal in the presence of foreground contamination, instrument noise and CMB anisotropies. This new avenue can probe the photon-ALP coupling over the ALP mass range from $10^{-13}$ eV to $10^{-12}$ eV with two orders of magnitude better accuracy from CMB-S4 than the current existing bounds.

Cosmological neutrinos have their greatest influence in voids: these are the regions with the highest neutrino to dark matter density ratios. The marked power spectrum can be used to emphasize low density regions over high density regions, and therefore is potentially much more sensitive than the power spectrum to the effects of neutrino masses. Using 22,000 N-body simulations from the Quijote suite, we quantify the information content in the marked power spectrum of the matter field, and show that it outperforms the standard power spectrum by setting constraints improved by a factor larger than 2 on all cosmological parameters. The combination of marked and standard power spectrum allows to place a 4.3{\sigma} constraint on the minimum sum of the neutrino masses with a volume equal to 1 (Gpc/h)^3 and without CMB priors. Combinations of different marked power spectra yield a 6{\sigma} constraint within the same conditions.

We present an in-depth exploration of the phenomenon of dynamical friction in a universe where the dark matter is composed entirely of so-called Fuzzy Dark Matter (FDM), ultralight bosons of mass $m\sim\mathcal{O}(10^{-22})\,$eV. We review the classical treatment of dynamical friction before presenting analytic results in the case of FDM for point masses, extended mass distributions, and FDM backgrounds with finite velocity dispersion. We then test these results against a large suite of fully non-linear simulations that allow us to assess the regime of applicability of the analytic results. We apply these results to a variety of astrophysical problems of interest, including infalling satellites in a galactic dark matter background, and determine that \emph{(1)}~for FDM masses $m\gtrsim 10^{-21}\, {\rm eV}\, c^{-2}$, the timing problem of the Fornax dwarf spheroidal's globular clusters is no longer solved and \emph{(2)}~the effects of FDM on the process of dynamical friction for satellites of total mass $M$ and relative velocity $v_{\rm rel}$ should require detailed numerical simulations for $\left(M/10^9~M_{\odot}\right) \left(m/10^{-22}~{\rm eV}\right)\left(100~{\rm km}~{\rm s}^{-1}/v_{\rm rel}\right) \sim 1$, parameters which would lie outside the validated range of applicability of any currently developed analytic theory, due to transient wave structures in the time-dependent regime.

Optimal analyses of many signals in the cosmic microwave background (CMB) require map-level extraction of individual components in the microwave sky, rather than measurements at the power spectrum level alone. To date, nearly all map-level component separation in CMB analyses has been performed exclusively using satellite data. In this paper, we implement a component separation method based on the internal linear combination (ILC) approach which we have designed to optimally account for the anisotropic noise (in the 2D Fourier domain) often found in ground-based CMB experiments. Using this method, we combine multi-frequency data from the Planck satellite and the Atacama Cosmology Telescope Polarimeter (ACTPol) to construct the first wide-area, arcminute-resolution component-separated maps (covering approximately 2100 sq. deg.) of the CMB temperature anisotropy and the thermal Sunyaev-Zel'dovich (tSZ) effect sourced by the inverse-Compton scattering of CMB photons off hot, ionized gas. Our ILC pipeline allows for explicit deprojection of various contaminating signals, including a modified blackbody approximation of the cosmic infrared background (CIB) spectral energy distribution. The cleaned CMB maps will be a useful resource for CMB lensing reconstruction, kinematic SZ cross-correlations, and primordial non-Gaussianity studies. The tSZ maps will be used to study the pressure profiles of galaxies, groups, and clusters through cross-correlations with halo catalogs, with dust contamination controlled via CIB deprojection. The data products described in this paper will be made available on LAMBDA.

We make use of the IllustrisTNG cosmological, hydrodynamical simulations to test fundamental assumptions of the mass-based Halo Occupation Distribution (HOD) approach to modelling the galaxy-halo connection. By comparing the clustering of galaxies measured in the 300 Mpc TNG box (TNG300) with that predicted by the standard ("basic") HOD model (bHOD), we find that, on average, the bHOD model underpredicts the real-space correlation function in the TNG300 box by $\sim$ 15\% on scales of $1 \ {\rm Mpc}/h < r < 20 \ {\rm Mpc}/h$, which is well beyond the target precision demanded of next-generation galaxy redshift surveys. We perform several tests to establish the robustness of our findings to systematic effects, including the effect of finite box size and the choice of halo finder. In our exploration of "secondary" parameters with which to augment the bHOD, we find that the local environment of the halo, the velocity dispersion anisotropy, $\beta$, and the product of the half-mass radius and the velocity dispersion, $\sigma^2 R_{\rm halfmass}$, are the three most effective measures of assembly bias that help reconcile the bHOD-predicted clustering with that in TNG300. In addition, we test other halo properties such as halo spin, formation epoch and halo concentration. We also find that at fixed halo mass, galaxies in one type of environment cluster differently from galaxies in another. We demonstrate that a more complete model of the galaxy-halo connection can be constructed if both the mass and information regarding the local environment in which a halo is embedded are combined.

Cosmic voids, the underdense regions in the universe, are particularly sensitive to diffuse density components such as cosmic neutrinos. This sensitivity is enhanced by the match between void sizes and the free-streaming scale of massive neutrinos. Using the massive neutrino simulations \texttt{MassiveNuS}, we investigate the effect of neutrino mass on dark matter halos as a function of environment. We find that the halo mass function depends strongly on neutrino mass and that this dependence is more pronounced in voids than in high-density environments. An observational program that measured the characteristic mass of the most massive halos in voids should be able to place novel constraints on the sum of the masses of neutrinos $\sum m_\nu$. The neutrino mass effect in the simulations is quite strong: In a 512$^3$ $h^{-3}$ Mpc$^3$ survey, the mean mass of the 1000 most massive halos in the void interiors is $(4.82 \pm 0.11) \times 10^{12} h^{-1}M_{\odot}$ for $\sum m_\nu = 0.6$ eV and $(8.21 \pm 0.13) \times 10^{12} h^{-1}M_{\odot}$ for $\sum m_\nu = 0.1$ eV. Subaru (SuMIRe), Euclid and WFIRST will have both spectroscopic and weak lensing surveys. Covering volumes at least 50 times larger than our simulations, they should be sensitive probes of neutrino mass through void substructure.

Understanding the properties of the galactic emission at millimetre wavelengths is important for studies of the cosmic microwave background (CMB). In this work we explore the bispectrum, the harmonic equivalent of the three point function, from galactic dust and synchrotron emission. We investigate these effects across a broad range of frequencies using the synchrotron dominated S-band Polarization All Sky Survey (SPASS) maps at 2.3 GHz, the Planck satellite maps (30-857 GHz) and dust dominated Infrared Astronomical Satellite (IRAS) maps at 3 THz. We measure bispectra of total intensity fields, T, as well as the gradient, E, and curl modes, B, of the polarization field. We find that the synchrotron and galactic dust have strong temperature bispectra with significant contributions in the squeezed limit, which probes the correlations between two small scale modes with a large scale mode. Additionally, we find that the dust also has strong polarised bispectra that also peak in the squeezed configuration. We explore parity odd bispectra, such as BTT bispectra, and find strong parity odd bispectra for the galactic dust notably in BTT, BTE and BEE configurations. After masking bright sources, we find no evidence for polarised synchrotron bispectra and no evidence for cross bispectra between the dust and synchrotron emission. The strong foreground bispectra discussed here need to be carefully controlled to avoid biasing measurements of primordial non-Gaussianity. Finally we use these bispectra tools to test for residual foregrounds in the component separated Planck maps and find no evidence of residual foregrounds. These tools will be useful for characterizing residual foregrounds in component separated maps, particularly for experiments with less frequency coverage than the Planck satellite.

Improved distance measurements to millisecond pulsars can enhance pulsar timing array (PTA) sensitivity to gravitational-waves, improve tests of general relativity with binary pulsars, improve constraints on fuzzy dark matter, and more. Here we report the parallax and distance measurements to seven Gaia DR2 objects associated with 7 International PTA pulsars: J0437-4715, J1012+5307, J1024-0719, J1732-5049, J1910+1256, J1843-1113, and J1949+3106. By multiplying the posteriors of the PTA-based parallax measurements with the Gaia parallax measurement to the pulsar's companion, we improve the distance measurements, and provide a tentative detection of a previously unknown binary companion to J1843-1113. We also find an order of magnitude improvement in the parallax measurement to J1949+3106. Finally, we recommend that future Gaia data releases use J0437-4715 as a calibration point, since its distance estimate in Gaia DR2 is relatively poor compared to pulsar timing measurements.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

We demonstrate an algorithm for learning a flexible color-magnitude diagram from noisy parallax and photometry measurements using a normalizing flow, a deep neural network capable of learning an arbitrary multi-dimensional probability distribution. We present a catalog of 640M photometric distance posteriors to nearby stars derived from this data-driven model using Gaia DR2 photometry and parallaxes. Dust estimation and dereddening is done iteratively inside the model and without prior distance information, using the Bayestar map. The signal-to-noise (precision) of distance measurements improves on average by more than 48% over the raw Gaia data, and we also demonstrate how the accuracy of distances have improved over other models, especially in the noisy-parallax regime. Applications are discussed, including significantly improved Milky Way disk separation and substructure detection. We conclude with a discussion of future work, which exploits the normalizing flow architecture to allow us to exactly marginalize over missing photometry, enabling the inclusion of many surveys without losing coverage.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

Do void statistics contain information beyond the tracer 2-point correlation function? Yes! As we vary the sum of the neutrino masses, we find void statistics contain information absent when using just tracer 2-point statistics. Massive neutrinos uniquely affect cosmic voids. We explore their impact on void clustering using both the DEMNUni and MassiveNuS simulations. For voids, neutrino effects depend on the observed void tracers. As the neutrino mass increases, the number of small voids traced by cold dark matter particles increases and the number of large voids decreases. Surprisingly, when massive, highly biased, halos are used as tracers, we find the opposite effect. The scale at which voids cluster, as well as the void correlation, is similarly sensitive to the sum of neutrino masses and the tracers. This scale dependent trend is not due to simulation volume or halo density. The interplay of these signatures in the void abundance and clustering leaves a distinct fingerprint that could be detected with observations and potentially help break degeneracies between different cosmological parameters. This paper paves the way to exploit cosmic voids in future surveys to constrain the mass of neutrinos.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

The Wide Field Infrared Survey Telescope (WFIRST) will monitor $\sim 2$ deg$^2$ toward the Galactic bulge in a wide ($\sim 1-2~\mu$m) W149 filter at 15-minute cadence with exposure times of $\sim$50s for 6 seasons of 72 days each, for a total $\sim$41,000 exposures taken over $\sim$432 days, spread over the 5-year prime mission. This will be one of the deepest exposures of the sky ever taken, reaching a photon-noise photometric precision of 0.01 mag per exposure and collecting a total of $\sim 10^9$ photons over the course of the survey for a W149$_{\rm AB}\sim 21$ star. Of order $4 \times 10^7$ stars will be monitored with W149$_{\rm AB}$<21, and 10$^8$ stars with W145$_{\rm AB}$<23. The WFIRST microlensing survey will detect $\sim$54,000 microlensing events, of which roughly 1% ($\sim$500) will be due to isolated black holes, and $\sim$3% ($\sim$1600) will be due to isolated neutron stars. It will be sensitive to (effectively) isolated compact objects with masses as low as the mass of Pluto, thereby enabling a measurement of the compact object mass function over 10 orders of magnitude. Assuming photon-noise limited precision, it will detect $\sim 10^5$ transiting planets with sizes as small as $\sim 2~R_\oplus$, perform asteroseismology of $\sim 10^6$ giant stars, measure the proper motions to $\sim 0.3\%$ and parallaxes to $\sim 10\%$ for the $\sim 6 \times 10^6$ disk and bulge stars in the survey area, and directly detect $\sim 5 \times 10^3$ Trans-Neptunian objects (TNOs) with diameters down to $\sim 10$ km, as well as detect $\sim 10^3$ occulations of stars by TNOs during the survey. All of this science will completely serendipitous, i.e., it will not require modifications of the WFIRST optimal microlensing survey design. Allowing for some minor deviation from the optimal design, such as monitoring the Galactic center, would enable an even broader range of transformational science.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

The nature of dark energy remains a profound mystery 20 years after the discovery of cosmic acceleration. A very high number density galaxy redshift survey over a wide area (HD GRS Wide) spanning the redshift range of 0.5<z<4 using the same tracer, carried out using massively parallel wide field multi-object slit spectroscopy from space, will provide definitive dark energy measurements with minimal observational systematics by design. The HD GRS Wide will illuminate the nature of dark energy, and lead to revolutionary advances in particle physics and cosmology. It will also trace the cosmic web of dark matter and provide key insight into large-scale structure in the Universe. The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

Do cool-core (CC) and noncool-core (NCC) clusters live in different environments? We make novel use of H$\alpha$ emission lines in the central galaxies of redMaPPer clusters as proxies to construct large (1,000's) samples of CC and NCC clusters, and measure their relative assembly bias using both clustering and weak lensing. We increase the statistical significance of the bias measurements from clustering by cross-correlating the clusters with an external galaxy redshift catalog from the Sloan Digital Sky Survey III, the LOWZ sample. Our cross-correlations can constrain assembly bias up to a statistical uncertainty of 6%. Given our H$\alpha$ criteria for CC and NCC, we find no significant differences in their clustering amplitude. Interpreting this difference as the absence of halo assembly bias, our results rule out the possibility of having different large-scale (tens of Mpc) environments as the source of diversity observed in cluster cores. Combined with recent observations of the overall mild evolution of CC and NCC properties, such as central density and CC fraction, this would suggest that either the cooling properties of the cluster core are determined early on solely by the local (<200 kpc) gas properties at formation or that local merging leads to stochastic CC relaxation and disruption in a periodic way, preserving the average population properties over time. Studying the small-scale clustering in clusters at high redshift would help shed light on the exact scenario.

Astrophysical transients have been observed for millennia and have shaped our most basic assumptions about the Universe. In the last century, systematic searches have grown from detecting handfuls of transients per year to over 7000 in 2018 alone. As these searches have matured, we have discovered both large samples of "normal" classes and new, rare classes. Recently, a transient was the first object observed in both gravitational waves and light. Ground-based observatories, including LSST, will discover thousands of transients in the optical, but these facilities will not provide the high-fidelity near-infrared (NIR) photometry and high-resolution imaging of a space-based observatory. WFIRST can fill this gap. With its survey designed to measure the expansion history of the Universe with Type Ia supernovae, WFIRST will also discover and monitor thousands of other transients in the NIR, revealing the physics for these high-energy events. Small-scale GO programs, either as a supplement to the planned survey or as specific target-of-opportunity observations, would significantly expand the scope of transient science that can be studied with WFIRST.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

Most secondary sources of cosmic microwave background anisotropy (radio sources, dusty galaxies -DSFGs-, thermal Sunyaev Zel'dovich -tSZ- distortions from hot gas, and gravitational lensing) are highly non-Gaussian. Statistics beyond the power spectrum are therefore potentially important sources of information about the physics of these processes. We combine data from the Atacama Cosmology Telescope and with data from the Planck satellite (only using Planck data in the overlapping region) to constrain the amplitudes of a set of theoretical bispectrum templates from the tSZ effect, DSFGs, gravitational lensing, and radio galaxies. We make a strong detection of radio galaxies (>5$\sigma$) and have hints of non-Gaussianity arising from the tSZ effect, DSFGs, from cross-correlations between the tSZ effect and DSFGs and from cross-correlations among the tSZ effect, DSFGs and radio galaxies. These results suggest that the same halos host radio sources, DSFGs, and have tSZ signal. We present a new method to calculate the non-Gaussian contributions to the template covariances. Using this method we find significant non-Gaussian contributions to the variance and covariance of our templates, with templates involving the tSZ effect most effected. Strong degeneracies exist between the various sources at the current noise levels. In light of these degeneracies, combined with theoretical uncertainty in the templates, these results are a demonstration of this technique. With these caveats, we demonstrate the utility of future bispectrum measurements by using the tSZ bispectrum measurement to constrain a combination of the amplitude of matter fluctuations and the matter density to be $\sigma_8 \Omega_m^{0.17}=0.65^{+0.05}_{-0.06}$. Improvements in signal to noise from upcoming observations will enable the separation of bispectrum components and robust constraints on cosmological parameters.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

We measure cosmic weak lensing shear power spectra with the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalog covering 137deg$^2$ of the sky. Thanks to the high effective galaxy number density of $\sim$17 arcmin$^{-2}$ even after conservative cuts such as magnitude cut of $i<24.5$ and photometric redshift cut of $0.3\leq z \leq 1.5$, we obtain a high significance measurement of the cosmic shear power spectra in 4 tomographic redshift bins, achieving a total signal-to-noise ratio of 16 in the multipole range $300 \leq \ell \leq 1900$. We carefully account for various uncertainties in our analysis including the intrinsic alignment of galaxies, scatters and biases in photometric redshifts, residual uncertainties in the shear measurement, and modeling of the matter power spectrum. The accuracy of our power spectrum measurement method as well as our analytic model of the covariance matrix are tested against realistic mock shear catalogs. For a flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, we find $S_8\equiv \sigma_8(\Omega_{\rm m}/0.3)^\alpha=0.800^{+0.029}_{-0.028}$ for $\alpha=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $\alpha=0.5$) from our HSC tomographic cosmic shear analysis alone. In comparison with Planck cosmic microwave background constraints, our results prefer slightly lower values of $S_8$, although metrics such as the Bayesian evidence ratio test do not show significant evidence for discordance between these results. We study the effect of possible additional systematic errors that are unaccounted in our fiducial cosmic shear analysis, and find that they can shift the best-fit values of $S_8$ by up to $\sim 0.6\sigma$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) is a high-contrast imager and integral field spectrograph that will enable the study of exoplanets and circumstellar disks at visible wavelengths. Ground-based high-contrast instrumentation has fundamentally limited performance at small working angles, even under optimistic assumptions for 30m-class telescopes. There is a strong scientific driver for better performance, particularly at visible wavelengths. Future flagship mission concepts aim to image Earth analogues with visible light flux ratios of more than 10^10. CGI is a critical intermediate step toward that goal, with a predicted 10^8-9 flux ratio capability in the visible. CGI achieves this through improvements over current ground and space systems in several areas: (i) Hardware: space-qualified (TRL9) deformable mirrors, detectors, and coronagraphs, (ii) Algorithms: wavefront sensing and control; post-processing of integral field spectrograph, polarimetric, and extended object data, and (iii) Validation of telescope and instrument models at high accuracy and precision. This white paper, submitted to the 2018 NAS Exoplanet Science Strategy call, describes the status of key CGI technologies and presents ways in which performance is likely to evolve as the CGI design matures.

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science.

We develop a unified model of darkgenesis and baryogenesis involving strongly interacting dark quarks, utilizing the gravitational anomaly of chiral gauge theories. In these models, both the visible and dark baryon asymmetries are generated by the gravitational anomaly induced by the presence of chiral primordial gravitational waves. We provide a concrete model of an SU(2) gauge theory with two massless quarks. In this model, the dark quarks condense and form a dark baryon charge superfluid (DBS), in which the Higgs-mode acts as cold dark matter. We elucidate the essential features of this dark matter scenario and discuss its phenomenological prospects.

We explore the effect of massive neutrinos on the weak lensing shear bispectrum using the Cosmological Massive Neutrino Simulations. We find that the primary effect of massive neutrinos is to suppress the amplitude of the bispectrum with limited effect on the bispectrum shape. The suppression of the bispectrum amplitude is a factor of two greater than the suppression of the small scale power-spectrum. For an LSST-like weak lensing survey that observes half of the sky with five tomographic redshift bins, we explore the constraining power of the bispectrum on three cosmological parameters: the sum of the neutrino mass $\sum m_\nu$, the matter density $\Omega_m$ and the amplitude of primordial fluctuations $A_s$. Bispectrum measurements alone provide similar constraints to the power spectrum measurements and combining the two probes leads to significant improvements than using the latter alone. We find that the joint constraints tighten the power spectrum $95\%$ constraints by $\sim 32\%$ for $\sum m_\nu$, $13\%$ for $\Omega_m$ and $57\%$ for $A_s$ .

The cosmic microwave background (CMB) fluctuations effectively measure the basic properties of the universe during the recombination epoch. CMB measurements fix the distance to the surface of last scatter, the sound horizon of the baryon-photon fluid and the fraction of the energy density in relativistic species. We show that the microwave background observations can also very effectively constrain the thickness of the last scattering surface, which is directly related to the ratio of the small-scale E-mode polarization signal to the small-scale temperature signal. The current cosmological data enables a 0.1\% measurement of the thickness of the surface of last scatter: $19 \pm 0.065$ Mpc. This constraint is relatively model-independent, so it can provide a new metric for systematic errors and an independent test of the $\Lambda {\rm CDM}$ model. On the other hand, it is sensitive to models which affect the reionization history of the universe such as models with annihilating dark matter and varying fundamental constants (e.g., the fine-structure constant, $\alpha_{\rm EM}$, and electron rest mass, $m_{\rm e}$) and as such can be used as a viable tool to constrain them.

The observation of GW170817 in both gravitational and electromagnetic waves provides a number of unique tests of general relativity. One question we can answer with this event is: Do large-wavelength gravitational waves and short-frequency photons experience the same number of spacetime dimensions? In models that include additional non-compact spacetime dimensions, as the gravitational waves propagate, they "leak" into the extra dimensions, leading to a reduction in the amplitude of the observed gravitational waves, and a commensurate systematic error in the inferred distance to the gravitational wave source. Electromagnetic waves would remain unaffected. We compare the inferred distance to GW170817 from the observation of gravitational waves, $d_L^\mathrm{GW}$, with the inferred distance to the electromagnetic counterpart NGC 4993, $d_L^\mathrm{EM}$. We constrain $d_L^\mathrm{GW} = (d_L^\mathrm{EM}/\mathrm{Mpc})^\mathrm{\gamma}$ with $\gamma = 1.01^{+0.04}_{-0.05}$ (for the SHoES value of $H_0$) or $\gamma = 0.99^{+0.03}_{-0.05}$ (for the Planck value of $H_0$), where all values are MAP and minimal 68% credible intervals. These constraints imply that gravitational waves propagate in $D=3+1$ spacetime dimensions, as expected in general relativity. In particular, we find that $D = 4.02^{+0.07}_{-0.10}$ (SHoES) and $D = 3.98^{+0.07}_{-0.09}$ (Planck). Furthermore, we place limits on the screening scale for theories with $D>4$ spacetime dimensions, finding that the screening scale must be greater than $\sim 20$ Mpc. We also place a lower limit on the lifetime of the graviton of $t > 4.50 \times 10^8$ yr.

Gravitational lensing deflects the paths of cosmic infrared background (CIB) photons, leaving a measurable imprint on CIB maps. The resulting statistical anisotropy can be used to reconstruct the matter distribution out to the redshifts of CIB sources. To this end, we generalize the CMB lensing quadratic estimator to any weakly non-Gaussian source field, by deriving the optimal lensing weights. We point out the additional noise and bias caused by the non-Gaussianity and the `self-lensing' of the source field. We propose methods to reduce, subtract or model these non-Gaussianities. We show that CIB lensing should be detectable with Planck data, and detectable at high significance for future CMB experiments like CCAT-Prime. The CIB thus constitutes a new source image for lensing studies, providing constraints on the amplitude of structure at intermediate redshifts between galaxies and the CMB. CIB lensing measurements will also give valuable information on the star formation history in the universe, constraining CIB halo models beyond the CIB power spectrum. By laying out a detailed treatment of lens reconstruction from a weakly non-Gaussian source field, this work constitutes a stepping stone towards lens reconstruction from continuum or line intensity mapping data, such as the Lyman-alpha emission, absorption, and the 21cm radiation.

The initial conditions of cosmological simulations are commonly drawn from a Gaussian ensemble. The limited number of modes inside a simulation volume gives rise to statistical fluctuations known as \textit{sample variance}, limiting the accuracy of simulation predictions. Fixed fields offer an alternative initialization strategy; they have the same power spectrum as standard Gaussian fields but without intrinsic amplitude scatter at linear order. Paired fixed fields consists of two fixed fields with opposite phases that cancel phase correlations which otherwise induce second-order scatter in the non-linear power spectrum. We study the statistical properties of those fields for 19 different quantities at different redshifts through a large set of 600 N-body and 506 state-of-the-art magneto-hydrodynamic simulations covering a wide range of scales, mass and spatial resolutions. We find that paired fixed simulations do not introduce a bias on any of the examined quantities. We quantify the statistical improvement brought by these simulations, over standard ones, on different power spectra such as matter, halos, CDM, gas, stars, black-holes and magnetic fields, finding that they can reduce their variance by factors as large as $10^6$. We quantify the improvement achieved by fixing and by pairing, showing that sample variance in some quantities can be highly suppressed by pairing after fixing. Paired fixed simulations do not change the scatter in quantities such as the probability distribution function of matter density, or the halo, void or stellar mass functions. We argue that procedures aiming at reducing the sample variance of those quantities are unlikely to work. Our results show that paired fixed simulations do not affect either mean relations or scatter of galaxy properties, and suggest that the information embedded in 1-pt statistics is highly complementary to that in clustering.

We study the signal from patchy reionization in view of the future high accuracy polarization measurements of the Cosmic Microwave Background (CMB). We implement an extraction procedure of the patchy reionization signal analogous to CMB lensing. We evaluate the signal to noise ratio (SNR) for the future Stage IV (S4) CMB experiment. The signal has a broad peak centered on the degree angular scales, with a long tail at higher multipoles. The CMB S4 experiment can effectively constrain the properties of reionization by measuring the signal on degree scales. The signal amplitude depends on the properties of the structure determining the reionization morphology. We describe bubbles having radii distributed log-normally. The expected S/N is sensitive to the mean bubble radius: $\bar{R}=5$ Mpc implies $S/N \approx 4$, $\bar{R}=10$ Mpc implies $S/N \approx 20$. The spread of the radii distribution strongly affects the integrated SNR, that changes by a factor of $10^2$ when $\sigma_{lnr}$ goes from $\ln 2$ to $\ln3$. Future CMB experiments will thus place important constraints on the physics of reionization.

[Abridged] We study the abundance and clustering properties of HI at redshifts $z\leqslant5$ using TNG100, a large state-of-the-art magneto-hydrodynamic simulation of a 75 Mpc/h box size. We show that most of the HI lies within dark matter halos and quantify the average HI mass hosted by halos of mass M at redshift z. We find that only halos with circular velocities larger than $\simeq$ 30 km/s contain HI. While the density profiles of HI exhibit a large halo-to-halo scatter, the mean profiles are universal across mass and redshift. The HI in low-mass halos is mostly located in the central galaxy, while in massive halos is concentrated in the satellites. We show that the HI and matter density probability distribution functions differ significantly. Our results point out that for small halos the HI bulk velocity goes in the same direction and has the same magnitude as the halo peculiar velocity, while in large halos differences show up. We find that halo HI velocity dispersion follows a power-law with halo mass. We find a complicated HI bias, with HI becoming non-linear already at $k=0.3$ h/Mpc at $z\gtrsim3$. Our simulation reproduces the DLAs bias value from observations. We find that the clustering of HI can be accurately reproduced by perturbative methods. We identify a new secondary bias, by showing that the clustering of halos depends not only on mass but also on HI content. We compute the amplitude of the HI shot-noise and find that it is small at all redshifts. We study the clustering of HI in redshift-space, and show that linear theory can explain the ratio between the monopoles in redshift- and real-space down to small scales at high redshift. We find that the amplitude of the Fingers-of-God effect is larger for HI than for matter. We point out that accurate 21 cm maps can be created from N-body or approximate simulations rather than full hydrodynamic simulations.

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.

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap."

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

We present the measurement of the kinematic Sunyaev-Zel'dovich (kSZ) effect in Fourier space, rather than in real space. We measure the density-weighted pairwise kSZ power spectrum, the first use of this promising approach, by cross-correlating a cleaned Cosmic Microwave Background (CMB) temperature map, which jointly uses both Planck Release 2 and Wilkinson Microwave Anisotropy Probe nine-year data, with the two galaxy samples, CMASS and LOWZ, derived fr om the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12. With the current data, we constrain the average optical depth $\tau$ multiplied by the ratio of the Hubble parameter at redshift $z$ and the present day, $E=H/H_0$; we find $\tau E = (3.95\pm1.62)\times10^{-5}$ for LOWZ and $\tau E = ( 1.25\pm 1.06)\times10^{-5}$ for CMASS, with the optimal angular radius of an aperture photometry filter to estimate the CMB temperature distortion associ ated with each galaxy. By repeating the pairwise kSZ power analysis for various aperture radii, we measure the optical depth as a function of aperture ra dii. While this analysis results in the kSZ signals with only evidence for a detection, ${\rm S/N}=2.54$ for LOWZ and $1.24$ for CMASS, the combination of future CMB and spectroscopic galaxy surveys should enable precision measurements. We estimate that the combination of CMB-S4 and data from DESI shoul d yield detections of the kSZ signal with ${\rm S/N}=70-100$, depending on the resolution of CMB-S4.

WFIRST will conduct a coronagraphic program of characterizing the atmospheres of planets around bright nearby stars. When observed with the WFIRST Wide Field Camera, these stars will saturate the detector and produce very strong diffraction spikes. In this paper, we forecast the astrometric precision that WFIRST can achieve by centering on the diffraction spikes of highly saturated stars. This measurement principle is strongly facilitated by the WFIRST H4RG detectors, which confine excess charges within the potential well of saturated pixels. By adopting a simplified analytical model of the diffraction spike caused by a single support strut obscuring the telescope aperture, integrated over the WFIRST pixel size, we predict the performance of this approach with the Fisher-matrix formalism. We discuss the validity of the model and find that 10 ${\mu}$as astrometric precision is achievable with a single 100 s exposure of a R = 6 or a J = 5 star. We discuss observational limitations from the optical distortion correction and pixel-level artifacts, which need to be calibrated at the level of 10 - 20 ${\mu}$as so as to not dominate the error budget. To suppress those systematics, we suggest a series of short exposures, dithered by at least several hundred pixels, to reach an effective per-visit astrometric precision of better than 10 ${\mu}$as. If this can be achieved, a dedicated WFIRST GO program will be able to detect Earth-mass exoplanets with orbital periods of 1 yr around stars within a few pc as well as Neptune-like planets with shorter periods or around more massive or distant stars. Such a program will also enable mass measurements of many anticipated direct-imaging exoplanet targets of the WFIRST coronagraph and a "starshade" occulter.

Calculations of the Cosmic Microwave Background lensing power implemented into the standard cosmological codes such as CAMB and CLASS usually treat the surface of last scatter as an infinitely thin screen. However, since the CMB anisotropies are smoothed out on scales smaller than the diffusion length due to the effect of Silk damping, the photons which carry information about the small-scale density distribution come from slightly earlier times than the standard recombination time. The dominant effect is the scale dependence of the mean redshift associated with the fluctuations during recombination. We find that fluctuations at $k = 0.01 {\rm \ Mpc^{-1}}$ come from a characteristic redshift of $z \approx 1090$, while fluctuations at $k = 0.3 {\rm \ Mpc^{-1}}$ come from a characteristic redshift of $z \approx 1130$. We then estimate the corrections to the lensing kernel and the related power spectra due to this effect. We conclude that neglecting it would result in a deviation from the true value of the lensing kernel at the half percent level at small CMB scales. For an all-sky, noise-free experiment, this corresponds to a $\sim 0.1 \sigma$ shift in the observed temperature power spectrum on small scales ($ 2500 \lesssim l \lesssim 4000$).

Using $\sim$140 deg$^2$ Subaru Hyper Suprime-Cam (HSC) survey data, we stack the weak lensing (WL) signal around five Planck clusters found within the footprint. This yields a 15$\sigma$ detection of the mean Planck cluster mass density profile. The five Planck clusters span a relatively wide mass range, $M_{\rm WL,500c} = (2-30)\times10^{14}\,M_\odot/h$ with a mean mass of $M_{\rm WL,500c} = (4.15\pm0.61)\times10^{14}\,M_\odot/h$. The ratio of the stacked Planck Sunyaev-Zel'dovich (SZ) mass to the stacked WL mass is $ \langle M_{\rm SZ}\rangle/\langle M_{\rm WL}\rangle = 1-b = 0.80\pm0.14$. This mass bias is consistent with previous WL mass calibrations of Planck clusters within the errors. We discuss the implications of our findings for the calibration of SZ cluster counts and the much discussed tension between Planck SZ cluster counts and Planck $\Lambda$CDM cosmology.

We present a catalog of 182 galaxy clusters detected through the Sunyaev-Zel'dovich effect by the Atacama Cosmology Telescope in a contiguous 987.5 deg$^{2}$ field. The clusters were detected as SZ decrements by applying a matched filter to 148 GHz maps that combine the original ACT equatorial survey with data from the first two observing seasons using the ACTPol receiver. Optical/IR confirmation and redshift measurements come from a combination of large public surveys and our own follow-up observations. Where necessary, we measured photometric redshifts for clusters using a pipeline that achieves accuracy $\Delta z/(1 + z)=0.015$ when tested on SDSS data. Under the assumption that clusters can be described by the so-called Universal Pressure Profile and its associated mass-scaling law, the full signal-to-noise > 4 sample spans the mass range $1.6 < M^{\rm UPP}_{\rm 500c}/10^{14}{\rm M}_{\odot}<9.1$, with median $M^{\rm UPP}_{\rm 500c}=3.1 \times 10^{14}$ M$_{\odot}$. The sample covers the redshift range $0.1 < z < 1.4$ (median $z = 0.49$) and 28 clusters are new discoveries (median $z = 0.80$). We compare our catalog with other overlapping cluster samples selected using the SZ, optical,and X-ray wavelengths. We find the ratio of the UPP-based SZ mass to richness-based weak-lensing mass is $\langle M^{\rm UPP}_{\rm 500c} \rangle / \langle M^{\rm \lambda WL}_{\rm 500c} \rangle = 0.68 \pm 0.11$. After applying this calibration, the mass distribution for clusters with $M_{\rm 500c} > 4 \times 10^{14}$ M$_{\odot}$ is consistent with the number of such clusters found in the South Pole Telescope SZ survey.

The non-zero mass of neutrinos suppresses the growth of cosmic structure on small scales. Since the level of suppression depends on the sum of the masses of the three active neutrino species, the evolution of large-scale structure is a promising tool to constrain the total mass of neutrinos and possibly shed light on the mass hierarchy. In this work, we investigate these effects via a large suite of N-body simulations that include massive neutrinos using an analytic linear-response approximation: the Cosmological Massive Neutrino Simulations (MassiveNuS). The simulations include the effects of radiation on the background expansion, as well as the clustering of neutrinos in response to the nonlinear dark matter evolution. We allow three cosmological parameters to vary: the neutrino mass sum M_nu in the range of 0-0.6 eV, the total matter density Omega_m, and the primordial power spectrum amplitude A_s. The rms density fluctuation in spheres of 8 comoving Mpc/h (sigma_8) is a derived parameter as a result. Our data products include N-body snapshots, halo catalogues, merger trees, ray- traced galaxy lensing convergence maps for four source redshift planes between z_s=1-2.5, and ray-traced cosmic microwave background lensing convergence maps. We describe the simulation procedures and code validation in this paper. The data are publicly available at http://columbialensing.org.

The brighter fatter effect has been postulated to arise due to the build up of a transverse electric field, produced as photo-charges accumulate in the pixels' potential wells. We investigate the brighter fatter effect in Hyper Suprime-Cam by examining flat fields and moments of stars. We observe deviations from the expected linear relation in the photon transfer curve, luminosity dependent correlations between pixels in flat field images and a luminosity dependent point spread function (PSF) in stellar observations. Under the key assumptions of translation invariance and Maxwell's equations in the quasi-static limit, we give a first-principles proof that the effect can be parametrized by a translationally invariant scalar kernel. We describe how this kernel can be estimated from flat fields and discuss how this kernel has been used to remove the brighter fatter distortions in Hyper Suprime-Cam images. We find that our correction restores the expected linear relation in the photon transfer curves and significantly reduces, but does not completely remove, the luminosity dependence of the PSF over a wide range of magnitudes.

We present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for four clusters in the current coverage area out of a sample of 22 MCXC clusters. We additionally analyze a non-MCXC cluster associated with one MCXC cluster. We show that H.E. masses for the MCXC clusters are correlated with cluster richness from the CAMIRA catalog (Oguri et al. 2017), while that for the non-MCXC cluster deviates from the scaling relation. The mass normalization of the relationship between the cluster richness and H.E. mass is compatible with one inferred by matching CAMIRA cluster abundance with a theoretical halo mass function. The mean gas mass fraction based on H.E. masses for the MCXC clusters is $\langle f_{\rm gas} \rangle = 0.125\pm0.012$ at spherical overdensity $\Delta=500$, which is $\sim80-90$ percent of the cosmic mean baryon fraction, $\Omega_b/\Omega_m$, measured by cosmic microwave background experiments. We find that the mean baryon fraction estimated from X-ray and HSC-SSP optical data is comparable to $\Omega_b/\Omega_m$. A weak-lensing shear catalog of background galaxies, combined with photometric redshifts, is currently available only for three clusters in our sample. Hydrostatic equilibrium masses roughly agree with weak-lensing masses, albeit with large uncertainty. This study demonstrates that further multiwavelength study for a large sample of clusters using X-ray, HSC-SSP optical and weak lensing data will enable us to understand cluster physics and utilize cluster-based cosmology.

The improving sensitivity of measurements of the kinetic Sunyaev-Zel'dovich (SZ) effect opens a new window into the thermodynamic properties of the baryons in halos. We propose a methodology to constrain these thermodynamic properties by combining the kinetic SZ, which is an unbiased probe of the free electron density, and the thermal SZ, which probes their thermal pressure. We forecast that our method constrains the average thermodynamic processes that govern the energetics of galaxy evolution like energetic feedback across all redshift ranges where viable halos sample are available. Current Stage-3 cosmic microwave background (CMB) experiments like AdvACT and SPT-3G can measure the kSZ and tSZ to greater than 100$\sigma$ if combined with a DESI-like spectroscopic survey. Such measurements translate into percent-level constraints on the baryonic density and pressure profiles and on the feedback and non-thermal pressure support parameters for a given ICM model. This in turn will provide critical thermodynamic tests for sub-grid models of feedback in cosmological simulations of galaxy formation. The high fidelity measurements promised by the next generation CMB experiment, CMB-S4, allow one to further sub-divide these constraints beyond redshift into other classifications, like stellar mass or galaxy type.

International efforts play a key role in driving all areas of astrophysics. Strategic planning is essential to explore possible partnerships and joint projects that otherwise could not be afforded, and to maximize their scientific return. For these reasons, the International Astronomical Union (IAU) established in 2016 the Working Group on Global Coordination of Ground and Space Astrophysics. The Working Group held a Kavli IAU Workshop on "Global Coordination: Future Space-Based Ultraviolet-Optical-Infrared Telescopes" in July 2017, bringing together scientific leaders from 17 countries to discuss compelling science drivers, technical requirements, political constraints, and opportunities for a future large-scale UV/optical/Infrared space mission. This report summarizes the workshop and recommends that astronomers worldwide intensify their activities to explore the possibilities for science with such a mission, including consideration of instrumentation and technology definition and development, construction, launch, data analysis, and complementary science. The workshop presentations are available at https://www.strw.leidenuniv.nl/KavliIAU2017.

We report and discuss the discovery of a comoving pair of bright solar-type stars, HD 240430 and HD 240429, with a significant difference in their chemical abundances. The two stars have an estimated 3D separation of $\approx 0.6$ pc ($\approx 0.01$ pc projected) at a distance of $r\approx 100$ pc with nearly identical three-dimensional velocities, as inferred from Gaia TGAS parallaxes and proper motions, and high-precision radial velocity measurements. Stellar parameters determined from high-resolution Keck HIRES spectra indicate that both stars are $\sim 4$ Gyr old. The more metal-rich of the two, HD 240430, shows an enhancement of refractory ($T_C>1200$ K) elements by $\approx 0.2$ dex and a marginal enhancement of (moderately) volatile elements ($T_C<1200$ K, C, N, O, Na, and Mn). This is the largest metallicity difference found in a wide binary pair yet. Additionally, HD 240430 shows an anomalously high surface lithium abundance ($A(\mathrm{Li})=2.75$), higher than its companion by $0.5$ dex. The proximity in phase-space and ages between the two stars suggests that they formed together with the same composition, at odds with the observed differences in metallicity and abundance patterns. We therefore suggest that the star HD~240430, "Kronos", accreted 15 $M_\oplus$ of rocky material after birth, selectively enhancing the refractory elements as well as lithium in its surface and convective envelope.

Members of widely-separated stellar binaries (> 0.1 pc) or loosely bound stellar associations and clusters are born with small spreads in kinematics and chemistry, but disrupt and dissociate on timescales comparable to their orbital time within the Milky Way. The kinematic properties of comoving star pairs---bound wide binaries or unbound but coeval stars---are therefore sensitive to both the smooth tidal field and gravitational perturbation spectrum of the Milky Way. In previous work, we identified >4,000 candidate comoving pairs using only astrometric data from the Tycho-Gaia Astrometric Solution (TGAS) catalog and find a large number of pairs with separations > 1 pc. Without radial velocity measurements for the vast majority of these pairs, the false-positive rate at these separations is significant (~50%). Here we present results from our own low-resolution radial-velocity survey of a random sample of the candidate pairs to identify and validate true comoving star pairs. Of the 311 observed comoving pairs, we confirm 127 comoving pairs with separations as large as ~10 pc, our original search limit. Between separations of 0.001 pc to ~0.5 pc, the number of confirmed pairs per separation decreases with increasing separation (number per log separation is uniform). From ~0.5 pc to 10 pc, the number per separation is approximately uniform (number per log separation increases). We confirm the discovery of a population of comoving star pairs at very large separations, suggesting that disrupted wide binaries and stellar associations remain approximately spatially coherent.

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.

(abridged) We investigate the signatures left by the cosmic neutrino background on the clustering of matter, CDM+baryons and halos in redshift-space using a set of more than 1000 N-body and hydrodynamical simulations with massless and massive neutrinos. We find that the effect neutrinos induce on the clustering of CDM+baryons in redshift-space on small scales is almost entirely due to the change in $\sigma_8$. Neutrinos imprint a characteristic signature in the quadrupole of the matter (CDM+baryons+neutrinos) field on small scales, that can be used to disentangle the effect of $\sigma_8$ and $M_\nu$. We show that the effect of neutrinos on the clustering of halos is very different, on all scales, to the one induced by $\sigma_8$. We find that the effects of neutrinos of the growth rate of CDM+baryons ranges from $\sim0.3\%$ to $2\%$ on scales $k\in[0.01, 0.5]~h{\rm Mpc}^{-1}$ for neutrinos with masses $M_\nu \leqslant 0.15$ eV. We compute the bias between the momentum of halos and the momentum of CDM+baryon and find it to be 1 on large scales for all models with massless and massive neutrinos considered. This point towards a velocity bias between halos and total matter on large scales that it is important to account for in order to extract unbiased neutrino information from velocity/momentum surveys such as kSZ observations. We show that baryonic effects can affect the clustering of matter and CDM+baryons in redshift-space by up to a few percent down to $k=0.5~h{\rm Mpc}^{-1}$. We find that hydrodynamics and astrophysical processes, as implemented in our simulations, only distort the relative effect that neutrinos induce on the anisotropic clustering of matter, CDM+baryons and halos in redshift-space by less than $1\%$. Thus, the effect of neutrinos in the fully non-linear regime can be written as a transfer function with very weak dependence on astrophysics.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

The primary sample of the {\it Gaia} Data Release 1 is the Tycho-Gaia Astrometric Solution (TGAS): $\approx$ 2 million Tycho-2 sources with improved parallaxes and proper motions relative to the initial catalog. This increased astrometric precision presents an opportunity to find new binary stars and moving groups. We search for high-confidence comoving pairs of stars in TGAS by identifying pairs of stars consistent with having the same 3D velocity using a marginalized likelihood ratio test to discriminate candidate comoving pairs from the field population. Although we perform some visualizations using (bias- corrected) inverse parallax as a point estimate of distance, the likelihood ratio is computed with a probabilistic model that includes the covariances of parallax and proper motions and marginalizes the (unknown) true distances and 3D velocities of the stars. We find 13,085 comoving star pairs among 10,606 unique stars with separations as large as 10 pc (our search limit). Some of these pairs form larger groups through mutual comoving neighbors: many of these pair networks correspond to known open clusters and OB associations, but we also report the discovery of several new comoving groups. Most surprisingly, we find a large number of very wide ($>1$ pc) separation comoving star pairs, the number of which increases with increasing separation and cannot be explained purely by false-positive contamination. Our key result is a catalog of high-confidence comoving pairs of stars in TGAS. We discuss the utility of this catalog for making dynamical inferences about the Galaxy, testing stellar atmosphere models, and validating chemical abundance measurements.

We present a new measurement of the kinematic Sunyaev-Zeldovich effect using data from the Atacama Cosmology Telescope (ACT) and the Baryon Oscillation Spectroscopic Survey (BOSS). Using 600 square degrees of overlapping sky area, we evaluate the mean pairwise baryon momentum associated with the positions of 50,000 bright galaxies in the BOSS DR11 Large Scale Structure catalog. A non-zero signal arises from the large-scale motions of halos containing the sample galaxies. The data fits an analytical signal model well, with the optical depth to microwave photon scattering as a free parameter determining the overall signal amplitude. We estimate the covariance matrix of the mean pairwise momentum as a function of galaxy separation, using microwave sky simulations, jackknife evaluation, and bootstrap estimates. The most conservative simulation-based errors give signal-to-noise estimates between 3.6 and 4.1 for varying galaxy luminosity cuts. We discuss how the other error determinations can lead to higher signal-to-noise values, and consider the impact of several possible systematic errors. Estimates of the optical depth from the average thermal Sunyaev-Zeldovich signal at the sample galaxy positions are broadly consistent with those obtained from the mean pairwise momentum signal.

We measure the gravitational lensing shear signal around dark matter halos hosting CMASS galaxies using light sources at $z\sim 1$ (background galaxies) and at the surface of last scattering at $z\sim 1100$ (the cosmic microwave background). The galaxy shear measurement uses data from the CFHTLenS survey, and the microwave background shear measurement uses data from the {\it Planck} satellite. The ratio of shears from these cross-correlations provides a purely geometric distance measurement across the longest possible cosmological lever arm. This is because the matter distribution around the halos, including uncertainties in galaxy bias and systematic errors such as miscentering, cancels in the ratio for halos in thin redshift slices. We measure this distance ratio in three different redshift slices of the CMASS sample, and combine them to obtain a $17\%$ measurement of the distance ratio, $r=0.390^{+0.070}_{-0.062}$ at an effective redshift of $z=0.53$. This is consistent with the predicted ratio from the {\it Planck} best-fit $\Lambda$CDM cosmology of $r=0.419$.

There are well-observed differences between cool-core (CC) and non-cool-core (NCC) clusters, but the origin of this distinction is still largely unknown. Competing theories can be divided into internal (inside-out), in which internal physical processes transform or maintain the NCC phase, and external (outside-in), in which the cluster type is determined by its initial conditions, which in turn lead to different formation histories (i.e., assembly bias). We propose a new method that uses the relative assembly bias of CC to NCC clusters, as determined via the two-point cluster-galaxy cross-correlation function (CCF), to test whether formation history plays a role in determining their nature. We apply our method to 48 ACCEPT clusters, which have well resolved central entropies, and cross-correlate with the SDSS-III/BOSS LOWZ galaxy catalog. We find that the relative bias of NCC over CC clusters is $b = 1.42 \pm 0.35$ ($1.6\sigma$ different from unity). Our measurement is limited by the small number of clusters with core entropy information within the BOSS footprint, 14 CC and 34 NCC. Future compilations of X-ray cluster samples, combined with deep all-sky redshift surveys, will be able to better constrain the relative assembly bias of CC and NCC clusters and determine the origin of the bimodality.

We present cosmological constraints from the combination of the full mission 9-year WMAP release and small-scale temperature data from the pre-Planck ACT and SPT generation of instruments. This is an update of the analysis presented in Calabrese et al. 2013 and highlights the impact on $\Lambda$CDM cosmology of a 0.06 eV massive neutrino - which was assumed in the Planck analysis but not in the ACT/SPT analyses - and a Planck-cleaned measurement of the optical depth to reionization. We show that cosmological constraints are now strong enough that small differences in assumptions about reionization and neutrino mass give systematic differences which are clearly detectable in the data. We recommend that these updated results be used when comparing cosmological constraints from WMAP, ACT and SPT with other surveys or with current and future full-mission Planck cosmology. Cosmological parameter chains are publicly available on the NASA's LAMBDA data archive.

Yes. Future CMB experiments such as Advanced ACTPol and CMB-S4 should achieve measurements with S/N of $> 0.1$ for the typical galaxies in redshift surveys. These measurements will provide complementary measurements of the growth rate of large scale structure $f$ and the expansion rate of the Universe $H$ to galaxy clustering measurements. This paper emphasizes that there is significant information in the anisotropy of the relative pairwise kSZ measurements. We expand the relative pairwise kSZ power spectrum in Legendre polynomials and consider up to its octopole. Assuming that the noise in the filtered maps is uncorrelated between the positions of galaxies in the survey, we derive a simple analytic form for the power spectrum covariance of the relative pairwise kSZ temperature in redshift space. While many previous studies have assumed optimistically that the optical depth of the galaxies $\tau_{\rm T}$ in the survey is known, we marginalize over $\tau_{\rm T}$, to compute constraints on the growth rate $f$ and the expansion rate $H$. For realistic sure parameters, we find that combining kSZ and galaxy redshift survey data reduces the marginalized $1$-$\sigma$ errors on $H$ and $f$ by $\sim50$-$70\%$ compared to the galaxy-only analysis.

We present the development of a photometrically selected massive galaxy catalog, targeting Luminous Red Galaxies (LRGs) and massive blue galaxies at redshifts $z\geq 0.55$. Massive galaxy candidates are selected using infrared/optical color-color cuts, with optical data from Sloan Digital Sky Survey (SDSS) and infrared data from "unWISE" forced photometry derived from the Wide-field Infrared Survey Explorer (WISE). The selection method is based on previously developed techniques to select LRGs with $z>0.5$, and is optimized using receiver operating characteristic (ROC) curves. The catalog contains 16,191,145 objects, selected over the full SDSS DR10 footprint. The redshift distribution of the resulting catalog is estimated using spectroscopic redshifts from the DEEP2 Galaxy Redshift Survey and photometric redshifts from COSMOS. Restframe $U-B$ colors from DEEP2 are used to estimate LRG selection efficiency. Using DEEP2, the resulting catalog has average redshift $z=0.65$, with standard deviation $\sigma = 2.0$, and average restframe $U-B=1.0$, with standard deviation $\sigma=0.27$. Using COSMOS, the resulting catalog has average redshift $z=0.60$, with standard deviation $\sigma = 1.8$. We estimate 34% of the catalog to be blue galaxies with $z\geq0.55$. An estimated 9.6% of selected objects are blue sources with redshift $z<0.55$. Stellar contamination is estimated to be 1.8%.

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.

Correlations between intrinsic galaxy shapes on large-scales arise due to the effect of the tidal field of the large-scale structure. Anisotropic primordial non-Gaussianity induces a distinct scale-dependent imprint in these tidal alignments on large scales. Motivated by the observational finding that the alignment strength of luminous red galaxies depends on how galaxy shapes are measured, we study the use of two different shape estimators as a multi-tracer probe of intrinsic alignments. We show, by means of a Fisher analysis, that this technique promises a significant improvement on anisotropic non-Gaussianity constraints over a single-tracer method. For future weak lensing surveys, the uncertainty in the anisotropic non-Gaussianity parameter, $A_2$, is forecast to be $\sigma(A_2)\approx 50$, $\sim 40\%$ smaller than currently available constraints from the bispectrum of the Cosmic Microwave Background. This corresponds to an improvement of a factor of $4-5$ over the uncertainty from a single-tracer analysis.

We present the first detection of a correlation between the Lyman-$\alpha$ forest and cosmic microwave background (CMB) lensing. For each Lyman-$\alpha$ forest in SDSS-III/BOSS DR12, we correlate the one-dimensional power spectrum with the CMB lensing convergence on the same line of sight from Planck. This measurement constitutes a position-dependent power spectrum, or a squeezed bispectrum, and quantifies the non-linear response of the Lyman-$\alpha$ forest power spectrum to a large-scale overdensity. The signal is measured at 5~$\sigma$ and is consistent with the $\Lambda$CDM expectation. We measure the linear bias of the Lyman-$\alpha$ forest with respect to the dark matter distribution, and constrain a combination of non-linear terms including the non-linear bias. This new observable provides a consistency check for the Lyman-$\alpha$ forest as a large-scale structure probe and tests our understanding of the relation between intergalactic gas and dark matter. In the future, it could be used to test hydrodynamical simulations and calibrate the relation between the Lyman-$\alpha$ forest and dark matter.

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

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

The kinematic Sunyaev-Zel'dovich (kSZ) effect --- the Doppler boosting of cosmic microwave background (CMB) photons due to Compton-scattering off free electrons with non-zero bulk velocity --- probes the abundance and distribution of baryons in the Universe. All kSZ measurements to date have explicitly required spectroscopic redshifts. Here, we implement a novel estimator for the kSZ -- large-scale structure cross-correlation based on projected fields: it does not require redshift estimates for individual objects, allowing kSZ measurements from large-scale imaging surveys. We apply this estimator to cleaned CMB temperature maps constructed from Planck and Wilkinson Microwave Anisotropy Probe data and a galaxy sample from the Wide-field Infrared Survey Explorer (WISE). We measure the kSZ effect at 3.8-4.5$\sigma$ significance, depending on the use of additional WISE galaxy bias constraints. We verify that our measurements are robust to possible dust emission from the WISE galaxies. Assuming the standard $\Lambda$CDM cosmology, we directly constrain $( {f_{b}}/{0.158} ) ( {f_{\rm free}}/{1.0} ) = 1.48 \pm 0.19$ (statistical error only) at redshift $z \approx 0.4$, where $f_{b}$ is the fraction of matter in baryonic form and $f_{\rm free}$ is the free electron fraction. This is the tightest kSZ-derived constraint reported to date on these parameters. The consistency between the $f_{b}$ value found here and the values inferred from analyses of the primordial CMB and Big Bang nucleosynthesis verifies that baryons approximately trace the dark matter distribution down to $\sim$Mpc scales. While our projected-field estimator is already competitive with other kSZ approaches when applied to current datasets (because we are able to use the full-sky WISE photometric survey), it will yield enormous signal-to-noise when applied to upcoming high-resolution, multi-frequency CMB surveys.

Mass calibration uncertainty is the largest systematic effect for using clusters of galaxies to constrain cosmological parameters. We present weak lensing mass measurements from the Canada-France-Hawaii Telescope Stripe 82 Survey for galaxy clusters selected through their high signal-to-noise thermal Sunyaev-Zeldovich (tSZ) signal measured with the Atacama Cosmology Telescope (ACT). For a sample of 9 ACT clusters with a tSZ signal-to-noise greater than five the average weak lensing mass is $\left(4.8\pm0.8\right)\,\times10^{14}\,\mathrm{M}_\odot$, consistent with the tSZ mass estimate of $\left(4.70\pm1.0\right)\,\times10^{14}\,\mathrm{M}_\odot$ which assumes a universal pressure profile for the cluster gas. Our results are consistent with previous weak-lensing measurements of tSZ-detected clusters from the Planck satellite. When comparing our results, we estimate the Eddington bias correction for the sample intersection of Planck and weak-lensing clusters which was previously excluded.

In recent years there have been significant improvements in the sensitivity and the angular resolution of the instruments dedicated to the observation of the Cosmic Microwave Background (CMB). ACTPol is the first polarization receiver for the Atacama Cosmology Telescope (ACT) and is observing the CMB sky with arcmin resolution over about 2000 sq. deg. Its upgrade, Advanced ACTPol (AdvACT), will observe the CMB in five frequency bands and over a larger area of the sky. We describe the optimization and implementation of the ACTPol and AdvACT surveys. The selection of the observed fields is driven mainly by the science goals, that is, small angular scale CMB measurements, B-mode measurements and cross-correlation studies. For the ACTPol survey we have observed patches of the southern galactic sky with low galactic foreground emissions which were also chosen to maximize the overlap with several galaxy surveys to allow unique cross-correlation studies. A wider field in the northern galactic cap ensured significant additional overlap with the BOSS spectroscopic survey. The exact shapes and footprints of the fields were optimized to achieve uniform coverage and to obtain cross-linked maps by observing the fields with different scan directions. We have maximized the efficiency of the survey by implementing a close to 24 hour observing strategy, switching between daytime and nighttime observing plans and minimizing the telescope idle time. We describe the challenges represented by the survey optimization for the significantly wider area observed by AdvACT, which will observe roughly half of the low-foreground sky. The survey strategies described here may prove useful for planning future ground-based CMB surveys, such as the Simons Observatory and CMB Stage IV surveys.

The next generation weak lensing surveys (i.e., LSST, Euclid and WFIRST) will require exquisite control over systematic effects. In this paper, we address shear calibration and present the most realistic forecast to date for LSST/Euclid/WFIRST and CMB lensing from a stage 4 CMB experiment (CMB S4). We use the CosmoLike code to simulate a joint analysis of all the two-point functions of galaxy density, galaxy shear and CMB lensing convergence. We include the full Gaussian and non-Gaussian covariances and explore the resulting joint likelihood with Monte Carlo Markov Chains. We constrain shear calibration biases while simultaneously varying cosmological parameters, galaxy biases and photometric redshift uncertainties. We find that CMB lensing from CMB S4 enables the calibration of the shear biases down to 0.2% - 3% in 10 tomographic bins for LSST (below the ~0.5% requirements in most tomographic bins), down to 0.4% - 2.4% in 10 bins for Euclid and 0.6% - 3.2% in 10 bins for WFIRST. For a given lensing survey, the method works best at high redshift where shear calibration is otherwise most challenging. This self-calibration is robust to Gaussian photometric redshift uncertainties and to a reasonable level of intrinsic alignment. It is also robust to changes in the beam and the effectiveness of the component separation of the CMB experiment, and slowly dependent on its depth, making it possible with third generation CMB experiments such as AdvACT and SPT-3G, as well as the Simons Observatory.

We show that the gravitational wave source counts distribution can test how gravitational radiation propagates on cosmological scales. This test does not require obtaining redshifts for the sources. If the signal-to-noise ratio (SNR, $\rho$) from a gravitational wave source is proportional to the strain then it falls as $R^{-1}$, thus we expect the source counts to follow $dN/d\rho \propto \rho^{-4}$. However, if gravitational waves decay as they propagate or propagate into other dimensions, then there can be deviations from this generic prediction. We consider the possibility that the strain falls as $R^{-\gamma}$, where $\gamma=1$ recovers the expected predictions in a Euclidean uniformly-filled universe, and forecast the sensitivity of future observations to deviations from standard General Relativity. We first consider the case of few objects, 7 sources, with a signal-to-noise from 8 to 24, and impose a lower limit on $\gamma$, finding $\gamma>0.33$ at $95\%$ confidence level. The distribution of our simulated sample is very consistent with the distribution of the trigger events reported by Advanced LIGO. Future measurements will improve these constraints: with 100 events, we estimate that $\gamma$ can be measured with an uncertainty of $15\%$. We generalize the formalism to account for a range of chirp masses and the possibility that the signal falls as $\exp(-R/R_0)/R^\gamma$.

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.

Cold Dark Matter (CDM) models struggle to match the observations at galactic scales. The tension can be reduced either by dramatic baryonic feedback effects or by modifying the particle physics of CDM. Here, we consider an ultra-light scalar field DM particle manifesting a wave nature below a DM particle mass-dependent Jeans scale. For DM mass $m\sim10^{-22}{\rm eV}$, this scenario delays galaxy formation and avoids cusps in the center of the dark matter haloes. We use new measurements of half-light mass in ultra-faint dwarf galaxies Draco II and Triangulum II to estimate the mass of the DM particle in this model. We find that if the stellar populations are within the core of the density profile then the data are in agreement with a wave dark matter model having a DM particle with $m\sim 3.7-5.6\times 10^{-22}{\rm eV}$. The presence of this extremely light particle will contribute to the formation of a central solitonic core replacing the cusp of a Navarro-Frenk-White profile and bringing predictions closer to observations of cored central density in dwarf galaxies.

The kinematic Sunyaev-Zel'dovich (kSZ) signal is a powerful probe of the cosmic baryon distribution. The kSZ signal is proportional to the integrated free electron momentum rather than the electron pressure (which sources the thermal SZ signal). Since velocities should be unbiased on large scales, the kSZ signal is an unbiased tracer of the large-scale electron distribution, and thus can be used to detect the "missing baryon" that evade most observational techniques. While most current methods for kSZ extraction rely on the availability of very accurate redshifts, we revisit a method that allows measurements even in the absence of redshift information for individual objects. It involves cross-correlating the square of an appropriately filtered cosmic microwave background (CMB) temperature map with a projected density map constructed from a sample of large-scale structure tracers. We show that this method will achieve high signal-to-noise when applied to the next generation of high-resolution CMB experiments, provided that component separation is sufficiently effective at removing foreground contamination. Considering statistical errors only, we forecast that this estimator can yield $S/N \approx $ 3, 120 and over 150 for Planck, Advanced ACTPol, and hypothetical Stage-IV CMB experiments, respectively, in combination with a galaxy catalog from WISE, and about 20% larger $S/N$ for a galaxy catalog from the proposed SPHEREx experiment. This work serves as a companion paper to the first kSZ measurement with this method, where we used CMB temperature maps constructed from Planck and WMAP data, together with galaxies from the WISE survey, to obtain a 3.8 - 4.5$\sigma$ detection of the kSZ$^2$ amplitude.

Observations of the kinetic Sunyaev-Zel'dovich (kSZ) effect measure the density-weighted velocity field, a potentially powerful cosmological probe. This paper presents an analytical method to predict the power spectrum and two-point correlation function of the density-weighted velocity in redshift space, the direct observables in kSZ surveys. We show a simple relation between the density power spectrum and the density-weighted velocity power spectrum that holds for both dark matter and halos. Using this relation, we can then extend familiar perturbation expansion techniques to the kSZ power spectrum. One of the most important features of density-weighted velocity statistics in redshift space is the change in sign of the cross-correlation between the density and density-weighted velocity at mildly small scales due to nonlinear redshift space distortions. Our model can explain this characteristic feature without any free parameters. As a result, our results can precisely predict the non-linear behavior of the density-weighted velocity field in redshift space up to $\sim30\ h^{-1} {\rm Mpc}$ for dark matter particles at the redshifts of $z=0.0$, $0.5$, and $1.0$.

We present significant evidence of halo assembly bias for SDSS redMaPPer galaxy clusters in the redshift range $[0.1, 0.33]$. By dividing the 8,648 clusters into two subsamples based on the average member galaxy separation from the cluster center, we first show that the two subsamples have very similar halo mass of $M_{\rm 200m}\simeq 1.9\times 10^{14}~h^{-1}M_\odot$ based on the weak lensing signals at small radii $R<\sim 10~h^{-1}{\rm Mpc}$. However, their halo bias inferred from both the large-scale weak lensing and the projected auto-correlation functions differs by a factor of $\sim$1.5, which is a signature of assembly bias. The same bias hypothesis for the two subsamples is excluded at 2.5$\sigma$ in the weak lensing and 4.4$\sigma$ in the auto-correlation data, respectively. This result could bring a significant impact on both galaxy evolution and precision cosmology.

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$\sigma$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $\Lambda$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies.

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.

The cosmic microwave background (CMB) energy spectrum is a near-perfect blackbody. The standard model of cosmology predicts small spectral distortions to this form, but no such distortion of the sky-averaged CMB spectrum has yet been measured. We calculate the largest expected distortion, which arises from the inverse Compton scattering of CMB photons off hot, free electrons, known as the thermal Sunyaev-Zel'dovich (tSZ) effect. We show that the predicted signal is roughly one order of magnitude below the current bound from the COBE-FIRAS experiment, but can be detected at enormous significance ($\gtrsim 1000\sigma$) by the proposed Primordial Inflation Explorer (PIXIE). Although cosmic variance reduces the effective signal-to-noise to $230\sigma$, this measurement will still yield a sub-percent constraint on the total thermal energy of electrons in the observable universe. Furthermore, we show that PIXIE can detect subtle relativistic effects in the sky-averaged tSZ signal at $30\sigma$, which directly probe moments of the optical depth-weighted intracluster medium electron temperature distribution. These effects break the degeneracy between the electron density and temperature in the mean tSZ signal, allowing a direct inference of the mean baryon density at low redshift. Future spectral distortion probes will thus determine the global thermodynamic properties of ionized gas in the universe with unprecedented precision. These measurements will impose a fundamental "integral constraint" on models of galaxy formation and the injection of feedback energy over cosmic time.

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.

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.

While observations of large-scale structure and the cosmic microwave background (CMB) provide strong constraints on the amplitude of the primordial power spectrum (PPS) on scales larger than 10~Mpc, the amplitude of the power spectrum on sub-galactic length scales is much more poorly constrained. We study early structure formation in a cosmological model with a blue-tilted PPS. We assume that the standard scale-invariant PPS is modified at small length scales as $P(k) \sim k^{m_{\rm s}}$ with $m_{\rm s} > 1$. We run a series of cosmological hydrodynamic simulations to examine the dependence of the formation epoch and the characteristic mass of primordial stars on the tilt of the PPS. In models with $m_{\rm s} > 1$, star-forming gas clouds are formed at $z > 100$, when formation of hydrogen molecules is inefficient because the intense CMB radiation destroys chemical intermediates. Without efficient coolant, the gas clouds gravitationally contract while keeping a high temperature. The protostars formed in such "hot" clouds grow very rapidly by accretion to become extremely massive stars that may leave massive black holes with a few hundred solar-masses at $z > 100$. The shape of the PPS critically affects the properties and the formation epoch of the first generation of stars. Future experiments of the CMB polarization and the spectrum distortion may provide important information on the nature of the first stars and their formation epoch, and hence on the shape of the small-scale power spectrum.

We use the Planck full mission temperature maps to examine the stacked thermal Sunyaev-Zel'dovich (tSZ) signal of 188042 "locally brightest galaxies'" (LBGs) selected from the Sloan Digital Sky Survey Data Release 7. Our LBG sample closely matches that of Planck Collaboration XI (2013, PCXI) but our analysis differs in several ways. We work directly in terms of physically observable quantities, requiring minimal assumptions about the gas pressure profile. We explicitly model the dust emission from each LBG and simultaneously measure both the stacked tSZ and dust signals as a function of stellar mass $M_*$. There is a small residual bias in stacked tSZ measurements; we measure this bias and subtract it from our results, finding that the effects are non-negligible at the lowest masses in the LBG sample. Finally, we compare our measurements with two pressure profile models, finding that the profile from Battaglia et al. (2012b) provides a better fit to the results than the Arnaud et al. (2010) "universal pressure profile." However, within the uncertainties, we find that the data are consistent with a self-similar scaling with mass --- more precise measurements are needed to detect the relatively small deviations from self-similarity predicted by these models. Consistent with PCXI, we measure the stacked tSZ signal from LBGs with stellar masses down to $\log_{10}(M_*/M_{\odot}) \sim 11.1-11.3$. For lower stellar masses, however, we do not see evidence for a stacked tSZ signal. We note that the stacked dust emission is comparable to, or larger than, the stacked tSZ signal for $\log_{10}(M_*/M_{\odot}) \lesssim 11.3$. Future tSZ analyses with larger samples and lower noise levels should be able to probe deviations from self-similarity and thus provide constraints on models of feedback and the evolution of hot halo gas over cosmic time.

The Prime Focus Spectrograph (PFS) is an optical/near-infrared multifiber spectrograph with 2394 science fibers distributed across a 1.3-deg diameter field of view at the Subaru 8.2-m telescope. The wide wavelength coverage from 0.38 {\mu}m to 1.26 {\mu}m, with a resolving power of 3000, simultaneously strengthens its ability to target three main survey programs: cosmology, galactic archaeology and galaxy/AGN evolution. A medium resolution mode with a resolving power of 5000 for 0.71 {\mu}m to 0.89 {\mu}m will also be available by simply exchanging dispersers. We highlight some of the technological aspects of the design. To transform the telescope focal ratio, a broad-band coated microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of the cable system, optimizing overall throughput; a fiber with low focal ratio degradation is selected for the fiber-positioner and fiber-slit components, minimizing the effects of fiber movements and fiber bending. Fiber positioning will be performed by a positioner consisting of two stages of piezo-electric rotary motors. The positions of these motors are measured by taking an image of artificially back-illuminated fibers with the metrology camera located in the Cassegrain container; the fibers are placed in the proper location by iteratively measuring and then adjusting the positions of the motors. Target light reaches one of the four identical fast-Schmidt spectrograph modules, each with three arms. The PFS project has passed several project-wide design reviews and is now in the construction phase.

Tidal debris from infalling satellites can leave observable structure in the phase-space distribution of the Galactic halo. Such substructure can be manifest in the spatial and/or velocity distributions of the stars in the halo. This paper focuses on a class of substructure that is purely kinematic in nature, with no accompanying spatial features. To study its properties, we use a simulated stellar halo created by dynamically populating the Via Lactea II high-resolution N-body simulation with stars. A significant fraction of the stars in the inner halo of Via Lactea share a common speed and metallicity, despite the fact that they are spatially diffuse. We argue that this kinematic substructure is a generic feature of tidal debris from older mergers and may explain the detection of radial-velocity substructure in the inner halo made by the Sloan Extension for Galactic Understanding and Exploration. The GAIA satellite, which will provide the proper motions of an unprecedented number of stars, should further characterize the kinematic substructure in the inner halo. Our study of the Via Lactea simulation suggests that the stellar halo can be used to map the speed distribution of the local dark-matter halo, which has important consequences for dark-matter direct-detection experiments.

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 perform a joint analysis of the abundance, the clustering and the galaxy-galaxy lensing signal of galaxies measured from Data Release 11 of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (SDSS III-BOSS) in our companion paper, Miyatake et al. (2014). The lensing signal was obtained by using the shape catalog of background galaxies from the Canada France Hawaii Telescope Legacy Survey, which was made publicly available by the CFHTLenS collaboration, with an area overlap of about 105 deg$^2$. We analyse the data in the framework of the halo model in order to fit halo occupation parameters and cosmological parameters ($\Omega_{\rm m}$ and $\sigma_8$) to these observables simultaneously, and thus break the degeneracy between galaxy bias and cosmology. Adopting a flat $\Lambda$CDM cosmology with priors on $\Omega_b h^2$, $n_{\rm s}$ and $h$ from the analysis of WMAP 9-year data, we obtain constraints on the stellar mass-halo mass relation of galaxies in our sample. Marginalizing over the halo occupation distribution parameters and a number of other nuisance parameters in our model, we obtain $\Omega_{\rm m}=0.310^{+0.019}_{-0.020}$ and $\sigma_8=0.785^{+0.044}_{-0.044}$ (68% confidence). We demonstrate the robustness of our results with respect to sample selection and a variety of systematics such as the halo off-centering effect and possible incompleteness in our sample. Our constraints are consistent, complementary and competitive with those obtained using other independent probes of these cosmological parameters. The cosmological analysis is the first of its kind to be performed at a redshift as high as $0.53$.

A joint analysis of the clustering of galaxies and their weak gravitational lensing signal is well-suited to simultaneously constrain the galaxy-halo connection as well as the cosmological parameters by breaking the degeneracy between galaxy bias and the amplitude of clustering signal. In a series of two papers, we perform such an analysis at the highest redshift ($z\sim0.53$) in the literature using CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Eleventh Data Release (SDSS-III/BOSS DR11) catalog spanning 8300~deg$^2$. In this paper, we present details of the clustering and weak lensing measurements of these galaxies. We define a subsample of 400,916 CMASS galaxies based on their redshifts and stellar mass estimates so that the galaxies constitute an approximately volume-limited and similar population over the redshift range $0.47\le z\le 0.59$. We obtain a signal-to-noise ratio $S/N\simeq 56$ for the galaxy clustering measurement. We also explore the redshift and stellar mass dependence of the clustering signal. For the weak lensing measurement, we use existing deeper imaging data from the CFHTLS with publicly available shape and photometric redshift catalogs from CFHTLenS, but only in a 105~deg$^2$ area which overlaps with BOSS. This restricts the lensing measurement to only 5,084 CMASS galaxies. After careful systematic tests, we find a highly significant detection of the CMASS weak lensing signal, with total $S/N\simeq 26$. These measurements form the basis of the halo occupation distribution and cosmology analysis presented in More et al. (Paper II).

This report describes the 2014 study by the Science Definition Team (SDT) of the Wide-Field Infrared Survey Telescope (WFIRST) mission. It is a space observatory that will address the most compelling scientific problems in dark energy, exoplanets and general astrophysics using a 2.4-m telescope with a wide-field infrared instrument and an optical coronagraph. The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Study Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality telescope assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012) and the previous WFIRST-2.4 DRM, reported by Spergel et. (2013). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The addition of an on-axis coronagraphic instrument to the baseline design enables imaging and spectroscopic studies of planets around nearby stars.

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 focus of this report is on the opportunities enabled by the combination of LSST, Euclid and WFIRST, the optical surveys that will be an essential part of the next decade's astronomy. The sum of these surveys has the potential to be significantly greater than the contributions of the individual parts. As is detailed in this report, the combination of these surveys should give us multi-wavelength high-resolution images of galaxies and broadband data covering much of the stellar energy spectrum. These stellar and galactic data have the potential of yielding new insights into topics ranging from the formation history of the Milky Way to the mass of the neutrino. However, enabling the astronomy community to fully exploit this multi-instrument data set is a challenging technical task: for much of the science, we will need to combine the photometry across multiple wavelengths with varying spectral and spatial resolution. We identify some of the key science enabled by the combined surveys and the key technical challenges in achieving the synergies.

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.

Weakly Interacting Massive Particles (WIMPs) may constitute most of the matter in the Universe. The ability to detect the directionality of recoil nuclei will considerably facilitate detection of WIMPs. In this paper we propose a novel type of dark matter detector: detectors made of DNA or RNA could provide nanometer resolution for tracking, an energy threshold of 0.5 keV, and can operate at room temperature. When a WIMP from the Galactic Halo elastically scatters off of a nucleus in the detector, the recoiling nucleus then traverses hundreds of strings of single stranded nucleic acids (ssNA) with known base sequences and severs ssNA strands along its trajectory. The location of the break can be identified by amplifying and identifying the segments of cut ssNA using techniques well known to biologists. Thus the path of the recoiling nucleus can be tracked to nanometer accuracy. In one such detector concept, the transducers are nanometer-thick Au-foils of 1m x 1m, and the direction of recoiling nuclei is measured by "NA Tracking Chamber" consisting of ordered array of ssNA strands. Polymerase Chain Reaction (PCR) and ssNA sequencing are used to read-out the detector. The proposed detector is smaller and cheaper than other alternatives: 1 kg of gold and 0.1 to 4 kg of ssNA (depending on length and strand density), packed into 0.01m$^3$, can be used to study 10 GeV WIMPs. A variety of other detector target elements could be used in this detector to optimize for different WIMP masses and to identify WIMP properties. By leveraging advances in molecular biology, we aim to achieve about 1,000-fold better spatial resolution than in conventional WIMP detectors at reasonable cost.

We consider the cosmological consequences if a small fraction ($f\lesssim 0.1$) of the dark matter is ultra-strongly self-interacting, with an elastic self-interaction cross-section per unit mass $\sigma\gg1\ \mathrm{cm^{2}/g}$. This possibility evades all current constraints that assume that the self-interacting component makes up the majority of the dark matter. Nevertheless, even a small fraction of ultra-strongly self-interacting dark matter (uSIDM) can have observable consequences on astrophysical scales. In particular, the uSIDM subcomponent can undergo gravothermal collapse and form seed black holes in the center of a halo. These seed black holes, which form within several hundred halo interaction times, contain a few percent of the total uSIDM mass in the halo. For reasonable values of $\sigma f$, these black holes can form at high enough redshifts to grow to $\sim10^9 M_\odot$ quasars by $z \gtrsim 6$, alleviating tension within the standard $\Lambda$CDM cosmology. The ubiquitous formation of central black holes in halos could also create cores in dwarf galaxies by ejecting matter during binary black hole mergers, potentially resolving the "too big to fail" problem.

We measure the cross-correlation of cosmic microwave background lensing convergence maps derived from Atacama Cosmology Telescope data with galaxy lensing convergence maps as measured by the Canada-France-Hawaii Telescope Stripe 82 Survey. The CMB-galaxy lensing cross power spectrum is measured for the first time with a significance of 4.2{\sigma}, which corresponds to a 12% constraint on the amplitude of density fluctuations at redshifts ~ 0.9. With upcoming improved lensing data, this novel type of measurement will become a powerful cosmological probe, providing a precise measurement of the mass distribution at intermediate redshifts and serving as a calibrator for systematic biases in weak lensing measurements.

The Ly$\alpha$ forest transmission probability distribution function (PDF) is an established probe of the intergalactic medium (IGM) astrophysics, especially the temperature-density relationship of the IGM. We measure the transmission PDF from 3393 Baryon Oscillations Spectroscopic Survey (BOSS) quasars from SDSS Data Release 9, and compare with mock spectra that include careful modeling of the noise, continuum, and astrophysical uncertainties. The BOSS transmission PDFs, measured at $\langle z \rangle = [2.3,2.6,3.0]$, are compared with PDFs created from mock spectra drawn from a suite of hydrodynamical simulations that sample the IGM temperature-density relationship, $\gamma$, and temperature at mean-density, $T_0$, where $T(\Delta) = T_0 \Delta^{\gamma-1}$. We find that a significant population of partial Lyman-limit systems with a column-density distribution slope of $\beta_\mathrm{pLLS} \sim -2$ are required to explain the data at the low-transmission end of transmission PDF, while uncertainties in the mean Ly$\alpha$ forest transmission affect the high-transmission end. After modelling the LLSs and marginalizing over mean-transmission uncertainties, we find that $\gamma=1.6$ best describes the data over our entire redshift range, although constraints on $T_0$ are affected by systematic uncertainties. Within our model framework, isothermal or inverted temperature-density relationships ($\gamma \leq 1$) are disfavored at a significance of over 4$\sigma$, although this could be somewhat weakened by cosmological and astrophysical uncertainties that we did not model.

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.

The Wide-Field InfraRed Survey Telescope (WFIRST) is a NASA space mission in study for launch in 2024. It has a 2.4 m telescope, wide-field IR instrument operating in the 0.7 - 2.0 micron range and an exoplanet imaging coronagraph instrument operating in the 400 - 1000 nm range. The observatory will perform galaxy surveys over thousands of square degrees to J=27 AB for dark energy weak lensing and baryon acoustic oscillation measurements and will monitor a few square degrees for dark energy SN Ia studies. It will perform microlensing observations of the galactic bulge for an exoplanet census and direct imaging observations of nearby exoplanets with a pathfinder coronagraph. The mission will have a robust and well-funded guest observer program for 25% of the observing time. WFIRST will be a powerful tool for time domain astronomy and for coordinated observations with gravitational wave experiments. Gravitational wave events produced by mergers of nearby binary neutron stars (LIGO-Virgo) or extragalactic supermassive black hole binaries (LISA) will produce electromagnetic radiation that WFIRST can observe.

We present a statistical analysis of the millimeter-wavelength properties of 1.4 GHz-selected sources and a detection of the Sunyaev-Zel'dovich (SZ) effect associated with the halos that host them. The Atacama Cosmology Telescope (ACT) has conducted a survey at 148 GHz, 218 GHz and 277 GHz along the celestial equator. Using samples of radio sources selected at 1.4 GHz from FIRST and NVSS, we measure the stacked 148, 218 and 277 GHz flux densities for sources with 1.4 GHz flux densities ranging from 5 to 200 mJy. At these flux densities, the radio source population is dominated by active galactic nuclei (AGN), with both steep and flat spectrum populations, which have combined radio-to-millimeter spectral indices ranging from 0.5 to 0.95, reflecting the prevalence of steep spectrum sources at high flux densities and the presence of flat spectrum sources at lower flux densities. The thermal SZ effect associated with the halos that host the AGN is detected at the 5$\sigma$ level through its spectral signature. When we compare the SZ effect with weak lensing measurements of radio galaxies, we find that the relation between the two is consistent with that measured by Planck for local bright galaxies. We present a detection of the SZ effect in some of the lowest mass halos (average $M_{200}\approx10^{13}$M$_{\odot}h_{70}^{-1}$) studied to date. This detection is particularly important in the context of galaxy evolution models, as it confirms that galaxies with radio AGN also typically support hot gaseous halos. With Herschel observations, we show that the SZ detection is not significantly contaminated by dust. We show that 5 mJy$<S_{1.4}<$200 mJy radio sources contribute $\ell(\ell+1)C_{\ell}/(2\pi)=0.37\pm0.03\mu$K$^2$ to the angular power spectrum at $\ell=3000$ at 148 GHz, after accounting for the SZ effect associated with their host halos.

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.

We apply our recently developed code to search for resonance features in the Planck CMB temperature data. We search both for log spaced oscillations or linear spaced oscillations and compare our findings with results of our WMAP9 analysis and the Planck team analysis. While there are hints of log spaced resonant features present in the WMAP9 data, the significance of these features weaken with more data. With more accurate small scale measurements, we also find that the best fit frequency has shifted and the amplitude has been reduced. We confirm the presence of a several low frequency peaks, earlier identified by the Planck team, but with a better improvement of fit (delta chi^2 ~ 12). We further investigate this improvement by allowing the lensing potential to vary as well, showing mild correlation between the amplitude of the oscillations and the lensing amplitude. We find that the improvement of the fit increases even more (delta chi^2 ~ 14) for the low frequencies that modify the spectrum in a way that mimics the lensing effect. Since these features were not present in the WMAP data, they are primarily due to better measurements of Planck at small angular scales. For linear spaced oscillations we find a maximum delta chi^2 ~ 13 scanning two orders of magnitude in frequency space, and the biggest improvements are at extremely high frequencies. We recover a best fit frequency very close to the one found in WMAP9, which confirms that the fit improvement is driven by low l. Further comparisons with WMAP9 show Planck contains many more features, both for linear and log space oscillations, but with a smaller improvement of fit. We discuss the improvement as a function of the number of modes and study the effect of the 217 GHz map, which appears to drive most of the improvement for log spaced oscillations. We conclude that none of the detected features are statistically significant.

A small deviation from scale invariance in the form of oscillations in the primordial correlation spectra has been predicted by various cosmological models. In this paper we review a recently developed method to search for these features in the data in a more effective way. By Taylor expanding the small features around the 'background' cosmology, we have shown we are able to improve the search for these features compared to previous analyses. In this short paper we will extend that work by combining this method with a multi nested sampler. We recover our previous findings and are able to do so in 192 CPU hours. We will also briefly discuss the possibility of a long wavelength feature in the data to alleviate some tension between CMB data and the LCDM+r concordance cosmology.

The Prime Focus Spectrograph (PFS) is an optical/near-infrared multi-fiber spectrograph with 2394 science fibers, which are distributed in 1.3 degree diameter field of view at Subaru 8.2-meter telescope. The simultaneous wide wavelength coverage from 0.38 um to 1.26 um, with the resolving power of 3000, strengthens its ability to target three main survey programs: cosmology, Galactic archaeology, and galaxy/AGN evolution. A medium resolution mode with resolving power of 5000 for 0.71 um to 0.89 um also will be available by simply exchanging dispersers. PFS takes the role for the spectroscopic part of the Subaru Measurement of Images and Redshifts project, while Hyper Suprime-Cam works on the imaging part. To transform the telescope plus WFC focal ratio, a 3-mm thick broad-band coated glass-molded microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of cable system, while one with a better FRD performance is selected for the fiber-positioner and fiber-slit components, given the more frequent fiber movements and tightly curved structure. Each Fiber positioner consists of two stages of piezo-electric rotary motors. Its engineering model has been produced and tested. Fiber positioning will be performed iteratively by taking an image of artificially back-illuminated fibers with the Metrology camera located in the Cassegrain container. The camera is carefully designed so that fiber position measurements are unaffected by small amounts of high special-frequency inaccuracies in WFC lens surface shapes. Target light carried through the fiber system reaches one of four identical fast-Schmidt spectrograph modules, each with three arms. Prototype VPH gratings have been optically tested. CCD production is complete, with standard fully-depleted CCDs for red arms and more-challenging thinner fully-depleted CCDs with blue-optimized coating for blue arms.

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.

BICEP2 has reported the detection of a degree-scale B-mode polarization pattern in the Cosmic Microwave Background (CMB) and has interpreted the measurement as evidence for primordial gravitational waves. Motivated by the profound importance of the discovery of gravitational waves from the early Universe, we examine to what extent a combination of Galactic foregrounds and lensed E-modes could be responsible for the signal. We reanalyze the BICEP2 results and show that the 100x150 GHz and 150x150 GHz data are consistent with a cosmology with r=0.2 and negligible foregrounds, but also with a cosmology with r=0 and a significant dust polarization signal. We give independent estimates of the dust polarization signal in the BICEP2 region using four different approaches. While these approaches are consistent with each other, the expected amplitude of the dust polarization power spectrum remains uncertain by about a factor of three. The lower end of the prediction leaves room for a primordial contribution, but at the higher end the dust in combination with the standard CMB lensing signal could account for the BICEP2 observations, without requiring the existence of primordial gravitational waves. By measuring the cross-correlations between the pre-Planck templates used in the BICEP2 analysis and between different versions of a data-based template, we emphasize that cross-correlations between models are very sensitive to noise in the polarization angles and that measured cross-correlations are likely underestimates of the contribution of foregrounds to the map. These results suggest that BICEP1 and BICEP2 data alone cannot distinguish between foregrounds and a primordial gravitational wave signal, and that future Keck Array observations at 100 GHz and Planck observations at higher frequencies will be crucial to determine whether the signal is of primordial origin. (abridged)

Naive estimates of the statistics of large scale structure and weak lensing power spectrum measurements that include only Gaussian errors exaggerate their scientific impact. Non-linear evolution and finite volume effects are both significant sources of non-Gaussian covariance that reduce the ability of power spectrum measurements to constrain cosmological parameters. Using a halo model formalism, we derive an intuitive understanding of the various contributions to the covariance and show that our analytical treatment agrees with simulations. This approach enables an approximate derivation of a joint likelihood for the cluster number counts, the weak lensing power spectrum and the bispectrum. We show that this likelihood is a good description of the ray-tracing simulation. Since all of these observables are sensitive to the same finite volume effects and contain information about the non-linear evolution, a combined analysis recovers much of the "lost" information and obviates the non-Gaussian covariance. For upcoming weak lensing surveys, we estimate that a joint analysis of power spectrum, number counts and bispectrum will produce an improvement of about $30-40\%$ in determinations of the matter density and the scalar amplitude. This improvement is equivalent to doubling the survey area.

While the precise relationship between the Milky Way disk and the symmetry planes of the dark matter halo remains somewhat uncertain, a time-varying disk orientation with respect to an inertial reference frame seems probable. Hierarchical structure formation models predict that the dark matter halo is triaxial and tumbles with a characteristic rate of ~2 rad/Hubble time (~30 muas/yr). These models also predict a time-dependent accretion of gas, such that the angular momentum vector of the disk should be misaligned with that of the halo. These effects, as well as tidal effects of the LMC, will result in the rotation of the angular momentum vector of the disk population with respect to the quasar reference frame. We assess the accuracy with which the positions and proper motions from Gaia can be referred to a kinematically non-rotating system, and show that the spin vector of the transformation from any rigid self-consistent catalog frame to the quasi-inertial system defined by quasars should be defined to better than 1 muas/yr. Determination of this inertial frame by Gaia will reveal any signature of the disk orientation varying with time, improve models of the potential and dynamics of the Milky Way, test theories of gravity, and provide new insights into the orbital evolution of the Sagittarius dwarf galaxy and the Magellanic Clouds.

Fluctuations in the intensity and polarization of the cosmic microwave background (CMB) and the large-scale distribution of matter in the universe each contain clues about the nature of the earliest moments of time. The next generation of CMB and large-scale structure (LSS) experiments are poised to test the leading paradigm for these earliest moments---the theory of cosmic inflation---and to detect the imprints of the inflationary epoch, thereby dramatically increasing our understanding of fundamental physics and the early universe. A future CMB experiment with sufficient angular resolution and frequency coverage that surveys at least 1% of the sky to a depth of 1 uK-arcmin can deliver a constraint on the tensor-to-scalar ratio that will either result in a 5-sigma measurement of the energy scale of inflation or rule out all large-field inflation models, even in the presence of foregrounds and the gravitational lensing B-mode signal. LSS experiments, particularly spectroscopic surveys such as the Dark Energy Spectroscopic Instrument, will complement the CMB effort by improving current constraints on running of the spectral index by up to a factor of four, improving constraints on curvature by a factor of ten, and providing non-Gaussianity constraints that are competitive with the current CMB bounds.

This is a report on the status and prospects of the quantification of neutrino properties through the cosmological neutrino background for the Cosmic Frontier of the Division of Particles and Fields Community Summer Study long-term planning exercise. Experiments planned and underway are prepared to study the cosmological neutrino background in detail via its influence on distance-redshift relations and the growth of structure. The program for the next decade described in this document, including upcoming spectroscopic galaxy surveys eBOSS and DESI and a new Stage-IV CMB polarization experiment CMB-S4, will achieve sigma(sum m_nu) = 16 meV and sigma(N_eff) = 0.020. Such a mass measurement will produce a high significance detection of non-zero sum m_nu, whose lower bound derived from atmospheric and solar neutrino oscillation data is about 58 meV. If neutrinos have a minimal normal mass hierarchy, this measurement will definitively rule out the inverted neutrino mass hierarchy, shedding light on one of the most puzzling aspects of the Standard Model of particle physics --- the origin of mass. This precise a measurement of N_eff will allow for high sensitivity to any light and dark degrees of freedom produced in the big bang and a precision test of the standard cosmological model prediction that N_eff = 3.046.

A coherent over- or under-density contrast across a finite survey volume causes an upward- or downward-fluctuation in the observed number of halos. This fluctuation in halo number adds a significant co-variant scatter in the observed amplitudes of weak lensing power spectrum at nonlinear, small scales -- the so-called super-sample variance or the halo sample variance. In this paper, we show that by measuring both the number counts of clusters and the power spectrum in the same survey region, we can mitigate this loss of information and significantly enhance the scientific return from the upcoming surveys. First, using the halo model approach, we derive the cross-correlation between the halo number counts and the weak lensing power spectrum, taking into account the super-sample covariance effect, which well matches the distributions measured from 1000 realizations for a Lambda-dominated cold dark matter model. Then we show that adding the observed number counts of massive halos with M> 10^14 Ms/h can significantly improve the information content of weak lensing power spectrum, almost recovering the Gaussian information up to lmax~1000, if the average mass profiles of the massive halos are known, which can be estimated from stacked lensing. When combined with the halo number counts for M > 3 or 1 x 10^14 Ms/h, the improvement is up to a factor of 1.4 or 2 at lmax~1000 --2000, equivalent to a factor 2 or 4 times larger survey volume, compared to the power spectrum measurement alone.

We present the temperature power spectrum of the Cosmic Microwave Background obtained by cross-correlating maps from the Atacama Cosmology Telescope (ACT) at 148 and 218 GHz with maps from the Planck satellite at 143 and 217 GHz, in two overlapping regions covering 592 square degrees. We find excellent agreement between the two datasets at both frequencies, quantified using the variance of the residuals between the ACT power spectra and the ACTxPlanck cross-spectra. We use these cross-correlations to calibrate the ACT data at 148 and 218 GHz, to 0.7% and 2% precision respectively. We find no evidence for anisotropy in the calibration parameter. We compare the Planck 353 GHz power spectrum with the measured amplitudes of dust and cosmic infrared background (CIB) of ACT data at 148 and 218 GHz. We also compare planet and point source measurements from the two experiments.

The nominal mission maps from the Planck satellite contain a wealth of information about secondary anisotropies in the cosmic microwave background (CMB), including those induced by the thermal Sunyaev-Zel'dovich (tSZ) effect and gravitational lensing. As both the tSZ and CMB lensing signals trace the large-scale matter density field, the anisotropies sourced by these processes are expected to be correlated. We report the first detection of this cross-correlation signal, which we measure at 6.2 sigma significance using the Planck data. We take advantage of Planck's multifrequency coverage to construct a tSZ map using internal linear combination techniques, which we subsequently cross-correlate with the publicly-released Planck CMB lensing potential map. [Abridged] We interpret the signal using halo model calculations, which indicate that the tSZ -- CMB lensing cross-correlation is a unique probe of the physics of intracluster gas in high-redshift, low-mass groups and clusters. Our results are consistent with extrapolations of existing gas physics models to this previously unexplored regime and show clear evidence for contributions from both the one- and two-halo terms, but no statistically significant evidence for contributions from diffuse, unbound gas outside of collapsed halos. We also show that the amplitude of the signal depends rather sensitively on the amplitude of fluctuations ($\sigma_8$) and the matter density ($\Omega_m$). We constrain the degenerate combination $\sigma_8 (\Omega_m/0.282)^{0.26} = 0.824 \pm 0.029$, a result that is in less tension with primordial CMB constraints than some recent tSZ analyses. Our detection is a direct confirmation that hot, ionized gas traces the dark matter distribution over a wide range of scales in the universe ($\sim 0.1$--$50 \, {\rm Mpc}/h$).

We study the non-linear behavior of the baryon acoustic oscillation in the power spectrum and the correlation function by decomposing the dark matter perturbations into the short- and long-wavelength modes. The evolution of the dark matter fluctuations can be described as a global coordinate transformation caused by the long-wavelength displacement vector acting on short-wavelength matter perturbation undergoing non-linear growth. Using this feature, we investigate the well known cancellation of the high-$k$ solutions in the standard perturbation theory. While the standard perturbation theory naturally satisfies the cancellation of the high-$k$ solutions, some of the recently proposed improved perturbation theories do not guarantee the cancellation. We show that this cancellation clarifies the success of the standard perturbation theory at the 2-loop order in describing the amplitude of the non-linear power spectrum even at high-$k$ regions.We propose an extension of the standard 2-loop level perturbation theory model of the non-linear power spectrum that more accurately models the non-linear evolution of the baryon acoustic oscillation than the standard perturbation theory. The model consists of simple and intuitive parts: the non-linear evolution of the smoothed power spectrum without the baryon acoustic oscillations and the non-linear evolution of the baryon acoustic oscillations due to the large-scale velocity of dark matter and due to the gravitational attraction between dark matter particles. Our extended model predicts the smoothing parameter of the baryon acoustic oscillation peak at $z=0.35$ as $\sim 7.7\ {\rm Mpc}/h$ and describes the small non-linear shift in the peak position due to the galaxy random motions.

We present the first set of maps and band-merged catalog from the Herschel Stripe 82 Survey (HerS). Observations at 250, 350, and 500 micron were taken with the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. HerS covers 79 deg$^2$ along the SDSS Stripe 82 to a depth of 13.0, 12.9, and 14.8 mJy beam$^{-1}$ (including confusion) at 250, 350, and 500 micron, respectively. HerS was designed to measure correlations with external tracers of the dark matter density field --- either point-like (i.e., galaxies selected from radio to X-ray) or extended (i.e., clusters and gravitational lensing) --- in order to measure the bias and redshift distribution of intensities of infrared-emitting dusty star-forming galaxies and AGN. By locating HeRS in Stripe 82, we maximize the overlap with available and upcoming cosmological surveys. The band-merged catalog contains 3.3x10$^4$ sources detected at a significance of >3 $\sigma$ (including confusion noise). The maps and catalog are available at http://www.astro.caltech.edu/hers/

The Integrated Sachs-Wolfe effect (ISW) measures the decay of the gravitational potential due to cosmic acceleration and is thus a direct probe of Dark Energy. In some of the earlier studies, the amplitude of the ISW effect was found to be in tension with the predictions of the standard $\Lambda$CDM model. We measure the cross-power of galaxies and AGN from the WISE mission with CMB temperature data from WMAP9 in order to provide an independent measurement of the ISW amplitude. Cross-correlations with the recently released Planck lensing potential maps are used to calibrate the bias and contamination fraction of the sources, thus avoiding systematic effects that could be present when using auto-spectra to measure bias. We find an amplitude of the cross-power of $\mathcal{A} = 1.24\pm 0.47$ from the galaxies and $\mathcal{A} = 0.88 \pm 0.74$ from the AGN, fully consistent with the $\Lambda$CDM prediction of $\mathcal{A} =1$. The ISW measurement signal-to-noise ratio is 2.7 and 1.2 respectively, giving a combined significance close to $3 \sigma$. Comparing the amplitudes of the galaxy and AGN cross-correlations, which arise from different redshifts, we find no evidence for redshift evolution in Dark Energy properties, consistent with a Cosmological Constant.

Yes! Upcoming galaxy shear surveys have the potential to significantly improve our understanding of dark energy and neutrino mass if lensing systematics can be sufficiently controlled. The cross-correlations between the weak lensing shear, galaxy number counts from a galaxy redshift survey, and the CMB lensing convergence can be used to calibrate the shear multiplicative bias, one of the most challenging systematics in lensing surveys. These cross-correlations can significantly reduce the deleterious effects of the uncertainties in multiplicative bias.

In this first of two papers, we present a new method for searching for oscillatory features in the primordial power spectrum. A wide variety of models predict these features in one of two different flavors: logarithmically spaced oscillations and linearly spaced oscillations. The proposed method treats the oscillations as perturbations on top of the scale-invariant power spectrum, allowing us to vary all cosmological parameters. This perturbative approach reduces the computational requirements for the search as the transfer functions and their derivatives can be precomputed. We show that the most significant degeneracy in the analysis is between the distance to last scattering and the overall amplitude at low frequencies. For models with logarithmic oscillations, this degeneracy leads to an uncertainty in the phase. For linear spaced oscillations, it affects the frequency of the oscillations. In this first of two papers, we test our code on simulated Planck-like data, and show we are able to recover fiducial input oscillations with an amplitude of a few times order 10^{-2}. We apply the code to WMAP9-year data and confirm the existence of two intriguing resonant frequencies for log spaced oscillations. For linear spaced oscillations we find a single resonance peak. We use numerical simulations to assess the significance of these features and conclude that the data do not provide compelling evidence for the existence of oscillatory features in the primordial spectrum.

We consider the cross-correlation between free electrons and neutral hydrogen during the epoch of reionization. The free electrons are traced by the optical depth to reionization tau, while the neutral hydrogen can be observed through 21 cm photon emission. As expected, this correlation is sensitive to the detailed physics of reionization. Foremost, if reionization occurs through the merger of relatively large halos hosting an ionizing source, the free electrons and neutral hydrogen are anti-correlated for most of the reionization history. A positive contribution to the correlation can occur when the halos that can form an ionizing source are small. A measurement of this sign change in the cross-correlation would directly measure the bubble bias and as, such, the halo mass. We estimate the signal-to-noise of the cross-correlation using the estimator for inhomogeneous reionization tau_{lm} proposed by Dvorkin and Smith (2009). We find that with upcoming radio interferometers and CMB experiments, the cross-correlation is measurable going up to multipoles l ~ 1000. We also derive parameter constraints and conclude that, despite the foregrounds, the cross-correlation proofs a complementary measurement of the EoR parameters to the 21 cm and CMB polarization auto-correlations expected to be observed in the coming decade.

We present constraints on cosmological and astrophysical parameters from high-resolution microwave background maps at 148 GHz and 218 GHz made by the Atacama Cosmology Telescope (ACT) in three seasons of observations from 2008 to 2010. A model of primary cosmological and secondary foreground parameters is fit to the map power spectra and lensing deflection power spectrum, including contributions from both the thermal Sunyaev-Zeldovich (tSZ) effect and the kinematic Sunyaev-Zeldovich (kSZ) effect, Poisson and correlated anisotropy from unresolved infrared sources, radio sources, and the correlation between the tSZ effect and infrared sources. The power ell^2 C_ell/2pi of the thermal SZ power spectrum at 148 GHz is measured to be 3.4 +\- 1.4 muK^2 at ell=3000, while the corresponding amplitude of the kinematic SZ power spectrum has a 95% confidence level upper limit of 8.6 muK^2. Combining ACT power spectra with the WMAP 7-year temperature and polarization power spectra, we find excellent consistency with the LCDM model. We constrain the number of effective relativistic degrees of freedom in the early universe to be Neff=2.79 +\- 0.56, in agreement with the canonical value of Neff=3.046 for three massless neutrinos. We constrain the sum of the neutrino masses to be Sigma m_nu < 0.39 eV at 95% confidence when combining ACT and WMAP 7-year data with BAO and Hubble constant measurements. We constrain the amount of primordial helium to be Yp = 0.225 +\- 0.034, and measure no variation in the fine structure constant alpha since recombination, with alpha/alpha0 = 1.004 +/- 0.005. We also find no evidence for any running of the scalar spectral index, dns/dlnk = -0.004 +\- 0.012.

The Atacama Cosmology Telescope has measured the angular power spectra of microwave fluctuations to arcminute scales at frequencies of 148 and 218 GHz, from three seasons of data. At small scales the fluctuations in the primordial Cosmic Microwave Background (CMB) become increasingly obscured by extragalactic foregounds and secondary CMB signals. We present results from a nine-parameter model describing these secondary effects, including the thermal and kinematic Sunyaev-Zel'dovich (tSZ and kSZ) power; the clustered and Poisson-like power from Cosmic Infrared Background (CIB) sources, and their frequency scaling; the tSZ-CIB correlation coefficient; the extragalactic radio source power; and thermal dust emission from Galactic cirrus in two different regions of the sky. In order to extract cosmological parameters, we describe a likelihood function for the ACT data, fitting this model to the multi-frequency spectra in the multipole range 500<ell<10000. We extend the likelihood to include spectra from the South Pole Telescope at frequencies of 95, 150, and 220 GHz. Accounting for different radio source levels and Galactic cirrus emission, the same model provides an excellent fit to both datasets simultaneously, with chi2/dof= 675/697 for ACT, and 96/107 for SPT. We then use the multi-frequency likelihood to estimate the CMB power spectrum from ACT in bandpowers, marginalizing over the secondary parameters. This provides a simplified `CMB-only' likelihood in the range 500<ell<3500 for use in cosmological parameter estimation.

We measure a significant correlation between the thermal Sunyaev-Zel'dovich effect in the Planck and WMAP maps and an X-ray cluster map based on ROSAT. We use the 100, 143 and 343 GHz Planck maps and the WMAP 94 GHz map to obtain this cluster cross spectrum. We check our measurements for contamination from dusty galaxies using the cross correlations with the 220, 545 and 843 GHz maps from Planck. Our measurement yields a direct characterization of the cluster power spectrum over a wide range of angular scales that is consistent with large cosmological simulations. The amplitude of this signal depends on cosmological parameters that determine the growth of structure (\sigma_8 and \Omega_M) and scales as \sigma_8^7.4 and \Omega_M^1.9 around the multipole (ell) ~ 1000. We constrain \sigma_8 and \Omega_M from the cross-power spectrum to be \sigma_8 (\Omega_M/0.30)^0.26 = 0.8 +/- 0.02. Since this cross spectrum produces a tight constraint in the \sigma_8 and \Omega_M plane the errors on a \sigma_8 constraint will be mostly limited by the uncertainties from external constraints. Future cluster catalogs, like those from eRosita and LSST, and pointed multi-wavelength observations of clusters will improve the constraining power of this cross spectrum measurement. In principle this analysis can be extended beyond \sigma_8 and \Omega_M to constrain dark energy or the sum of the neutrino masses.

We present a first measurement of the stellar mass component of galaxy clusters selected via the Sunyaev-Zel'dovich (SZ) effect, using 3.6 um and 4.5 um photometry from the Spitzer Space Telescope. Our sample consists of 14 clusters detected by the Atacama Cosmology Telescope (ACT), which span the redshift range 0.27 < z < 1.07 (median z = 0.50), and have dynamical mass measurements, accurate to about 30 per cent, with median M500 = 6.9 x 10^{14} MSun. We measure the 3.6 um and 4.5 um galaxy luminosity functions, finding the characteristic magnitude (m*) and faint-end slope (alpha) to be similar to those for IR-selected cluster samples. We perform the first measurements of the scaling of SZ-observables (Y500 and y0) with both brightest cluster galaxy (BCG) stellar mass and total cluster stellar mass (M500star). We find a significant correlation between BCG stellar mass and Y500 (E(z)^{-2/3} DA^2 Y500 ~ M*^{1.2 +/- 0.6}), although we are not able to obtain a strong constraint on the slope of the relation due to the small sample size. Additionally, we obtain E(z)^{-2/3} DA^2 Y500 ~ M500star^{1.0 +/- 0.6} for the scaling with total stellar mass. The mass fraction in stars spans the range 0.006-0.034, with the second ranked cluster in terms of dynamical mass (ACT-CL J0237-4939) having an unusually low total stellar mass and the lowest stellar mass fraction. For the five clusters with gas mass measurements available in the literature, we see no evidence for a shortfall of baryons relative to the cosmic mean value.

Nonlinear redshift-space distortions, the Finger-of-God (FoG) effect, can complicate the interpretation of the galaxy power spectrum. Here, we demonstrate the method proposed by Hikage et al. (2012) to use complimentary observations to directly constrain this effect on the data. We use catalogs of Luminous Red Galaxies (LRGs) and photometric galaxies from the SDSS DR7 to measure the redshift-space power spectrum of LRGs, the cross-correlation of LRGs with the shapes of background photometric galaxies (galaxy-galaxy weak lensing), and the projected cross-correlation of LRGs with photometric galaxies having similar photometric redshifts to the LRG spectroscopic redshift. All of these measurements use a reconstructed halo field. While we use the position of each LRG for single LRG systems, we compare the measurements using different halo-center proxies for multiple-LRG systems (4.5 per cent of all the halos): the brightest LRG position (BLRG), the faintest LRG position (FLRG) and their arithmetical mean position (Mean), respectively, in each system. We find significant differences in the measured correlations of different centers, showing consistent off-centering effects in the three observables. By comparing the measurements with a halo model that treats the satellite photometric galaxies as being distributed according to a generalized NFW profile, we find that about 40 (70) per cent of BLRGs (FLRGs) are off-centered satellite galaxies in the multiple-LRG systems. The satellite LRGs have typical off-centering radius of about 400 kpc/h, and velocity dispersion of about 500 km/s in host halos with a mean mass of 1.6x10^14 Ms/h. We show that, if LRGs in the single LRG systems have similar offsets, the residual FoG contamination in the LRG power spectrum can be significant at k>0.1 h/Mpc, which may cause a bias in cosmological parameters such as the neutrino mass.

The Subaru Prime Focus Spectrograph (PFS) is a massively-multiplexed fiber-fed optical and near-infrared 3-arm spectrograph (N_fiber=2400, 380<lambda<1260nm, 1.3 degree diameter FoV), offering unique opportunities in survey astronomy. Here we summarize the science case feasible for a survey of Subaru 300 nights. We describe plans to constrain the nature of dark energy via a survey of emission line galaxies spanning a comoving volume of 9.3 (Gpc/h)^3 in the redshift range 0.8<z<2.4. In each of 6 redshift bins, the cosmological distances will be measured to 3% precision via BAO, and redshift-space distortions will be used to constrain structure growth to 6% precision. In the GA program, radial velocities and chemical abundances of stars in the Milky Way and M31 will be used to infer the past assembly histories of spiral galaxies and the structure of their dark matter halos. Data will be secured for 10^6 stars in the Galactic thick-disk, halo and tidal streams as faint as V~22, including stars with V < 20 to complement the goals of the Gaia mission. A medium-resolution mode with R = 5000 to be implemented in the red arm will allow the measurement of multiple alpha-element abundances and more precise velocities for Galactic stars, elucidating the detailed chemo-dynamical structure and evolution of each of the main stellar components of the Milky Way Galaxy and of its dwarf spheroidal galaxies. For the extragalactic program, our simulations suggest the wide avelength range will be powerful in probing the galaxy population and its clustering over a wide redshift range. We propose to conduct a color-selected survey of 1<z<2 galaxies and AGN over 16 deg^2 to J~23.4, yielding a fair sample of galaxies with stellar masses above ~10^{10}Ms at z~2. A two-tiered survey of higher redshift LBGs and LAEs will quantify the properties of early systems close to the reionization epoch.

[Abridged] We present a catalog of 68 galaxy clusters, of which 19 are new discoveries, detected via the Sunyaev-Zel'dovich effect (SZ) at 148 GHz in the Atacama Cosmology Telescope (ACT) survey of 504 square degrees on the celestial equator. A subsample of 48 clusters within the 270 square degree region overlapping SDSS Stripe 82 is estimated to be 90% complete for M_500c > 4.5e14 Msun and 0.15 < z < 0.8. While matched filters are used to detect the clusters, the sample is studied further through a "Profile Based Amplitude Analysis" using a single filter at a fixed \theta_500 = 5.9' angular scale. This new approach takes advantage of the "Universal Pressure Profile" (UPP) to fix the relationship between the cluster characteristic size (R_500) and the integrated Compton parameter (Y_500). The UPP scalings are found to be nearly identical to an adiabatic model, while a model incorporating non-thermal pressure better matches dynamical mass measurements and masses from the South Pole Telescope. A high signal to noise ratio subsample of 15 ACT clusters is used to obtain cosmological constraints. We first confirm that constraints from SZ data are limited by uncertainty in the scaling relation parameters rather than sample size or measurement uncertainty. We next add in seven clusters from the ACT Southern survey, including their dynamical mass measurements based on galaxy velocity dispersions. In combination with WMAP7 these data simultaneously constrain the scaling relation and cosmological parameters, yielding \sigma_8 = 0.829 \pm 0.024 and \Omega_m = 0.292 \pm 0.025. The results include marginalization over a 15% bias in dynamical mass relative to the true halo mass. In an extension to LCDM that incorporates non-zero neutrino mass density, we combine our data with WMAP7+BAO+Hubble constant measurements to constrain \Sigma m_\nu < 0.29 eV (95% C. L.).

We present the first dynamical mass estimates and scaling relations for a sample of Sunyaev-Zel'dovich effect (SZE) selected galaxy clusters. The sample consists of 16 massive clusters detected with the Atacama Cosmology Telescope (ACT) over a 455 sq. deg. area of the southern sky. Deep multi-object spectroscopic observations were taken to secure intermediate-resolution (R~700-800) spectra and redshifts for ~60 member galaxies on average per cluster. The dynamical masses M_200c of the clusters have been calculated using simulation-based scaling relations between velocity dispersion and mass. The sample has a median redshift z=0.50 and a median mass M_200c~12e14 Msun/h70 with a lower limit M_200c~6e14 Msun/h70, consistent with the expectations for the ACT southern sky survey. These masses are compared to the ACT SZE properties of the sample, specifically, the match-filtered central SZE amplitude y, the central Compton parameter y0, and the integrated Compton signal Y_200c, which we use to derive SZE-Mass scaling relations. All SZE estimators correlate with dynamical mass with low intrinsic scatter (<~20%), in agreement with numerical simulations. We explore the effects of various systematic effects on these scaling relations, including the correlation between observables and the influence of dynamically disturbed clusters. Using the 3-dimensional information available, we divide the sample into relaxed and disturbed clusters and find that ~50% of the clusters are disturbed. There are hints that disturbed systems might bias the scaling relations but given the current sample sizes these differences are not significant; further studies including more clusters are required to assess the impact of these clusters on the scaling relations.

We develop a novel abundance matching method to construct a mock catalog of luminous red galaxies (LRGs) in SDSS, using catalogs of halos and subhalos in N-body simulations for a LCDM model. Motivated by observations suggesting that LRGs are passively-evolving, massive early-type galaxies with a typical age >5Gyr, we assume that simulated halos at z=2 (z2-halo) are progenitors for LRG-host subhalos observed today, and we label the most tightly bound particles in each progenitor z2-halo as LRG ``stars''. We then identify the subhalos containing these stars to z=0.3 (SDSS redshift) in descending order of the masses of z2-halos until the comoving number density of the matched subhalos becomes comparable to the measured number density of SDSS LRGs, n=10^{-4} (h/Mpc)^3. Once the above prescription is determined, our only free parameter is the number density of halos identified at z=2 and this parameter is fixed to match the observed number density at z = 0.3. By tracing subsequent merging and assembly histories of each progenitor z2-halo, we can directly compute, from the mock catalog, the distributions of central and satellite LRGs and their internal motions in each host halo at z=0.3. While the SDSS LRGs are galaxies selected by the magnitude and color cuts from the SDSS images and are not necessarily a stellar-mass-selected sample, our mock catalog reproduces a host of SDSS measurements: the halo occupation distribution for central and satellite LRGs, the projected auto-correlation function of LRGs, the cross-correlation of LRGs with shapes of background galaxies (LRG-galaxy weak lensing), and the nonlinear redshift-space distortion effect, the Finger-of-God effect, in the angle-averaged redshift-space power spectrum.

We present cosmological parameter constraints based on the final nine-year WMAP data, in conjunction with additional cosmological data sets. The WMAP data alone, and in combination, continue to be remarkably well fit by a six-parameter LCDM model. When WMAP data are combined with measurements of the high-l CMB anisotropy, the BAO scale, and the Hubble constant, the densities, Omegabh2, Omegach2, and Omega_L, are each determined to a precision of ~1.5%. The amplitude of the primordial spectrum is measured to within 3%, and there is now evidence for a tilt in the primordial spectrum at the 5sigma level, confirming the first detection of tilt based on the five-year WMAP data. At the end of the WMAP mission, the nine-year data decrease the allowable volume of the six-dimensional LCDM parameter space by a factor of 68,000 relative to pre-WMAP measurements. We investigate a number of data combinations and show that their LCDM parameter fits are consistent. New limits on deviations from the six-parameter model are presented, for example: the fractional contribution of tensor modes is limited to r<0.13 (95% CL); the spatial curvature parameter is limited to -0.0027 (+0.0039/-0.0038); the summed mass of neutrinos is <0.44 eV (95% CL); and the number of relativistic species is found to be 3.84+/-0.40 when the full data are analyzed. The joint constraint on Neff and the primordial helium abundance agrees with the prediction of standard Big Bang nucleosynthesis. We compare recent PLANCK measurements of the Sunyaev-Zel'dovich effect with our seven-year measurements, and show their mutual agreement. Our analysis of the polarization pattern around temperature extrema is updated. This confirms a fundamental prediction of the standard cosmological model and provides a striking illustration of acoustic oscillations and adiabatic initial conditions in the early universe.

We present the final nine-year maps and basic results from the WMAP mission. We provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate CMB anisotropy from foreground emission, and both types of signals are analyzed in detail. The WMAP mission has resulted in a highly constrained LCDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that Big Bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (3.84+/-0.40). The model fit also implies that the age of the universe is 13.772+/-0.059 Gyr, and the fit Hubble constant is H0 = 69.32+/-0.80 km/s/Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity reported earlier by WMAP now has high statistical significance (n_s = 0.9608+/-0.0080); and the universe is close to flat/Euclidean, Omega_k = -0.0027 (+0.0039/-0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter LCDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further. With no significant anomalies and an adequate goodness-of-fit, the inflationary flat LCDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

We present a model of inflation based on the interaction between a homogeneous and isotropic configuration of a U(1) gauge field and fermionic charge density $\mathcal{J}_{0}$. The regulated fermionic charge density is generated from a Bunch-Davies vacuum state using the methods of Koksma and Prokopec \cite{Koksma:2009tc}, and is found to redshift as $1/a(\eta)$. The time-like component of gauge field is sourced by the fermionic charge leading to a growth in the gauge field $A(\eta)_{0}\sim a(\eta)$. As a result inflation is dominated by the energy density contained in the gauge field and fermionic charge interaction, $A_{0}\, \mathcal{J}^{0}$, which remains constant during inflation. We also provide a mechanism to generate a net lepton asymmetry. The coupling of a pseudo scalar to the Chern-Simons term converts the gauge field fluctuations into lepton number and all three Sahkarov conditions are satisfied during inflation. Finally, the rapid oscillation of the pseudo scalar field near its minimum thermalizes the gauge field and ends inflation. We provide the necessary initial condition on the gauge field and fermionic charge to simultaneously generate enough e-folds and baryon asymmetry index.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. The report of the WFIRST-AFTA Science Definition Team (SDT) presents a Design Reference Mission for WFIRST that employs one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the smaller aperture designs previously considered for WFIRST, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This short article, produced as a companion to the SDT report, summarizes the key points of the WFIRST-2.4 DRM. It highlights the remarkable opportunity that the 2.4-m telescope affords for advances in many fields of astrophysics and cosmology, including dark energy, the demographics and characterization of exoplanets, the evolution of galaxies and quasars, and the stellar populations of the Milky Way and its neighbors.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Project Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This document presents the final report of the SDT.

We present the measured Sunyaev-Zel'dovich (SZ) flux from 474 optically-selected MaxBCG clusters that fall within the Atacama Cosmology Telescope (ACT) Equatorial survey region. The ACT Equatorial region used in this analysis covers 510 square degrees and overlaps Stripe 82 of the Sloan Digital Sky Survey. We also present the measured SZ flux stacked on 52 X-ray-selected MCXC clusters that fall within the ACT Equatorial region and an ACT Southern survey region covering 455 square degrees. We find that the measured SZ flux from the X-ray-selected clusters is consistent with expectations. However, we find that the measured SZ flux from the optically-selected clusters is both significantly lower than expectations and lower than the recovered SZ flux measured by the Planck satellite. Since we find a lower recovered SZ signal than Planck, we investigate the possibility that there is a significant offset between the optically-selected brightest cluster galaxies (BCGs) and the SZ centers, to which ACT is more sensitive due to its finer resolution. Such offsets can arise due to either an intrinsic physical separation between the BCG and the center of the gas concentration or from misidentification of the cluster BCG. We find that the entire discrepancy for both ACT and Planck can be explained by assuming that the BCGs are offset from the SZ maxima with a uniform random distribution between 0 and 1.5 Mpc. Such large offsets between gas peaks and BCGs for optically-selected cluster samples seem unlikely given that we find the physical separation between BCGs and X-ray peaks for an X-ray-selected subsample of MaxBCG clusters to have a much narrower distribution that peaks within 0.2 Mpc. It is possible that other effects are lowering the ACT and Planck signals by the same amount, with offsets between BCGs and SZ peaks explaining the remaining difference between measurements. (Abridged)

Recent data from the WMAP, ACT and SPT experiments provide precise measurements of the cosmic microwave background temperature power spectrum over a wide range of angular scales. The combination of these observations is well fit by the standard, spatially flat LCDM cosmological model, constraining six free parameters to within a few percent. The scalar spectral index, n_s = 0.9690 +/- 0.0089, is less than unity at the 3.6 sigma level, consistent with simple models of inflation. The damping tail of the power spectrum at high resolution, combined with the amplitude of gravitational lensing measured by ACT and SPT, constrains the effective number of relativistic species to be N_eff = 3.28 +/- 0.40, in agreement with the standard model's three species of light neutrinos.

We present the temperature power spectra of the cosmic microwave background (CMB) derived from the three seasons of data from the Atacama Cosmology Telescope (ACT) at 148 GHz and 218 GHz, as well as the cross-frequency spectrum between the two channels. We detect and correct for contamination due to the Galactic cirrus in our equatorial maps. We present the results of a number of tests for possible systematic error and conclude that any effects are not significant compared to the statistical errors we quote. Where they overlap, we cross-correlate the ACT and the South Pole Telescope (SPT) maps and show they are consistent. The measurements of higher-order peaks in the CMB power spectrum provide an additional test of the Lambda CDM cosmological model, and help constrain extensions beyond the standard model. The small angular scale power spectrum also provides constraining power on the Sunyaev-Zel'dovich effects and extragalactic foregrounds. We also present a measurement of the CMB gravitational lensing convergence power spectrum at 4.6-sigma detection significance.

We present the optical and X-ray properties of 68 galaxy clusters selected via the Sunyaev-Zel'dovich Effect at 148 GHz by the Atacama Cosmology Telescope (ACT). Our sample, from an area of 504 square degrees centered on the celestial equator, is divided into two regions. The main region uses 270 square degrees of the ACT survey that overlaps with the co-added ugriz imaging from the Sloan Digital Sky Survey (SDSS) over Stripe 82 plus additional near-infrared pointed observations with the Apache Point Observatory 3.5-meter telescope. We confirm a total of 49 clusters to z~1.3, of which 22 (all at z>0.55) are new discoveries. For the second region the regular-depth SDSS imaging allows us to confirm 19 more clusters up to z~0.7, of which 10 systems are new. We present the optical richness, photometric redshifts, and separation between the SZ position and the brightest cluster galaxy (BCG). We find no significant offset between the cluster SZ centroid and BCG location and a weak correlation between optical richness and SZ-derived mass. We also present X-ray fluxes and luminosities from the ROSAT All Sky Survey which confirm that this is a massive sample. One of the newly discovered clusters, ACT-CL J0044.4+0113 at z=1.1 (photometric), has an integrated XMM-Newton X-ray temperature of kT_x=7.9+/-1.0 keV and combined mass of M_200a=8.2(-2.5,+3.3)x10^14 M_sun/h70 placing it among the most massive and X-ray-hot clusters known at redshifts beyond z=1. We also highlight the optically-rich cluster ACT-CL J2327.4-0204 (RCS2 2327) at z=0.705 (spectroscopic) as the most significant detection of the whole equatorial sample with a Chandra-derived mass of M_200a=1.9(-0.4,+0.6)x10^15 M_sun/h70, comparable to some of the most massive known clusters like "El Gordo" and the Bullet Cluster.

We show that subhalos falling into the Milky Way create a flow of tidally-stripped debris particles near the galactic center with characteristic velocity behavior. In the Via Lactea-II N-body simulation, this unvirialized component constitutes a few percent of the local density and has velocities peaked at 340 km/s in the solar neighborhood. Such velocity substructure has important implications for surveys of low-metallicity stars, as well as direct detection experiments sensitive to dark matter with large scattering thresholds.

We present the BOSS Lyman-alpha (Lya) Forest Sample from SDSS Data Release 9, comprising 54,468 quasar spectra with zqso > 2.15 suitable for Lya forest analysis. This data set probes the intergalactic medium with absorption redshifts 2.0 < z_alpha < 5.7 over an area of 3275 square degrees, and encompasses an approximate comoving volume of 20 h^-3 Gpc^3. With each spectrum, we have included several products designed to aid in Lya forest analysis: improved sky masks that flag pixels where data may be unreliable, corrections for known biases in the pipeline estimated noise, masks for the cores of damped Lya systems and corrections for their wings, and estimates of the unabsorbed continua so that the observed flux can be converted to a fractional transmission. The continua are derived using a principal component fit to the quasar spectrum redwards of restframe Lya (lambda > 1216 Ang), extrapolated into the forest region and normalized by a linear function to fit the expected evolution of the Lya forest mean-flux. The estimated continuum errors are ~5% rms. We also discuss possible systematics arising from uncertain spectrophotometry and artifacts in the flux calibration; global corrections for the latter are provided. Our sample provides a convenient starting point for users to analyze clustering in BOSS Lya forest data, and it provides a fiducial data set that can be used to compare results from different analyses of baryon acoustic oscillations in the Lya forest. The full data set is available from the SDSS-III DR9 web site.

We present a detection of the unnormalized skewness <T^3> induced by the thermal Sunyaev-Zel'dovich (tSZ) effect in filtered Atacama Cosmology Telescope (ACT) 148 GHz cosmic microwave background temperature maps. Contamination due to infrared and radio sources is minimized by template subtraction of resolved sources and by constructing a mask using outlying values in the 218 GHz (tSZ-null) ACT maps. We measure <T^3>= -31 +- 6 \mu K^3 (measurement error only) or +- 14 \mu K^3 (including cosmic variance error) in the filtered ACT data, a 5-sigma detection. We show that the skewness is a sensitive probe of sigma_8, and use analytic calculations and tSZ simulations to obtain cosmological constraints from this measurement. From this signal alone we infer a value of sigma_8= 0.79 +0.03 -0.03 (68 % C.L.) +0.06 -0.06 (95 % C.L.). Our results demonstrate that measurements of non-Gaussianity can be a useful method for characterizing the tSZ effect and extracting the underlying cosmological information.

We measure the cross-correlation of Atacama Cosmology Telescope CMB lensing convergence maps with quasar maps made from the Sloan Digital Sky Survey DR8 SDSS-XDQSO photometric catalog. The CMB lensing-quasar cross-power spectrum is detected for the first time at a significance of 3.8 sigma, which directly confirms that the quasar distribution traces the mass distribution at high redshifts z>1. Our detection passes a number of null tests and systematic checks. Using this cross-power spectrum, we measure the amplitude of the linear quasar bias assuming a template for its redshift dependence, and find the amplitude to be consistent with an earlier measurement from clustering; at redshift z ~ 1.4, the peak of the distribution of quasars in our maps, our measurement corresponds to a bias of b = 2.5 +/- 0.6. With the signal-to-noise ratio on CMB lensing measurements likely to improve by an order of magnitude over the next few years, our results demonstrate the potential of CMB lensing cross-correlations to probe astrophysics at high redshifts.

We show how the correlation between the thermal Sunyaev Zel'dovich effect (tSZ) from galaxy clusters and dust emission from cosmic infrared background (CIB) sources can be calculated in a halo model framework. Using recent tSZ and CIB models, we find that the size of the tSZ x CIB cross-correlation is approximately 20 per cent at 150 GHz. The contribution to the total angular power spectrum is of order -2 \mu K^2 at ell=3000, however, this value is uncertain by a factor of two to three, primarily because of CIB source modelling uncertainties. We expect the large uncertainty in this component to degrade upper limits on the kinematic Sunyaev Zel'dovich effect (kSZ), due to similarity in the frequency dependence of the tSZ x CIB and kSZ across the frequency range probed by current Cosmic Microwave Background missions. We also find that the degree of tSZ x CIB correlation is higher for mm x sub-mm spectra than mm x mm, because more of the sub-mm CIB originates at lower redshifts (z<2), where most tSZ clusters are found.

We discuss scientific, technical and programmatic issues related to the use of an NRO 2.4m telescope for the WFIRST initiative of the 2010 Decadal Survey. We show that this implementation of WFIRST, which we call "NEW WFIRST," would achieve the goals of the NWNH Decadal Survey for the WFIRST core programs of Dark Energy and Microlensing Planet Finding, with the crucial benefit of deeper and/or wider near-IR surveys for GO science and a potentially Hubble-like Guest Observer program. NEW WFIRST could also include a coronagraphic imager for direct detection of dust disks and planets around neighboring stars, a high-priority science and technology precursor for future ambitious programs to image Earth-like planets around neighboring stars.

We describe the conceptual design of the camera cryostats, detectors, and detector readout electronics for the SuMIRe Prime Focus Spectrograph (PFS) being developed for the Subaru telescope. The SuMIRe PFS will consist of four identical spectrographs, each receiving 600 fibers from a 2400 fiber robotic positioner at the prime focus. Each spectrograph will have three channels covering wavelength ranges 3800 {\AA} - 6700 {\AA}, 6500 {\AA} - 10000 {\AA}, and 9700 {\AA} - 13000 {\AA}, with the dispersed light being imaged in each channel by a f/1.10 vacuum Schmidt camera. In the blue and red channels a pair of Hamamatsu 2K x 4K edge-buttable CCDs with 15 um pixels are used to form a 4K x 4K array. For the IR channel, the new Teledyne 4K x 4K, 15 um pixel, mercury-cadmium-telluride sensor with substrate removed for short-wavelength response and a 1.7 um cutoff will be used. Identical detector geometry and a nearly identical optical design allow for a common cryostat design with the only notable difference being the need for a cold radiation shield in the IR camera to mitigate thermal background. This paper describes the details of the cryostat design and cooling scheme, relevant thermal considerations and analysis, and discusses the detectors and detector readout electronics.

The Prime Focus Spectrograph (PFS) of the Subaru Measurement of Images and Redshifts (SuMIRe) project has been endorsed by Japanese community as one of the main future instruments of the Subaru 8.2-meter telescope at Mauna Kea, Hawaii. This optical/near-infrared multi-fiber spectrograph targets cosmology with galaxy surveys, Galactic archaeology, and studies of galaxy/AGN evolution. Taking advantage of Subaru's wide field of view, which is further extended with the recently completed Wide Field Corrector, PFS will enable us to carry out multi-fiber spectroscopy of 2400 targets within 1.3 degree diameter. A microlens is attached at each fiber entrance for F-ratio transformation into a larger one so that difficulties of spectrograph design are eased. Fibers are accurately placed onto target positions by positioners, each of which consists of two stages of piezo-electric rotary motors, through iterations by using back-illuminated fiber position measurements with a wide-field metrology camera. Fibers then carry light to a set of four identical fast-Schmidt spectrographs with three color arms each: the wavelength ranges from 0.38 {\mu}m to 1.3 {\mu}m will be simultaneously observed with an average resolving power of 3000. Before and during the era of extremely large telescopes, PFS will provide the unique capability of obtaining spectra of 2400 cosmological/astrophysical targets simultaneously with an 8-10 meter class telescope. The PFS collaboration, led by IPMU, consists of USP/LNA in Brazil, Caltech/JPL, Princeton, & JHU in USA, LAM in France, ASIAA in Taiwan, and NAOJ/Subaru.

We present a Subaru weak lensing measurement of ACT-CL J0022.2-0036, one of the most luminous, high-redshift (z=0.81) Sunyaev-Zel'dovich (SZ) clusters discovered in the 268 deg^2 equatorial region survey of the Atacama Cosmology Telescope. For the weak lensing analysis using i'-band images, we use a model-fitting (Gauss-Laguerre shapelet) method to measure shapes of galaxy images, where we fit galaxy images in different exposures simultaneously to obtain best-fit ellipticities taking into account the different PSFs in each exposure. We also take into account the astrometric distortion effect on galaxy images by performing the model fitting in the world coordinate system. To select background galaxies behind the cluster at z=0.81, we use photometric redshift (photo-z) estimates for every galaxy derived from the co-added images of multi-passband Br'i'z'Y, with PSF matching/homogenization. After a photo-z cut for background galaxy selection, we detect the tangential weak lensing distortion signal with a total signal-to-noise ratio of about 3.7. By fitting a Navarro-Frenk-White model to the measured shear profile, we find the cluster mass to be M_200\bar{\rho}_m = [7.5^+3.2_-2.8(stat.)^+1.3_-0.6(sys.)] x 10^14 M_\odot/h. The weak lensing-derived mass is consistent with previous mass estimates based on the SZ observation, with assumptions of hydrostatic equilibrium and virial theorem, as well as with scaling relations between SZ signal and mass derived from weak lensing, X-ray, and velocity dispersion, within the measurement errors. We also show that the existence of ACT-CL J0022.2-0036 at z=0.81 is consistent with the cluster abundance prediction of the \Lambda-dominated cold dark matter structure formation model. We thus demonstrate the capability of Subaru-type ground-based images for studying weak lensing of high-redshift clusters.

Tidal stripping of dark matter from subhalos falling into the Milky Way produces narrow, cold tidal streams as well as more spatially extended "debris flows" in the form of shells, sheets, and plumes. Here we focus on the debris flow in the Via Lactea II simulation, and show that this incompletely phase-mixed material exhibits distinctive high velocity behavior. Unlike tidal streams, which may not necessarily intersect the Earth's location, debris flow is spatially uniform at 8 kpc and thus guaranteed to be present in the dark matter flux incident on direct detection experiments. At Earth-frame speeds greater than 450 km/s, debris flow comprises more than half of the dark matter at the Sun's location, and up to 80% at even higher speeds. Therefore, debris flow is most important for experiments that are particularly sensitive to the high speed tail of the dark matter distribution, such as searches for light or inelastic dark matter or experiments with directional sensitivity. We show that debris flow yields a distinctive recoil energy spectrum and a broadening of the distribution of incidence direction.

We present a description of the data reduction and mapmaking pipeline used for the 2008 observing season of the Atacama Cosmology Telescope (ACT). The data presented here at 148 GHz represent 12% of the 90 TB collected by ACT from 2007 to 2010. In 2008 we observed for 136 days, producing a total of 1423 hours of data (11 TB for the 148 GHz band only), with a daily average of 10.5 hours of observation. From these, 1085 hours were devoted to a 850 deg^2 stripe (11.2 hours by 9.1 deg) centered on a declination of -52.7 deg, while 175 hours were devoted to a 280 deg^2 stripe (4.5 hours by 4.8 deg) centered at the celestial equator. We discuss sources of statistical and systematic noise, calibration, telescope pointing, and data selection. Out of 1260 survey hours and 1024 detectors per array, 816 hours and 593 effective detectors remain after data selection for this frequency band, yielding a 38% survey efficiency. The total sensitivity in 2008, determined from the noise level between 5 Hz and 20 Hz in the time-ordered data stream (TOD), is 32 micro-Kelvin sqrt{s} in CMB units. Atmospheric brightness fluctuations constitute the main contaminant in the data and dominate the detector noise covariance at low frequencies in the TOD. The maps were made by solving the least-squares problem using the Preconditioned Conjugate Gradient method, incorporating the details of the detector and noise correlations. Cross-correlation with WMAP sky maps, as well as analysis from simulations, reveal that our maps are unbiased at multipoles ell > 300. This paper accompanies the public release of the 148 GHz southern stripe maps from 2008. The techniques described here will be applied to future maps and data releases.

Using high-resolution microwave sky maps made by the Atacama Cosmology Telescope, we for the first time present strong evidence for motions of galaxy clusters and groups via microwave background temperature distortions due to the kinematic Sunyaev-Zel'dovich effect. Galaxy clusters are identified by their constituent luminous galaxies observed by the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. We measure the mean pairwise momentum of clusters, with a probability of the signal being due to random errors of 0.002, and the signal is consistent with the growth of cosmic structure in the standard model of cosmology.

We perform a combined fit to angular power spectra of unresolved infrared (IR) point sources from the Planck satellite (at 217, 353, 545 and 857 GHz, over angular scales 100 < l < 2200), the Balloon-borne Large-Aperture Submillimeter Telescope (BLAST; 250, 350 and 500 um; 1000 < l < 9000), and from correlating BLAST and Atacama Cosmology Telescope (ACT; 148 and 218 GHz) maps. We find that the clustered power over the range of angular scales and frequencies considered is well fit by a simple power law of the form C_l \propto l^-n with n = 1.25 +/- 0.06. While the IR sources are understood to lie at a range of redshifts, with a variety of dust properties, we find that the frequency dependence of the clustering power can be described by the square of a modified blackbody, nu^beta B(nu,T_eff), with a single emissivity index beta = 2.20 +/- 0.07 and effective temperature T_eff = 9.7 K. Our predictions for the clustering amplitude are consistent with existing ACT and South Pole Telescope results at around 150 and 220 GHz, as is our prediction for the effective dust spectral index, which we find to be alpha_150-220 = 3.68 +/- 0.07 between 150 and 220 GHz. Our constraints on the clustering shape and frequency dependence can be used to model the IR clustering as a contaminant in Cosmic Microwave Background anisotropy measurements. The combined Planck and BLAST data also rule out a linear bias clustering model.

We present an update to the search for a non-trivial topology of the universe by searching for matching circle pairs in the cosmic microwave background using the WMAP 7 year data release. We extend the exisiting bounds to encompass a wider range of possible topologies by searching for matching circle pairs with opening angles 10 degree < \alpha < 90 degree and separation angles 11 degree < \theta < 180 degree. The extended search reveal two small anomalous regions in the CMB sky. Numerous pairs of well-matched circles are found where both circles pass through one or the other of those regions. As this is not the signature of any known manifold, but is a likely consequence of contamination in those sky regions, we repeat the search excluding circle pairs where both pass through either of the two regions. We then find no statistically significant pairs of matched circles, and so no hints of a non-trivial topology. The absence of matched circles increases the lower limit on the length of the shortest closed null geodesic that self-intersects at our location in the universe (equivalently the injectivity radius at our location) to 98.5% of the diameter of the last scattering surface or approximately 26 Gpc. It extends the limit to any manifolds in which the intersecting arcs of said geodesic form an angle greater than 10^o.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

[ABRIDGED] Inhomogeneous reionization acts as a source of arcminute-scale anisotropies in the cosmic microwave background (CMB), the most important of which is the kinetic Sunyaev-Zel'dovich (kSZ) effect. Observational efforts with the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT) are poised to detect this signal for the first time. Indeed, recent SPT measurements place a bound on the dimensionless kSZ power spectrum around a multipole of l=3000 of P_tot < 2.8 (6) micro K^2 at 95% C.L., by ignoring (allowing) correlations between the thermal Sunyaev-Zel'dovich (tSZ) effect and the cosmic infrared background (CIB). To interpret these and upcoming observations, we compute the kSZ signal from a suite of ~ 100 reionization models using the publicly available code 21cmFAST. Our physically motivated reionization models are parameterized by the ionizing efficiency of high-redshift galaxies, the minimum virial temperature of halos capable of hosting stars, and the ionizing photon mean free path. We predict the contribution of patchy reionization to be P_patchy = 1.5-3.5 micro K^2. Therefore, even when conservatively adopting a low estimate of the post-reionization signal, P_OV ~ 2 micro K^2, none of our models are consistent with the aggressive 2sigma SPT bound that does not include correlations. This implies that either: (i) the early stages of reionization occurred in a much more homogeneous manner than suggested by the stellar-driven scenarios we explore, such as would be the case if, e.g., very high energy X-rays or exotic particles contributed significantly; and/or (ii) that there is a significant correlation between the CIB and the tSZ. On the other hand, the conservative SPT bound is compatible with all of our models, and is on the boarder of constraining reionization.

We present a detection of the Sunyaev-Zel'dovich (SZ) decrement associated with the Luminous Red Galaxy (LRG) sample of the Sloan Digital Sky Survey. The SZ data come from 148 GHz maps of the equatorial region made by the Atacama Cosmology Telescope (ACT). The LRG sample is divided by luminosity into four bins, and estimates for the central Sunyaev-Zel'dovich temperature decrement are calculated through a stacking process. We detect and account for a bias of the SZ signal due to weak radio sources. We use numerical simulations to relate the observed decrement to Y200 and clustering properties to relate the galaxy luminosity bins to mass. We also use a relation between brightest cluster galaxy luminosity and cluster mass based on stacked gravitational lensing measurements to estimate the characteristic halo masses. The masses are found to be around 1e14 M_sun.

We present a detailed analysis from new multi-wavelength observations of the exceptional galaxy cluster ACT-CL J0102-4915 "El Gordo," likely the most massive, hottest, most X-ray luminous and brightest Sunyaev-Zeldovich (SZ) effect cluster known at z>0.6. The Atacama Cosmology Telescope collaboration discovered El Gordo as the most significant SZ decrement in a sky survey area of 755 deg^2. Our VLT/FORS2 spectra of 89 member galaxies yield a cluster redshift, z=0.870, and velocity dispersion, s=1321+/-106 km/s. Our Chandra observations reveal a hot and X-ray luminous system with an integrated temperature of Tx=14.5+/-1.0 keV and 0.5-2.0 keV band luminosity of Lx=(2.19+/-0.11)x10^45 h70^-2 erg/s. We obtain several statistically consistent cluster mass estimates; using mass scaling relations with velocity dispersion, X-ray Yx, and integrated SZ, we estimate a cluster mass of M200a=(2.16+/-0.32)x10^15 M_sun/h70. The Chandra and VLT/FORS2 optical data also reveal that El Gordo is undergoing a major merger between components with a mass ratio of approximately 2 to 1. The X-ray data show significant temperature variations from a low of 6.6+/-0.7 keV at the merging low-entropy, high-metallicity, cool core to a high of 22+/-6 keV. We also see a wake in the X-ray surface brightness caused by the passage of one cluster through the other. Archival radio data at 843 MHz reveal diffuse radio emission that, if associated with the cluster, indicates the presence of an intense double radio relic, hosted by the highest redshift cluster yet. El Gordo is possibly a high-redshift analog of the famous Bullet Cluster. Such a massive cluster at this redshift is rare, although consistent with the standard L-CDM cosmology in the lower part of its allowed mass range. Massive, high-redshift mergers like El Gordo are unlikely to be reproduced in the current generation of numerical N-body cosmological simulations.

For decades, cosmologists have been using galaxies to trace the large-scale distribution of matter. At present, the largest source of systematic uncertainty in this analysis is the challenge of modeling the complex relationship between galaxy redshift and the distribution of dark matter. If all galaxies sat in the centers of halos, there would be minimal Finger-of-God (FoG) effects and a simple relationship between the galaxy and matter distributions. However, many galaxies, even some of the luminous red galaxies (LRGs), do not lie in the centers of halos. Because the galaxy-galaxy lensing is also sensitive to the off-centered galaxies, we show that we can use the lensing measurements to determine the amplitude of this effect and to determine the expected amplitude of FoG effects. We develop an approach for using the lensing data to model how the FoG suppresses the power spectrum amplitudes and show that the current data implies a 30% suppression at wavenumber k=0.2h/Mpc. Our analysis implies that it is important to complement a spectroscopic survey with an imaging survey with sufficient depth and wide field coverage. Joint imaging and spectroscopic surveys allow a robust, unbiased use of the power spectrum amplitude information: it improves the marginalized error of growth rate fg=dln D/dln a by up to a factor of 2 over a wide range of redshifts z<1.4. We also find that the dark energy equation-of-state parameter, w0, and the neutrino mass, fnu, can be unbiasedly constrained by combining the lensing information, with an improvement of 10--25% compared to a spectroscopic survey without lensing calibration.

We introduce to astrophysics the threshold probability functions S_2, C_2, and D_2 first derived by \citet{torq+88}, which effectively samples the flux probability distribution (PDF) of the Lya forest at different spatial scales. These statistics are tested on mock Lya forest spectra based on various toy models for HeII reionization, with homogeneous models with various temperature-density relations as well as models with temperature inhomogeneities. These mock samples have systematics and noise added to simulate the latest Sloan Digital Sky Survey Data Release 7 (SDSS DR7) data. We find that the flux PDF from SDSS DR7 can be used to constrain the temperature-density relation $\gamma$ (where $T \propto (1 + \Delta)^{\gamma-1}$) of the intergalactic medium (IGM) at z=2.5 to a precision of $\Delta \gamma = 0.2$ at $\sim 4\sigma$ confidence. The flux PDF is degenerate to temperature inhomogeneities in the IGM arising from HeII reionization, but we find S_2 can detect these inhomogeneities at $\sim 3 \sigma$, with the assumption that the flux continuum of the Lya forest can be determined to 9% accuracy, approximately the error from current fitting methods. If the quasar continuum can be determined to 3% accuracy, then S_2 is capable of constraining the characteristic scale of temperature inhomogeneities, with $\sim 4 \sigma$ differentiation between toy models with hot bubble radii of 50 Mpc/h and 25 Mpc/h, comoving.

Continuum fitting is an important aspect of Lyman-alpha forest science, since errors in the estimated optical depths scale with the fractional continuum error. However, traditional methods of estimating continua in noisy and moderate-resolution spectra (S/N < 10 pixel^-1 and R ~ 2000, respectively, such as SDSS) using power-law extrapolation or the mean spectrum, achieve no better than ~ 10-15% RMS accuracy. To improve on this, we introduce mean-flux regulated/principal component analysis (MF-PCA) continuum fitting. In this technique, PCA fitting is carried out redwards of the quasar Lyman-alpha line in order to provide a prediction for the shape of the Lyman-alpha forest continuum. The slope and amplitude of this continuum prediction is then corrected using external constraints for the Lyman-alpha forest mean-flux. From tests on mock spectra, we find that MF-PCA reduces the errors to 8% RMS in S/N ~ 2 spectra, and < 5% RMS in spectra with S/N > 5. The residual Fourier power in the continuum is decreased by a factor of a few in comparison with dividing by the mean continuum, enabling Lyman-alpha flux power spectrum measurements to be extended to ~2x larger scales. Using this new technique, we make available continuum fits for 12,069 z>2.3 Lyman-alpha forest spectra from SDSS DR7 for use by the community. This technique is also applicable to future releases of the ongoing BOSS survey, which is obtaining spectra for ~ 150,000 Lyman-alpha forest spectra at low signal-to-noise (S/N ~ 2).

We present follow-up observations with the Sunyaev-Zel'dovich Array (SZA) of optically-confirmed galaxy clusters found in the equatorial survey region of the Atacama Cosmology Telescope (ACT): ACT-CL J0022-0036, ACT-CL J2051+0057, and ACT-CL J2337+0016. ACT-CL J0022-0036 is a newly-discovered, massive (10^15 Msun), high-redshift (z=0.81) cluster revealed by ACT through the Sunyaev-Zel'dovich effect (SZE). Deep, targeted observations with the SZA allow us to probe a broader range of cluster spatial scales, better disentangle cluster decrements from radio point source emission, and derive more robust integrated SZE flux and mass estimates than we can with ACT data alone. For the two clusters we detect with the SZA we compute integrated SZE signal and derive masses from the SZA data only. ACT-CL J2337+0016, also known as Abell 2631, has archival Chandra data that allow an additional X-ray-based mass estimate. Optical richness is also used to estimate cluster masses and shows good agreement with the SZE and X-ray-based estimates. Based on the point sources detected by the SZA in these three cluster fields and an extrapolation to ACT's frequency, we estimate that point sources could be contaminating the SZE decrement at the <= 20% level for some fraction of clusters.

The conjecture that the ancient globular clusters (GCs) formed at the center of their own dark matter halos was first proposed by Peebles (1984), and has recently been revived to explain the puzzling abundance patterns observed within many GCs. In this paper we demonstrate that the outer stellar density profile of isolated GCs is very sensitive to the presence of an extended dark halo. The GCs NGC 2419, located at 90 kpc from the center of our Galaxy, and MGC1, located at ~200 kpc from the center of M31, are ideal laboratories for testing the scenario that GCs formed at the centers of massive dark halos. Comparing analytic models to observations of these GCs, we conclude that these GCs cannot be embedded within dark halos with a virial mass greater than 10^6 Msun, or, equivalently, the dark matter halo mass-to-stellar mass ratio must be Mdm/M_*<1. If these GCs have indeed orbited within weak tidal fields throughout their lifetimes, then these limits imply that these GCs did not form within their own dark halos. Recent observations of an extended stellar halo in the GC NGC 1851 are also interpreted in the context of our analytic models. Implications of these results for the formation of GCs are briefly discussed.

For the first time, measurements of the cosmic microwave background radiation (CMB) alone favor cosmologies with $w=-1$ dark energy over models without dark energy at a 3.2-sigma level. We demonstrate this by combining the CMB lensing deflection power spectrum from the Atacama Cosmology Telescope with temperature and polarization power spectra from the Wilkinson Microwave Anisotropy Probe. The lensing data break the geometric degeneracy of different cosmological models with similar CMB temperature power spectra. Our CMB-only measurement of the dark energy density $\Omega_\Lambda$ confirms other measurements from supernovae, galaxy clusters and baryon acoustic oscillations, and demonstrates the power of CMB lensing as a new cosmological tool.

We report the first detection of the gravitational lensing of the cosmic microwave background through a measurement of the four-point correlation function in the temperature maps made by the Atacama Cosmology Telescope. We verify our detection by calculating the levels of potential contaminants and performing a number of null tests. The resulting convergence power spectrum at 2-degree angular scales measures the amplitude of matter density fluctuations on comoving length scales of around 100 Mpc at redshifts around 0.5 to 3. The measured amplitude of the signal agrees with Lambda Cold Dark Matter cosmology predictions. Since the amplitude of the convergence power spectrum scales as the square of the amplitude of the density fluctuations, the 4-sigma detection of the lensing signal measures the amplitude of density fluctuations to 12%.

We report on twenty-three clusters detected blindly as Sunyaev-Zel'dovich (SZ) decrements in a 148 GHz, 455 square-degree map of the southern sky made with data from the Atacama Cosmology Telescope 2008 observing season. All SZ detections announced in this work have confirmed optical counterparts. Ten of the clusters are new discoveries. One newly discovered cluster, ACT-CL J0102-4915, with a redshift of 0.75 (photometric), has an SZ decrement comparable to the most massive systems at lower redshifts. Simulations of the cluster recovery method reproduce the sample purity measured by optical follow-up. In particular, for clusters detected with a signal-to-noise ratio greater than six, simulations are consistent with optical follow-up that demonstrated this subsample is 100% pure. The simulations further imply that the total sample is 80% complete for clusters with mass in excess of 6x10^14 solar masses referenced to the cluster volume characterized by five hundred times the critical density. The Compton y -- X-ray luminosity mass comparison for the eleven best detected clusters visually agrees with both self-similar and non-adiabatic, simulation-derived scaling laws.

We report on extragalactic sources detected in a 455 square-degree map of the southern sky made with data at a frequency of 148 GHz from the Atacama Cosmology Telescope 2008 observing season. We provide a catalog of 157 sources with flux densities spanning two orders of magnitude: from 15 to 1500 mJy. Comparison to other catalogs shows that 98% of the ACT detections correspond to sources detected at lower radio frequencies. Three of the sources appear to be associated with the brightest cluster galaxies of low redshift X-ray selected galaxy clusters. Estimates of the radio to mm-wave spectral indices and differential counts of the sources further bolster the hypothesis that they are nearly all radio sources, and that their emission is not dominated by re-emission from warm dust. In a bright (>50 mJy) 148 GHz-selected sample with complete cross-identifications from the Australia Telescope 20 GHz survey, we observe an average steepening of the spectra between 5, 20, and 148 GHz with median spectral indices of $\alpha_{\rm 5-20} = -0.07 \pm 0.06$, $\alpha_{\rm 20-148} = -0.39 \pm0.04$, and $\alpha_{\rm 5-148} = -0.20 \pm 0.03$. When the measured spectral indices are taken into account, the 148 GHz differential source counts are consistent with previous measurements at 30 GHz in the context of a source count model dominated by radio sources. Extrapolating with an appropriately rescaled model for the radio source counts, the Poisson contribution to the spatial power spectrum from synchrotron-dominated sources with flux density less than 20 mJy is $C^{\rm Sync} = (2.8 \pm 0.3) \times 10^{-6} \micro\kelvin^2$.

We present constraints on the primordial power spectrum of adiabatic fluctuations using data from the 2008 Southern Survey of the Atacama Cosmology Telescope (ACT). The angular resolution of ACT provides sensitivity to scales beyond \ell = 1000 for resolution of multiple peaks in the primordial temperature power spectrum, which enables us to probe the primordial power spectrum of adiabatic scalar perturbations with wavenumbers up to k \simeq 0.2 Mpc^{-1}. We find no evidence for deviation from power-law fluctuations over two decades in scale. Matter fluctuations inferred from the primordial temperature power spectrum evolve over cosmic time and can be used to predict the matter power spectrum at late times; we illustrate the overlap of the matter power inferred from CMB measurements (which probe the power spectrum in the linear regime) with existing probes of galaxy clustering, cluster abundances and weak lensing constraints on the primordial power. This highlights the range of scales probed by current measurements of the matter power spectrum.

The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission to measure the gravity-wave signature of primordial inflation through its distinctive imprint on the linear polarization of the cosmic microwave background. The instrument consists of a polarizing Michelson interferometer configured as a nulling polarimeter to measure the difference spectrum between orthogonal linear polarizations from two co-aligned beams. Either input can view the sky or a temperature-controlled absolute reference blackbody calibrator. PIXIE will map the absolute intensity and linear polarization (Stokes I, Q, and U parameters) over the full sky in 400 spectral channels spanning 2.5 decades in frequency from 30 GHz to 6 THz (1 cm to 50 um wavelength). Multi-moded optics provide background-limited sensitivity using only 4 detectors, while the highly symmetric design and multiple signal modulations provide robust rejection of potential systematic errors. The principal science goal is the detection and characterization of linear polarization from an inflationary epoch in the early universe, with tensor-to-scalar ratio r < 10^{-3} at 5 standard deviations. The rich PIXIE data set will also constrain physical processes ranging from Big Bang cosmology to the nature of the first stars to physical conditions within the interstellar medium of the Galaxy.

We present measurements of the auto- and cross-frequency correlation power spectra of the cosmic (sub)millimeter background at: 250, 350, and 500 um (1200, 860, and 600 GHz) from observations made with the Balloon-borne Large Aperture Submillimeter Telescope, BLAST; and at 1380 and 2030 um (218 and 148 GHz) from observations made with the Atacama Cosmology Telescope, ACT. The overlapping observations cover 8.6 deg^2 in an area relatively free of Galactic dust near the south ecliptic pole (SEP). The ACT bands are sensitive to radiation from the CMB, the Sunyaev-Zel'dovich (SZ) effect from galaxy clusters, and to emission by radio and dusty star-forming galaxies (DSFGs), while the dominant contribution to the BLAST bands is from DSFGs. We confirm and extend the BLAST analysis of clustering with an independent pipeline, and also detect correlations between the ACT and BLAST maps at over 25sigma significance, which we interpret as a detection of the DSFGs in the ACT maps. In addition to a Poisson component in the cross-frequency power spectra, we detect a clustered signal at >4sigma, and using a model for the DSFG evolution and number counts, we successfully fit all our spectra with a linear clustering model and a bias that depends only on redshift and not on scale. Finally, the data are compared to, and generally agree with, phenomenological models for the DSFG population. This study represents a first of its kind, and demonstrates the constraining power of the cross-frequency correlation technique to constrain models for the DSFGs. Similar analyses with more data will impose tight constraints on future models.

(Abridged) A simple six-parameter LCDM model provides a successful fit to WMAP data, both when the data are analyzed alone and in combination with other cosmological data. Even so, it is appropriate to search for any hints of deviations from the now standard model of cosmology, which includes inflation, dark energy, dark matter, baryons, and neutrinos. The cosmological community has subjected the WMAP data to extensive and varied analyses. While there is widespread agreement as to the overall success of the six-parameter LCDM model, various "anomalies" have been reported relative to that model. In this paper we examine potential anomalies and present analyses and assessments of their significance. In most cases we find that claimed anomalies depend on posterior selection of some aspect or subset of the data. Compared with sky simulations based on the best fit model, one can select for low probability features of the WMAP data. Low probability features are expected, but it is not usually straightforward to determine whether any particular low probability feature is the result of the a posteriori selection or of non-standard cosmology. We examine in detail the properties of the power spectrum with respect to the LCDM model. We examine several potential or previously claimed anomalies in the sky maps and power spectra, including cold spots, low quadrupole power, quadropole-octupole alignment, hemispherical or dipole power asymmetry, and quadrupole power asymmetry. We conclude that there is no compelling evidence for deviations from the LCDM model, which is generally an acceptable statistical fit to WMAP and other cosmological data.

(Abridged) We present WMAP seven-year observations of bright sources which are often used as calibrators at microwave frequencies. Ten objects are studied in five frequency bands (23 - 94 GHz): the outer planets (Mars, Jupiter, Saturn, Uranus and Neptune) and five fixed celestial sources (Cas A, Tau A, Cyg A, 3C274 and 3C58). The seven-year analysis of Jupiter provides temperatures which are within 1-sigma of the previously published WMAP five-year values, with slightly tighter constraints on variability with orbital phase, and limits (but no detections) on linear polarization. Scaling factors are provided which, when multiplied by the Wright Mars thermal model predictions at 350 micron, reproduce WMAP seasonally averaged observations of Mars within ~2%. An empirical model is described which fits brightness variations of Saturn due to geometrical effects and can be used to predict the WMAP observations to within 3%. Seven-year mean temperatures for Uranus and Neptune are also tabulated. Uncertainties in Uranus temperatures are 3%-4% in the 41, 61 and 94 GHz bands; the smallest uncertainty for Neptune is ~8% for the 94 GHz band. Intriguingly, the spectrum of Uranus appears to show a dip at ~30 GHz of unidentified origin, although the feature is not of high statistical significance. Flux densities for the five selected fixed celestial sources are derived from the seven-year WMAP sky maps, and are tabulated for Stokes I, Q and U, along with polarization fraction and position angle. Fractional uncertainties for the Stokes I fluxes are typically 1% to 3%. Source variability over the seven-year baseline is also estimated. Significant secular decrease is seen for Cas A and Tau A: our results are consistent with a frequency independent decrease of about 0.53% per year for Cas A and 0.22% per year for Tau A.