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.

We explore in this paper the use of neural networks designed for point-clouds and sets on a new meta-learning task. We present experiments on the astronomical challenge of characterizing the stellar population of stellar streams. Stellar streams are elongated structures of stars in the outskirts of the Milky Way that form when a (small) galaxy breaks up under the Milky Way's gravitational force. We consider that we obtain, for each stream, a small 'support set' of stars that belongs to this stream. We aim to predict if the other stars in that region of the sky are from that stream or not, similar to one-class classification. Each "stream task" could also be transformed into a binary classification problem in a highly imbalanced regime (or supervised anomaly detection) by using the much bigger set of "other" stars and considering them as noisy negative examples. We propose to study the problem in the meta-learning regime: we expect that we can learn general information on characterizing a stream's stellar population by meta-learning across several streams in a fully supervised regime, and transfer it to new streams using only positive supervision. We present a novel use of Deep Sets, a model developed for point-cloud and sets, trained in a meta-learning fully supervised regime, and evaluated in a one-class classification setting. We compare it against Random Forests (with and without self-labeling) in the classic setting of binary classification, retrained for each task. We show that our method outperforms the Random-Forests even though the Deep Sets is not retrained on the new tasks, and accesses only a small part of the data compared to the Random Forest. We also show that the model performs well on a real-life stream when including additional fine-tuning.

We combine hierarchical Bayesian modeling with a flow-based deep generative network, in order to demonstrate that one can efficiently constraint numerical gravitational wave (GW) population models at a previously intractable complexity. Existing techniques for comparing data to simulation,such as discrete model selection and Gaussian process regression, can only be applied efficiently to moderate-dimension data. This limits the number of observable (e.g. chirp mass, spins.) and hyper-parameters (e.g. common envelope efficiency) one can use in a population inference. In this study, we train a network to emulate a phenomenological model with 6 observables and 4 hyper-parameters, use it to infer the properties of a simulated catalogue and compare the results to the phenomenological model. We find that a 10-layer network can emulate the phenomenological model accurately and efficiently. Our machine enables simulation-based GW population inferences to take on data at a new complexity level.

Observationally, it has been reported that the densest stellar system in the Universe does not exceed a maximum stellar surface density, $\Sigma^{\max}_{*}$ = $3\times10^5$M$_{\odot}$pc$^{-2}$, throughout a wide physical scale ranging from star cluster to galaxy. This suggests there exists a fundamental physics which regulates the star formation and stellar density. However, factors that determine this maximum limit are not clear. In this study, we show that $\Sigma^{\max}_{*}$ of galaxies is not a constant as previous work reported, but actually depends on the stellar mass. We select galaxy sample from the Sloan Digital Sky Survey Data Release 12 at $z=0.01-0.5$. In contrast to a constant maximum predicted by theoretical models, $\Sigma^{\max}_{*}$ strongly depends on stellar mass especially for less massive galaxies with $\sim10^{10}$M$_{\odot}$. We also found that a majority of high-$\Sigma_{*}$ galaxies show red colours and low star-formation rates. These galaxies probably reach the $\Sigma^{\max}_{*}$ as a consequence of the galaxy evolution from blue star forming to red quiescent by quenching star formation. One possible explanation of the stellar-mass dependency of $\Sigma^{\max}_{*}$ is a mass dependent efficiency of stellar feedback. The stellar feedback could be relatively more efficient in a shallower gravitational potential, which terminates star formation quickly before the stellar system reaches a high stellar density.

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.

We present the physical properties of AKARI sources without optical counterparts in optical images from the Hyper Suprime-Cam (HSC) on the Subaru telescope. Using the AKARI infrared (IR) source catalog and HSC optical catalog, we select 583 objects that do not have HSC counterparts in the AKARI North Ecliptic Pole (NEP) wide survey field ($\sim 5$ deg$^{2}$). Because the HSC limiting magnitude is deep ($g_{\rm AB}$ $\sim 28.6$), these are good candidates for extremely red star-forming galaxies (SFGs) and/or active galactic nuclei (AGNs), possibly at high redshifts. We compile multi-wavelength data out to 500 $\mu$m and use it for Spectral Energy Distribution (SED) fitting with CIGALE to investigate the physical properties of AKARI galaxies without optical counterparts. We also compare their physical quantities with AKARI mid-IR selected galaxies with HSC counterparts. The estimated redshifts of AKARI objects without HSC counterparts range up to $z\sim 4$, significantly higher than that of AKARI objects with HSC counterparts. We find that: (i) 3.6 $-$ 4.5 $\mu$m color, (ii) AGN luminosity, (iii) stellar mass, (iv) star formation rate, and (v) $V$-band dust attenuation in the interstellar medium of AKARI objects without HSC counterparts are systematically larger than those of AKARI objects with counterparts. These results suggest that our sample includes luminous, heavily dust-obscured SFGs/AGNs at $z\sim 1-4$ that are missed by previous optical surveys, providing very interesting targets for the coming James Webb Space Telescope era.

This paper introduces the Wavelet Phase Harmonics (WPH) statistics. They are interpretable low-dimensional statistics which describe 2D non-Gaussian density fields. These statistics are built from WPH moments, which have been recently introduced in data science and machine learning community. In this paper, we applied WPH statistics to projected matter density fields from the Quijote N-body simulations. We find by computing the Fisher information matrix, that the WPH statistics can place more stringent constraints on 5 cosmological parameters when compared to the combination of power-spectrum and bi-spectrum. We also use the WPH statistics to successfully generate from a maximum entropy model new 2D density fields that reproduce the PDF, mean, power-spectrum, bispectrum and the Minkowski functionals of the input density fields. While separate methods have been proven very efficient for parameter estimations and statistical syntheses for for large scale structure, WPH statistics are the first statistics that can both achieve a more stringent cosmological parameter constraint, and produce a sufficiently accurate simulation of the Universe while being interpretable.

The Tri-Agency Cosmological Simulations (TACS) Task Force was formed when Program Managers from the Department of Energy (DOE), the National Aeronautics and Space Administration (NASA), and the National Science Foundation (NSF) expressed an interest in receiving input into the cosmological simulations landscape related to the upcoming DOE/NSF Vera Rubin Observatory (Rubin), NASA/ESA's Euclid, and NASA's Wide Field Infrared Survey Telescope (WFIRST). The Co-Chairs of TACS, Katrin Heitmann and Alina Kiessling, invited community scientists from the USA and Europe who are each subject matter experts and are also members of one or more of the surveys to contribute. The following report represents the input from TACS that was delivered to the Agencies in December 2018.

Fast radio bursts (FRBs) are mysterious radio bursts with a time scale of approximately milliseconds. Two populations of FRB, namely repeating and non-repeating FRBs, are observationally identified. However, the differences between these two and their origins are still cloaked in mystery. Here we show the time-integrated luminosity-duration ($L_{\nu}$-$w_{\rm int,rest}$) relations and luminosity functions (LFs) of repeating and non-repeating FRBs in the FRB Catalogue project. These two populations are obviously separated in the $L_{\nu}$-$w_{\rm int,rest}$ plane with distinct LFs, i.e., repeating FRBs have relatively fainter $L_{\nu}$ and longer $w_{\rm int,rest}$ with a much lower LF. In contrast with non-repeating FRBs, repeating FRBs do not show any clear correlation between $L_{\nu}$ and $w_{\rm int,rest}$. These results suggest essentially different physical origins of the two. The faint ends of the LFs of repeating and non-repeating FRBs are higher than volumetric occurrence rates of neutron-star mergers and accretion-induced collapse (AIC) of white dwarfs, and are consistent with those of soft gamma-ray repeaters (SGRs), type Ia supernovae, magnetars, and white-dwarf mergers. This indicates two possibilities: either (i) faint non-repeating FRBs originate in neutron-star mergers or AIC and are actually repeating during the lifetime of the progenitor, or (ii) faint non-repeating FRBs originate in any of SGRs, type Ia supernovae, magnetars, and white-dwarf mergers. The bright ends of LFs of repeating and non-repeating FRBs are lower than any candidates of progenitors, suggesting that bright FRBs are produced from a very small fraction of the progenitors regardless of the repetition. Otherwise, they might originate in unknown progenitors.

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 Simons Observatory (SO) is an upcoming polarization-sensitive Cosmic Microwave Background (CMB) experiment on the Cerro Toco Plateau (Chile) with large overlap with other optical and infrared surveys (e.g., DESI, LSST, HSC). To enable the readout of \bigO(10,000) detectors in each of the four telescopes of SO, we will employ the microwave SQUID multiplexing technology. With a targeted multiplexing factor of \bigO{(1,000)}, microwave SQUID multiplexing has never been deployed on the scale needed for SO. Here we present the design of the cryogenic coaxial cable and RF component chain that connects room temperature readout electronics to superconducting resonators that are coupled to Transition Edge Sensor bolometers operating at sub-Kelvin temperatures. We describe design considerations including cryogenic RF component selection, system linearity, noise, and thermal power dissipation.

One of the most important checks for systematic errors in CMB studies is the cross correlation of maps made by independent experiments. In this paper we report on the cross correlation between maps from the Atacama B-mode Search (ABS) and Atacama Cosmology Telescope (ACT) experiments in both temperature and polarization. These completely different measurements have a clear correlation with each other and with the Planck satellite in both the EE and TE spectra at $\ell<400$ over the roughly $1100$ deg$^2$ common to all three. The TB, EB, and BB cross spectra are consistent with noise. Exploiting such cross-correlations will be important for future experiments operating in Chile that aim to probe the $30<\ell<8,000$ range.

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.

The Simons Observatory (SO) is an upcoming cosmic microwave background (CMB) experiment located on Cerro Toco, Chile, that will map the microwave sky in temperature and polarization in six frequency bands spanning 27 to 285 GHz. SO will consist of one 6-meter Large Aperture Telescope (LAT) fielding $\sim$30,000 detectors and an array of three 0.42-meter Small Aperture Telescopes (SATs) fielding an additional 30,000 detectors. This synergy will allow for the extremely sensitive characterization of the CMB over angular scales ranging from an arcmin to tens of degrees, enabling a wide range of scientific output. Here we focus on the SATs targeting degree angular scales with successive dichroic instruments observing at Mid-Frequency (MF: 93/145 GHz), Ultra-High-Frequency (UHF: 225/285 GHz), and Low-Frequency (LF: 27/39 GHz). The three SATs will be able to map $\sim$10% of the sky to a noise level of 2 $\mu$K-arcmin when combining 93 and 145 GHz. The multiple frequency bands will allow the CMB to be separated from galactic foregrounds (primarily synchrotron and dust), with the primary science goal of characterizing the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma \left(r\right) \approx 0.003$.

The Simons Observatory is building both large (6 m) and small (0.5 m) aperture telescopes in the Atacama desert in Chile to observe the cosmic microwave background (CMB) radiation with unprecedented sensitivity. Simons Observatory telescopes in total will use over 60,000 transition edge sensor (TES) detectors spanning center frequencies between 27 and 285 GHz and operating near 100 mK. TES devices have been fabricated for the Simons Observatory by NIST, Berkeley, and HYPRES/SeeQC corporation. Iterations of these devices have been tested cryogenically in order to inform the fabrication of further devices, which will culminate in the final TES designs to be deployed in the field. The detailed design specifications have been independently iterated at each fabrication facility for particular detector frequencies. We present test results for prototype devices, with emphasis on NIST high frequency detectors. A dilution refrigerator was used to achieve the required temperatures. Measurements were made both with 4-lead resistance measurements and with a time domain Superconducting Quantum Interference Device (SQUID) multiplexer system. The SQUID readout measurements include analysis of current vs voltage (IV) curves at various temperatures, square wave bias step measurements, and detector noise measurements. Normal resistance, superconducting critical temperature, saturation power, thermal and natural time constants, and thermal properties of the devices are extracted from these measurements.

The Simons Observatory will consist of a single large (6 m diameter) telescope and a number of smaller (0.5 m diameter) refracting telescopes designed to measure the polarization of the Cosmic Microwave Background to unprecedented accuracy. The large aperture telescope is the same design as the CCAT-prime telescope, a modified Crossed Dragone design with a field-of-view of over 7.8 degrees diameter at 90 GHz. This paper presents an overview of the cold reimaging optics for this telescope and what drove our choice of 350-400 mm diameter silicon lenses in a 2.4 m cryostat over other possibilities. We will also consider the future expandability of this design to CMB Stage-4 and beyond.

Advanced ACTPol is the second generation polarization-sensitive upgrade to the $6\, \rm m$ aperture Atacama Cosmology Telescope (ACT), which increased detector count and frequency coverage compared to the previous ACTPol receiver. Advanced ACTPol utilizes a new two-stage time-division multiplexing readout architecture based on superconducting quantum interference devices (SQUIDs) to achieve a multiplexing factor as high as 64 (rows), fielding a 2,012 detector camera at 150/220 GHz and two 90/150 GHz cameras containing 1,716 detectors each. In a time domain system, aliasing introduces noise to the readout. In this work we present a figure of merit to measure this noise contribution and present measurements of the aliased noise fraction of the Advanced ACTPol receiver as deployed.

We constrain gas inflow speeds in star-forming galaxies with color gradients consistent with inside-out disk growth. Our method combines new measurements of disk orientation with previously described circumgalactic absorption in background quasar spectra. Two quantities, a position angle and an axis ratio, describe the projected shape of each galactic disk on the sky, leaving an ambiguity about which side of the minor axis is tipped toward the observer. This degeneracy regarding the 3D orientation of disks has compromised previous efforts to measure gas inflow speeds. We present HST and Keck/LGSAO imaging that resolves the spiral structure in five galaxies at redshift $z\approx0.2$. We determine the sign of the disk inclination for four galaxies, under the assumption that spiral arms trail the rotation. We project models for both radial infall in the disk plane and circular orbits onto each quasar sightline. We compare the resulting line-of-sight velocities to the observed velocity range of Mg II absorption in spectra of background quasars, which intersect the disk plane at radii between 69 and 115 kpc. For two sightlines, we constrain the maximum radial inflow speeds as 30-40 km s$^{-1}$. We also rule out a velocity component from radial inflow in one sightline, suggesting that the structures feeding gas to these growing disks do not have unity covering factor. We recommend appropriate selection criteria for building larger samples of galaxy--quasar pairs that produce orientations sensitive to constraining inflow properties.

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 introduce an approach for imposing physically motivated inductive biases on graph networks to learn interpretable representations and improved zero-shot generalization. Our experiments show that our graph network models, which implement this inductive bias, can learn message representations equivalent to the true force vector when trained on n-body gravitational and spring-like simulations. We use symbolic regression to fit explicit algebraic equations to our trained model's message function and recover the symbolic form of Newton's law of gravitation without prior knowledge. We also show that our model generalizes better at inference time to systems with more bodies than had been experienced during training. Our approach is extensible, in principle, to any unknown interaction law learned by a graph network, and offers a valuable technique for interpreting and inferring explicit causal theories about the world from implicit knowledge captured by deep learning.

The Wide-Field InfraRed Space Telescope (WFIRST) will be capable of delivering precise astrometry for faint sources over the enormous field of view of its main camera, the Wide-Field Imager (WFI). This unprecedented combination will be transformative for the many scientific questions that require precise positions, distances, and velocities of stars. We describe the expectations for the astrometric precision of the WFIRST WFI in different scenarios, illustrate how a broad range of science cases will see significant advances with such data, and identify aspects of WFIRST's design where small adjustments could greatly improve its power as an astrometric instrument.

