Published July 22, 2022 | Version Submitted
Discussion Paper Open

Disruption of Dark Matter Minihaloes in the Milky Way environment: Implications for Axion Miniclusters and Early Matter Domination

Abstract

Many theories of dark matter beyond the Weakly Interacting Massive Particles (WIMP) paradigm feature an enhanced matter power spectrum on sub-parsec scales, leading to the formation of dense dark matter minihaloes. While these minihaloes are currently weakly constrained, future local observations, through a variety of techniques, may strongly constrain such substructures. The survival probability of those dense minihaloes in the Milky Way environment is crucial for interpreting local observations. In this work, we investigate two disruption effects: stellar disruption and (smooth) tidal disruption. These two mechanisms are studied using semi-analytic models and idealized N-body simulations. For stellar disruption, we perform a series of N-body simulations of isolated minihalo-star encounters to test and calibrate analytic models of stellar encounters, and apply the model to the realistic Milky Way disk environment. For tidal disruption, we also perform N-body simulations to confirm the effectiveness of the analytic treatment. Finally, we propose a framework to combine the stellar and tidal disruption of minihaloes with an orbit model, using it to make predictions for the overall survival probability of minihaloes in the Milky Way. We find the survival fraction for dense dark matter minihaloes, e.g. for axion miniclusters and minihaloes from Early Matter Domination, is ∼70%, with the relatively low-mass, compact population surviving. The survival fraction is insensitive to the detailed model parameters. We discuss various implications for their mass functions and future detection prospects.

Additional Information

We thank Andrea Mitridate for useful discussions on the mass-concentration relation and Gabriel Aguiar for the collaboration in the early stages of this work. HX is supported in part by the United States Department of Energy (DOE) under grant number DE-SC0011637. Support for XS and PFH was provided by the National Science Foundation (NSF) Research Grants 1911233, 20009234, 2108318, NSF CAREER grant 1455342, the National Aeronautics and Space Administration (NASA) grants 80NSSC18K0562, HST-AR-15800. Numerical calculations were run on the Caltech computing cluster "Wheeler", allocations AST21010 and AST20016 supported by the NSF and the Texas Advanced Computing Center (TACC), and NASA HEC SMD-16-7592. DATA AVAILABILITY. The simulation data of this work was generated and stored on the computing system "Wheeler" at California Institute of Technology. The code for the semi-analytic model and the summary of simulation results are available at project repository. The raw data of the idealized simulations will be shared on reasonable request to the corresponding author.

Attached Files

Submitted - 2207.11276.pdf

Files

2207.11276.pdf

Files (4.3 MB)

Name Size Download all
md5:2ce51dbdec332e3a71dd743ff19b9783
4.3 MB Preview Download

Additional details

Identifiers

Eprint ID
116317
Resolver ID
CaltechAUTHORS:20220816-192431512

Funding

Department of Energy (DOE)
DE-SC0011637
NSF
AST-1911233
NSF
AST-20009234
NSF
AST-2108318
NSF
AST-1455342
NASA
80NSSC18K0562
NASA
HST-AR-15800.001-A
NSF
AST-21010
NSF
AST-20016
Texas Advanced Computing Center (TACC)
NASA
SMD-16-7592

Dates

Created
2022-08-16
Created from EPrint's datestamp field
Updated
2023-06-02
Created from EPrint's last_modified field

Caltech Custom Metadata

Caltech groups
TAPIR, Walter Burke Institute for Theoretical Physics, Astronomy Department