We present an algorithm enabling computation of the anisotropic redshift-space galaxy threepoint correlation function (3PCF) scaling as N2, with N the number of galaxies. Our previous work showed how to compute the isotropic 3PCF with this scaling by expanding the radially binned density field around each galaxy in the survey into spherical harmonics and combining these coefficients to form multipole moments. The N2 scaling occurred because this approach never explicitly required the relative angle between a galaxy pair about the primary galaxy. Here, we generalize this work, demonstrating that in the presence of azimuthally symmetric anisotropy produced by redshift-space distortions (RSD), the 3PCF can be described by two triangle side lengths, two independent total angular momenta, and a spin. This basis for the anisotropic 3PCF allows its computation with negligible additional work over the isotropic 3PCF. We also present the covariance matrix of the anisotropic 3PCF measured in this basis. Our algorithm tracks the full 5D redshift-space 3PCF, uses an accurate line of sight to each triplet, is exact in angle, and easily handles edge correction. It will enable use of the anisotropic large-scale 3PCF as a probe of RSD in current and upcoming large-scale redshift surveys.

We present a configuration-space model of the large-scale galaxy 3-point correlation function (3PCF) based on leading-order perturbation theory and including redshift-space distortions (RSD). This model should be useful in extracting distance-scale information from the 3PCF via the baryon acoustic oscillation method. We include the first redshift-space treatment of biasing by the baryon-dark matter relative velocity. Overall, on large scales the effect of RSD is primarily a renormalization of the 3PCF that is roughly independent of both physical scale and triangle opening angle; for our adoptedm and bias values, the rescaling is a factor of ∼1.8. We also present an efficient scheme for computing 3PCF predictions from our model, important for allowing fast exploration of the space of cosmological parameters in future analyses.

A measurement of a primordial non-Gaussianity (PNG) signal through late- or early-Universe probes has the potential to transform our understanding of the physics of the primordial Universe. While large-scale structure observables in principle contain vital information, interpreting these measurements is challenging due to poorly understood astrophysical effects. Luckily, N-body simulations, such as the abacuspng set presented in this study, consisting of 9 boxes, each of size Lbox=2 Gpc/h and particle mass of 1.01×1010M⊙/h, provide a viable path forward. As validation, we find good agreement between the simulations and our expectations from one-loop perturbation theory (PT) and the "separate universe"method for the matter bispectrum, matter power spectrum and the halo bias parameter associated with PNG, bφ. As a science application, we investigate the link between halo assembly bias and bφ for halo properties known to play a vital role in accurately predicting galaxy clustering: concentration, shear (environment), and accretion rate. We find a strong response for all three parameters, suggesting that the connection between bφ and the assembly history of halos needs to be taken into account by future PNG analyses. We further perform the first study of the bφ parameter from fits to early DESI data of the luminous red galaxy (LRG) and quasistellar object (QSO) samples and comment on the effect on fNL constraints for the allowed galaxy-halo models (note that σ[fNL]∝σ[bφ]bφ). We find that the error on fNL is 15, 8, 7 for the LRGs at z=0.5 and z=0.8 and QSOs at z=1.4, respectively, suggesting that a thorough understanding of galaxy assembly bias is warranted so as to perform robust high-precision analysis of local-type PNG with future surveys.

We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical haonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $ {O}(N_math{g}2)$ for Ng particles, or $ {O}(N_math{FFT}log N_math{FFT})$ when using a Fast Fourier Transfo with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical haonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in $100$ CPU-hours, including the corrections necessary for non-unifo survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

We test for galactic conformity at 0.2 < z < 1.0 to a projected distance of 5Mpc using spectroscopic redshifts from the PRism MUlti-object Survey (PRIMUS). Our sample consists of ∼60,000 galaxies in five separate fields covering a total of ∼5.5 square degrees, which allows us to account for cosmic variance. We identify star-forming and quiescent "isolated primary" (i.e., central) galaxies using isolation criteria and cuts in specific star formation rate. We match the redshift and stellar mass distributions of these samples to control for correlations between quiescent fraction and redshift and stellar mass. We detect a significant (σ) one-halo conformity signal, or an excess of star-forming neighbors around star-forming central galaxies, of ∼5% on scales of 0-1 Mpc and a 2.5σ two-halo signal of ∼1% on scales of 1-3 Mpc. These signals are weaker than those detected in the Sloan Digital Sky Survey and are consistent with galactic conformity being the result of large-scale tidal fields and reflecting assembly bias. We also measure the star-forming fraction of central galaxies at fixed stellar mass as a function of large-scale environment and find that central galaxies are more likely to be quenched in overdense environments, independent of stellar mass. However, we find that environment does not affect the star formation efficiency of central galaxies, as long as they are forming stars. We test for redshift and stellar mass dependence of the conformity signal within our sample and show that large volumes and multiple fields are required at intermediate redshift to adequately account for cosmic variance.

