Engineering topics

Slone, Oren

Publications and source records attributed to Slone, Oren.

8 records

DOE OSTI2024

Numerical challenges in modeling gravothermal collapse in Self-Interacting Dark Matter halos

When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature. To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m = 50 cm 2 /g). When simulating these halos, we also vary the mass resolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses. The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10 6 particles per halo), we find that variations of order 10% in collapse time are still present. Here, the results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

A semi-analytic study of self-interacting dark-matter haloes with baryons

ABSTRACT We combine the isothermal Jeans model and the model of adiabatic halo contraction into a semi-analytic procedure for computing the density profile of self-interacting dark-matter (SIDM) haloes with the gravitational influence from the inhabitant galaxies. The model agrees well with cosmological SIDM simulations over the entire core-forming stage up to the onset of gravothermal core-collapse. Using this model, we show that the halo response to baryons is more diverse in SIDM than in CDM and depends sensitively on galaxy size, a desirable feature in the context of the structural diversity of bright dwarfs. The fast speed of the method facilitates analyses that would be challenging for numerical simulations – notably, we quantify the SIDM halo response as functions of the baryonic properties, on a fine mesh grid spanned by the baryon-to-total-mass ratio, Mb/Mvir, and galaxy compactness, r1/2/Rvir; we show with high statistical precision that for typical Milky-Way-like systems, the SIDM profiles are similar to their CDM counterparts; and we delineate the regime of core-collapse in the Mb/Mvir − r1/2/Rvir space, for a given cross section and concentration. Finally, we compare the isothermal Jeans model with the more sophisticated gravothermal fluid model, and show that the former yields faster core formation and agrees better with cosmological simulations. We attribute the difference to whether the target CDM halo is used as a boundary condition or as the initial condition for the gravothermal evolution, and thus comment on possible improvements of the fluid model. We have made our model publicly available at https://github.com/JiangFangzhou/SIDM.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Orbital evolution of satellite galaxies in self-interacting dark matter models

Dark matter self-interactions can leave distinctive signatures on the properties of satellite galaxies around Milky Way–like hosts through their impact on tidal stripping, ram pressure, and gravothermal collapse. We delineate the regions of self-interacting dark matter parameter space—specified by interaction cross section and a velocity scale—where each of these effects dominates and show how the relative mass loss depends on the satellite’s initial mass, density profile, and orbit. We obtain novel, conservative constraints in this parameter space using Milky Way satellite galaxies with notably high central densities and small pericenter distances. Furthermore, our results for self-interacting dark matter models, in combination with constraints from clusters of galaxies, favor either velocity-dependent cross sections that lead to gravothermal core collapse in the densest satellites or small cross sections which more closely resemble cold and collisionless dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

A preference for cold dark matter over Superfluid Dark Matter in local Milky Way data

There are many well-known correlations between dark matter and baryons that exist on galactic scales. These correlations can essentially be encompassed by a simple scaling relation between observed and baryonic accelerations, historically known as the Mass Discrepancy Acceleration Relation (MDAR). The existence of such a relation has prompted many theories that attempt to explain the correlations by invoking additional fundamental forces on baryons. The standard lore has been that a theory that reduces to the MDAR on galaxy scales but behaves like cold dark matter (CDM) on larger scales provides an excellent fit to data, since CDM is desirable on scales of clusters and above. However, this statement should be revised in light of recent results showing that a fundamental force that reproduces the MDAR is challenged by local Milky Way dynamics and rotation curve data between 5–18 kpc. In this study, we test this claim on the example of Superfluid Dark Matter. We find that a standard CDM model is preferred over a static superfluid profile assuming a steady-state Galactic disk and discuss the robustness of this conclusion to disequilibrium effects. This preference is due to the fact that the superfluid model over-predicts vertical accelerations, even while reproducing galactic rotation curves. Our results establish an important criterion that any dark matter model must satisfy within the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Dynamics of stellar disc tilting from satellite mergers