Cosmological simulations play an important role in the interpretation of astronomical data, in particular in comparing observed data to our theoretical expectations. However, to compare data with these simulations, the simulations in principle need to include gravity, magneto-hydrodyanmics, radiative transfer, etc. These ideal large-volume simulations (gravo-magneto-hydrodynamical) are incredibly computationally expensive which can cost tens of millions of CPU hours to run. In this paper, we propose a deep learning approach to map from the dark-matter-only simulation (computationally cheaper) to the galaxy distribution (from the much costlier cosmological simulation). The main challenge of this task is the high sparsity in the target galaxy distribution: space is mainly empty. We propose a cascade architecture composed of a classification filter followed by a regression procedure. We show that our result outperforms a state-of-the-art model used in the astronomical community, and provides a good trade-off between computational cost and prediction accuracy.

Measuring the sum of the three active neutrino masses, $M_\nu$, is one of the most important challenges in modern cosmology. Massive neutrinos imprint characteristic signatures on several cosmological observables in particular on the large-scale structure of the Universe. In order to maximize the information that can be retrieved from galaxy surveys, accurate theoretical predictions in the non-linear regime are needed. Currently, one way to achieve those predictions is by running cosmological numerical simulations. Unfortunately, producing those simulations requires high computational resources -- seven hundred CPU hours for each neutrino mass case. In this work, we propose a new method, based on a deep learning network (U-Net), to quickly generate simulations with massive neutrinos from standard $\Lambda$CDM simulations without neutrinos. We computed multiple relevant statistical measures of deep-learning generated simulations, and conclude that our method accurately reproduces the 3-dimensional spatial distribution of matter down to non-linear scales: $k < 0.7$ h/Mpc. Finally, our method allows us to generate massive neutrino simulations 10,000 times faster than the traditional methods.

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.

Hidden monopole is a plausible dark matter candidate due to its stability, but its direct experimental search is extremely difficult due to feeble interactions with the standard model particles in the minimal form. Then, we introduce an axion, $a$, connecting the hidden monopole and the standard model particles and examine the current limits and future prospects of direct dark matter searches and beam-dump experiments. We find two parameter regions around $m_a = {\cal O}(10)$ MeV, $f_a = {\cal O}(10^{5})$ GeV and $m_a = {\cal O}(100)$ MeV, $f_a = {\cal O}(10^{4})$ GeV where monopole dark matter and the axion are respectively within the reach of the future experiments such as PICO-500 and SHiP. We also note that the hidden photons mainly produced by the axion decay contribute to dark radiation with $\Delta N_{\rm eff} \simeq 0.6$ which can relax the $H_0$ tension.

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.

Matter evolved under influence of gravity from minuscule density fluctuations. Non-perturbative structure formed hierarchically over all scales, and developed non-Gaussian features in the Universe, known as the Cosmic Web. To fully understand the structure formation of the Universe is one of the holy grails of modern astrophysics. Astrophysicists survey large volumes of the Universe and employ a large ensemble of computer simulations to compare with the observed data in order to extract the full information of our own Universe. However, to evolve trillions of galaxies over billions of years even with the simplest physics is a daunting task. We build a deep neural network, the Deep Density Displacement Model (hereafter D$^3$M), to predict the non-linear structure formation of the Universe from simple linear perturbation theory. Our extensive analysis, demonstrates that D$^3$M outperforms the second order perturbation theory (hereafter 2LPT), the commonly used fast approximate simulation method, in point-wise comparison, 2-point correlation, and 3-point correlation. We also show that D$^3$M is able to accurately extrapolate far beyond its training data, and predict structure formation for significantly different cosmological parameters. Our study proves, for the first time, that deep learning is a practical and accurate alternative to approximate simulations of the gravitational structure formation of the Universe.

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.

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.

We present a detection of the splashback feature around galaxy clusters selected using their Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, $r_{\rm sp}$, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that $r_{\rm sp}$ inferred from the observed radial distribution of galaxies is affected by selection effects related to the optical cluster-finding algorithms. We test this possibility by measuring the splashback feature in clusters selected via the SZ effect in data from the South Pole Telescope SZ survey and the Atacama Cosmology Telescope Polarimeter survey. The measurement is accomplished by correlating these clusters with galaxies detected in the Dark Energy Survey Year 3 data. The SZ observable used to select clusters in this analysis is expected to have a tighter correlation with halo mass and to be more immune to projection effects and aperture-induced biases than optically selected clusters. We find that the measured $r_{\rm sp}$ for SZ-selected clusters is consistent with the expectations from simulations, although the small number of SZ-selected clusters makes a precise comparison difficult. In agreement with previous work, when using optically selected redMaPPer clusters, $r_{\rm sp}$ is $\sim$ $2\sigma$ smaller than in the simulations. These results motivate detailed investigations of selection biases in optically selected cluster catalogs and exploration of the splashback feature around larger samples of SZ-selected clusters. Additionally, we investigate trends in the galaxy profile and splashback feature as a function of galaxy color, finding that blue galaxies have profiles close to a power law with no discernible splashback feature, which is consistent with them being on their first infall into the cluster.

Halos and galaxies acquire their angular momentum during the collapse of surrounding large-scale structure. This process imprints alignments between galaxy spins and nearby filaments and sheets. Low mass halos grow by accretion onto filaments, aligning their spins with the filaments, whereas high mass halos grow by mergers along filaments, generating spins perpendicular to the filament. We search for this alignment signal using filaments identified with the "Cosmic Web Reconstruction" algorithm applied to the Sloan Digital Sky Survey Main Galaxy Sample and galaxy spins from the MaNGA integral-field unit survey. MaNGA produces a map of the galaxy's rotational velocity, allowing direct measurement of the galaxy's spin direction, or unit angular momentum vector projected onto the sky. We find no evidence for alignment between galaxy spins and filament directions. We do find hints of a mass-dependent alignment signal, which is in 2-3$\sigma$ tension with the mass-dependent alignment signal in the MassiveBlack-II and Illustris hydrodynamical simulations. However, the tension vanishes when galaxy spin is measured using the H$\alpha$ emission line velocity rather than stellar velocity. Finally, in simulations we find that the mass-dependent transition from aligned to anti-aligned dark matter halo spins is not necessarily present in stellar spins: we find a stellar spin transition in Illustris but not in MassiveBlack-II, highlighting the sensitivity of spin-filament alignments to feedback prescriptions and subgrid physics.

One of the most promising ways to observe the Universe is by detecting the 21cm emission from cosmic neutral hydrogen (HI) through radio-telescopes. Those observations can shed light on fundamental astrophysical questions only if accurate theoretical predictions are available. In order to maximize the scientific return of these surveys, those predictions need to include different observables and be precise on non-linear scales. Currently, one of the best ways to achieve this is via cosmological hydrodynamic simulations; however, the computational cost of these simulations is high -- tens of millions of CPU hours. In this work, we use Wasserstein Generative Adversarial Networks (WGANs) to generate new high-resolution ($35~h^{-1}{\rm kpc}$) 3D realizations of cosmic HI at $z=5$. We do so by sampling from a 100-dimension manifold, learned by the generator, that characterizes the fully non-linear abundance and clustering of cosmic HI from the state-of-the-art simulation IllustrisTNG. We show that different statistical properties of the produced samples -- 1D PDF, power spectrum, bispectrum, and void size function -- match very well those of IllustrisTNG, and outperform state-of-the-art models such as Halo Occupation Distributions (HODs). Our WGAN samples reproduce the abundance of HI across 9 orders of magnitude, from the Ly$\alpha$ forest to Damped Lyman Absorbers. WGAN can produce new samples orders of magnitude faster than hydrodynamic simulations.

We present observations of 50 pairs of redshift z ~ 0.2 star-forming galaxies and background quasars. These sightlines probe the circumgalactic medium (CGM) out to half the virial radius, and we describe the circumgalactic gas kinematics relative to the reference frame defined by the galactic disks. We detect halo gas in MgII absorption, measure the equivalent-width-weighted Doppler shifts relative to each galaxy, and find that the CGM has a component of angular momentum that is aligned with the galactic disk. No net counter-rotation of the CGM is detected within 45 degrees of the major axis at any impact parameter. The velocity offset of the circumgalactic gas correlates with the projected rotation speed in the disk plane out to disk radii of roughly 70 kpc. We confirm previous claims that the MgII absorption becomes stronger near the galactic minor axis and show that the equivalent width correlates with the velocity range of the absorption. We cannot directly measure the location of any absorber along the sightline, but we explore the hypothesis that individual velocity components can be associated with gas orbiting in the disk plane or flowing radially outward in a conical outflow. We conclude that centrifugal forces partially support the low-ionization gas and galactic outflows kinematically disturb the CGM producing excess absorption. Our results firmly rule out schema for the inner CGM that lack rotation and suggest that angular momentum as well as galactic winds should be included in any viable model for the low-redshift CGM.

We show that the cosmological abundance of string axions is much smaller than naive estimates if the Hubble scale of inflation, $H_{\rm inf}$, is sufficiently low (but can still be much higher than the axion masses) and if the inflation lasts sufficiently long. The reason is that the initial misalignment angles of the string axions follow the Bunch-Davies distribution peaked at the potential minima. As a result, the cosmological moduli problem induced by the string axions can be significantly relaxed by low-scale inflation, and astrophysical and cosmological bounds are satisfied over a wide range of the mass without any fine-tuning of the initial misalignment angles. Specifically, the axion with its decay constant $f_\phi = 10^{16}$\,GeV satisfies the bounds over $10^{-18}{\rm \, eV} \lesssim m_\phi \lesssim 10{\rm\,TeV}$ for $H_{\rm inf} \lesssim 10{\rm\,keV}- 10^{6}$\,{\rm GeV}. We also discuss cases with multiple axions and the QCD axion.

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.

Cosmological surveys aim at answering fundamental questions about our Universe, including the nature of dark matter or the reason of unexpected accelerated expansion of the Universe. In order to answer these questions, two important ingredients are needed: 1) data from observations and 2) a theoretical model that allows fast comparison between observation and theory. Most of the cosmological surveys observe galaxies, which are very difficult to model theoretically due to the complicated physics involved in their formation and evolution; modeling realistic galaxies over cosmological volumes requires running computationally expensive hydrodynamic simulations that can cost millions of CPU hours. In this paper, we propose to use deep learning to establish a mapping between the 3D galaxy distribution in hydrodynamic simulations and its underlying dark matter distribution. One of the major challenges in this pursuit is the very high sparsity in the predicted galaxy distribution. To this end, we develop a two-phase convolutional neural network architecture to generate fast galaxy catalogues, and compare our results against a standard cosmological technique. We find that our proposed approach either outperforms or is competitive with traditional cosmological techniques. Compared to the common methods used in cosmology, our approach also provides a nice trade-off between time-consumption (comparable to fastest benchmark in the literature) and the quality and accuracy of the predicted simulation. In combination with current and upcoming data from cosmological observations, our method has the potential to answer fundamental questions about our Universe with the highest accuracy.

Galaxy redshift surveys, such as 2dF, SDSS, 6df, GAMA and VIPERS, have shown that the spatial distribution of matter forms a rich web, known as the cosmic web. The majority of galaxy survey analyses measure the amplitude of galaxy clustering as a function of scale, ignoring information beyond a small number of summary statistics. Since the matter density field becomes highly non-Gaussian as structure evolves under gravity, we expect other statistical descriptions of the field to provide us with additional information. One way to study the non-Gaussianity is to study filaments, which evolve non-linearly from the initial density fluctuations produced in the primordial Universe. In our study, we report the first detection of CMB (Cosmic Microwave Background) lensing by filaments and we apply a null test to confirm our detection. Furthermore, we propose a phenomenological model to interpret the detected signal and we measure how filaments trace the matter distribution on large scales through filament bias, which we measure to be around 1.5. Our study provides a new scope to understand the environmental dependence of galaxy formation. In the future, the joint analysis of lensing and Sunyaev-Zel'dovich observations might reveal the properties of `missing baryons', the vast majority of the gas which resides in the intergalactic medium and has so far evaded most observations.

Numerous observations indicate that galaxies need a continuous gas supply to fuel star formation and explain the star formation history. However, direct observational evidence of gas accretion remains rare. Using the EAGLE cosmological hydrodynamic simulation suite, we study cold gas accretion onto galaxies and the observational signatures of the cold gas kinematics. For EAGLE galaxies at z=0.27, we find that cold gas accretes onto galaxies anisotropically with typical inflow speeds between 20 km s$^{-1}$ and 60 km s$^{-1}$. Most of these galaxies have comparable mass inflow rates and star formation rates, implying that the cold inflowing gas plausibly accounts for sustaining the star-forming activities of the galaxies. As motivation for future work to compare the cold gas kinematics with measurements from quasar sightline observations, we select an EAGLE galaxy with an extended cold gas disk, and we probe the cold gas using mock quasar sightlines. We demonstrate that by viewing the disk edge-on, sightlines at azimuthal angles below 10$^{\circ}$ and impact parameters out to 60 pkpc can detect cold gas that corotates with the galaxy disk. This example suggests cold gas disks that extend beyond the optical disks possibly explain the sightline observations that detect corotating cold gas near galaxy major axes.

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.

The outer regions of globular clusters can enable us to answer many fundamental questions concerning issues ranging from the formation and evolution of clusters and their multiple stellar populations to the study of stars near and beyond the hydrogen-burning limit and to the dynamics of the Milky Way. The outskirts of globular clusters are still uncharted territories observationally. A very efficient way to explore them is through high-precision proper motions of low-mass stars over a large field of view. The Wide Field InfraRed Survey Telescope (WFIRST) combines all these characteristics in a single telescope, making it the best observational tool to uncover the wealth of information contained in the clusters' outermost regions.

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.

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

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

We study the cross-correlation between the Planck CMB lensing convergence map and the eBOSS quasar overdensity obtained from the Sloan Digital Sky Survey (SDSS) IV, in the redshift range $0.9 < z < 2.2$. We detect the CMB lensing convergence-quasar cross power spectrum at $5.4 \sigma$ significance. The cross power spectrum provides a quasar clustering bias measurement that is expected to be particularly robust against systematic effects. The redshift distribution of the quasar sample has a median redshift $z \approx 1.55$, and an effective redshift about $1.51$. The best fit bias of the quasar sample is $b_q = 2.43 \pm 0.45$, corresponding to a host halo mass of $\log_{10}\left( \frac{M}{h^{-1} M_\odot} \right) = 12.54^{+0.25}_{-0.36}$. This is broadly consistent with the previous literature on quasars with a similar redshift range and selection. Since our constraint on the bias comes from the cross-correlation between quasars and CMB lensing, we expect it to be robust to a wide range of possible systematic effects that may contaminate the auto correlation of quasars. We checked for a number of systematic effects from both CMB lensing and quasar overdensity, and found that all systematics are consistent with null within $2 \sigma$. The data is not sensitive to a possible scale dependence of the bias at present, but we expect that as the number of quasars increases (in future surveys such as DESI), it is likely that strong constraints on the scale dependence of the bias can be obtained.

Previous studies have shown the filamentary structures in the cosmic web influence the alignments of nearby galaxies. We study this effect in the LOWZ sample of the Sloan Digital Sky Survey using the "Cosmic Web Reconstruction" filament catalogue. We find that LOWZ galaxies exhibit a small but statistically significant alignment in the direction parallel to the orientation of nearby filaments. This effect is detectable even in the absence of nearby galaxy clusters, which suggests it is an effect from the matter distribution in the filament. A nonparametric regression model suggests that the alignment effect with filaments extends over separations of 30-40 Mpc. We find that galaxies that are bright and early-forming align more strongly with the directions of nearby filaments than those that are faint and late-forming; however, trends with stellar mass are less statistically significant, within the narrow range of stellar mass of this sample.

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.

