We use self-similarity in N-body simulations of scale-free models to test for resolution dependence in the mass function and two-point correlation functions of dark matter halos. We use 1024$^3$ particle simulations performed with the ABACUS N-body code and define halos with two different algorithms, "friends of friends" (FOF) and ROCKSTAR. The FOF mass functions show a systematic deviation from self-similarity which is explained by a resolution dependence of the FOF mass assignment previously reported in the literature. Evidence for convergence is observed only starting from halos of several thousand particles, and mass functions are overestimated by as much as 20-25 percent for halos of 50 particles. The mass function of ROCKSTAR halos, on the other hand, shows good convergence from of order 50 to 100 particles per halo, with no detectable evidence at the one percent level of any systematic dependence for larger particle number. Tests show that the mass unbinding procedure in ROCKSTAR is the key factor in obtaining this much improved resolution. Applying the same analysis to the halo-halo two point correlation function, we find again strong evidence for convergence only for ROCKSTAR halos, at separations sufficiently large so that halos do not overlap. At these separations we can exclude dependence on resolution at the few percent level once halos have of order 50 to 100 particles. At smaller separations results are not converged even at significantly larger particle number, and bigger simulations would be required to establish the resolution required for convergence.

We demonstrate that testing for self-similarity in scale-free simulations provides an excellent tool to quantify the resolution at small scales of cosmological N-body simulations. Analysing two-point correlation functions measured in simulations using ABACUS, we show how observed deviations from self-similarity reveal the range of time and distance scales in which convergence is obtained. While the well-converged scales show accuracy below 1 percent, our results show that, with a small force softening length, the spatial resolution is essentially determined by the mass resolution. At later times the lower cut-off scale on convergence evolves in comoving units as $a^{-1/2}$ ($a$ being the scale factor), consistent with a hypothesis that it is set by two-body collisionality. A corollary of our results is that N-body simulations, particularly at high red-shift, contain a significant spatial range in which clustering appears converged with respect to the time-stepping and force softening but has not actually converged to the physical continuum result. The method developed can be applied to determine the resolution of any clustering statistic and extended to infer resolution limits for non-scale-free simulations.

We determined the rotational periods of 3122 Florence, 3830 Trelleborg, and (131077) 2000 YH105 with the Harvard Clay Telescope and KeplerCam at the Fred L. Whipple Observatory. We found the rotational periods to be 2.3580 $\pm$ 0.0015 h, 17.059 $\pm$ 0.017 h, and 1.813 $\pm$ 0.00003 h, respectively. Our measurement of 3122 Florence's period agrees with Warner (2016), who reported 2.3580 $\pm$ 0.0002 h.

We describe our nonlinear emulation (i.e., interpolation) framework that combines the halo occupation distribution (HOD) galaxy bias model with $N$-body simulations of nonlinear structure formation, designed to accurately predict the projected clustering and galaxy-galaxy lensing signals from luminous red galaxies (LRGs) in the redshift range $0.16 < z < 0.36$ on comoving scales $0.6 < r_p < 30$ \hMpc. The interpolation accuracy is $\lesssim 1-2$ per cent across the entire physically plausible range of parameters for all scales considered. We correctly recover the true value of the cosmological parameter $S_8 = ({\sigma_8}/{0.8228}) ({\Omega_{\text{m}}}/{0.3107})^{0.6}$ from mock measurements produced via subhalo abundance matching (SHAM)-based lightcones designed to approximately match the properties of the SDSS LOWZ galaxy sample. Applying our model to Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 14 (DR14) LOWZ galaxy clustering and galaxy-shear cross-correlation measurements made with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) imaging, we perform a prototype cosmological analysis marginalizing over $w$CDM cosmological parameters and galaxy HOD parameters. We obtain a 4.4 per cent measurement of $S_8 = 0.847 \pm 0.037$, in $3.5\sigma$ tension with the Planck cosmological results of $1.00 \pm 0.02$. We discuss the possibility of underestimated systematic uncertainties or astrophysical effects that could explain this discrepancy.

Correlation functions are widely used in extra-galactic astrophysics to extract insights into how galaxies occupy dark matter halos and in cosmology to place stringent constraints on cosmological parameters. A correlation function fundamentally requires computing pair-wise separations between two sets of points and then computing a histogram of the separations. Corrfunc is an existing open-source, high-performance software package for efficiently computing a multitude of correlation functions. In this paper, we will discuss the SIMD AVX512F kernels within Corrfunc, capable of processing 16 floats or 8 doubles at a time. The latest manually implemented Corrfunc AVX512F kernels show a speedup of up to $\sim 4\times$ relative to compiler-generated code for double-precision calculations. The AVX512F kernels show $\sim 1.6\times$ speedup relative to the AVX kernels and compare favorably to a theoretical maximum of $2\times$. In addition, by pruning pairs with too large of a minimum possible separation, we achieve a $\sim 5-10\%$ speedup across all the SIMD kernels. Such speedups highlight the importance of programming explicitly with SIMD vector intrinsics for complex calculations that can not be efficiently vectorized by compilers. Corrfunc is publicly available at https://github.com/manodeep/Corrfunc/.

