Effects of New Forces on Scalar Dark Matter Solitons

Basic Information

• Paper ID: 2511.15916
• Title: Effects of New Forces on Scalar Dark Matter Solitons
• Authors: Alize Sucsuzer (Tufts University), Mark P. Hertzberg (Tufts University), Michiru Uwabo-Niibo (IBS, Korea)
• Classification: astro-ph.CO, astro-ph.GA, gr-qc, hep-ph, hep-th
• Submission Date: November 19, 2025
• Paper Link: https://arxiv.org/abs/2511.15916

Abstract

This paper investigates the effects of new long-range forces on scalar dark matter solitons (boson stars). The authors introduce a light but non-zero mass mediator field that transmits new interactions between dark matter particles beyond gravity. The choice of mediator mass allows it to affect small-scale galactic behavior without impacting current cosmological observations. Through numerical simulations, the authors determine the density profiles of solitons in the presence of new forces and extend the analysis to multi-mediator scenarios. The results show that new forces modify the relationship between galactic core density and core radius, improving agreement with observational data to some extent, though improvements are limited for couplings of gravitational strength magnitude.

Research Background and Motivation

Core Problems

1. Galactic Core Problem: Observational data reveal a relationship between galactic core density ρc and core radius Rc given by ρc ∝ 1/R^q_c, where q ≈ 1.3. However, ultralight axion dark matter interacting purely gravitationally predicts q = 4, which severely contradicts observations.
2. Possibility of Dark Matter Self-Interactions: While new long-range forces on ordinary matter are tightly constrained, observational constraints on dark matter are relatively weak, providing space to explore non-gravitational interactions between dark matter particles.

Research Significance

• Exploring Dark Matter Nature: The nature of dark matter remains a fundamental question in modern physics; exploring possible self-interactions is crucial for understanding dark matter properties
• Small-Scale Structure Problems: The standard cold dark matter model faces potential issues at galactic scales; new interaction mechanisms may provide solutions
• Theoretical Extension: Offers new phenomenological possibilities for dark matter theory

Limitations of Existing Approaches

• Ultralight axion models considering only gravitational interactions cannot explain the observed core density-radius relationship
• Even with local scalar potentials V(φ), the problem persists
• New physical mechanisms are needed to improve consistency between theoretical predictions and observations

Core Contributions

1. Theoretical Framework Construction: First systematic study of scalar dark matter soliton theory with new force mediators, establishing a complete non-relativistic weak-field framework
2. Numerical Solution Methods: Developed iterative numerical methods successfully solving soliton configuration equations with new forces
3. Density-Radius Relation Correction: Demonstrated that new forces modify the core density-radius relationship, with different effective gravitational constants in regions Rc ≪ m^(-1)_χ and Rc ≫ m^(-1)_χ
4. Multi-Mediator Extension: Extended analysis to multiple force mediators, demonstrating richer physical behavior
5. Parameter Space Exploration: Systematically investigated effects of dimensionless coupling constant gχ and dimensionless mediator mass β' on soliton properties

Methodology Details

Task Definition

Inputs:

• Dark matter field φ mass: mφ
• Mediator field χ mass: mχ
• Coupling constants: c1, c2 (or combination C = c1 + c2/2)

Outputs:

• Density profile f(r) of static spherically symmetric soliton solution
• Newtonian potential φN(r)
• Mediator field χ(r)
• Relationship between core density ρc and core radius Rc

Constraints:

• Weak-field approximation: φN ≪ 1, ψ ≪ 1
• Non-relativistic limit: characteristic wavelength ≫ Compton wavelength
• Small-field approximation: potential dominated by mass term

Theoretical Architecture

1. Action Construction

The complete action is: $S = \int d^4x \sqrt{-g} [L_\phi + L_\chi + L_{int} + L_{grav}]$

where each term is:

• Dark Matter Field: $L_\phi = \frac{1}{2}g^{\mu\nu}\partial_\mu\phi\partial_\nu\phi - V(\phi)$
• Mediator Field: $L_\chi = \frac{1}{2}g^{\mu\nu}\partial_\mu\chi\partial_\nu\chi - \frac{1}{2}m^2_\chi\chi^2$
• Interaction: $L_{int} = -\chi[c_1g^{\mu\nu}\partial_\mu\phi\partial_\nu\phi + c_2\tilde{V}(\phi)]$
• Gravity: $L_{grav} = -\frac{R}{16\pi G}$