The gravitational potential well of an $M>10^{13}$ $M_\odot$ galaxy will lead to a gravitational redshift differential of order 1 km/s between the galaxy core and its outskirts. Current surveys of massive galaxies with spatially resolved spectroscopy have reached a size which makes feasible attempts to detect gravitational redshifts within galaxies. We use spectra from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) experiment to attempt a measurement of the averaged stellar redshift profile of large elliptical galaxies. We find that systematic effects (possibly related to charge transfer or wavelength calibration accuracy) make the standard MaNGA data pipeline unsuitable for measuring the relevant sub km/s wavelength shifts. We therefore develop a cross-correlation technique to mitigate these effects, but caution that we are working beyond the design accuracy of the MaNGA experiment. With a sample of $272$ galaxies in halos with $\log (M/M_\odot)>13$, we attempt a measurement of the gravitational redshift profile, achieving $1 \sigma$ errors of size $\sim 0.5$ km/s, but are unable to make a significant detection of the signal. Even without a detection, our measurement can be used to limit the dark matter mass within the half light radius of elliptical galaxies to $1.2$ times the stellar mass, at the 68% confidence level. We also perform weighting scheme tests and split sample tests, and address target selection issues and other relativistic effects, including the transverse Doppler effect and relativistic beaming of stars. Future detections could lead to new constraints on the galaxy mass distribution that are different from kinematic and lensing determinations and open a window on galaxy properties and tests of gravity.

Gravitational redshifts and other relativistic effects are beginning to be studied in the context of galaxy clustering. Distortions consistent with those expected in General Relativity have been measured in galaxy cluster redshift profiles by Wojtak et al. and others and in the the cross-correlation function of galaxy populations by Alam et al. On scales below ~20 Mpc/$h$ simulations have shown that gravitational redshift dominates over other effects. However, this signal is related to the shape and depth of gravitational potentials, and therefore the matter density in galaxies and galaxy clusters that is responsible for them. We investigate the effects of baryonic physics on the gravitational redshift profiles of massive (group and cluster-sized) halos. We compare the profiles of different components in halos taken from the MassiveBlack-II cosmological hydrodynamic simulation and a dark matter-only version of the same simulation. We find that inclusion of baryons, cooling, star formation and feedback significantly alters the relevant inner density profiles. These baryonic effects lead to overall increases in both gravitational redshifts and the transverse relativistic Doppler effects by up to ~50% for group size halos. We show how modified Navarro Frenk White halo profiles can be used to parametrize these differences, and provide relevant halo profile fits.

