-
Confronting the Diversity Problem: The Limits of Galaxy Rotation Curves as a tool to Understand Dark Matter Profiles
Authors:
Isabel S. Sands,
Philip F. Hopkins,
Xuejian Shen,
Michael Boylan-Kolchin,
James Bullock,
Claude-Andre Faucher-Giguere,
Francisco J. Mercado,
Jorge Moreno,
Lina Necib,
Xiaowei Ou,
Sarah Wellons,
Andrew Wetzel
Abstract:
While galaxy rotation curves provide one of the most powerful methods for measuring dark matter profiles in the inner regions of rotation-supported galaxies, at the dwarf scale there are factors that can complicate this analysis. Given the expectation of a universal profile in dark matter-only simulations, the diversity of observed rotation curves has become an often-discussed issue in Lambda Cold…
▽ More
While galaxy rotation curves provide one of the most powerful methods for measuring dark matter profiles in the inner regions of rotation-supported galaxies, at the dwarf scale there are factors that can complicate this analysis. Given the expectation of a universal profile in dark matter-only simulations, the diversity of observed rotation curves has become an often-discussed issue in Lambda Cold Dark Matter cosmology on galactic scales. We analyze a suite of Feedback in Realistic Environments (FIRE) simulations of $10^{10}-10^{12}$ $M_\odot$ halos with standard cold dark matter, and compare the true circular velocity to rotation curve reconstructions. We find that, for galaxies with well-ordered gaseous disks, the measured rotation curve may deviate from true circular velocity by at most 10% within the radius of the disk. However, non-equilibrium behavior, non-circular motions, and non-thermal and non-kinetic stresses may cause much larger discrepancies of 50% or more. Most rotation curve reconstructions underestimate the true circular velocity, while some reconstructions transiently over-estimate it in the central few kiloparsecs due to dynamical phenomena. We further demonstrate that the features that contribute to these failures are not always visibly obvious in HI observations. If such dwarf galaxies are included in galaxy catalogs, they may give rise to the appearance of "artificial" rotation curve diversity that does not reflect the true variation in underlying dark matter profiles.
△ Less
Submitted 24 April, 2024;
originally announced April 2024.
-
Wide-Angle Effects in the Power Spectrum Multipoles in Next-Generation Redshift Surveys
Authors:
Joshua N. Benabou,
Isabel Sands,
Henry S. Grasshorn Gebhardt,
Chen Heinrich,
Olivier Doré
Abstract:
As galaxy redshift surveys expand to larger areas on the sky, effects coming from the curved nature of the sky become important, introducing wide-angle (WA) corrections to the power spectrum multipoles at large galaxy-pair separations. These corrections particularly impact the measurement of physical effects that are predominantly detected on large scales, such the local primordial non-Gaussianiti…
▽ More
As galaxy redshift surveys expand to larger areas on the sky, effects coming from the curved nature of the sky become important, introducing wide-angle (WA) corrections to the power spectrum multipoles at large galaxy-pair separations. These corrections particularly impact the measurement of physical effects that are predominantly detected on large scales, such the local primordial non-Gaussianities. In this paper, we examine the validity of the perturbative approach to modeling WA effects in the power spectrum multipoles for upcoming surveys by comparing to measurements on simulated galaxy catalogs using the Yamamoto estimator. We show that on the scales $k \lesssim 2π/χ$, where $χ$ is the comoving distance to the galaxies, the estimated power spectrum monopole differs by up to $5\%$ from the second-order perturbative result, with similar absolute deviations for higher multipoles. To enable precision comparison, we pioneer an improved treatment of the $μ$-leakage effects in the Yamamoto estimator. Additionally, we devise a solution to include $f_{\rm NL}$ in the perturbative WA calculations, avoiding divergences in the original framework through the integral constraint. This allows us to conclude that WA effects can mimic a $f_{\rm NL}\sim5$ signal in the lowest SPHEREx redshift bin. We recommend using non-perturbative methods to model large scale power spectrum multipoles for $f_{\rm NL}$ measurements. A companion paper, Wen et al. 2024, addresses this by introducing a new non-perturbative method going through the spherical Fourier-Bessel basis.
△ Less
Submitted 7 April, 2024;
originally announced April 2024.
-
Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations
Authors:
Philip F. Hopkins,
Ethan O. Nadler,
Michael Y. Grudic,
Xuejian Shen,
Isabel Sands,
Fangzhou Jiang
Abstract:
Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve adaptively, in any high-dynamic range simulation, but this poses unique challenges of consistency, conservation, and accuracy, especially in multi-physics simula…
▽ More
Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve adaptively, in any high-dynamic range simulation, but this poses unique challenges of consistency, conservation, and accuracy, especially in multi-physics simulations where species with different softening laws may interact. We therefore derive a generalized form of the energy-and-momentum conserving gravitational equations of motion, applicable to arbitrary rules used to determine the force softening, together with consistent associated timestep criteria, interaction terms between species with different softening laws, and arbitrary maximum/minimum softenings. We also derive new methods to maintain better accuracy and conservation when symmetrizing forces between particles. We review and extend previously-discussed adaptive softening schemes based on the local neighbor particle density, and present several new schemes for scaling the softening with properties of the gravitational field, i.e. the potential or acceleration or tidal tensor. We show that the tidal softening scheme not only represents a physically-motivated, translation and Galilean invariant and equivalence-principle respecting (and therefore conservative) method, but imposes negligible timestep or other computational penalties, ensures that pairwise two-body scattering is small compared to smooth background forces, and can resolve outstanding challenges in properly capturing tidal disruption of substructures (minimizing artificial destruction) while also avoiding excessive N-body heating. We make all of this public in the GIZMO code.
△ Less
Submitted 28 August, 2023; v1 submitted 13 December, 2022;
originally announced December 2022.
-
Schrödinger-Poisson Solitons: Perturbation Theory
Authors:
J. Luna Zagorac,
Isabel Sands,
Nikhil Padmanabhan,
Richard Easther
Abstract:
Self-gravitating quantum matter may exist in a wide range of cosmological and astrophysical settings from the very early universe through to present-day boson stars. Such quantum matter arises in a number of different theories, including the Peccei-Quinn axion and UltraLight (ULDM) or Fuzzy (FDM) dark matter scenarios. We consider the dynamical evolution of perturbations to the spherically symmetr…
▽ More
Self-gravitating quantum matter may exist in a wide range of cosmological and astrophysical settings from the very early universe through to present-day boson stars. Such quantum matter arises in a number of different theories, including the Peccei-Quinn axion and UltraLight (ULDM) or Fuzzy (FDM) dark matter scenarios. We consider the dynamical evolution of perturbations to the spherically symmetric soliton, the ground state solution to the Schrödinger-Poisson system common to all these scenarios. We construct the eigenstates of the Schrödinger equation, holding the gravitational potential fixed to its ground state value. We see that the eigenstates qualitatively capture the properties seen in full ULDM simulations, including the soliton "breathing" mode, the random walk of the soliton center, and quadrupolar distortions of the soliton. We then show that the time-evolution of the gravitational potential and its impact on the perturbations can be well described within the framework of time-dependent perturbation theory. Applying our formalism to a synthetic ULDM halo reveals considerable mixing of eigenstates, even though the overall density profile is relatively stable. Our results provide a new analytic approach to understanding the evolution of these systems as well as possibilities for faster approximate simulations.
△ Less
Submitted 15 April, 2022; v1 submitted 4 September, 2021;
originally announced September 2021.