2. Non-Relativistic Reduction

Through field decomposition: $\phi(t,\vec{x}) = \frac{1}{\sqrt{2m_\phi}}[e^{-im_\phi t}F(t,\vec{x}) + e^{im_\phi t}F^*(t,\vec{x})]$

After time averaging, with the special choice c1 = c2/2 (making χ couple primarily to energy density), the simplified Lagrangian is: $L = \frac{i}{2}(\dot{F}F^* - F\dot{F}^*) - \frac{1}{2m_\phi}|\nabla F|^2 - \phi_N\rho_\phi - \frac{1}{8\pi G}(\nabla\phi_N)^2 - \frac{1}{2}(\nabla\chi)^2 - \frac{1}{2}m^2_\chi\chi^2 - C\chi\rho_\phi$

3. Static Spherically Symmetric Solitons

For condensate solutions $F(t,\vec{x}) = e^{-i\mu t}f(\vec{x})$, the equation system is:

Dark Matter Field Equation: $\mu f = -\frac{1}{2m_\phi}\left(f'' + \frac{2}{r}f'\right) + m_\phi(\phi_N + C\chi)f$

Poisson Equation: $\nabla^2\phi_N = 4\pi G\rho_\phi, \quad \rho_\phi = m_\phi f^2$

Screened Poisson Equation: $\nabla^2\chi - m^2_\chi\chi = C\rho_\phi$

4. Dimensionless Formulation

Introducing dimensionless variables:

• $f(r) = \alpha\tilde{f}(r)$, where $\alpha = f(0)$
• $r = \beta\tilde{r}$, where $\beta = (Gm^3_\phi\alpha^2)^{-1/4}$

The system is controlled by two dimensionless parameters:

• Dimensionless Coupling: $g_\chi \equiv \frac{C^2}{4\pi G}$ (strength of new force relative to gravity)
• Dimensionless Mass: $\beta' \equiv m_\chi\beta$ (soliton scale relative to mediator Compton wavelength)

Technical Innovations

1. Special Coupling Choice: By selecting c1 = c2/2, the mediator field χ couples primarily to energy density rather than pressure, avoiding significant χ waves and simplifying analysis
2. Effective Gravitational Constant Concept:
   • $\beta' \ll 1$ (long-range force): $G_{eff} = (1+g_\chi)G$
   • $\beta' \gg 1$ (short-range force): $G_{eff} = G$
   • Intermediate region requires numerical solution
3. Iterative Solution Algorithm:
   • Initially set χ̃ = 0
   • Solve {f̃, φ̃N} system, adjusting chemical potential μ̃ via bisection to ensure solution vanishes at large radius
   • Substitute f̃ into integral representation to obtain χ̃
   • Iterate until convergence
4. Core Radius Definition: Using full width at half maximum definition $f(R_c)^2 = \frac{1}{2}f(0)^2$, obtaining: $\rho_c = \frac{\zeta}{G_{eff}m^2_\phi R^4_c}, \quad \zeta \approx 0.226$

Experimental Setup

Numerical Implementation Details

Boundary Conditions:

• $\tilde{f}(0) = 1$, $\tilde{f}'(0) = 0$
• $\tilde{\phi}_N(0) = 0$, $\tilde{\phi}'_N(0) = 0$
• $\chi \to 0$ as $r \to \infty$

Numerical Methods:

• Bisection method for rapid chemical potential convergence
• χ field solved via radial integration of Yukawa-type Green's function
• Convergence criterion: relative field value changes < threshold

Parameter Ranges:

• gχ = 1, 2 (gravitational strength magnitude)
• β' ~ 0.3, 1.3, 5.7 (covering different scale regions)
• Multi-mediator cases: mass ratio mχ2/mχ1 = 3, 30

Evaluation Metrics

1. Density Profile: Shape and decay behavior of f̃(r̃)
2. Core Density-Radius Relation: Power-law exponent of ρc vs Rc
3. Rescaled Density: Behavior of ρcR^4_c, reflecting effective gravitational constant changes
4. Fitting Quality: Agreement between numerical solutions and theoretical asymptotic behavior

Experimental Results

Main Results

1. Field Profile Characteristics (Figure 1)

• Dark Matter Field f̃: Monotonically decreasing in all cases, approaching zero at r̃ ~ 6-8
• Newtonian Potential φ̃N: Monotonically increasing, reflecting gravitational potential
• Mediator Field χ̃:
  • β' ~ 0.3 (long-range): χ̃ significant and shape similar to φ̃N
  • β' ~ 5.7 (short-range): χ̃ exponentially decaying, effect negligible
  • χ̃ amplitude at gχ = 2 approximately twice that at gχ = 1

2. Density-Radius Relation (Figure 2)

Upper Panel (Log-Log Plot):

• Clearly shows ρc ∝ R^(-4)_c power-law relationship
• Curves for different β' values converge to same asymptote at large Rc (pure gravity case)
• Curves deviate at small Rc, reflecting new force contribution

Lower Panel (Rescaled Density):

• gχ = 1: ρcR^4_c rises from ~0.113 at small Rc to ~0.226 at large Rc, ratio 1/2 = 1/(1+gχ)
• gχ = 2: ρcR^4_c rises from ~0.075 to ~0.226, ratio 1/3 = 1/(1+gχ)
• Transition occurs near Rc ~ m^(-1)_χ

Fitting Function: $(\rho_c R^4_c)_{fit} = \frac{\zeta}{(1+g_\chi)Gm^2_\phi}\left(1 + \frac{g_\chi}{1+a(R_c m_\chi)^{-p}}\right)$ Shows excellent agreement with numerical results.

3. Multi-Mediator Cases (Figures 3-4)

Two Mediators (mχ2 = 30mχ1, gχ1 = gχ2 = 2):

• Exhibits three regions:
  1. Rc ≪ m^(-1)_χ2: Geff = (1+gχ1+gχ2)G
  2. m^(-1)_χ2 ≪ Rc ≪ m^(-1)_χ1: Geff = (1+gχ1)G
  3. Rc ≫ m^(-1)_χ1: Geff = G
• ρcR^4_c curve shows step-like rise, with "steps" appearing each time mediator mass threshold is crossed
• With mass ratio 30, intermediate region "bulge" is noticeable but not dramatic

Ablation Studies

While not explicitly labeled as ablation studies, systematic parameter variations provide equivalent analysis:

1. Effects of gχ:
   • gχ = 0: Degenerates to pure gravity case
   • gχ = 1, 2: New force contribution gradually strengthens
   • Convergence slower for larger gχ (not fully explored)
2. Effects of β':
   • β' ≪ 1: New force long-range, significantly modifies soliton structure
   • β' ≫ 1: New force exponentially suppressed, effect negligible
   • β' ~ 1: Transition region, requires numerical solution
3. Multi-Mediator Superposition Effects:
   • Single vs. dual mediator: latter provides richer scale-dependent behavior
   • Larger mass hierarchies produce more prominent intermediate plateaus

Key Findings

1. Limited Slope Improvement: For gχ ~ O(1), while new forces modify the ρc-Rc relation, they cannot significantly reduce power-law index from q = 4 to observed q ≈ 1.3
2. Scale Dependence: New force effects strongly depend on ratio of core radius to mediator Compton wavelength
3. Multi-Mediator Potential: Multiple mass-scale mediators can produce gentler ρcR^4_c rise, theoretically potentially closer to observations
4. Numerical Accuracy: High consistency between fitting function and numerical solutions validates numerical method reliability

Related Work

Dark Matter Self-Interactions

1. Elastic Scattering Self-Interactions:
   • Spergel & Steinhardt (2000) first proposed self-interacting dark matter
   • Bullet Cluster observations constrain: σ/m ≲ cm²/g
   • Primarily study dynamical friction effects
2. Differences in This Work: Studies long-range or intermediate-range forces rather than short-range scattering

Ultralight Axion Dark Matter

1. Foundational Theory:
   • QCD axions: mφ ~ 10^(-5) - 10^(-6) eV
   • Ultralight axions: mφ ~ 10^(-22) - 10^(-19) eV
   • de Broglie wavelengths reach astronomical scales
2. Observational Constraints:
   • Lyman-α forest (Irš ič et al. 2017, Armengaud et al. 2017)
   • Dwarf galaxy dynamics (Dalal & Kravtsov 2022)
   • Gravitational lensing (Powell et al. 2023)
3. Multiple Scenarios:
   • Multi-axion models (Svrček & Witten 2006)
   • Numerical simulations (Huang et al. 2023, Gosenca et al. 2023)
   • This work extends to multi-mediator new forces

Boson Stars/Solitons

1. Theoretical Studies:
   • Chavanis (2011): Analytical mass-radius relation results
   • Liebling & Palenzuela (2023): Dynamical boson stars review
   • Visinelli et al. (2018): Dilute and dense axion stars
2. Formation Mechanisms:
   • Numerical simulations show formation at galactic centers (Schwabe et al. 2016, Schive et al. 2014)
   • Formation from mini-clusters and halos (Eggemeier & Niemeyer 2019)
   • Semi-analytical work (Levkov et al. 2018)
3. This Work's Contribution: First systematic study of new force effects on solitons

Core-Cusp Problem

1. Observations: Rodrigues et al. (2017) found ρc ∝ R^(-1.3)_c
2. Theoretical Attempts:
   • Hu et al. (2000): Fuzzy dark matter proposal
   • Harko (2011), Robles & Matos (2012): Bose-Einstein condensation
   • Deng et al. (2018): Proved pure gravity case q = 4 cannot explain observations
3. This Work's Position: Explores new forces as possible solution, acknowledging limited improvements for gχ ~ O(1)

Conclusions and Discussion

Main Conclusions

1. New Forces Modify Soliton Properties: Introducing mediator field χ with mass mχ significantly changes scalar dark matter soliton structure, particularly near Rc ~ m^(-1)_χ scale
2. Effective Gravitational Constant: New force effects can be understood via scale-dependent effective gravitational constant Geff, enhanced to (1+gχ)G at small scales, returning to G at large scales
3. Density-Radius Relation Correction: While maintaining ρc ∝ R^(-4)_c asymptotic behavior, proportionality coefficient changes make overall relation less steep
4. Multi-Mediator Richness: Multiple mediators produce multiple mass thresholds, enabling more complex scale-dependent behavior
5. Limited Improvement: For gχ ~ O(1) (gravitational strength magnitude), cannot significantly reduce power-law index from q = 4 to observed q ≈ 1.3

Limitations

1. Coupling Strength Constraints:
   • Numerical methods converge slowly for gχ ≫ 1
   • Strong coupling region (potentially showing more significant effects) insufficiently explored
2. Special Coupling Choice:
   • Only considers c1 = c2/2 case
   • Other couplings (e.g., c1 = -c2/4, coupling to stress-energy tensor trace) may produce χ waves and new phenomena
3. Idealized Assumptions:
   • Neglects baryonic matter
   • Static spherically symmetric configurations
   • Small-field approximation (V(φ) ≈ m²φφ²/2)
4. Cosmological Effects Not Considered:
   • Focuses on galactic scales
   • CMB-era potential effects not detailed
5. Missing Detailed Observational Comparison:
   • No systematic galaxy sample fitting
   • Statistical significance of improvements not quantified

Future Directions

1. Strong Coupling Region: Develop more efficient numerical methods exploring gχ ≫ 1
2. General Couplings: Study χ wave production and astrophysical effects when c1 ≠ c2/2
3. Cosmological Constraints:
   • CMB-scale effects (at t ~ 380,000 years, Hubble scale ~ 180 kpc)
   • Systematic parameterization using effective field theory of large-scale structure (EFTofLSS)
4. Dynamical Studies:
   • Soliton formation processes
   • Soliton mergers
   • Non-spherically symmetric configurations
5. Observational Comparison:
   • Detailed galaxy rotation curve fitting
   • Comparison with dwarf and low surface brightness galaxy data
   • Statistical analysis of improvement magnitude
6. Equivalence Principle Violation: Explore observational consequences and constraints of dark matter violating equivalence principle at < m^(-1)_χ scales

In-Depth Evaluation

Strengths

1. Theoretical Rigor:
   • Systematic derivation of weak-field non-relativistic limit from complete action
   • Clear mathematical derivations and explicit physical pictures
   • Dimensionless treatment clarifies parameter space structure
2. Numerical Method Innovation:
   • Simple and effective iterative solution scheme
   • Integral representation for χ field avoids direct differential equation solving difficulties
   • Bisection method rapidly converges chemical potential
3. Deep Physical Insights:
   • Effective gravitational constant concept elegantly summarizes new force effects
   • Fitting function (Equation 54) successfully captures transition between asymptotic regions
   • Multi-mediator extension demonstrates rich physical possibilities
4. Clear Result Presentation:
   • Well-designed figures with high information content
   • Log-log and rescaled plots complement each other, showing results from different perspectives
   • Numerical-theoretical prediction comparisons enhance persuasiveness
5. Clear Research Motivation:
   • Addresses core-cusp problem in ultralight axion dark matter models
   • Systematic exploration of dark matter new physics
   • Establishes foundation for subsequent research

Weaknesses

1. Limited Practical Improvement:
   • For gravitational-strength couplings, cannot resolve core-cusp problem
   • Significant improvements require gχ ≫ 1, insufficiently explored
2. Insufficient Observational Verification:
   • Lacks quantitative comparison with actual galaxy data
   • No discussion of constraining parameters mχ and gχ from observations
   • No assessment of this model's advantages relative to alternatives
3. Unclear Physical Mechanisms:
   • No discussion of mediator field χ physical origin
   • Lack of fundamental theory support for specific coupling choice
   • Unclear connection to particle physics standard model
4. Incomplete Parameter Space Exploration:
   • gχ only considers values 1 and 2
   • Multi-mediator cases only consider two mediators
   • Insufficient systematic mass ratio scanning
5. Missing Dynamical Effects:
   • Only considers static solutions
   • Soliton stability not discussed
   • Formation processes and evolution not addressed
6. Insufficient Equivalence Principle Violation Discussion:
   • New forces on dark matter imply equivalence principle violation
   • Existing constraints and observable effects insufficiently discussed

Impact Assessment

Contribution to Field:

• Moderate to High: Provides new theoretical framework for dark matter self-interaction research
• Systematic numerical methods available for subsequent studies
• Multi-mediator extension opens new parameter space

Practical Value:

• Limited but Promising: Current parameters cannot resolve core problem, but indicate possible directions
• If future evidence for gχ ≫ 1 emerges, this work provides theoretical foundation
• Numerical code applicable to broader parameter exploration

Reproducibility:

• High: Equations and boundary conditions clearly described
• Numerical method steps explicit
• Lacks code release and specific numerical data tables

Limitations:

• Strongly depends on specific coupling choice (c1 = c2/2)
• Further work needed for meaningful observational comparison

Applicable Scenarios

1. Theoretical Research:
   • Theoretical work exploring dark matter new physics
   • Systematic boson star/soliton property studies
   • Modified gravity or new interaction phenomenology
2. Numerical Simulations:
   • Subgrid physics in large-scale structure simulations
   • Dark matter core modeling in galaxy formation simulations
3. Observational Constraints:
   • Provides theoretical templates for observational data analysis
   • Helps design observation strategies targeting new forces
4. Inapplicable Scenarios:
   • Standard ΛCDM cosmology (new forces suppressed at large scales)
   • Short-range scattering self-interactions (this work focuses on long/intermediate-range forces)
   • Strong relativistic regimes (this work uses non-relativistic approximation)

Key References

1. Deng et al. (2018) 51: Proved ultralight axion with pure gravity cannot explain core-cusp problem, direct motivation for this work
2. Rodrigues et al. (2017) 50: Provides observational data ρc ∝ R^(-1.3)_c, target for theoretical explanation
3. Chavanis (2011) 28, 29: Foundational work on boson star mass-radius relations
4. Hu, Barkana & Gruzinov (2000) 10: Pioneering fuzzy dark matter work
5. Spergel & Steinhardt (2000) 1: Early self-interacting dark matter proposal

Overall Assessment: This is a theoretically rigorous and methodologically innovative paper providing new perspectives on dark matter self-interaction research. While practical effectiveness in resolving the core-cusp problem is limited, it establishes systematic theoretical frameworks and numerical methods, laying groundwork for subsequent research. The paper's main value lies in exploring new physical possibilities rather than providing definitive observational explanations. Future extension to strong coupling regimes and detailed observational data comparison would significantly enhance its impact.