Marked correlation functions, which are sensitive to the clustering of galaxies in different environments, have been proposed as constraints on modified gravity models. We present measurements of the marked correlation functions of galaxies in redshift space using 361,761 LOWZ ($z_{\rm eff} = 0.32$) galaxies from the Sloan Digital Sky Survey III (SDSS III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) and compare them to $\Lambda$CDM+General Relativity simulations. We apply mass cuts to find the best match between the redshift space autocorrelation function of subhaloes in the simulation and in the observations. We then compare the marked correlation functions, finding no significant evidence for deviations of the marked correlation functions of LOWZ galaxies from $\Lambda$CDM on scales $6 \ h^{-1}$Mpc $ \leq s \leq$ $69 \ h^{-1}$Mpc. The constraining power of marked correlation functions in our analysis is limited by our ability to model the autocorrelation function of galaxies on small scales including the effect of redshift distortions. The statistical errors are well below the differences seen between marked correlation functions of $f$(R) gravity models and $\Lambda$CDM in recent publications (Armijo et al., Hern\'{a}ndez-Aguayo et al.) indicating that improved future theoretical analyses should be able to rule out some models definitively.

One of the most powerful cosmological datasets when it comes to constraining neutrino masses is represented by galaxy power spectrum measurements, $P_{gg}(k)$. The constraining power of $P_{gg}(k)$ is however severely limited by uncertainties in the modeling of the scale-dependent galaxy bias $b(k)$. In this Letter we present a new method to constrain $b(k)$ by using the cross-correlation between the Cosmic Microwave Background (CMB) lensing signal and galaxy maps ($C_\ell^{\rm \kappa g}$) using a simple but theoretically well-motivated parametrization for $b(k)$. We apply the method using $C_\ell^{\rm \kappa g}$ measured by cross-correlating Planck lensing maps and the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 (DR11) CMASS galaxy sample, and $P_{gg}(k)$ measured from the BOSS DR12 CMASS sample. We detect a non-zero scale-dependence at moderate significance, which suggests that a proper modeling of $b(k)$ is necessary in order to reduce the impact of non-linearities and minimize the corresponding systematics. The accomplished increase in constraining power of $P_{gg}(k)$ is demonstrated by determining a 95% C.L. upper bound on the sum of the three active neutrino masses $M_{\nu}$ of $M_{\nu}<0.19\, {\rm eV}$. This limit represents a significant improvement over previous bounds with comparable datasets. Our method will prove especially powerful and important as future large-scale structure surveys will overlap more significantly with the CMB lensing kernel providing a large cross-correlation signal.

We present measurements of $E_G$, a probe of gravity from large-scale structure, using BOSS LOWZ and CMASS spectroscopic samples, with lensing measurements from SDSS (galaxy lensing) and Planck (CMB lensing). Using SDSS lensing and the BOSS LOWZ sample, we measure $\langle{E_G}\rangle=0.40^{+0.05}_{-0.04}$ (stat), $\pm 0.026$ (systematic), consistent with the predicted value from the Planck $\Lambda$CDM model, $E_G=0.46$. Using CMB lensing, we measure $\langle{E_G}\rangle=0.46^{+0.08}_{-0.09}$ (stat) for LOWZ (statistically consistent with galaxy lensing and Planck predictions) and $\langle{E_G}\rangle=0.39^{+0.05}_{-0.05}$ (stat) for the CMASS sample, consistent with the Planck prediction of $E_G=0.40$ given the higher redshift of the sample. We also study the redshift evolution of $E_G$ by splitting the LOWZ sample into two samples based on redshift, with results being consistent with model predictions. We estimate systematic uncertainties on the above $\langle{E_G}\rangle$ numbers to be $\sim 6$% (when using galaxy-galaxy lensing) or $\sim 3$% (when using CMB lensing), subdominant to the quoted statistical errors. These systematic error budgets are dominated by observational systematics in galaxy-galaxy lensing and by theoretical modeling uncertainties, respectively. We do not estimate observational systematics in galaxy-CMB lensing cross correlations.

We report on measurements of the polarization of extragalactic sources at 148 GHz made during the first two seasons of the Atacama Cosmology Telescope Polarization (ACTPol) survey. The survey covered 680 deg$^{2}$ of the sky on the celestial equator. Polarization measurements of 169 intensity-selected sources brighter than 30 mJy, that are predominantly Active Galactic Nuclei, are presented. Above a total flux of 215 mJy where the noise bias removal in the polarization measurement is reliable, we detect 26 sources, 14 of which have a detection of linear polarization at greater than 3$\sigma_{p}$ significance. The distribution of the fractional polarization as a function of total source intensity is analyzed. Our result is consistent with the scenario that the fractional polarization of our measured radio source population is independent of total intensity down to the limits of our measurements and well described by a Gaussian distribution with a mean fractional polarization $p=0.028\pm$0.005 and standard deviation $\sigma_{\mathrm{p}}=0.054$, truncated at $p=0$. Extrapolating this model for the distribution of source polarization below the ACTPol detection threshold, we predict that one could get a clean measure of the E-mode polarization power spectrum of the microwave background out to $\ell=6000$ with 1 $\mu$K-arcminute maps over 10$\%$ of the sky from a future survey. We also study the spectral energy distribution of the total and polarized source flux densities by cross-matching with low radio frequency catalogs. We do not find any correlation between the spectral indices for total flux and polarized flux.

Deep learning is a promising tool to determine the physical model that describes our universe. To handle the considerable computational cost of this problem, we present CosmoFlow: a highly scalable deep learning application built on top of the TensorFlow framework. CosmoFlow uses efficient implementations of 3D convolution and pooling primitives, together with improvements in threading for many element-wise operations, to improve training performance on Intel(C) Xeon Phi(TM) processors. We also utilize the Cray PE Machine Learning Plugin for efficient scaling to multiple nodes. We demonstrate fully synchronous data-parallel training on 8192 nodes of Cori with 77% parallel efficiency, achieving 3.5 Pflop/s sustained performance. To our knowledge, this is the first large-scale science application of the TensorFlow framework at supercomputer scale with fully-synchronous training. These enhancements enable us to process large 3D dark matter distribution and predict the cosmological parameters $\Omega_M$, $\sigma_8$ and n$_s$ with unprecedented accuracy.

We revisit the adiabatic conversion between the QCD axion and axion-like particle (ALP) at level crossing, which can occur in the early universe as a result of the existence of a hypothetical mass mixing. This is similar to the Mikheyev-Smirnov-Wolfenstein effect in neutrino oscillations. After refining the conditions for the adiabatic conversion to occur, we focus on a scenario where the ALP produced by the adiabatic conversion of the QCD axion explains the observed dark matter abundance. Interestingly, we find that the ALP decay constant can be much smaller than the ordinary case in which the ALP is produced by the realignment mechanism. As a consequence, the ALP-photon coupling is enhanced by a few orders of magnitude, which is advantageous for the future ALP and axion-search experiments using the ALP-photon coupling.

The Simons Observatory (SO) will measure the cosmic microwave background (CMB) in both temperature and polarization over a wide range of angular scales and frequencies from 27-270 GHz with unprecedented sensitivity. One technology for coupling light onto the $\sim$50 detector wafers that SO will field is spline-profiled feedhorns, which offer tunability between coupling efficiency and control of beam polarization leakage effects. We will present efforts to scale up feedhorn production for SO and their viability for future CMB experiments, including direct-machining metal feedhorn arrays and laser machining stacked Si arrays.

In this proceeding, we present studies of instrumental systematic effects for the Simons Obsevatory (SO) that are associated with the detector system and its interaction with the full SO experimental systems. SO will measure the Cosmic Microwave Background (CMB) temperature and polarization anisotropies over a wide range of angular scales in six bands with bandcenters spanning from 27 GHz to 270 GHz. We explore effects including intensity-to-polarization leakage due to coupling optics, bolometer nonlinearity, uncalibrated gain variations of bolometers, and readout crosstalk. We model the level of signal contamination, discuss proposed mitigation schemes, and present instrument requirements to inform the design of SO and future CMB projects.

The Simons Observatory (SO) will make precision temperature and polarization measurements of the cosmic microwave background (CMB) using a series of telescopes which will cover angular scales between one arcminute and tens of degrees and sample frequencies between 27 and 270 GHz. Here we present the current design of the large aperture telescope receiver (LATR), a 2.4 m diameter cryostat that will be mounted on the SO 6 m telescope and will be the largest CMB receiver to date. The cryostat size was chosen to take advantage of the large focal plane area having high Strehl ratios, which is inherent to the Cross-Dragone telescope design. The LATR will be able to accommodate thirteen optics tubes, each having a 36 cm diameter aperture and illuminating several thousand transition-edge sensor (TES) bolometers. This set of equipment will provide an opportunity to make measurements with unparalleled sensitivity. However, the size and complexity of the LATR also pose numerous technical challenges. In the following paper, we present the design of the LATR and include how we address these challenges. The solutions we develop in the process of designing the LATR will be informative for the general CMB community, and for future CMB experiments like CMB-S4.

We present a search for CII emission over cosmological scales at high-redshifts. The CII line is a prime candidate to be a tracer of star formation over large-scale structure since it is one of the brightest emission lines from galaxies. Redshifted CII emission appears in the submillimeter regime, meaning it could potentially be present in the higher frequency intensity data from the Planck satellite used to measure the cosmic infrared background (CIB). We search for CII emission over redshifts z=2-3.2 in the Planck 545 GHz intensity map by cross-correlating the 3 highest frequency Planck maps with spectroscopic quasars and CMASS galaxies from the Sloan Digital Sky Survey III (SDSS-III), which we then use to jointly fit for CII intensity, CIB parameters, and thermal Sunyaev-Zeldovich (SZ) emission. We report a measurement of an anomalous emission $\mathrm{I_\nu}=6.6^{+5.0}_{-4.8}\times10^4$ Jy/sr at 95% confidence, which could be explained by CII emission, favoring collisional excitation models of CII emission that tend to be more optimistic than models based on CII luminosity scaling relations from local measurements; however, a comparison of Bayesian information criteria reveal that this model and the CIB & SZ only model are equally plausible. Thus, more sensitive measurements will be needed to confirm the existence of large-scale CII emission at high redshifts. Finally, we forecast that intensity maps from Planck cross-correlated with quasars from the Dark Energy Spectroscopic Instrument (DESI) would increase our sensitivity to CII emission by a factor of 5, while the proposed Primordial Inflation Explorer (PIXIE) could increase the sensitivity further.

The Simons Observatory (SO) will make precise temperature and polarization measurements of the cosmic microwave background (CMB) using a set of telescopes which will cover angular scales between 1 arcminute and tens of degrees, contain over 60,000 detectors, and observe at frequencies between 27 and 270 GHz. SO will consist of a 6 m aperture telescope coupled to over 30,000 transition-edge sensor bolometers along with three 42 cm aperture refractive telescopes, coupled to an additional 30,000+ detectors, all of which will be located in the Atacama Desert at an altitude of 5190 m. The powerful combination of large and small apertures in a CMB observatory will allow us to sample a wide range of angular scales over a common survey area. SO will measure fundamental cosmological parameters of our universe, constrain primordial fluctuations, find high redshift clusters via the Sunyaev-Zel`dovich effect, constrain properties of neutrinos, and trace the density and velocity of the matter in the universe over cosmic time. The complex set of technical and science requirements for this experiment has led to innovative instrumentation solutions which we will discuss. The large aperture telescope will couple to a cryogenic receiver that is 2.4 m in diameter and nearly 3 m long, creating a number of technical challenges. Concurrently, we are designing the array of cryogenic receivers housing the 42 cm aperture telescopes. We will discuss the sensor technology SO will use and we will give an overview of the drivers for and designs of the SO telescopes and receivers, with their cold optical components and detector arrays.

The Advanced ACTPol (AdvACT) upgrade to the Atacama Cosmology Telescope (ACT) features arrays of aluminum manganese transition-edge sensors (TESes) optimized for ground-based observations of the Cosmic Microwave Background (CMB). Array testing shows highly responsive detectors with anticipated in-band noise performance under optical loading. We report on TES parameters measured with impedance data taken on a subset of TESes. We then compare modeled noise spectral densities to measurements. We find excess noise at frequencies around 100 Hz, nearly outside of the signal band of CMB measurements. In addition, we describe full-array noise measurements in the laboratory and in the field for two new AdvACT mid-frequency arrays, sensitive at bands centered on 90 and 150 GHz, and data for the high-frequency array (150/230 GHz) as deployed.

The Advanced ACTPol upgrade on the Atacama Cosmology Telescope aims to improve the measurement of the cosmic microwave background anisotropies and polarization, using four new dichroic detector arrays fabricated on 150-mm silicon wafers. These bolometric cameras use AlMn transition-edge sensors, coupled to feedhorns with orthomode transducers for polarization sensitivity. The first deployed camera is sensitive to both 150 GHz and 230 GHz. Here we present the lab characterization of the thermal parameters and optical efficiencies for the two newest fielded arrays, each sensitive to both 90 GHz and 150 GHz. We provide assessments of the parameter uniformity across each array with evaluation of systematic uncertainties. Lastly, we show the arrays' initial performance in the field.

We develop an automated technique for detecting damped Lyman-$\alpha$ absorbers (DLAs) along spectroscopic lines of sight to quasi-stellar objects (QSOs or quasars). The detection of DLAs in large-scale spectroscopic surveys such as SDSS-III sheds light on galaxy formation at high redshift, showing the nucleation of galaxies from diffuse gas. We use nearly 50 000 QSO spectra to learn a novel tailored Gaussian process model for quasar emission spectra, which we apply to the DLA detection problem via Bayesian model selection. We propose models for identifying an arbitrary number of DLAs along a given line of sight. We demonstrate our method's effectiveness using a large-scale validation experiment, with excellent performance. We also provide a catalog of our results applied to 162 858 spectra from SDSS-III data release 12.

Advanced ACTPol (AdvACT) is a third generation polarization upgrade to the Atacama Cosmology Telescope, designed to observe the cosmic microwave background (CMB). AdvACT expands on the 90 and 150 GHz transition edge sensor (TES) bolometer arrays of the ACT Polarimeter (ACTPol), adding both high frequency (HF, 150/230 GHz) and low frequency (LF, 27/39 GHz) multichroic arrays. The addition of the high and low frequency detectors allows for the characterization of synchrotron and spinning dust emission at the low frequencies and foreground emission from galactic dust and dusty star forming galaxies at the high frequencies. The increased spectral coverage of AdvACT will enable a wide range of CMB science, such as improving constraints on dark energy, the sum of the neutrino masses, and the existence of primordial gravitational waves. The LF array will be the final AdvACT array, replacing one of the MF arrays for a single season. Prior to the fabrication of the final LF detector array, we designed and characterized prototype TES bolometers. Detector geometries in these prototypes are varied in order to inform and optimize the bolometer designs for the LF array, which requires significantly lower noise levels and saturation powers (as low as ${\sim}1$ pW) than the higher frequency detectors. Here we present results from tests of the first LF prototype TES detectors for AdvACT, including measurements of the saturation power, critical temperature, thermal conductance and time constants. We also describe the modifications to the time-division SQUID readout architecture compared to the MF and HF arrays.

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.

Standard present day large-scale structure (LSS) analyses make a major assumption in their Bayesian parameter inference --- that the likelihood has a Gaussian form. For summary statistics currently used in LSS, this assumption, even if the underlying density field is Gaussian, cannot be correct in detail. We investigate the impact of this assumption on two recent LSS analyses: the Beutler et al. (2017) power spectrum multipole ($P_\ell$) analysis and the Sinha et al. (2017) group multiplicity function ($\zeta$) analysis. Using non-parametric divergence estimators on mock catalogs originally constructed for covariance matrix estimation, we identify significant non-Gaussianity in both the $P_\ell$ and $\zeta$ likelihoods. We then use Gaussian mixture density estimation and Independent Component Analysis on the same mocks to construct likelihood estimates that approximate the true likelihood better than the Gaussian $pseudo$-likelihood. Using these likelihood estimates, we accurately estimate the true posterior probability distribution of the Beutler et al. (2017) and Sinha et al. (2017) parameters. Likelihood non-Gaussianity shifts the $f\sigma_8$ constraint by $-0.44\sigma$, but otherwise, does not significantly impact the overall parameter constraints of Beutler et al. (2017). For the $\zeta$ analysis, using the pseudo-likelihood significantly underestimates the uncertainties and biases the constraints of Sinha et al. (2017) halo occupation parameters. For $\log M_1$ and $\alpha$, the posteriors are shifted by $+0.43\sigma$ and $-0.51\sigma$ and broadened by $42\%$ and $66\%$, respectively. The divergence and likelihood estimation methods we present provide a straightforward framework for quantifying the impact of likelihood non-Gaussianity and deriving more accurate parameter constraints.

We investigate the potential sources of theoretical systematics in the anisotropic Baryon Acoustic Oscillation (BAO) distance scale measurements from the clustering of galaxies in configuration space using the final Data Release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS). We perform a detailed study of the impact on BAO measurements from choices in the methodology such as fiducial cosmology, clustering estimators, random catalogues, fitting templates, and covariance matrices. The theoretical systematic uncertainties in BAO parameters are found to be 0.002 in the isotropic dilation $\alpha$ and 0.003 in the quadrupolar dilation $\epsilon$. The leading source of systematic uncertainty is related to the reconstruction techniques. Theoretical uncertainties are sub-dominant compared with the statistical uncertainties for BOSS survey, accounting $0.2\sigma_{stat}$ for $\alpha$ and $0.25\sigma_{stat}$ for $\epsilon$ ($\sigma_{\alpha,stat} \sim$0.010 and $\sigma_{\epsilon,stat}\sim$ 0.012 respectively). We also present BAO-only distance scale constraints from the anisotropic analysis of the correlation function. Our constraints on the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$, including both statistical and theoretical systematic uncertainties, are 1.5\% and 2.8\% at $z_{\rm eff}=0.38$, 1.4\% and 2.4\% at $z_{\rm eff}=0.51$, and 1.7\% and 2.6\% at $z_{\rm eff}=0.61$. This paper is part of a set that analyzes the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are cross-checked with other BAO analysis in \citet{Acacia16}. The systematic error budget concerning the methodology on post-reconstruction BAO analysis presented here is used in \citet{Acacia16} to produce the final cosmological constraints from BOSS.

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 Atacama B-mode Search (ABS) is an experiment designed to measure cosmic microwave background (CMB) polarization at large angular scales ($\ell>40$). It operated from the ACT site at 5190~m elevation in northern Chile at 145 GHz with a net sensitivity (NEQ) of 41 $\mu$K$\sqrt{\rm s}$. It employed an ambient-temperature sapphire half-wave plate rotating at 2.55 Hz to modulate the incident polarization signal and reduce systematic effects. We report here on the analysis of data from a 2400 deg$^2$ patch of sky centered at declination $-42^\circ$ and right ascension $25^\circ$. We perform a blind analysis. After unblinding, we find agreement with the Planck TE and EE measurements on the same region of sky. We marginally detect polarized dust emission and give an upper limit on the tensor-to-scalar ratio of $r<2.3$ (95% cl) with the equivalent of 100 on-sky days of observation. We also present a new measurement of the polarization of Tau A and introduce new methods associated with HWP-based observations.

Using some of the latest cosmological datasets publicly available, we derive the strongest bounds in the literature on the sum of the three active neutrino masses, $M_\nu$, within the assumption of a background flat $\Lambda$CDM cosmology. In the most conservative scheme, combining Planck cosmic microwave background (CMB) temperature anisotropies and baryon acoustic oscillations (BAO) data, as well as the up-to-date constraint on the optical depth to reionization ($\tau$), the tightest $95\%$ confidence level (C.L.) upper bound we find is $M_\nu<0.151$~eV. The addition of Planck high-$\ell$ polarization data, which however might still be contaminated by systematics, further tightens the bound to $M_\nu<0.118$~eV. A proper model comparison treatment shows that the two aforementioned combinations disfavor the IH at $\sim 64\%$~C.L. and $\sim 71\%$~C.L. respectively. In addition, we compare the constraining power of measurements of the full-shape galaxy power spectrum versus the BAO signature, from the BOSS survey. Even though the latest BOSS full shape measurements cover a larger volume and benefit from smaller error bars compared to previous similar measurements, the analysis method commonly adopted results in their constraining power still being less powerful than that of the extracted BAO signal. Our work uses only cosmological data; imposing the constraint $M_\nu>0.06\,{\rm eV}$ from oscillations data would raise the quoted upper bounds by ${\cal O}(0.1\sigma)$ and would not affect our conclusions.

A grand challenge of the 21st century cosmology is to accurately estimate the cosmological parameters of our Universe. A major approach to estimating the cosmological parameters is to use the large-scale matter distribution of the Universe. Galaxy surveys provide the means to map out cosmic large-scale structure in three dimensions. Information about galaxy locations is typically summarized in a "single" function of scale, such as the galaxy correlation function or power-spectrum. We show that it is possible to estimate these cosmological parameters directly from the distribution of matter. This paper presents the application of deep 3D convolutional networks to volumetric representation of dark-matter simulations as well as the results obtained using a recently proposed distribution regression framework, showing that machine learning techniques are comparable to, and can sometimes outperform, maximum-likelihood point estimates using "cosmological models". This opens the way to estimating the parameters of our Universe with higher accuracy.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

Large redshift surveys of galaxies and clusters are providing the first opportunities to search for distortions in the observed pattern of large-scale structure due to such effects as gravitational redshift. We focus on non-linear scales and apply a quasi-Newtonian approach using N-body simulations to predict the small asymmetries in the cross-correlation function of two galaxy different populations. Following recent work by Bonvin et al., Zhao and Peacock and Kaiser on galaxy clusters, we include effects which enter at the same order as gravitational redshift: the transverse Doppler effect, light-cone effects, relativistic beaming, luminosity distance perturbation and wide-angle effects. We find that all these effects cause asymmetries in the cross-correlation functions. Quantifying these asymmetries, we find that the total effect is dominated by the gravitational redshift and luminosity distance perturbation at small and large scales, respectively. By adding additional subresolution modelling of galaxy structure to the large-scale structure information, we find that the signal is significantly increased, indicating that structure on the smallest scales is important and should be included. We report on comparison of our simulation results with measurements from the SDSS/BOSS galaxy redshift survey in a companion paper.

General relativistic effects have long been predicted to subtly influence the observed large-scale structure of the universe. The current generation of galaxy redshift surveys have reached a size where detection of such effects is becoming feasible. In this paper, we report the first detection of the redshift asymmetry from the cross-correlation function of two galaxy populations which is consistent with relativistic effects. The dataset is taken from the Sloan Digital Sky Survey DR12 CMASS galaxy sample, and we detect the asymmetry at the $2.7\sigma$ level by applying a shell-averaged estimator to the cross-correlation function. Our measurement dominates at scales around $10$ h$^{-1}$Mpc, larger than those over which the gravitational redshift profile has been recently measured in galaxy clusters, but smaller than scales for which linear perturbation theory is likely to be accurate. The detection significance varies by 0.5$\sigma$ with the details of our measurement and tests for systematic effects. We have also devised two null tests to check for various survey systematics and show that both results are consistent with the null hypothesis. We measure the dipole moment of the cross-correlation function, and from this the asymmetry is also detected, at the $2.8 \sigma$ level. The amplitude and scale-dependence of the clustering asymmetries are approximately consistent with the expectations of General Relativity and a biased galaxy population, within large uncertainties. We explore theoretical predictions using numerical simulations in a companion paper.

We study the asymmetry in the two-point cross-correlation function of two populations of galaxies focusing in particular on the relativistic effects that include the gravitational redshift. We derive the cross-correlation function on small and large scales using two different approaches: General Relativistic and Newtonian perturbation theory. Following recent work by Bonvin et al., Gaztanaga et al. and Croft, we calculate the dipole and the shell estimator with the two procedures and we compare our results. We find that while General Relativistic Perturbation Theory (GRPT) is able to make predictions of relativistic effects on very large, obviously linear scales (r > 50 Mpc/h), the presence of non-linearities physically occurring on much smaller scales (down to those describing galactic potential wells) can strongly affect the asymmetry estimators. These can lead to cancellations of the relativistic terms, and sign changes in the estimators on scales up to r ~ 50 Mpc/h. On the other hand, with an appropriate non-linear gravitational potential, the results obtained using Newtonian theory can successfully describe the asymmetry on smaller, non-linear scales (r < 20 Mpc/h) where gravitational redshift is the dominant term. On larger scales the asymmetry is much smaller in magnitude, and measurement is not within reach of current observations. This is in agreement with the observational results obtained by Gaztnaga et al. and the first detection of relativistic effects (on (r < 20 Mpc/h) scales) by Alam et al.

In a galaxy redshift survey the objects to be targeted for spectra are selected from a photometrically observed sample. The observed magnitudes and colours of galaxies in this parent sample will be affected by their peculiar velocities, through relativistic Doppler and relativistic beaming effects. In this paper we compute the resulting expected changes in galaxy photometry. The magnitudes of the relativistic effects are a function of redshift, stellar mass, galaxy velocity and velocity direction. We focus on the CMASS sample from the Sloan Digital Sky Survey (SDSS), Baryon Oscillation Spectroscopic Survey (BOSS), which is selected on the basis of colour and magnitude. We find that 0.10\% of the sample ($\sim 585$ galaxies) has been scattered into the targeted region of colour-magnitude space by relativistic effects, and conversely 0.09\% of the sample ($\sim 532$ galaxies) has been scattered out. Observational consequences of these effects include an asymmetry in clustering statistics, which we explore in a companion paper. Here we compute a set of weights which can be used to remove the effect of modulations introduced into the density field inferred from a galaxy sample. We conclude by investigating the possible effects of these relativistic modulation on large scale clustering of the galaxy sample.

CMB-S4 is a proposed experiment to map the polarization of the Cosmic Microwave Background (CMB) to nearly the cosmic variance limit for angular scales that are accessible from the ground. The science goals and capabilities of CMB-S4 in illuminating cosmic inflation, measuring the sum of neutrino masses, searching for relativistic relics in the early universe, characterizing dark energy and dark matter, and mapping the matter distribution in the universe have been described in the CMB-S4 Science Book. This Technology Book is a companion volume to the Science Book. The ambitious science goals of CMB-S4, a "Stage-4" experiment, require a step forward in experimental capability from the current Stage=II experiments. To guide this process, we summarize the current state of CMB instrumentation technology, and identify R&D efforts necessary to advance it for use in CMB-S4. The book focuses on technical challenges in four broad areas: Telescope Design; Receiver Optics; Focal-Plane Optical Coupling; and Focal-Plane Sensor and Readout.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We present a measurement of the linear growth rate of structure, \textit{f} from the Sloan Digital Sky Survey III (SDSS III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) using Convolution Lagrangian Perturbation Theory (CLPT) with Gaussian Streaming Redshift-Space Distortions (GSRSD) to model the two point statistics of BOSS galaxies in DR12. The BOSS-DR12 dataset includes 1,198,006 massive galaxies spread over the redshift range $0.2 < z < 0.75$. These galaxy samples are categorized in three redshift bins. Using CLPT-GSRSD in our analysis of the combined sample of the three redshift bins, we report measurements of $f \sigma_8$ for the three redshift bins. We find $f \sigma_8 = 0.430 \pm 0.054$ at $z_{\rm eff} = 0.38$, $f \sigma_8 = 0.452 \pm 0.057$ at $z_{\rm eff} = 0.51$ and $f \sigma_8 = 0.457 \pm 0.052$ at $z_{\rm eff} = 0.61$. Our results are consistent with the predictions of Planck $\Lambda$CDM-GR. Our constraints on the growth rates of structure in the Universe at different redshifts serve as a useful probe, which can help distinguish between a model of the Universe based on dark energy and models based on modified theories of gravity. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. 2016 to produce the final cosmological constraints from BOSS.

We perform a tomographic baryon acoustic oscillations analysis using the two-point galaxy correlation function measured from the combined sample of BOSS DR12, which covers the redshift range of $0.2<z<0.75$. Splitting the sample into multiple overlapping redshift slices to extract the redshift information of galaxy clustering, we obtain a measurement of $D_A(z)/r_d$ and $H(z)r_d$ at nine effective redshifts with the full covariance matrix calibrated using MultiDark-Patchy mock catalogues. Using the reconstructed galaxy catalogues, we obtain the precision of $1.3\%-2.2\%$ for $D_A(z)/r_d$ and $2.1\%-6.0\%$ for $H(z)r_d$. To quantify the gain from the tomographic information, we compare the constraints on the cosmological parameters using our 9-bin BAO measurements, the consensus 3-bin BAO and RSD measurements at three effective redshifts in \citet{Alam2016}, and the non-tomographic (1-bin) BAO measurement at a single effective redshift. Comparing the 9-bin with 1-bin constraint result, it can improve the dark energy Figure of Merit by a factor of 1.24 for the Chevallier-Polarski-Linder parametrisation for equation of state parameter $w_{\rm DE}$. The errors of $w_0$ and $w_a$ from 9-bin constraints are slightly improved when compared to the 3-bin constraint result.

We describe the kinematics of circumgalactic gas near the galactic plane, combining new measurements of galaxy rotation curves and spectroscopy of background quasars. The sightlines pass within 19--93 kpc of the target galaxy and generally detect Mg II absorption. The Mg II Doppler shifts have the same sign as the galactic rotation, so the cold gas co-rotates with the galaxy. Because the absorption spans a broader velocity range than disk rotation can explain, we explore simple models for the circumgalactic kinematics. Gas spiraling inwards (near the disk plane) offers a successful description of the observations. An Appendix describes the addition of tangential and radial gas flows and illustrates how the sign of the disk inclination produces testable differences in the projected line-of-sight velocity range. This inflow interpretation implies that cold flow disks remain common down to redshift $z \approx 0.2$ and prolong star formation by supplying gas to the disk.

We present a Bayesian phase-space reconstruction of the cosmic large-scale matter density and velocity fields from the SDSS-III Baryon Oscillations Spectroscopic Survey Data Release 12 (BOSS DR12) CMASS galaxy clustering catalogue. We rely on a given $\Lambda$CDM cosmology, a mesh resolution in the range of 6-10 $h^{-1}$ Mpc, and a lognormal-Poisson model with a redshift dependent nonlinear bias. The bias parameters are derived from the data and a general renormalised perturbation theory approach. We use combined Gibbs and Hamiltonian sampling, implemented in the \textsc{argo} code, to iteratively reconstruct the dark matter density field and the coherent peculiar velocities of individual galaxies, correcting hereby for coherent redshift space distortions (RSD). Our tests relying on accurate $N$-body based mock galaxy catalogues, show unbiased real space power spectra of the nonlinear density field up to $k\sim0.2\, h$ Mpc$^{-1}$, and vanishing quadrupoles down to $r\sim20\,h^{-1}$ Mpc. We also demonstrate that the nonlinear cosmic web can be obtained from the tidal field tensor based on the Gaussian component of the reconstructed density field. We find that the reconstructed velocities have a statistical correlation coefficient compared to the true velocities of each individual lightcone mock galaxy of $r\sim0.68$ including about 10% of satellite galaxies with virial motions (about $r=0.75$ without satellites). The power spectra of the velocity divergence agree well with theoretical predictions up to $k\sim0.2\,h\,{\rm Mpc}^{-1}$. This work will be especially useful to improve, e.g. BAO reconstructions, kinematic Sunyaev-Zeldovich (kSZ), integrated Sachs-Wolfe (ISW) measurements, or environmental studies.

We study the effects of filaments on galaxy properties in the Sloan Digital Sky Survey (SDSS) Data Release 12 using filaments from the `Cosmic Web Reconstruction' catalogue (Chen et al. 2016), a publicly available filament catalogue for SDSS. Since filaments are tracers of medium-to-high density regions, we expect that galaxy properties associated with the environment are dependent on the distance to the nearest filament. Our analysis demonstrates that a red galaxy or a high-mass galaxy tend to reside closer to filaments than a blue or low-mass galaxy. After adjusting the effect from stellar mass, on average, early-forming galaxies or large galaxies have a shorter distance to filaments than late-forming galaxies or small galaxies. For the Main galaxy sample (MGS), all signals are very significant ($>6\sigma$). For the LOWZ and CMASS sample, the stellar mass and size are significant ($>2 \sigma$). The filament effects we observe persist until $z = 0.7$ (the edge of the CMASS sample). Comparing our results to those using the galaxy distances from redMaPPer galaxy clusters as a reference, we find a similar result between filaments and clusters. Moreover, we find that the effect of clusters on the stellar mass of nearby galaxies depends on the galaxy's filamentary environment. Our findings illustrate the strong correlation of galaxy properties with proximity to density ridges, strongly supporting the claim that density ridges are good tracers of filaments.

We present new estimates for the statistical properties of damped Lyman-$\alpha$ absorbers (DLAs). We compute the column density distribution function at $z>2$, the line density, $\mathrm{d}N/\mathrm{d}X$, and the neutral hydrogen density, $\Omega_\mathrm{DLA}$. Our estimates are derived from the DLA catalogue of Garnett 2016, which uses the SDSS-III DR12 quasar spectroscopic survey. This catalogue provides a probability that a given spectrum contains a DLA, allowing us to use even the noisiest data without biasing our results and thus substantially increase our sample size. We measure a non-zero column density distribution function at $95\%$ confidence for all column densities $N_\mathrm{HI} < 5\times 10^{22}$ cm$^{-2}$. We make the first measurements from SDSS of $\mathrm{d}N/\mathrm{d}X$ and $\Omega_\mathrm{DLA}$ at $z>4$. We show that our results are insensitive to the signal-to-noise ratio of the spectra, but that there is a residual dependence on quasar redshift for $z<2.5$, which may be due to remaining systematics in our analysis.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

We measure a combination of gravitational lensing, galaxy clustering, and redshift-space distortions called $E_G$. The quantity $E_G$ probes both parts of metric potential and is insensitive to galaxy bias and $\sigma_8$. These properties make it an attractive statistic to test $\Lambda$CDM, General Relativity and its alternate theories. We have combined CMASS DR11 with CFHTLenS and recent measurements of $\beta$ from RSD analysis, and find $E_G(z = 0.57) = 0.42 \pm 0.056$, an 13\% measurement in agreement with the prediction of general relativity $E_G(z = 0.57) = 0.396 \pm 0.011$ using the Planck 2015 cosmological parameters. We have corrected our measurement for various observational and theoretical systematics. Our measurement is consistent with the first measurement of $E_G$ using CMB lensing in place of galaxy lensing (Pullen et. al. 2015a) at small scales, but shows 2.8$\sigma$ tension when compared with their final results including large scales. This analysis with future surveys will provide improved statistical error and better control over systematics to test General Relativity and its alternate theories.

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

The reconstruction algorithm introduced by \cite{Eis07}, which is widely used in clustering analysis, is based on the inference of the first order Lagrangian displacement field from the Gaussian smoothed galaxy density field in redshift space. The \modif2{smoothing scale} applied to the density field affects the inferred displacement field that is used to move {the galaxies}, and partially \modif2{erases} the nonlinear evolution {of the density field}. In this article, we explore this crucial step \modif2{in} the reconstruction algorithm. We study the performance of the reconstruction technique using two metrics: first, we study the performance using the anisotropic clustering, extending previous studies focused on isotropic clustering; second, we study its effect on the displacement field. We find that smoothing has a strong effect in the quadrupole of the correlation function and affects the accuracy and precision \modif2{with} which we can measure $D_A (z)$ and $H(z)$. We find that the optimal smoothing scale to use in the reconstruction algorithm applied to BOSS-CMASS is between 5-10 $h^{-1}$Mpc. Varying from the "usual" 15$h^{-1}$Mpc to $5 h^{-1}$Mpc \modif2{shows} $\sim$ 0.3\% variations in $D_A(z)$ and $\sim$ 0.4\% $H(z)$ and uncertainties are also reduced by 40\% and 30\% respectively. We also find that the accuracy of velocity field reconstruction depends strongly on the smoothing scale used for the density field. We measure the bias and uncertainties associated with different choices of smoothing length.

We study the impact of assumptions about neutrino properties on the estimation of inflationary parameters from cosmological data, with a specific focus on the allowed contours in the $n_s/r$ plane. We study the following neutrino properties: (i) the total neutrino mass $ M_\nu =\sum_i m_i$; (ii) the number of relativistic degrees of freedom $N_{eff}$; and (iii) the neutrino hierarchy: whereas previous literature assumed 3 degenerate neutrino masses or two massless neutrino species (that do not match neutrino oscillation data), we study the cases of normal and inverted hierarchy. Our basic result is that these three neutrino properties induce $< 1 \sigma$ shift of the probability contours in the $n_s/r$ plane with both current or upcoming data. We find that the choice of neutrino hierarchy has a negligible impact. However, the minimal cutoff on the total neutrino mass $M_{\nu,{min}}=0 $ that accompanies previous works using the degenerate hierarchy does introduce biases in the $n_s/r$ plane and should be replaced by $M_{\nu,min}= 0.059$ eV as required by oscillation data. Using current CMB data from Planck and Bicep/Keck, marginalizing over $ M_\nu$ and over $r$ can lead to a shift in the mean value of $n_s$ of $\sim0.3\sigma$ towards lower values. However, once BAO measurements are included, the standard contours in the $n_s/r$ plane are basically reproduced. Larger shifts of the contours in the $n_s/r$ plane (up to 0.8$\sigma$) arise for nonstandard values of $N_{eff}$. We also provide forecasts for the future CMB experiments COrE and Stage-IV and show that the incomplete knowledge of neutrino properties, taken into account by a marginalization over $M_\nu$, could induce a shift of $\sim0.4\sigma$ towards lower values in the determination of $n_s$ (or a $\sim 0.8\sigma$ shift if one marginalizes over $N_{eff}$). Comparison to specific inflationary models is shown.

We present baryon acoustic oscillation (BAO) scale measurements determined from the clustering of 1.2 million massive galaxies with redshifts 0.2 < z < 0.75 distributed over 9300 square degrees, as quantified by their redshift-space correlation function. In order to facilitate these measurements, we define, describe, and motivate the selection function for galaxies in the final data release (DR12) of the SDSS III Baryon Oscillation Spectroscopic Survey (BOSS). This includes the observational footprint, masks for image quality and Galactic extinction, and weights to account for density relationships intrinsic to the imaging and spectroscopic portions of the survey. We simulate the observed systematic trends in mock galaxy samples and demonstrate that they impart no bias on baryon acoustic oscillation (BAO) scale measurements and have a minor impact on the recovered statistical uncertainty. We measure transverse and radial BAO distance measurements in 0.2 < z < 0.5, 0.5 < z < 0.75, and (overlapping) 0.4 < z < 0.6 redshift bins. In each redshift bin, we obtain a precision that is 2.7 per cent or better on the radial distance and 1.6 per cent or better on the transverse distance. The combination of the redshift bins represents 1.8 per cent precision on the radial distance and 1.1 per cent precision on the transverse distance. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

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

We present an evaluation of systematic effects associated with a continuously-rotating, ambient-temperature half-wave plate (HWP) based on two seasons of data from the Atacama B-Mode Search (ABS) experiment located in the Atacama Desert of Chile. The ABS experiment is a microwave telescope sensitive at 145 GHz. Here we present our in-field evaluation of celestial (CMB plus galactic foreground) temperature-to-polarization leakage. We decompose the leakage into scalar, dipole, and quadrupole leakage terms. We report a scalar leakage of ~0.01%, consistent with model expectations and an order of magnitude smaller than other CMB experiments have reported. No significant dipole or quadrupole terms are detected; we constrain each to be <0.07% (95% confidence), limited by statistical uncertainty in our measurement. Dipole and quadrupole leakage at this level lead to systematic error on r<0.01 before any mitigation due to scan cross-linking or boresight rotation. The measured scalar leakage and the theoretical level of dipole and quadrupole leakage produce systematic error of r<0.001 for the ABS survey and focal-plane layout before any data correction such as so-called deprojection. This demonstrates that ABS achieves significant beam systematic error mitigation from its HWP and shows the promise of continuously-rotating HWPs for future experiments.

Post-starburst galaxies are in the transitional stage between blue, star-forming galaxies and red, quiescent galaxies, and therefore hold important clues for our understanding of galaxy evolution. In this paper, we systematically searched for and identified a large sample of post-starburst galaxies from the spectroscopic dataset of the Sloan Digital Sky Survey (SDSS) Data Release 9. In total, we found more than 6000 objects with redshifts between $z\sim0.05$ and $z\sim1.3$, making this the largest sample of post-starburst galaxies in the literature. We calculated the luminosity function of the post-starburst galaxies using two uniformly selected subsamples: the SDSS Main Galaxy Sample and the Baryon Oscillation Spectroscopic Survey CMASS Sample. The luminosity functions are reasonably fit by half-Gaussian functions. The peak magnitudes shift as a function of redshift from $M\sim-23.5$ at $z\sim0.8$ to $M\sim-20.3$ at $z\sim0.1$. This is consistent with the downsizing trend, whereby more massive galaxies form earlier than low-mass galaxies. We compared the mass of the post-starburst stellar population found in our sample to the decline of the global star-formation rate and found that only a small amount ($\sim1\%$) of all star-formation quenching in the redshift range $z=0.2-0.7$ results in post-starburst galaxies in the luminosity range our sample is sensitive to. Therefore, luminous post-starburst galaxies are not the place where most of the decline in star-formation rate of the universe is happening.

We demonstrate a new method to constrain gravity on the largest cosmological scales by combining measurements of cosmic microwave background (CMB) lensing and the galaxy velocity field. $E_G$ is a statistic, constructed from a gravitational lensing tracer and a measure of velocities such as redshift-space distortions (RSD), that can discriminate between gravity models while being independent of clustering bias and $\sigma_8$. While traditionally, the lensing field for $E_G$ has been probed through galaxy lensing, CMB lensing has been proposed as a more robust tracer of the lensing field for $E_G$ at higher redshifts while avoiding intrinsic alignments. We perform the largest-scale measurement of $E_G$ ever, up to 150 Mpc/$h$, by cross-correlating the Planck CMB lensing map with the Sloan Digital Sky Survey III (SDSS-III) CMASS galaxy sample and combining this with our measurement of the CMASS auto-power spectrum and the RSD parameter $\beta$. We report $E_G(z=0.57)=0.243\pm0.060$ (stat) $\pm0.013$ (sys), a measurement in tension with the general relativity prediction at a level of 2.6$\sigma$. Note that our $E_G$ measurement deviates from GR only at scales greater than 80 Mpc/$h$, scales which have not been probed by previous $E_G$ tests. Upcoming surveys, which will provide an order-of-magnitude reduction in statistical errors, can significantly constrain alternative gravity models when combined with better control of systematics.

The Atacama Cosmology Telescope Polarimeter (ACTPol) is a polarization sensitive upgrade to the Atacama Cosmology Telescope. Located at an elevation of 5190 m, ACTPol measures the Cosmic Microwave Background (CMB) temperature and polarization with arcminute-scale angular resolution. Calibration of the detector angles is a critical step in producing maps of the CMB polarization. Polarization angle offsets in the detector calibration can cause leakage in polarization from E to B modes and induce a spurious signal in the EB and TB cross correlations, which eliminates our ability to measure potential cosmological sources of EB and TB signals, such as cosmic birefringence. We present our optical modeling and measurements associated with calibrating the detector angles in ACTPol.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is a co-planar interferometer array operating at a wavelength of 3mm to measure the Sunyaev-Zeldovich effect (SZE) of galaxy clusters. In the first phase of operation -- with a compact 7-element array with 0.6m antennas (AMiBA-7) -- we observed six clusters at angular scales from 5\arcmin to 23\arcmin. Here, we describe the expansion of AMiBA to a 13-element array with 1.2m antennas (AMiBA-13), its subsequent commissioning, and our cluster SZE observing program. The most important changes compared to AMiBA-7 are (1) array re-configuration with baselines ranging from 1.4m to 4.8m covering angular scales from 2\arcmin to 11.5\arcmin, (2) thirteen new lightweight carbon-fiber-reinforced plastic (CFRP) 1.2m reflectors, and (3) additional correlators and six new receivers. From the AMiBA-13 SZE observing program, we present here maps of a subset of twelve clusters. In highlights, we combine AMiBA-7 and AMiBA-13 observations of Abell 1689 and perform a joint fitting assuming a generalized NFW pressure profile. Our cylindrically integrated Compton-y values for this cluster are consistent with the BIMA/OVRA, SZA, and Planck results. We report the first targeted SZE detection towards the optically selected galaxy cluster RCS J1447+0828, and we demonstrate the ability of AMiBA SZE data to serve as a proxy for the total cluster mass. Finally, we show that our AMiBA-SZE derived cluster masses are consistent with recent lensing mass measurements in the literature.

We search for a galaxy clustering bias due to a modulation of galaxy number with the baryon-dark matter relative velocity resulting from recombination-era physics. We find no detected signal and place the constraint $b_v < 0.01$ on the relative velocity bias for the CMASS galaxies. This bias is an important potential systematic of Baryon Acoustic Oscillation (BAO) method measurements of the cosmic distance scale using the 2-point clustering. Our limit on the relative velocity bias indicates a systematic shift of no more than $0.3\%$ rms in the distance scale inferred from the BAO feature in the BOSS 2-point clustering, well below the $1\%$ statistical error of this measurement. This constraint is the most stringent currently available and has important implications for the ability of upcoming large-scale structure surveys such as DESI to self-protect against the relative velocity as a possible systematic.

We present the large-scale 3-point correlation function (3PCF) of the SDSS DR12 CMASS sample of $777,202$ Luminous Red Galaxies, the largest-ever sample used for a 3PCF or bispectrum measurement. We make the first high-significance ($4.5\sigma$) detection of Baryon Acoustic Oscillations (BAO) in the 3PCF. Using these acoustic features in the 3PCF as a standard ruler, we measure the distance to $z=0.57$ to $1.7\%$ precision (statistical plus systematic). We find $D_{\rm V}= 2024\pm29\;{\rm Mpc\;(stat)}\pm20\;{\rm Mpc\;(sys)}$ for our fiducial cosmology (consistent with Planck 2015) and bias model. This measurement extends the use of the BAO technique from the 2-point correlation function (2PCF) and power spectrum to the 3PCF and opens an avenue for deriving additional cosmological distance information from future large-scale structure redshift surveys such as DESI. Our measured distance scale from the 3PCF is fairly independent from that derived from the pre-reconstruction 2PCF and is equivalent to increasing the length of BOSS by roughly 10\%; reconstruction appears to lower the independence of the distance measurements. Fitting a model including tidal tensor bias yields a moderate significance ($2.6\sigma)$ detection of this bias with a value in agreement with the prediction from local Lagrangian biasing.

Advanced ACTPol is an instrument upgrade for the six-meter Atacama Cosmology Telescope (ACT) designed to measure the cosmic microwave background (CMB) temperature and polarization with arcminute-scale angular resolution. To achieve its science goals, Advanced ACTPol utilizes a larger readout multiplexing factor than any previous CMB experiment to measure detector arrays with approximately two thousand transition-edge sensor (TES) bolometers in each 150 mm detector wafer. We present the implementation and testing of the Advanced ACTPol time-division multiplexing readout architecture with a 64-row multiplexing factor. This includes testing of individual multichroic detector pixels and superconducting quantum interference device (SQUID) multiplexing chips as well as testing and optimizing of the integrated readout electronics. In particular, we describe the new automated multiplexing SQUID tuning procedure developed to select and optimize the thousands of SQUID parameters required to readout each Advanced ACTPol array. The multichroic detector pixels in each array use separate channels for each polarization and each of the two frequencies, such that four TESes must be read out per pixel. Challenges addressed include doubling the number of detectors per multiplexed readout channel compared to ACTPol and optimizing the Nyquist inductance to minimize detector and SQUID noise aliasing.

The next generation Advanced ACTPol (AdvACT) experiment is currently underway and will consist of four Transition Edge Sensor (TES) bolometer arrays, with three operating together, totaling ~5800 detectors on the sky. Building on experience gained with the ACTPol detector arrays, AdvACT will utilize various new technologies, including 150mm detector wafers equipped with multichroic pixels, allowing for a more densely packed focal plane. Each set of detectors includes a feedhorn array of stacked silicon wafers which form a spline profile leading to each pixel. This is then followed by a waveguide interface plate, detector wafer, back short cavity plate, and backshort cap. Each array is housed in a custom designed structure manufactured from high purity copper and then gold plated. In addition to the detector array assembly, the array package also encloses cryogenic readout electronics. We present the full mechanical design of the AdvACT high frequency (HF) detector array package along with a detailed look at the detector array stack assemblies. This experiment will also make use of extensive hardware and software previously developed for ACT, which will be modified to incorporate the new AdvACT instruments. Therefore, we discuss the integration of all AdvACT arrays with pre-existing ACTPol infrastructure.

We analyse the broad-range shape of the monopole and quadrupole correlation functions of the BOSS Data Release 12 (DR12) CMASS and LOWZ galaxy sample to obtain constraints on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalised growth rate $f(z)\sigma_8(z)$, and the physical matter density $\Omega_mh^2$. We adopt wide and flat priors on all model parameters in order to ensure the results are those of a `single-probe' galaxy clustering analysis. We also marginalise over three nuisance terms that account for potential observational systematics affecting the measured monopole. However, such Monte Carlo Markov Chain analysis is computationally expensive for advanced theoretical models, thus we develop a new methodology to speed up our analysis. We obtain $\{D_A(z)r_{s,fid}/r_s$Mpc, $H(z)r_s/r_{s,fid}$kms$^{-1}$Mpc$^{-1}$, $f(z)\sigma_8(z)$, $\Omega_m h^2\}$ = $\{956\pm28$ , $75.0\pm4.0$ , $0.397 \pm 0.073$, $0.143\pm0.017\}$ at $z=0.32$ and $\{1421\pm23$, $96.7\pm2.7$ , $0.497 \pm 0.058$, $0.137\pm0.015\}$ at $z=0.59$ where $r_s$ is the comoving sound horizon at the drag epoch and $r_{s,fid}=147.66$Mpc for the fiducial cosmology in this study. In addition, we divide the galaxy sample into four redshift bins to increase the sensitivity of redshift evolution. However, we do not find improvements in terms of constraining dark energy model parameters. Combining our measurements with Planck data, we obtain $\Omega_m=0.306\pm0.009$, $H_0=67.9\pm0.7$kms$^{-1}$Mpc$^{-1}$, and $\sigma_8=0.815\pm0.009$ assuming $\Lambda$CDM; $\Omega_k=0.000\pm0.003$ assuming oCDM; $w=-1.01\pm0.06$ assuming $w$CDM; and $w_0=-0.95\pm0.22$ and $w_a=-0.22\pm0.63$ assuming $w_0w_a$CDM. Our results show no tension with the flat $\Lambda$CDM cosmological paradigm. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

We develop a new methodology called double-probe analysis with the aim of minimizing informative priors in the estimation of cosmological parameters. We extract the dark-energy-model-independent cosmological constraints from the joint data sets of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy sample and Planck cosmic microwave background (CMB) measurement. We measure the mean values and covariance matrix of $\{R$, $l_a$, $\Omega_b h^2$, $n_s$, $log(A_s)$, $\Omega_k$, $H(z)$, $D_A(z)$, $f(z)\sigma_8(z)\}$, which give an efficient summary of Planck data and 2-point statistics from BOSS galaxy sample, where $R=\sqrt{\Omega_m H_0^2}\,r(z_*)$, and $l_a=\pi r(z_*)/r_s(z_*)$, $z_*$ is the redshift at the last scattering surface, and $r(z_*)$ and $r_s(z_*)$ denote our comoving distance to $z_*$ and sound horizon at $z_*$ respectively. The advantage of this method is that we do not need to put informative priors on the cosmological parameters that galaxy clustering is not able to constrain well, i.e. $\Omega_b h^2$ and $n_s$. Using our double-probe results, we obtain $\Omega_m=0.304\pm0.009$, $H_0=68.2\pm0.7$, and $\sigma_8=0.806\pm0.014$ assuming $\Lambda$CDM; and $\Omega_k=0.002\pm0.003$ and $w=-1.00\pm0.07$ assuming o$w$CDM. The results show no tension with the flat $\Lambda$CDM cosmological paradigm. By comparing with the full-likelihood analyses with fixed dark energy models, we demonstrate that the double-probe method provides robust cosmological parameter constraints which can be conveniently used to study dark energy models. We extend our study to measure the sum of neutrino mass and obtain $\Sigma m_\nu<0.10/0.22$ (68\%/95\%) assuming $\Lambda$CDM and $\Sigma m_\nu<0.26/0.52$ (68\%/95\%) assuming $w$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

We present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg^2 and volume of 18.7 Gpc^3, divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51, and 0.61. We measure the angular diameter distance DM and Hubble parameter H from the baryon acoustic oscillation (BAO) method after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product DM*H from the Alcock-Paczynski (AP) effect and the growth of structure, quantified by f{\sigma}8(z), from redshift-space distortions (RSD). We combine measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method. Combined with Planck 2015 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature {\Omega}_K =0.0003+/-0.0026 and a dark energy equation of state parameter w = -1.01+/-0.06, in strong affirmation of the spatially flat cold dark matter model with a cosmological constant ({\Lambda}CDM). Our RSD measurements of f{\sigma}_8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H0 = 67.3+/-1.0 km/s/Mpc even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H0 = 67.8+/-1.2 km/s/Mpc. Assuming flat {\Lambda}CDM we find {\Omega}_m = 0.310+/-0.005 and H0 = 67.6+/-0.5 km/s/Mpc, and we find a 95% upper limit of 0.16 eV/c^2 on the neutrino mass sum.

With the largest spectroscopic galaxy survey volume drawn from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), we can extract cosmological constraints from the measurements of redshift and geometric distortions at quasi-linear scales (e.g. above 50 $h^{-1}$Mpc). We analyze the broad-range shape of the monopole and quadrupole correlation functions of the BOSS Data Release 12 (DR12) CMASS galaxy sample, at the effective redshift $z=0.59$, to obtain constraints on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalized growth rate $f(z)\sigma_8(z)$, and the physical matter density $\Omega_mh^2$. We obtain robust measurements by including a polynomial as the model for the systematic errors, and find it works very well against the systematic effects, e.g., ones induced by stars and seeing. We provide accurate measurements $\{D_A(0.59)r_{s,fid}/r_s$ $\rm Mpc$, $H(0.59)r_s/r_{s,fid}$ $km s^{-1} Mpc^{-1}$, $f(0.59)\sigma_8(0.59)$, $\Omega_m h^2\}$ = $\{1427\pm26$, $97.3\pm3.3$, $0.488 \pm 0.060$, $0.135\pm0.016\}$, where $r_s$ is the comoving sound horizon at the drag epoch and $r_{s,fid}=147.66$ Mpc is the sound scale of the fiducial cosmology used in this study. The parameters which are not well constrained by our galaxy clustering analysis are marginalized over with wide flat priors. Since no priors from other data sets, e.g., cosmic microwave background (CMB), are adopted and no dark energy models are assumed, our results from BOSS CMASS galaxy clustering alone may be combined with other data sets, i.e., CMB, SNe, lensing or other galaxy clustering data to constrain the parameters of a given cosmological model. The uncertainty on the dark energy equation of state parameter, $w$, from CMB+CMASS is about 8 per cent. The uncertainty on the curvature fraction, $\Omega_k$, is 0.3 per cent. We do not find deviation from flat $\Lambda$CDM.

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.

We present Intercut, a Python-based program that applies secondary line identification and photometric cuts to mock galaxy surveys, in order to simulate interloper identification. This program can be used to optimize the removal of interloper contamination in upcoming surveys. Intercut reads a mock galaxy survey and an emission line sensitivity and simulates interloper removal through secondary line identification and broad-band photometry. This program is designed to use the COSMOS mock catalog, although the program can be modified for an alternative mock catalog. The output of the program returns an interloper fraction for each emission line, as well as the total fraction over all lines, as a function of redshift. We test Intercut by predicting interloper rates for the WFIRST emission line sensitivity, finding agreement with previous work. This program is publically available on Github

[abridged] We present an anisotropic analysis of the baryonic acoustic oscillation (BAO) scale in the twelfth and final data release of the Baryonic Oscillation Spectroscopic Survey (BOSS). We independently analyse the LOWZ and CMASS galaxy samples: the LOWZ sample contains contains 361 762 galaxies with an effective redshift of $z_{\rm LOWZ}=0.32$; the CMASS sample consists of 777 202 galaxies with an effective redshift of $z_{\rm CMASS}=0.57$. We extract the BAO peak position from the monopole power spectrum moment, $\alpha_0$, and from the $\mu^2$ moment, $\alpha_2$, where $\mu$ is the cosine of the angle to the line-of-sight. The $\mu^2$-moment provides equivalent information to that available in the quadrupole but is simpler to analyse. After applying a reconstruction algorithm to reduce the BAO suppression by bulk motions, we measure the BAO peak position in the monopole and $\mu^2$-moment, which are related to radial and angular shifts in scale. We report $H(z_{\rm LOWZ})r_s(z_d)=(11.60\pm0.60)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.66\pm0.16$ with a cross-correlation coefficient of $r_{HD_A}=0.41$, for the LOWZ sample; and $H(z_{\rm CMASS})r_s(z_d)=(14.56\pm0.37)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.42\pm0.13$ with a cross-correlation coefficient of $r_{HD_A}=0.47$, for the CMASS sample. We combine these results with the measurements of the BAO peak position in the monopole and quadrupole correlation function of the same dataset \citep[][companion paper]{Cuestaetal2015} and report the consensus values: $H(z_{\rm LOWZ})r_s(z_d)=(11.63\pm0.69)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.67\pm0.15$ with $r_{HD_A}=0.35$ for the LOWZ sample; $H(z_{\rm CMASS})r_s(z_d)=(14.67\pm0.42)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.47\pm0.12$ with $r_{HD_A}=0.52$ for the CMASS sample.

We measure and analyse the clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) relative to the line-of-sight (LOS), for LOWZ and CMASS galaxy samples drawn from the final Data Release 12 (DR12). The LOWZ sample contains 361\,762 galaxies with an effective redshift of $z_{\rm lowz}=0.32$, and the CMASS sample 777\,202 galaxies with an effective redshift of $z_{\rm cmass}=0.57$. From the power spectrum monopole and quadrupole moments around the LOS, we measure the growth of structure parameter $f$ times the amplitude of dark matter density fluctuations $\sigma_8$ by modeling the Redshift-Space Distortion signal. When the geometrical Alcock-Paczynski effect is also constrained from the same data, we find joint constraints on $f\sigma_8$, the product of the Hubble constant and the comoving sound horizon at the baryon drag epoch $H(z)r_s(z_d)$, and the angular distance parameter divided by the sound horizon $D_A(z)/r_s(z_d)$. We find $f(z_{\rm lowz})\sigma_8(z_{\rm lowz})=0.394\pm0.062$, $D_A(z_{\rm lowz})/r_s(z_d)=6.35\pm0.19$, $H(z_{\rm lowz})r_s(z_d)=(11.41\pm 0.56)\,{10^3\rm km}s^{-1}$ for the LOWZ sample, and $f(z_{\rm cmass})\sigma_8(z_{\rm cmass})=0.444\pm0.038$, $D_A(z_{\rm cmass})/r_s(z_d)=9.42\pm0.15$, $H(z_{\rm cmass})r_s(z_d)=(13.92 \pm 0.44)\, {10^3\rm km}s^{-1}$ for the CMASS sample. We find general agreement with previous BOSS DR11 measurements. Assuming the Hubble parameter and angular distance parameter are fixed at fiducial $\Lambda$CDM values, we find $f(z_{\rm lowz})\sigma_8(z_{\rm lowz})=0.485\pm0.044$ and $f(z_{\rm cmass})\sigma_8(z_{\rm cmass})=0.436\pm0.022$ for the LOWZ and CMASS samples, respectively.

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

The most recent measurements of the temperature and low-multipole polarization anisotropies of the Cosmic Microwave Background (CMB) from the Planck satellite, when combined with galaxy clustering data from the Baryon Oscillation Spectroscopic Survey (BOSS) in the form of the full shape of the power spectrum, and with Baryon Acoustic Oscillation measurements, provide a $95\%$ confidence level (CL) upper bound on the sum of the three active neutrinos $\sum m _\nu< 0.183$ eV, among the tightest neutrino mass bounds in the literature, to date, when the same datasets are taken into account. This very same data combination is able to set, at $\sim70\%$ CL, an upper limit on $\sum m _\nu$ of $0.0968$ eV, a value that approximately corresponds to the minimal mass expected in the inverted neutrino mass hierarchy scenario. If high-multipole polarization data from Planck is also considered, the $95\%$ CL upper bound is tightened to $\sum m _\nu< 0.176$ eV. Further improvements are obtained by considering recent measurements of the Hubble parameter. These limits are obtained assuming a specific non-degenerate neutrino mass spectrum; they slightly worsen when considering other degenerate neutrino mass schemes. Current cosmological data, therefore, start to be mildly sensitive to the neutrino mass ordering. Low-redshift quantities, such as the Hubble constant or the reionization optical depth, play a very important role when setting the neutrino mass constraints. We also comment on the eventual shifts in the cosmological bounds on $\sum m_\nu$ when possible variations in the former two quantities are addressed.

We present distance scale measurements from the baryon acoustic oscillation signal in the CMASS and LOWZ samples from the Data Release 12 of the Baryon Oscillation Spectroscopic Survey (BOSS). The total volume probed is 14.5 Gpc$^3$, a 10 per cent increment from Data Release 11. From an analysis of the spherically averaged correlation function, we infer a distance to $z=0.57$ of $D_V(z)r^{\rm fid}_{\rm d}/r_ {\rm d}=2028\pm21$ Mpc and a distance to $z=0.32$ of $D_V(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1264\pm22$ Mpc assuming a cosmology in which $r^{\rm fid}_{\rm d}=147.10$ Mpc. From the anisotropic analysis, we find an angular diameter distance to $z=0.57$ of $D_{\rm A}(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1401\pm21$ Mpc and a distance to $z=0.32$ of $981\pm20$ Mpc, a 1.5 per cent and 2.0 per cent measurement respectively. The Hubble parameter at $z=0.57$ is $H(z)r_{\rm d}/r^{\rm fid}_{\rm d}=100.3\pm3.7$ km s$^{-1}$ Mpc$^{-1}$ and its value at $z=0.32$ is $79.2\pm5.6$ km s$^{-1}$ Mpc$^{-1}$, a 3.7 per cent and 7.1 per cent measurement respectively. These cosmic distance scale constraints are in excellent agreement with a $\Lambda$CDM model with cosmological parameters released by the recent Planck 2015 results.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

Modern cosmology predicts that a galaxy overdensity is associated to a large reservoir of the intergalactic gas, which can be traced by the Ly$\alpha$ forest absorption. We have undertaken a systematic study of the relation between Coherently Strong intergalactic Ly$\alpha$ Absorption systems (CoSLAs), which have highest optical depth ($\tau$) in $\tau$ distribution, and mass overdensities on the scales of $\sim$ 10 - 20 $h^{-1}$ comoving Mpc. On such large scales, our cosmological simulations show a strong correlation between the effective optical depth ($\tau_{\rm{eff}}$) of the CoSLAs and the 3-D mass overdensities. In moderate signal-to-noise spectra, however, the profiles of CoSLAs can be confused with high column density absorbers. For $z>2.6$, where the corresponding Ly$\beta$ is redshifted to the optical, we have developed the technique to differentiate between these two alternatives. We have applied this technique to SDSS-III quasar survey at $z = 2.6$ - 3.3, and we present a sample of five CoSLA candidates with $\tau_{\rm{eff}}$ on 15 $h^{-1}$ Mpc greater than $4.5\times$ the mean optical depth. At lower redshifts of $z < 2.6$, where the background quasar density is higher, the overdensity can be traced by intergalactic absorption groups using multiple sight lines. Our overdensity searches fully utilize the current and next generation of Ly$\alpha$ forest surveys which cover a survey volume of $> (1\ h^{-1}$ Gpc)$^3$. In addition, systems traced by CoSLAs will build a uniform sample of the most massive overdensities at $z > 2$ to constrain the models of structure formation, and offer a unique laboratory to study the interactions between galaxy overdensities and the intergalactic medium.

We report a measurement of the large-scale 3-point correlation function of galaxies using the largest dataset for this purpose to date, 777, 202 Luminous Red Galaxies in the Sloan Digital Sky Survey Baryon Acoustic Oscillation Spectroscopic Survey (SDSS BOSS) DR12 CMASS sample. This work exploits the novel algorithm of Slepian & Eisenstein (2015b) to compute the multipole moments of the 3PCF in $\mathcal{O}(N^2)$ time, with $N$ the number of galaxies. Leading-order perturbation theory models the data well in a compressed basis where one triangle side is integrated out. We also present an accurate and computationally efficient means of estimating the covariance matrix. With these techniques the redshift-space linear and non-linear bias are measured, with 2.6% precision on the former if $\sigma_8$ is fixed. The data also indicates a $2.8\sigma$ preference for the BAO, confirming the presence of BAO in the 3-point function.

We describe mock data-sets generated to simulate the high-redshift quasar sample in Data Release 11 (DR11) of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). The mock spectra contain Ly{\alpha} forest correlations useful for studying the 3D correlation function including Baryon Acoustic Oscillations (BAO). They also include astrophysical effects such as quasar continuum diversity and high-density absorbers, instrumental effects such as noise and spectral resolution, as well as imperfections introduced by the SDSS pipeline treatment of the raw data. The Ly{\alpha} forest BAO analysis of the BOSS collaboration, described in Delubac et al. 2014, has used these mock data-sets to develop and cross-check analysis procedures prior to performing the BAO analysis on real data, and for continued systematic cross checks. Tests presented here show that the simulations reproduce sufficiently well important characteristics of real spectra. These mock data-sets will be made available together with the data at the time of the Data Release 11.

The Atacama B-Mode Search (ABS) instrument is a cryogenic ($\sim$10 K) crossed-Dragone telescope located at an elevation of 5190 m in the Atacama Desert in Chile that observed for three seasons between February 2012 and October 2014. ABS observed the Cosmic Microwave Background (CMB) at large angular scales ($40<\ell<500$) to limit the B-mode polarization spectrum around the primordial B-mode peak from inflationary gravity waves at $\ell \sim100$. The ABS focal plane consists of 480 transition-edge sensor (TES) bolometers. They are coupled to orthogonal polarizations from a planar ortho-mode transducer (OMT) and observe at 145 GHz. ABS employs an ambient-temperature, rapidly rotating half-wave plate (HWP) to mitigate systematic effects and move the signal band away from atmospheric $1/f$ noise, allowing for the recovery of large angular scales. We discuss how the signal at the second harmonic of the HWP rotation frequency can be used for data selection and for monitoring the detector responsivities.

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

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.

The Baryon Oscillation Spectroscopic Survey (BOSS), part of the Sloan Digital Sky Survey (SDSS) III project, has provided the largest survey of galaxy redshifts available to date, in terms of both the number of galaxy redshifts measured by a single survey, and the effective cosmological volume covered. Key to analysing the clustering of these data to provide cosmological measurements is understanding the detailed properties of this sample. Potential issues include variations in the target catalogue caused by changes either in the targeting algorithm or properties of the data used, the pattern of spectroscopic observations, the spatial distribution of targets for which redshifts were not obtained, and variations in the target sky density due to observational systematics. We document here the target selection algorithms used to create the galaxy samples that comprise BOSS. We also present the algorithms used to create large scale structure catalogues for the final Data Release (DR12) samples and the associated random catalogues that quantify the survey mask. The algorithms are an evolution of those used by the BOSS team to construct catalogues from earlier data, and have been designed to accurately quantify the galaxy sample. The code used, designated MKSAMPLE, is released with this paper.

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

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

We introduce a new method for estimating the covariance matrix for the galaxy correlation function in surveys of large-scale structure. Our method combines simple theoretical results with a realistic characterization of the survey to dramatically reduce noise in the covariance matrix. For example, with an investment of only ~1,000 CPU hours we can produce a model covariance matrix with noise levels that would otherwise require ~35,000 mocks. Non-Gaussian contributions to the model are calibrated against mock catalogs, after which the model covariance is found to be in impressive agreement with the mock covariance matrix. Since calibration of this method requires fewer mocks than brute force approaches, we believe that it could dramatically reduce the number of mocks required to analyse future surveys.

We use measurements from the Planck satellite mission and galaxy redshift surveys over the last decade to test three of the basic assumptions of the standard model of cosmology, $\Lambda$CDM: the spatial curvature of the universe, the nature of dark energy and the laws of gravity on large scales. We obtain improved constraints on several scenarios that violate one or more of these assumptions. We measure $w_0=-0.94\pm0.17$ (18\% measurement) and $1+w_a=1.16\pm0.36$ (31\% measurement) for models with a time-dependent equation of state, which is an improvement over current best constraints \citep{Aubourg2014}. In the context of modified gravity, we consider popular scalar tensor models as well as a parametrization of the growth factor. In the case of one-parameter $f(R)$ gravity models with a $\Lambda$CDM background, we constrain $B_0 < 1.36 \times 10^{-5} $ (1$\sigma$ C.L.), which is an improvement by a factor of 4 on the current best \citep{XU2015}. We provide the very first constraint on the coupling parameters of general scalar-tensor theory and stringent constraint on the only free coupling parameter of Chameleon models. We also derive constraints on extended Chameleon models, improving the constraint on the coupling by a factor of 6 on the current best \citep{Hojjati2011} . We also measure $\gamma = 0.612 \pm 0.072$ (11.7\% measurement) for growth index parametrization. We improve all the current constraints by combining results from various galaxy redshift surveys in a coherent way, which includes a careful treatment of scale-dependence introduced by modified gravity.

We construct a catalogue for filaments using a novel approach called SCMS (subspace constrained mean shift; Ozertem & Erdogmus 2011; Chen et al. 2015). SCMS is a gradient-based method that detects filaments through density ridges (smooth curves tracing high-density regions). A great advantage of SCMS is its uncertainty measure, which allows an evaluation of the errors for the detected filaments. To detect filaments, we use data from the Sloan Digital Sky Survey, which consist of three galaxy samples: the NYU main galaxy sample (MGS), the LOWZ sample and the CMASS sample. Each of the three dataset covers different redshift regions so that the combined sample allows detection of filaments up to z = 0.7. Our filament catalogue consists of a sequence of two-dimensional filament maps at different redshifts that provide several useful statistics on the evolution cosmic web. To construct the maps, we select spectroscopically confirmed galaxies within 0.050 < z < 0.700 and partition them into 130 bins. For each bin, we ignore the redshift, treating the galaxy observations as a 2-D data and detect filaments using SCMS. The filament catalogue consists of 130 individual 2-D filament maps, and each map comprises points on the detected filaments that describe the filamentary structures at a particular redshift. We also apply our filament catalogue to investigate galaxy luminosity and its relation with distance to filament. Using a volume-limited sample, we find strong evidence (6.1$\sigma$ - 12.3$\sigma$) that galaxies close to filaments are generally brighter than those at significant distance from filaments.

The measured redshift ($z$) of an astronomical object is a combination of Hubble recession, gravitational redshift and peculiar velocity. The line of sight distance to a galaxy inferred from redshift is affected by the peculiar velocity component of galaxy redshift, which is observed as an anisotropy in the correlation function. This anisotropy allows us to measure the linear growth rate of matter ($f\sigma_8$). We measure the $f\sigma_8$ at $z=0.57$ using the CMASS sample from Data Release 11 of Sloan Digital Sky Survey III (SDSS III) Baryon Oscillations Spectroscopic Survey (BOSS). The galaxy sample consists of 690,826 massive galaxies in the redshift range 0.43-0.7 covering 8498 deg$^2$. Here we report the first simultaneous measurement of $f\sigma_8$ and background cosmological parameters using Convolution Lagrangian Perturbation Theory (CLPT) with Gaussian streaming model (GSRSD). We arrive at a constraint of $f\sigma_8=0.462\pm0.041$ (9\% accuracy) at effective redshift ($\bar{z}=0.57$) when we include Planck CMB likelihood while marginalizing over all other cosmological parameters. We also measure $b\sigma_8=1.19\pm0.03$, $H(z=0.57)=89.2\pm3.6$ km s$^{-1}$ Mpc$^{-1}$ and $D_A(z=0.57)=1401\pm23$ Mpc. Our analysis also improves the constraint on $\Omega_c h^2=0.1196\pm0.0009$ by a factor of 3 when compared to the Planck only measurement($\Omega_c h^2=0.1196 \pm 0.0031$). Our results are consistent with Planck $\Lambda$CDM-GR prediction and all other CMASS measurements, even though our theoretical models are fairly different. This consistency suggests that measurement of $f\sigma_8$ from Redshift space distortions at multiple redshifts will be a sensitive probe of the theory of gravity that is largely model independent, allowing us to place model-independent constraints on alternative models of gravity.

The detection and characterization of filamentary structures in the cosmic web allows cosmologists to constrain parameters that dictates the evolution of the Universe. While many filament estimators have been proposed, they generally lack estimates of uncertainty, reducing their inferential power. In this paper, we demonstrate how one may apply the Subspace Constrained Mean Shift (SCMS) algorithm (Ozertem and Erdogmus (2011); Genovese et al. (2012)) to uncover filamentary structure in galaxy data. The SCMS algorithm is a gradient ascent method that models filaments as density ridges, one-dimensional smooth curves that trace high-density regions within the point cloud. We also demonstrate how augmenting the SCMS algorithm with bootstrap-based methods of uncertainty estimation allows one to place uncertainty bands around putative filaments. We apply the SCMS method to datasets sampled from the P3M N-body simulation, with galaxy number densities consistent with SDSS and WFIRST-AFTA and to LOWZ and CMASS data from the Baryon Oscillation Spectroscopic Survey (BOSS). To further assess the efficacy of SCMS, we compare the relative locations of BOSS filaments with galaxy clusters in the redMaPPer catalog, and find that redMaPPer clusters are significantly closer (with p-values $< 10^{-9}$) to SCMS-detected filaments than to randomly selected galaxies.

We examine the luminosity function (LF) of [OII] emission-line galaxies in the high-resolution cosmological simulation MassiveBlack-II (MBII). From the spectral energy distribution of each galaxy, we select a sub-sample of star-forming galaxies at $0.06 \le z \le 3.0$ using the [OII] emission line luminosity L([OII]). We confirm that the specific star formation rate matches that in the GAMA survey. We show that the [OII] LF at z=1.0 from the MBII shows a good agreement with the LFs from several surveys below L([OII])=$10^{43.0}$ erg/s while the low redshifts ($z \le 0.3$) show an excess in the prediction of bright [OII] galaxies, but still displaying a good match with observations below L([OII])=$10^{41.6}$ erg/s. Based on the validity in reproducing the properties of [OII] galaxies at low redshift ($z \le 1$), we forecast the evolution of the [OII] LF at high redshift ($z \le 3$), which can be tested by upcoming surveys such as the HETDEX and DESI. The slopes of the LFs at bright and faint ends range from -3 to -2 showing minima at z=2. The slope of the bright end evolves approximately as 1/(z+1) at z=2 while the faint end evolves as ~3/(z+1) at $0.6 \le z \le 2$. In addition, a similar analysis is applied for the evolution of [OIII] LFs, which is to be explored in the forthcoming survey WFIRST-AFTA. Finally, we show that the auto-correlation function of [OII] and [OIII] emitting galaxies shows a rapid evolution from z=2 to 1.

We describe the algorithm used to select the Luminous Red Galaxy (LRG) sample for the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV) using photometric data from both the SDSS and the Wide-Field Infrared Survey Explorer (WISE). LRG targets are required to meet a set of color selection criteria and have z-band and i-band MODEL magnitudes z < 19.95 and 19.9 < i < 21.8, respectively. Our algorithm selects roughly 50 LRG targets per square degree, the great majority of which lie in the redshift range 0.6 < z < 1.0 (median redshift 0.71). We demonstrate that our methods are highly effective at eliminating stellar contamination and lower-redshift galaxies. We perform a number of tests using spectroscopic data from SDSS-III/BOSS to determine the redshift reliability of our target selection and its ability to meet the science requirements of eBOSS. The SDSS spectra are of high enough signal-to-noise ratio that at least 89% of the target sample yields secure redshift measurements. We also present tests of the uniformity and homogeneity of the sample, demonstrating that it should be clean enough for studies of the large-scale structure of the universe at higher redshifts than SDSS-III/BOSS LRGs reached.

In this paper, we study the filamentary structures and the galaxy alignment along filaments at redshift $z=0.06$ in the MassiveBlack-II simulation, a state-of-the-art, high-resolution hydrodynamical cosmological simulation which includes stellar and AGN feedback in a volume of (100 Mpc$/h$)$^3$. The filaments are constructed using the subspace constrained mean shift (SCMS; Ozertem & Erdogmus (2011) and Chen et al. (2015a)). First, we show that reconstructed filaments using galaxies and reconstructed filaments using dark matter particles are similar to each other; over $50\%$ of the points on the galaxy filaments have a corresponding point on the dark matter filaments within distance $0.13$ Mpc$/h$ (and vice versa) and this distance is even smaller at high-density regions. Second, we observe the alignment of the major principal axis of a galaxy with respect to the orientation of its nearest filament and detect a $2.5$ Mpc$/h$ critical radius for filament's influence on the alignment when the subhalo mass of this galaxy is between $10^9M_\odot/h$ and $10^{12}M_\odot/h$. Moreover, we find the alignment signal to increase significantly with the subhalo mass. Third, when a galaxy is close to filaments (less than $0.25$ Mpc$/h$), the galaxy alignment toward the nearest galaxy group depends on the galaxy subhalo mass. Finally, we find that galaxies close to filaments or groups tend to be rounder than those away from filaments or groups.

We describe a methodology to probe gravity with the cosmic microwave background (CMB) lensing convergence $\kappa$, specifically by measuring $E_G$, the ratio of the Laplacian of the gravitational scalar potential difference with the velocity divergence. Using CMB lensing instead of galaxy-galaxy lensing avoids intrinsic alignments while also lacking a hard limit on the lens redshift and significant uncertainties in the source plane. We model $E_G$ for general relativity and modified gravity, finding that $E_G$ for $f(R)$ gravity should be scale-dependent due to the scale-dependence of the growth rate $f$. Next, we construct an estimator for $E_G$ in terms of the galaxy-CMB lensing and galaxy clustering angular power spectra, along with the RSD parameter $\beta$. We also forecast statistical errors for $E_G$ from the current Planck CMB lensing map and the spectroscopic galaxy and quasar samples from the Sloan Digital Sky Survey Data Release 11, being 9% with galaxies and 8% when quasars are included. We also find that upcoming spectroscopic and photometric surveys, combined with the final Planck lensing map, can measure precisely the redshift- and scale-dependence of $E_G$ out to redshifts $z=2$ and higher, with photometric surveys having an advantage due to their high number densities. Advanced ACTPol's lensing map will increase the $E_G$ sensitivity even further. Finally, we find that Advanced ACTPol cross-correlated with spectroscopic (photometric) surveys can differentiate between general relativity and $f(R)$ gravity at the level of $3\sigma$ ($13\sigma$). Performing a $<1$% measurement of $E_G$ requires a 10% precision in $\beta$ from Euclid or LSST, currently achievable with a spectroscopic survey but difficult with only a photometric survey.

Many astronomers now participate in large international collaborations, and it is important to examine whether these structures foster a scientific climate that is inclusive and diverse. The Committee on the Participation of Women in the Sloan Digital Sky Survey (CPWS) was formed to evaluate the demographics and gender climate within SDSS-IV, one of the largest and most geographically distributed astronomical collaborations. In April 2014, the CPWS administered a demographic survey to establish a baseline for the incipient SDSS-IV. We received responses from 250 participants (46% of the active membership). Half of the survey respondents were located in the US or Canada and 30% were based in Europe. Eleven percent of survey respondents considered themselves to be an ethnic minority at their current institution. Twenty-five percent of the SDSS-IV collaboration members are women, a fraction that is consistent with the US astronomical community, but substantially higher than the fraction of women in the IAU (16%). Approximately equal fractions of men and women report holding positions of leadership. When binned by academic age and career level, men and women also assume leadership roles at approximately equal rates, in a way that increases steadily for both genders with increasing seniority. In this sense, SDSS-IV has been successful in recruiting leaders that are representative of the collaboration. Yet, more progress needs to be made towards achieving gender balance and increasing diversity in the field of astronomy, and there is still room for improvement in the membership and leadership of SDSS-IV. For example, at the highest level of SDSS-IV leadership, women disproportionately assume roles related to education and public outreach. The CPWS plans to use these initial data to establish a baseline for tracking demographics over time as we work to assess and improve the climate of SDSS-IV.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We investigate the utility and robustness of a new statistic, $\omega_{\ell}\left(r_{c}\right)$, for analyzing Baryon Acoustic Oscillations (BAO). We apply $\omega_{\ell}\left(r_{c}\right)$, introduced in Xu et al. (2010), to mocks and data from the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) included in the SDSS Data Release Eleven (DR11). We fit the anisotropic clustering using the monopole and quadrupole of the $\omega_{\ell}\left(r_{c}\right)$ statistic in a manner similar to conventional multipole fitting methods using the correlation function as detailed in (Xu et al. 2012). To test the performance of the $\omega_{\ell}\left(r_{c}\right)$ statistic we compare our results to those obtained using the multipoles. The results are in agreement. We also conduct a brief investigation into some of the possible advantages of using the $\omega_{\ell}\left(r_{c}\right)$ statistic for BAO analysis. The $\omega_{\ell}\left(r_{c}\right)$ analysis matches the stability of the multipoles analysis in response to artificially introduced distortions in the data, without using extra nuisance parameters to improve the fit. When applied to data with systematics, the $\omega_{\ell}\left(r_{c}\right)$ statistic again matches the performance of fitting the multipoles without using nuisance parameters. In all the analyzed circumstances, we find that fitting the $\omega_{\ell}\left(r_{c}\right)$ statistic removes the requirement for extra nuisance parameters.

The Sloan Digital Sky Survey has surveyed 14,555 square degrees of the sky, and delivered over a trillion pixels of imaging data. We present the large-scale clustering of 1.6 million quasars between z = 0.5 and z = 2.5 that have been classified from this imaging, representing the highest density of quasars ever studied for clustering measurements. This data set spans ~11,000 square degrees and probes a volume of 80(Gpc/h)^3. In principle, such a large volume and medium density of tracers should facilitate high-precision cosmological constraints. We measure the angular clustering of photometrically classified quasars using an optimal quadratic estimator in four redshift slices with an accuracy of ~25% over a bin width of l ~10 - 15 on scales corresponding to matter-radiation equality and larger (l ~ 2 - 30). Observational systematics can strongly bias clustering measurements on large scales, which can mimic cosmologically relevant signals such as deviations from Gaussianity in the spectrum of primordial perturbations. We account for systematics by employing a new method recently proposed by Agarwal et al. (2014) to the clustering of photometrically classified quasars. We carefully apply our methodology to mitigate known observational systematics and further remove angular bins that are contaminated by unknown systematics. Combining quasar data with the photometric luminous red galaxy (LRG) sample of Ross et al. (2011) and Ho et al. (2012), and marginalizing over all bias and shot noise-like parameters, we obtain a constraint on local primordial non-Gaussianity of fNL = -113+/-154 (1\sigma error). [Abridged]

Reconstructing an estimate of linear Baryon Acoustic Oscillations (BAO) from an evolved galaxy field has become a standard technique in recent analyses. By partially removing non-linear damping caused by bulk motions, the real-space BAO peak in the correlation function is sharpened, and oscillations in the power spectrum are visible to smaller scales. In turn these lead to stronger measurements of the BAO scale. Future surveys are being designed assuming that this improvement has been applied, and this technique is therefore of critical importance for future BAO measurements. A number of reconstruction techniques are available, but the most widely used is a simple algorithm that decorrelates large-scale and small-scale modes approximately removing the bulk-flow displacements by moving the overdensity field (Eisenstein et al. 2007; Padmanabhan, White & Cohn 2009). We consider the practical implementation of this algorithm, looking at the efficiency of reconstruction as a function of the assumptions made for the bulk-flow scale, the shot noise level in a random catalogue used to quantify the mask and the method used to estimate the bulk-flow shifts. We also examine the efficiency of reconstruction against external factors including galaxy density, volume and edge effects, and consider their impact for future surveys. Throughout we make use of the mocks catalogues created for the Baryon Oscillation Spectroscopic Survey (BOSS) Date Release 11 samples covering 0.43 < z < 0.7 (CMASS) and 0.15 < z < 0.43 (LOWZ), to empirically test these changes.

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.

On the arcminute angular scales probed by Planck, the CMB anisotropies are gently perturbed by gravitational lensing. Here we present a detailed study of this effect, detecting lensing independently in the 100, 143, and 217GHz frequency bands with an overall significance of greater than 25sigma. We use the temperature-gradient correlations induced by lensing to reconstruct a (noisy) map of the CMB lensing potential, which provides an integrated measure of the mass distribution back to the CMB last-scattering surface. Our lensing potential map is significantly correlated with other tracers of mass, a fact which we demonstrate using several representative tracers of large-scale structure. We estimate the power spectrum of the lensing potential, finding generally good agreement with expectations from the best-fitting LCDM model for the Planck temperature power spectrum, showing that this measurement at z=1100 correctly predicts the properties of the lower-redshift, later-time structures which source the lensing potential. When combined with the temperature power spectrum, our measurement provides degeneracy-breaking power for parameter constraints; it improves CMB-alone constraints on curvature by a factor of two and also partly breaks the degeneracy between the amplitude of the primordial perturbation power spectrum and the optical depth to reionization, allowing a measurement of the optical depth to reionization which is independent of large-scale polarization data. Discarding scale information, our measurement corresponds to a 4% constraint on the amplitude of the lensing potential power spectrum, or a 2% constraint on the RMS amplitude of matter fluctuations at z~2.

This document presents the results from the Distances subgroup of the Cosmic Frontier Community Planning Study (Snowmass 2013). We summarize the current state of the field as well as future prospects and challenges. In addition to the established probes using Type IA supernovae and baryon acoustic oscillations, we also consider prospective methods based on clusters, active galactic nuclei, gravitational wave sirens and strong lensing time delays.

We present the first observation of a galaxy (z=0.2) that exhibits metal-line absorption back-illuminated by the galaxy ("down-the-barrel") and transversely by a background quasar at a projected distance of 58 kpc. Both absorption systems, traced by MgII, are blueshifted relative to the galaxy systemic velocity. The quasar sight-line, which resides almost directly along the projected minor axis of the galaxy, probes MgI and MgII absorption obtained from Keck/LRIS and Lya, SiII and SiIII absorption obtained from HST/COS. For the first time, we combine two independent models used to quantify the outflow properties for down-the-barrel and transverse absorption. We find that the modeled down-the-barrel deprojected outflow velocities range between $V_{dtb}=45-255$ km/s. The transverse bi-conical outflow model, assuming constant-velocity flows perpendicular to the disk, requires wind velocities $V_{outflow}=40-80$ km/s to reproduce the transverse MgII absorption kinematics, which is consistent with the range of $V_{dtb}$. The galaxy has a metallicity, derived from H$\alpha$ and NII, of $[{\rm O/H}]=-0.21\pm0.08$, whereas the transverse absorption has $[{\rm X/H}]=-1.12\pm0.02$. The galaxy star-formation rate is constrained between $4.6-15$ M$_{\odot}$/yr while the estimated outflow rate ranges between $1.6-4.2$ M$_{\odot}$/yr and yields a wind loading factor ranging between $0.1-0.9$. The galaxy and gas metallicities, the galaxy-quasar sight-line geometry, and the down-the-barrel and transverse modeled outflow velocities collectively suggest that the transverse gas originates from ongoing outflowing material from the galaxy. The $\sim$1 dex decrease in metallicity from the base of the outflow to the outer halo suggests metal dilution of the gas by the time it reached 58 kpc.

We analyse the anisotropic clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS Data Release 11 (DR11) sample, which consists of $690\,827$ galaxies in the redshift range $0.43 < z < 0.7$ and has a sky coverage of $8\,498\,\text{deg}^2$. We perform our analysis in Fourier space using a power spectrum estimator suggested by Yamamoto et al. (2006). We measure the multipole power spectra in a self-consistent manner for the first time in the sense that we provide a proper way to treat the survey window function and the integral constraint, without the commonly used assumption of an isotropic power spectrum and without the need to split the survey into sub-regions. The main cosmological signals exploited in our analysis are the Baryon Acoustic Oscillations and the signal of redshift space distortions, both of which are distorted by the Alcock-Paczynski effect. Together, these signals allow us to constrain the distance ratio $D_V(z_{\rm eff})/r_s(z_d) = 13.89\pm 0.18$, the Alcock-Paczynski parameter $F_{\rm AP}(z_{\rm eff}) = 0.679\pm0.031$ and the growth rate of structure $f(z_{\rm eff})\sigma_8(z_{\rm eff}) = 0.419\pm0.044$ at the effective redshift $z_{\rm eff}=0.57$. We did not find significant systematic uncertainties for $D_V/r_s$ or $F_{\rm AP}$ but include a systematic error for $f\sigma_8$ of $3.1\%$. Combining our dataset with Planck to test General Relativity (GR) through the simple $\gamma$-parameterisation, reveals a $\sim 2\sigma$ tension between the data and the prediction by GR. The tension between our result and GR can be traced back to a tension in the clustering amplitude $\sigma_8$ between CMASS and Planck.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $\Lambda$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,\sigma$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.