Constraining Ultra-Light Dark Matter mass with Dwarf Galaxy Rotation Curves

Basic Information

• Paper ID: 2510.12848
• Title: Constraining Ultra-Light Dark Matter mass with Dwarf Galaxy Rotation Curves
• Authors: Tian-Yao Fang, Ming-Chung Chu (Chinese University of Hong Kong)
• Classification: astro-ph.GA (Astrophysics - Galaxies)
• Publication Date: October 16, 2025
• Paper Link: https://arxiv.org/abs/2510.12848

Abstract

Although ultra-light bosonic dark matter (ULDM) in a Bose-Einstein condensate (BEC) state can naturally explain the central core structures of certain galaxies and resolve the core-cusp problem, the dark matter density distribution in the outer regions of galaxies remains insufficiently explored. This paper proposes a trial wave function to model the ULDM distribution beyond the BEC core, derives the corresponding rotation velocity curves, and demonstrates excellent agreement with observational data from 12 dwarf spheroidal galaxies. The optimal-fit ULDM particle mass for each dwarf galaxy falls within a remarkably narrow range: $m=(1.8-3.2)\times 10^{-23}\text{ eV}$.

Research Background and Motivation

Core Issues

1. Core-Cusp Problem: While the ΛCDM model successfully describes cosmic structure formation on large scales, it faces persistent challenges at galactic scales (≤10 kpc). Numerical simulations predict dense cusps at galaxy centers, but observations reveal flat core structures.
2. ULDM Theoretical Gaps: Existing ULDM research falls into two categories: (1) numerical solutions of the Schrödinger-Poisson equation with high computational costs; (2) phenomenological approaches relying on empirical relations rather than fundamental wave function derivations.

Research Significance

• The macroscopic de Broglie wavelength of ultra-light dark matter (mass ~10^-22 eV) can reach kpc scales, introducing quantum pressure that suppresses small-scale structures and naturally produces flat core profiles
• Different astrophysical probes provide contradictory constraints on ULDM mass, necessitating a unified theoretical framework

Research Motivation

This work aims to establish an analytical theoretical framework intermediate between computationally precise but expensive numerical methods and phenomenological approaches, constraining ULDM particle mass through dwarf galaxy rotation curves.

Core Contributions

1. Proposes Novel Analytical Trial Wave Function: First analytical model for ULDM in low-density regions beyond the BEC core
2. Derives Rotation Curve Formula: Derives corresponding rotation velocity curves from the trial wave function
3. Observational Data Validation: Achieves excellent fits with rotation curve data from 12 dwarf spheroidal galaxies
4. Constrains ULDM Mass: Obtains highly consistent particle mass range $m=(1.8-3.2)\times 10^{-23}\text{ eV}$
5. Reveals Universality: Discovers that rotation curves of different galaxies exhibit the same universal shape

Methodology Details

Task Definition

Input: Observational rotation velocity data of dwarf galaxies Output: Constraints on ULDM particle mass Objective: Determine optimal ULDM particle mass parameters through rotation curve fitting

Theoretical Framework

Schrödinger-Poisson Equation

Under spherical symmetry assumption, the Hamiltonian of the ULDM system is: $H = \sum_i \left(\frac{\vec{p}_i^2}{2m} - \frac{GM_b m}{|\vec{r}_i|}\right) - \sum_{i<j} \frac{Gm^2}{|\vec{r}_i - \vec{r}_j|}$