The nature of dark energy and the complete theory of gravity are two central questions currently facing cosmology. A vital tool for addressing them is the 3-point correlation function (3PCF), which probes deviations from a spatially random distribution of galaxies. However, the 3PCF's formidable computational expense has prevented its application to astronomical surveys comprising millions to billions of galaxies. We present Galactos, a high-performance implementation of a novel, O(N2) algorithm that uses a load-balanced k-d tree and spherical harmonic expansions to compute the anisotropic 3PCF. Our implementation is optimized for the Intel Xeon Phi architecture, exploiting SIMD parallelism, instruction and thread concurrency, and significant L1 and L2 cache reuse, reaching 39% of peak performance on a single node. Galactos scales to the full Cori system, achieving 9.8 PF (peak) and 5.06 PF (sustained) across 9636 nodes, making the 3PCF easily computable for all galaxies in the observable universe.

We present an observationally motivated model to connect the active galactic nucleus (AGN) and galaxy populations at 0.2 < z < 1.0 and predict the AGN X-ray luminosity function (XLF). We start with measurements of the stellar mass function of galaxies (from the Prism Multi-object Survey) and populate galaxies with AGNs using models for the probability of a galaxy hosting an AGN as a function of specific accretion rate. Our model is based on measurements indicating that the specific accretion rate distribution is a universal function across a wide range of host stellar masses with slope γ1 -0.65 and an overall normalization that evolves with redshift. We test several simple assumptions to extend this model to high specific accretion rates (beyond the measurements) and compare the predictions for the XLF with the observed data. We find good agreement with a model that allows for a break in the specific accretion rate distribution at a point corresponding to the Eddington limit, a steep power-law tail to super-Eddington ratios with slope , and a scatter of 0.38 dex in the scaling between black hole and host stellar mass. Our results show that samples of low luminosity AGNs are dominated by moderately massive galaxies () growing with a wide range of accretion rates due to the shape of the galaxy stellar mass function rather than a preference for AGN activity at a particular stellar mass. Luminous AGNs may be a severely skewed population with elevated black hole masses relative to their host galaxies and in rare phases of rapid accretion. © 2013. The American Astronomical Society. All rights reserved.

We report small-scale clustering measurements from the PRIsm MUlti-object Survey (PRIMUS) spectroscopic redshift survey as a function of color and luminosity. We measure the real-space cross-correlations between 62,106 primary galaxies with PRIMUS redshifts and a tracer population of ∼545,000 photometric galaxies over redshifts from z = 0.2 to z = 1. We separately fit a power-law model in redshift and luminosity to each of three independent color-selected samples of galaxies. We report clustering amplitudes at fiducial values of z = 0.5 and The clustering of the red galaxies is times as strong as that of the blue galaxies and as strong as that of the green galaxies. We also find that the luminosity dependence of the clustering is strongly dependent on physical scale, with greater luminosity dependence being found between and , compared to the to range. Moreover, over a range of two orders of magnitude in luminosity, a single power-law fit to the luminosity dependence is not sufficient to explain the increase in clustering at both the bright and faint ends at the smaller scales. We argue that luminosity-dependent clustering at small scales is a necessary component of galaxy-halo occupation models for blue, star-forming galaxies as well as for red, quenched galaxies.

We utilize ΛCDM halo occupation models of galaxy clustering to investigate the evolving stellar mass dependent clustering of galaxies in the PRIsm MUlti-object Survey (PRIMUS) and DEEP2 Redshift Survey over the past eight billion years of cosmic time, between 0.2 < z < 1.2. These clustering measurements provide new constraints on the connections between dark matter halo properties and galaxy properties in the context of the evolving largescale structure of the universe. Using both an analytic model and a set of mock galaxy catalogs, we find a strong correlation between central galaxy stellar mass and dark matter halo mass over the range Mhalo ∼ 10111013 h-1 M⊙, approximately consistent with previous observations and theoretical predictions. However, the stellarto-halo mass relation and the mass scale where star formation efficiency reaches a maximum appear to evolve more strongly than predicted by other models, including models based primarily on abundance-matching constraints. We find that the fraction of satellite galaxies in halos of a given mass decreases significantly from z ∼ 0.5 to z ∼ 0.9, partly due to the fact that halos at fixed mass are rarer at higher redshift and have lower abundances. We also find that the M1/Mmin ratio, a model parameter that quantifies the critical mass above which halos host at least one satellite, decreases from ≈20 at z ∼ 0 to ≈13 at z ∼ 0.9. Considering the evolution of the subhalo mass function vis-à-vis satellite abundances, this trend has implications for relations between satellite galaxies and halo substructures and for intracluster mass, which we argue has grown due to stripped and disrupted satellites between z ∼ 0.9 and z ∼ 0.5.