The two-point correlation function (2PCF) is the most widely used tool for quantifying the spatial distribution of galaxies. Since the distribution of galaxies is determined by galaxy formation physics as well as the underlying cosmology, fitting an observed correlation function yields valuable insights into both. The calculation for a 2PCF involves computing pair-wise separations and consequently, the computing time scales quadratically with the number of galaxies. The next-generation galaxy surveys are slated to observe many millions of galaxies, and computing the 2PCF for such surveys would be prohibitively time-consuming. Additionally, modern modelling techniques require the 2PCF to be calculated thousands of times on simulated galaxy catalogues of {\em at least} equal size to the data and would be completely unfeasible for the next generation surveys. Thus, calculating the 2PCF forms a substantial bottleneck in improving our understanding of the fundamental physics of the universe, and we need high-performance software to compute the correlation function. In this paper, we present Corrfunc --- a suite of highly optimised, OpenMP parallel clustering codes. The improved performance of Corrfunc arises from both efficient algorithms as well as software design that suits the underlying hardware of modern CPUs. Corrfunc can compute a wide range of 2-D and 3-D correlation functions in either simulation (Cartesian) space or on-sky coordinates. Corrfunc runs efficiently in both single- and multi-threaded modes and can compute a typical 2-point projected correlation function ($w_p(r_p)$) for ~1 million galaxies within a few seconds on a single thread. Corrfunc is designed to be both user-friendly and fast and is publicly available at https://github.com/manodeep/Corrfunc.

We present a deep machine learning (ML)-based technique for accurately determining $\sigma_8$ and $\Omega_m$ from mock 3D galaxy surveys. The mock surveys are built from the AbacusCosmos suite of $N$-body simulations, which comprises 40 cosmological volume simulations spanning a range of cosmological models, and we account for uncertainties in galaxy formation scenarios through the use of generalized halo occupation distributions (HODs). We explore a trio of ML models: a 3D convolutional neural network (CNN), a power-spectrum-based fully connected network, and a hybrid approach that merges the two to combine physically motivated summary statistics with flexible CNNs. We describe best practices for training a deep model on a suite of matched-phase simulations and we test our model on a completely independent sample that uses previously unseen initial conditions, cosmological parameters, and HOD parameters. Despite the fact that the mock observations are quite small ($\sim0.07h^{-3}\,\mathrm{Gpc}^3$) and the training data span a large parameter space (6 cosmological and 6 HOD parameters), the CNN and hybrid CNN can constrain $\sigma_8$ and $\Omega_m$ to $\sim3\%$ and $\sim4\%$, respectively.

We present the discovery of KELT-24 b, a massive hot Jupiter orbiting a bright (V=8.3 mag, K=7.2 mag) young F-star with a period of 5.6 days. The host star, KELT-24 (HD 93148), has a $T_{\rm eff}$ =$6509^{+50}_{-49}$ K, a mass of $M_{*}$ = $1.460^{+0.055}_{-0.059}$ $M_{\odot}$, radius of $R_{*}$ = $1.506\pm0.022$ $R_{\odot}$, and an age of $0.78^{+0.61}_{-0.42}$ Gyr. Its planetary companion (KELT-24 b) has a radius of $R_{\rm P}$ = $1.272\pm0.021$ $R_{\rm J}$, a mass of $M_{\rm P}$ = $5.18^{+0.21}_{-0.22}$ $M_{\rm J}$, and from Doppler tomographic observations, we find that the planet's orbit is well-aligned to its host star's projected spin axis ($\lambda$ = $2.6^{+5.1}_{-3.6}$). The young age estimated for KELT-24 suggests that it only recently started to evolve from the zero-age main sequence. KELT-24 is the brightest star known to host a transiting giant planet with a period between 5 and 10 days. Although the circularization timescale is much longer than the age of the system, we do not detect a large eccentricity or significant misalignment that is expected from dynamical migration. The brightness of its host star and its moderate surface gravity make KELT-24b an intriguing target for detailed atmospheric characterization through spectroscopic emission measurements since it would bridge the current literature results that have primarily focused on lower mass hot Jupiters and a few brown dwarfs.

