Engineering topics

Hearin, Andrew P. (ORCID:0000000322196852)

Publications and source records attributed to Hearin, Andrew P. (ORCID:0000000322196852).

4 records

DOE OSTI2023

A differentiable model of the evolution of dark matter halo concentration

ABSTRACT We introduce a new model of the evolution of the concentration of dark matter haloes, c(t). For individual haloes, our model approximates c(t) as a power law with a time-dependent index, such that at early times, concentration has a nearly constant value of c ≈ 3–4, and as cosmic time progresses, c(t) smoothly increases. Using large samples of halo merger trees taken from the Bolshoi–Planck and MultiDark Planck 2 cosmological simulations, we demonstrate that our three-parameter model can approximate the evolution of the concentration of individual haloes with a typical accuracy of 0.1 dex for $t\gtrsim 2\, {\rm Gyr}$ for all Bolshoi–Planck and MultiDark Planck 2 haloes of present-day peak mass $M_{0}\gtrsim 10^{11.5}\, {\rm M}_{\odot }$. We additionally present a new model of the evolution of the concentration of halo populations, which we show faithfully reproduces both average concentration growth and the diversity of smooth trajectories of c(t), including capturing correlations with halo mass and halo assembly history. Our publicly available source code, diffprof, can be used to generate Monte Carlo realizations of the concentration histories of cosmologically representative halo populations. diffprof is differentiable due to its implementation in the jax autodiff library, which facilitates the incorporation of our model into existing analytical halo model frameworks.

Stevanovich, Dash↗

DOE OSTI2023

Constraints on S 8 from a full-scale and full-shape analysis of redshift-space clustering and galaxy–galaxy lensing in BOSS

ABSTRACT We present a novel simulation-based cosmological analysis of galaxy–galaxy lensing and galaxy redshift-space clustering. Compared to analysis methods based on perturbation theory, our simulation-based approach allows us to probe a much wider range of scales, $0.4 \, h^{-1} \, \mathrm{Mpc}$ to $63 \, h^{-1} \, \mathrm{Mpc}$, including highly non-linear scales, and marginalizes over astrophysical effects such as assembly bias. We apply this framework to data from the Baryon Oscillation Spectroscopic Survey LOWZ sample cross-correlated with state-of-the-art gravitational lensing catalogues from the Kilo Degree Survey and the Dark Energy Survey. We show that gravitational lensing and redshift-space clustering when analysed over a large range of scales place tight constraints on the growth-of-structure parameter $S_8 = \sigma _8 \sqrt{\Omega _{\rm m} / 0.3}$. Overall, we infer S8 = 0.792 ± 0.022 when analysing the combination of galaxy–galaxy lensing and projected galaxy clustering and S8 = 0.771 ± 0.027 for galaxy redshift-space clustering. These findings highlight the potential constraining power of full-scale studies over studies analysing only large scales and also showcase the benefits of analysing multiple large-scale structure surveys jointly. Our inferred values for S8 fall below the value inferred from the CMB, S8 = 0.834 ± 0.016. While this difference is not statistically significant by itself, our results mirror other findings in the literature whereby low-redshift large-scale structure probes infer lower values for S8 than the CMB, the so-called S8-tension.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

DSPS: Differentiable stellar population synthesis

ABSTRACT Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for SPS. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, ii) support for the Diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modelling the colours of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speed-up over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy–halo connection that fully forward model galaxy spectra and photometry.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

The LSST DESC DC2 Simulated Sky Survey

Here, we describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep area of approximately 300 deg 2 , as well as a deep drilling field of approximately 1 deg 2 . We simulate 5 yr of the planned 10 yr survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the data set to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic test bed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time domain cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