Dr Harry Desmond

Abstract:

<jats:p><jats:italic>Context.</jats:italic> Current and future large-scale structure surveys aim to constrain the neutrino mass and the equation of state of dark energy. To do this efficiently, rapid yet accurate evaluation of the matter power spectrum in the presence of these effects is essential.</jats:p> <jats:p><jats:italic>Aims.</jats:italic> We aim to construct accurate and interpretable symbolic approximations of the linear and nonlinear matter power spectra as a function of cosmological parameters in extended ΛCDM models that contain massive neutrinos and nonconstant equations of state for dark energy. This constitutes an extension of the S<jats:sc>YREN-HALOFIT</jats:sc> emulators to incorporate these two effects, which we call S<jats:sc>YREN-NEW</jats:sc> (SYmbolic-Regression-ENhanced power spectrum emulator with NEutrinos and <jats:italic>W</jats:italic><jats:sub>0</jats:sub>−<jats:italic>w</jats:italic><jats:sub><jats:italic>a</jats:italic></jats:sub>). We also wish to obtain a simple approximation of the derived parameter, <jats:italic>σ</jats:italic><jats:sub>8</jats:sub>, as a function of the cosmological parameters for these models.</jats:p> <jats:p><jats:italic>Methods.</jats:italic> We utilizedd symbolic regression to efficiently search through candidate analytic expressions to approximate the various quantities of interest. Our results for the linear power spectrum are designed to emulate C<jats:sc>LASS</jats:sc>, whereas for the nonlinear case we aim to match the results of E<jats:sc>UCLIDEMULATOR</jats:sc>2. We compared our results to existing emulators and <jats:italic>N</jats:italic>-body simulations.</jats:p> <jats:p><jats:italic>Results.</jats:italic> Our analytic emulators for <jats:italic>σ</jats:italic><jats:sub>8</jats:sub>, and the linear and nonlinear power spectra achieve root mean squared errors of 0.1%, 0.3%, and 1.3%, respectively, across a wide range of cosmological parameters, redshifts and wavenumbers. The error on the nonlinear power spectrum is reduced by approximately a factor of 2 when considering observationally plausible dark energy models and neutrino masses. We verify that emulator-related discrepancies are subdominant compared to observational errors and other modeling uncertainties when computing shear power spectra for LSST-like surveys. Our expressions have similar accuracy to existing (numerical) emulators, but are at least an order of magnitude faster, both on a CPU and a GPU.</jats:p> <jats:p><jats:italic>Conclusions.</jats:italic> Our work greatly improves the accuracy, speed, and applicability range of current symbolic approximations of the linear and nonlinear matter power spectra. These now cover the same range of cosmological models as many numerical emulators with similar accuracy, but are much faster and more interpretable. We provide publicly available code for all symbolic approximations found.</jats:p>

Abstract:

Resolved rotation curves (RCs) are the gold-standard measurements for inferring dark matter distributions in Lambda cold dark matter and testing alternative theories of dynamics in galaxies. However, they are expensive to obtain, making them prohibitive for large galaxy samples and at higher redshift. Spatially integrated flux profiles are more accessible and present the information in a different form, but – except in a highly compressed form, as linewidths – have not so far been compared in detail with RCs or employed for dynamical inferences. Here, we study the consistency of RCs and surface density profiles from SPARC with spatially integrated flux profiles from ALFALFA, by combining the resolved properties in a forward model for the flux profile. We define a new metric for asymmetry in the flux profiles, enabling us to cleanly identify those unsuitable for our axisymmetric method. Among all SPARC galaxies the agreement between RCs and flux profiles is satisfactory within the limitations of the data – with most galaxies having an uncertainty-normalized mean squared error (MSE) below 10 – whilst no galaxy deemed symmetric has a MSE above 1.2. Most cases of good agreement prefer an gas dispersion of 13 km s, consistent with resolved studies of gas dispersion from the literature. These results open the door for spatially integrated flux profiles to be used as proxies for spatially resolved dynamics, including a robust appraisal of the degree of asymmetry.