The corresponding Schrödinger-Poisson equation: $\left(-\frac{\hbar^2}{2m}\nabla^2 - \frac{GM_b m}{r}\right)\phi(r) - \int_{|\vec{r}'|>r_{min}} d^3r' \frac{Gm^2}{|\vec{r} - \vec{r}'|} \rho(r')\phi(r) = \mu\phi(r)$

where $\rho(\vec{r}) = \phi(\vec{r})^2$ is the particle number density.

Trial Wave Function Design

Outer Region ($r > r_0$): $\phi(r) = n_1 e^{-\eta r}$

Inner Region ($r < r_0$): $\phi(r) = n_2 r^{-1} e^{-f(r)}$

where the key function: $f''= \frac{2a}{r} + \frac{2b}{r+c}$

Unified Form: $\phi(r) = r^{-\beta(r)-1}$

$\beta(r) = \frac{1-\tanh[g(r-r_0)]}{2}[c_1 + c_2 r + 2ar(\ln r - 1) + 2br(\ln(r+c) - 1 + \frac{c}{r}\ln(r+c))] + \frac{\tanh[g(r-r_0)] + 1}{2}(\eta r - \frac{n_1/n_2}{\ln r} - 1)$

Technical Innovations

1. Piecewise Continuous Wave Function: Smoothly connects inner and outer regions via tanh function, ensuring third-derivative continuity
2. Parameterization Scheme: Introduces dimensionless parameter $t = r_0/r_{min}$ to control dark matter halo size and system energy
3. Deviation Minimization: Defines deviation function $\delta = \frac{H(r)\phi(r)}{(-\eta^2/2)\phi(r)} - 1$ and minimizes it
4. Rotation Curve Derivation: Calculates linear velocity directly from total potential energy

Experimental Setup

Dataset

• Dwarf Galaxy Sample: Rotation curve observational data from 12 dwarf spheroidal galaxies
• Large Galaxy Verification: Outer rotation curves of the Milky Way (MW) and Andromeda (M31)
• Data Sources: Observational studies including Karukes & Salucci (2016)

Fitting Method

• Fitting Tool: Python's curve_fit package
• Fitting Range: Primarily in the $r < 10r_{min}$ region
• Free Parameters: ULDM particle mass $m$ as the primary fitting parameter

Model Parameters

Optimal fitting parameters for different $t$ values:

Experimental Results

Main Results

Dwarf Galaxy Fitting

• Fitting Quality: Excellent fits achieved for rotation curves of all 12 dwarf galaxies
• Mass Constraints: ULDM particle mass range of $(1.8-3.2) \times 10^{-23}$ eV
• Consistency: Highly consistent mass estimates across different galaxies

Specific Fitting Results Examples

• DDO168: $m = 0.274 \pm 0.011 \times 10^{-22}$ eV
• UGC1281: $m = 0.264 \pm 0.007 \times 10^{-22}$ eV
• DDO133: $m = 0.319 \pm 0.011 \times 10^{-22}$ eV

Model Verification

Wave Function Deviation Analysis

• In the $R > R_0$ region, deviation $\delta$ rapidly decays to near zero
• Between $R_{min}$ and $R_0$, deviation remains within 5%
• High precision achieved for rotation curves in large $R$ regions

Large Galaxy Verification

• Milky Way: Consistent with observational data within 500 kpc range
• M31: Outer rotation curve shape matches model predictions
• Universality: Rotation curves for different $t$ values exhibit identical shapes after appropriate scaling

Experimental Findings

1. Universal Rotation Curves: All galaxies' outer rotation curves exhibit the same universal shape
2. Mass Consistency: ULDM mass estimates from 12 dwarf galaxies show remarkable agreement
3. Scale Independence: Rotation curve shape independent of system radius and mass

Related Work

Main Research Directions

1. Numerical Methods: Kaup (1968), Guzmán & Ureña-López (2004) and others solving Schrödinger-Poisson equations numerically
2. Phenomenological Constraints: Iršič et al. (2017) using Lyman-α data, Bar et al. (2018) employing multiple astrophysical probes
3. Theoretical Analysis: Hui et al. (2017) comprehensive review, Chavanis (2011) theoretical framework

Advantages of This Work

• Provides analytical framework intermediate between precise numerical methods and phenomenological approaches
• High computational efficiency suitable for large galaxy surveys
• Based on fundamental physics principles rather than empirical relations

Conclusions and Discussion

Main Conclusions

1. Successfully establishes analytical model for ULDM in BEC core outer regions
2. Constrains ULDM particle mass to $m=(1.8-3.2)\times 10^{-23}$ eV through 12 dwarf galaxy rotation curves
3. Discovers universal properties of galaxy rotation curves, supporting ULDM theoretical paradigm

Limitations

1. Simplified Assumptions: Neglects ULDM particle self-interactions and baryonic gravitational potential
2. BEC Boundary: Assumes BEC-outer matter non-interference; boundary treatment is coarse
3. Small Radius Deviations: Trial wave function shows larger deviations at smaller $r$
4. Baryonic Contribution: Does not account for visible matter contribution to rotation curves

Future Directions

1. Improve trial wave function to reduce small-radius deviations
2. Include effects of baryonic matter and self-interactions
3. Extend to larger galaxy samples for statistical analysis
4. Combine with other astrophysical observations to further constrain ULDM properties

In-Depth Evaluation

Strengths

Theoretical Innovation:

• First analytical model for ULDM distribution beyond BEC core
• Clever piecewise wave function design with smooth connection method
• Rigorous derivation from fundamental Schrödinger-Poisson equation

Experimental Verification:

• Excellent fits with observational data from 12 dwarf galaxies
• Consistent mass constraints across different galaxies
• Large galaxy verification enhances model credibility

Methodological Contribution:

• Balances computational precision and efficiency
• Provides practical tools for large-scale galaxy research

Shortcomings

Physical Assumptions:

• Neglecting baryonic contribution may lead to mass underestimation
• Assumption of BEC-outer matter non-interference is oversimplified
• Spherical symmetry assumption limits model applicability

Technical Limitations:

• Larger wave function deviations at small radii
• Parameter fitting exhibits certain model dependence
• Lacks systematic comparison with other ULDM constraint methods

Impact

Academic Contributions:

• Provides new observational support for ULDM theory
• Bridges quantum physics and astrophysics
• Advances precise constraint studies of dark matter candidates

Practical Value:

• Provides analytical tools for future large survey projects
• Extensible to research on other galaxy types
• Helps distinguish between different dark matter models

Applicable Scenarios

1. Dwarf Galaxy Research: Dark matter-dominated low-mass galaxies
2. Large Galaxy Outer Regions: Regions where baryonic contribution is negligible
3. Cosmological Simulations: Provides analytical benchmarks for numerical simulations
4. Multi-Messenger Astronomy: Combines with other observational methods like gravitational waves

References

1. Hui, L., Ostriker, J. P., Tremaine, S., & Witten, E. (2017). Ultralight scalars as cosmological dark matter. Physical Review D, 95(4), 043541.
2. Karukes, E. V., & Salucci, P. (2017). The universal rotation curve of dwarf disk galaxies. Monthly Notices of the Royal Astronomical Society, 465(4), 4703-4722.
3. Iršič, V., et al. (2017). New constraints on the free-streaming of warm dark matter from intermediate and small scale Lyman-α forest data. Physical Review Letters, 119(3), 031302.

Overall Assessment: This is a high-quality theoretical astrophysics paper that makes important contributions to ULDM research. Despite certain limitations in physical assumptions, its theoretical innovation, comprehensive experimental verification, and result consistency demonstrate significant progress in the field. The paper provides a new observational window for understanding dark matter's nature, with important scientific value and practical applications.