We present a high-fidelity realization of the cosmological $N$-body simulation from the Schneider et al. (2016) code comparison project. The simulation was performed with our Abacus $N$-body code, which offers high force accuracy, high performance, and minimal particle integration errors. The simulation consists of $2048^3$ particles in a $500\ h^{-1}\mathrm{Mpc}$ box, for a particle mass of $1.2\times 10^9\ h^{-1}\mathrm{M}_\odot$ with $10\ h^{-1}\mathrm{kpc}$ spline softening. Abacus executed 1052 global time steps to $z=0$ in 107 hours on one dual-Xeon, dual-GPU node, for a mean rate of 23 million particles per second per step. We find Abacus is in good agreement with Ramses and Pkdgrav3 and less so with Gadget3. We validate our choice of time step by halving the step size and find sub-percent differences in the power spectrum and 2PCF at nearly all measured scales, with $<0.3\%$ errors at $k<10\ \mathrm{Mpc}^{-1}h$. On large scales, Abacus reproduces linear theory better than $0.01\%$. Simulation snapshots are available at http://nbody.rc.fas.harvard.edu/public/S2016 .

We present a method for generating suites of dark-matter halo catalogs with only a few $N$-body simulations, focusing on making small changes to the underlying cosmology of a simulation with high precision. In the context of blind challenges, this allows us to reuse a simulation by giving it a new cosmology after the original cosmology is revealed. Starting with full $N$-body realizations of an original cosmology and a target cosmology, we fit a transfer function that displaces halos in the original so that the galaxy/HOD power spectrum matches that of the target cosmology. This measured transfer function can then be applied to a new realization of the original cosmology to create a new realization of the target cosmology. For a 1% change in $\sigma_8$, we achieve 0.1% accuracy to $k = 1h\,\mathrm{Mpc}^{-1}$ in the real-space power spectrum; this degrades to 0.3% when the transfer function is applied to a new realization. We achieve similar accuracy in the redshift-space monopole and quadrupole. In all cases, the result is better than the sample variance of our $1.1h^{-1}\,\mathrm{Gpc}$ simulation boxes.

We present the GeneRalized ANd Differentiable Halo Occupation Distribution (GRAND-HOD) routine that generalizes the standard 5 parameter halo occupation distribution model (HOD) with various halo-scale physics and assembly bias. We describe the methodology of 4 different generalizations: satellite distribution generalization, velocity bias, closest approach distance generalization, and assembly bias. We showcase the signatures of these generalizations in the 2-point correlation function (2PCF) and the squeezed 3-point correlation function (squeezed 3PCF). We identify generalized HOD prescriptions that are nearly degenerate in the projected 2PCF and demonstrate that these degeneracies are broken in the redshift-space anisotropic 2PCF and the squeezed 3PCF. We also discuss the possibility of identifying degeneracies in the anisotropic 2PCF and further demonstrate the extra constraining power of the squeezed 3PCF on galaxy-halo connection models. We find that within our current HOD framework, the anisotropic 2PCF can predict the squeezed 3PCF better than its statistical error. This implies that a discordant squeezed 3PCF measurement could falsify the particular HOD model space. Alternatively, it is possible that further generalizations of the HOD model would open opportunities for the squeezed 3PCF to provide novel parameter measurements. The GRAND-HOD Python package is publicly available at https://github.com/SandyYuan/GRAND-HOD .

We present a public data release of halo catalogs from a suite of 125 cosmological $N$-body simulations from the Abacus project. The simulations span 40 $w$CDM cosmologies centered on the Planck 2015 cosmology at two mass resolutions, $4\times 10^{10}\;h^{-1}M_\odot$ and $1\times 10^{10}\;h^{-1}M_\odot$, in $1.1\;h^{-1}\mathrm{Gpc}$ and $720\;h^{-1}\mathrm{Mpc}$ boxes, respectively. The boxes are phase-matched to suppress sample variance and isolate cosmology dependence. Additional volume is available via 16 boxes of fixed cosmology and varied phase; a few boxes of single-parameter excursions from Planck 2015 are also provided. Catalogs spanning $z=1.5$ to $0.1$ are available for friends-of-friends and Rockstar halo finders and include particle subsamples. All data products are available at https://lgarrison.github.io/AbacusCosmos

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

In preparation for deep extragalactic imaging with the James Webb Space Telescope, we explore the clustering of massive halos at $z=8$ and $10$ using a large N-body simulation. We find that halos with masses $10^9$ to $10^{11}$ $h^{-1}\;M_\odot$, which are those expected to host galaxies detectable with JWST, are highly clustered with bias factors ranging from 5 and 30 depending strongly on mass, as well as on redshift and scale. This results in correlation lengths of 5--10$h^{-1}\;{\rm Mpc}$, similar to that of today's galaxies. Our results are based on a simulation of 130 billion particles in a box of $250h^{-1}\;{\rm Mpc}$ size using our new high-accuracy ABACUS simulation code, the corrections to cosmological initial conditions of (Garrison et al. 2016, 2016MNRAS.461.4125G), and the Planck 2015 cosmology. We use variations between sub-volumes to estimate the detectability of the clustering. Because of the very strong inter-halo clustering, we find that surveys of order 25$h^{-1}\;{\rm Mpc}$ comoving transverse size may be able to detect the clustering of $z=8$--$10$ galaxies with only 500-1000 survey objects if the galaxies indeed occupy the most massive dark matter halos.