The Milky Way’s stellar disc can tilt in response to torques that result from infalling satellite galaxies and their associated tidal debris. In this work, we explore the dynamics of disc tilting by running N-body simulations of mergers in an isolated, isotropic Milky Way-like host galaxy, varying over satellite virial mass, initial position, and orbit. We develop and validate a first-principles understanding of the dynamics that govern how the host galaxy’s stellar disc responds to the satellite’s dark matter (DM) debris. We find that the degree of disc tilting can be large for cosmologically motivated merger histories. In particular, our results suggest that the Galactic disc may still be tilting in response to Gaia-Sausage-Enceladus, one of the most significant recent mergers in the Milky Way’s history. In conclusion, these findings have implications for terrestrial direct detection experiments as disc tilting changes the relative location of the Sun with respect to DM substructure left behind by a merging galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

A critical assessment of solutions to the galaxy diversity problem

Galactic rotation curves exhibit a diverse range of inner slopes. Observational data indicates that explaining this diversity may require a mechanism that correlates a galaxy's surface brightness with the central-most region of its dark matter halo. Here, in this work, we compare several concrete models that capture the relevant physics required to explain the galaxy diversity problem. We focus specifically on a Self-Interacting Dark Matter (SIDM) model with an isothermal core and two Cold Dark Matter (CDM) models with/without baryonic feedback. In contrast to the CDM case, the SIDM model can lead to the formation of an isothermal core in the halo, and is also mostly insensitive to baryonic feedback processes, which act on longer time-scales. Using rotation curves from 90 galaxies in the Spitzer Photometry & Accurate Rotation Curves (SPARC) catalog, we perform a comprehensive model comparison that addresses issues of statistical methodology from prior works. The best-fit halo models that we recover are consistent with standard CDM concentration-mass and abundance matching relations. We find that both the SIDM and feedback-affected CDM models are better than a CDM model with no feedback in explaining the rotation curves of low and high surface brightness galaxies in the sample. However, when compared to each other, there is no strong statistical preference for either the SIDM or the feedback-affected CDM halo model as the source of galaxy diversity in the SPARC catalog.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

A measurement of stellar surface gravity hidden in radial velocity differences of comoving stars

The gravitational redshift induced by stellar surface gravity is notoriously difficult to measure for non-degenerate stars, since its amplitude is small in comparison with the typical Doppler shift induced by stellar radial velocity. In this study, we make use of the large observational data set of the Gaia mission to achieve a significant reduction of noise caused by these random stellar motions. By measuring the differences in velocities between the components of the pairs of comoving stars and wide binaries, in this work we are able to statistically measure the combined effects of gravitational redshift and convective blueshifting of spectral lines, and nullify the effect of the peculiar motions of the stars. For the subset of stars considered in this study, we find a positive correlation between the observed differences in Gaia radial velocities and the differences in surface gravity and convective blueshift inferred from effective temperature and luminosity measurements. The results rule out a null signal at the 5σ level for our full data set. Additionally, we study the subdominant effects of binary motion, and possible systematic errors in radial velocity measurements within Gaia. Results from the technique presented in this study are expected to improve significantly with data from the next Gaia data release. Such improvements could be used to constrain the mass–luminosity relation and stellar models that predict the magnitude of convective blueshift.

(stars:) binaries: general↗

DOE OSTI2021

Dark Higgs dark matter

A new U(1) dark gauge group coupled to the Standard Model (SM) via the kinetic mixing portal provides a natural dark matter candidate in the form of the Higgs field, hd, responsible for generating the mass of the dark photon, γ d . We show that the condition m h d ≤m γ d , together with smallness of the kinetic mixing parameter, ϵ, and/or dark gauge coupling, g d , leads the dark Higgs to be sufficiently metastable to constitute dark matter. We analyze the Universe's thermal history and show that both freeze-in, SM→{γ d ,h d }, and freeze-out, {γ d ,h d }→SM, processes can lead to viable dark Higgs dark matter with a sub-GeV mass and a kinetic mixing parameter in the range 10 –13 ≲ϵ≲10 –6 . Observable signals in astrophysics and cosmology include modifications to primordial elemental abundances, altered energetics of supernovae explosions, dark Higgs decays in the late Universe, and dark matter self-interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