The combination of galaxy-galaxy lensing (GGL) with galaxy clustering is one of the most promising routes to determining the amplitude of matter clustering at low redshifts. We show that extending clustering+GGL analyses from the linear regime down to $\sim 0.5 \, h^{-1}$ Mpc scales increases their constraining power considerably, even after marginalizing over a flexible model of non-linear galaxy bias. Using a grid of cosmological N-body simulations, we construct a Taylor-expansion emulator that predicts the galaxy autocorrelation $\xi_{\text{gg}}(r)$ and galaxy-matter cross-correlation $\xi_{\text{gm}}(r)$ as a function of $\sigma_8$, $\Omega_m$, and halo occupation distribution (HOD) parameters, which are allowed to vary with large scale environment to represent possible effects of galaxy assembly bias. We present forecasts for a fiducial case that corresponds to BOSS LOWZ galaxy clustering and SDSS-depth weak lensing (effective source density $\sim 0.3$ arcmin$^{-2}$). Using tangential shear and projected correlation function measurements over $0.5 \leq r_p \leq 30 \, h^{-1}$ Mpc yields a 1.8% constraint on the parameter combination $\sigma_8\Omega_m^{0.58}$, a factor of two better than a constraint that excludes non-linear scales ($r_p > 2 \, h^{-1}$ Mpc, $4 \, h^{-1}$ Mpc for $\gamma_t,w_p$). Much of this improvement comes from the non-linear clustering information, which breaks degeneracies among HOD parameters that would otherwise degrade the inference of matter clustering from GGL. Increasing the effective source density to $3$ arcmin$^{-2}$ sharpens the constraint on $\sigma_8\Omega_m^{0.58}$ by a further factor of two. With robust modeling into the non-linear regime, low-redshift measurements of matter clustering at the 1-percent level with clustering+GGL alone are well within reach of current data sets such as those provided by the Dark Energy Survey.

We investigate the three-point correlation function (3PCF) in the squeezed limit by considering galaxy pairs as discrete objects and cross-correlating them with the galaxy field. We develop an efficient algorithm using Fast Fourier Transforms to compute such cross-correlations and their associated pair-galaxy bias bpg and the squeezed 3PCF coefficient Qeff. We implement our method using N-body cosmological simulations and a fiducial Halo Occupation Distribution (HOD) and present the results in both the real space and redshift space. In real space, we observe a peak in bpg and Qeff at pair separation of ~ 2 Mpc, attributed to the fact that galaxy pairs at 2 Mpc separation trace the most massive dark matter halos. We also see strong anisotropy in the bpg and Qeff signals that track the large-scale filamentary structure. In redshift space, both the 2 Mpc peak and the anisotropy are significantly smeared out along the line-of-sight due to Finger-of-God effect. In both the real space and redshift space, the squeezed 3PCF shows a factor of 2 variation, contradicting the hierarchical ansatz but offering rich information on the galaxy-halo connection. Thus, we explore the possibility of using the squeezed 3PCF to constrain the HOD. When we compare two simple HOD models that are closely matched in their projected two-point correlation function (2PCF), we do not yet see a strong variation in the 3PCF that is clearly disentangled from variations in the projected 2PCF. Nevertheless, we propose that more complicated HOD models, e.g. those incorporating assembly bias, can break degeneracies in the 2PCF and show a distinguishable squeezed 3PCF signal.

In cosmological $N$-body simulations, the representation of dark matter as discrete "macroparticles" suppresses the growth of structure, such that simulations no longer reproduce linear theory on small scales near $k_{\rm Nyquist}$. Marcos et al. demonstrate that this is due to sparse sampling of modes near $k_{\rm Nyquist}$ and that the often-assumed continuum growing modes are not proper growing modes of the particle system. We develop initial conditions that respect the particle linear theory growing modes and then rescale the mode amplitudes to account for growth suppression. These ICs also allow us to take advantage of our very accurate $N$-body code Abacus to implement 2LPT in configuration space. The combination of 2LPT and rescaling improves the accuracy of the late-time power spectra, halo mass functions, and halo clustering. In particular, we achieve 1% accuracy in the power spectrum down to $k_{\rm Nyquist}$, versus $k_{\rm Nyquist}/4$ without rescaling or $k_{\rm Nyquist}/13$ without 2LPT, relative to an oversampled reference simulation. We anticipate that our 2LPT will be useful for large simulations where FFTs are expensive and that rescaling will be useful for suites of medium-resolution simulations used in cosmic emulators and galaxy survey mock catalogs. Code to generate initial conditions is available at https://github.com/lgarrison/zeldovich-PLT