Effects of Self-Interaction and of an Ideal Gas in Binary Mergers of Bosonic Dark Matter Cores

This study uses 3D simulations to demonstrate that while binary mergers of bosonic dark matter cores consistently produce stable solitonic remnants, the final mass retention is non-universal and critically depends on whether repulsive or attractive self-interactions are present, whereas the inclusion of an ideal gas component modifies the gravitational background without altering the intrinsic scaling of the bosonic core.

The Gist

Imagine the universe is filled with a mysterious, invisible substance called Dark Matter. For a long time, scientists thought this stuff was just cold, clumpy dust. But a newer, cooler theory suggests it might actually be a giant, cosmic wave made of ultra-light particles. Think of it less like sand and more like a giant, invisible ocean of "fuzzy" energy.

In this paper, the authors act like cosmic weather forecasters. They wanted to see what happens when two massive "storms" of this fuzzy dark matter crash into each other. Do they merge into one giant storm? Do they bounce off? Do they lose a lot of energy?

To find out, they ran super-computer simulations of these collisions, but they added two new ingredients to the mix to see how they change the outcome:

1. Self-Interaction: What if the dark matter particles can push or pull on each other?
2. Ideal Gas: What if there's a cloud of normal gas (like the air in a balloon) mixed in with the dark matter?

Here is the breakdown of their findings using some everyday analogies:

1. The Basic Collision: The "Magic 60%" Rule

First, they looked at the "pure" version of this dark matter (no extra forces, no gas). When two fuzzy dark matter cores crash, they don't just stick together perfectly. It's like two waves crashing into each other; some of the water splashes out.

• The Result: No matter how big the two waves were, the final merged wave always kept about 60% to 63% of the original mass. The rest was "ejected" into a diffuse halo around the center.
• The Analogy: Imagine two people trying to merge their piles of sand into one big pile. As they push them together, some sand flies off into the air. Surprisingly, they always end up with roughly the same percentage of sand left in the main pile, regardless of how hard they pushed. The authors call this a "magic fraction."

2. Adding "Self-Interaction": The Rubber Band vs. The Magnet

Next, they asked: What if the dark matter particles have a personality? What if they can push each other away or pull each other closer?

• Repulsive Interaction (The Rubber Band): Imagine the particles are like people wearing bouncy rubber bands. If they get too close, the bands push them apart.

  • What happened: When the authors simulated this, the final core was heavier and held onto more mass. The "rubber bands" acted like extra pressure, keeping the core tight and preventing sand from flying off.
  • The Takeaway: Repulsion helps the core stay strong and keep its mass.

• Attractive Interaction (The Magnet): Now imagine the particles are magnets. They want to stick together tightly.

  • What happened: When they simulated this, the final core was lighter. The strong pull made the collision more violent, causing more mass to be flung out into space.
  • The Takeaway: Being too clingy (attractive) actually makes the merger messier and causes you to lose more of your stuff.

3. Adding the "Ideal Gas": The Ghost in the Machine

Finally, they added a cloud of normal gas (like the air in a room) to the dark matter ocean. This is like asking: "What happens if a fuzzy dark matter storm crashes into a cloud of regular air?"

• The Result: The dark matter didn't care much. Even if the gas was 10 times heavier than the dark matter, the dark matter still formed a tight, compact core and kept that same ~60% of its mass.
• The Analogy: Imagine a heavy, dense bowling ball (the dark matter core) rolling through a room full of fluffy cotton candy (the gas).
  • The cotton candy gets pushed around, swirls, and spreads out. It doesn't form a tight ball on its own.
  • But the bowling ball? It just keeps rolling, ignoring the cotton candy. The cotton candy just acts like a soft background that the bowling ball moves through.
• The Takeaway: The dark matter core is incredibly robust. It forms a solid "heart" regardless of whether it's surrounded by empty space or a thick fog of gas. The gas just follows the dark matter's lead, never taking control.

The Big Picture: Why Does This Matter?

The authors explain these results using "energy math."

• No Gas/No Interaction: The energy rules are strict, leading to that universal 60% rule.
• With Repulsion: The rules change slightly, allowing the core to keep more mass (like a tighter knot).
• With Gas: The gas is too "loose" to change the rules of the dark matter. It's just a passenger, not the driver.

In simple terms:
This paper tells us that the "fuzzy" dark matter cores are like tough, resilient seeds. Whether they are pushed apart by rubber bands or pulled together by magnets, they change how much mass they keep. But if you put them in a room full of gas, they don't care; they will still form a solid, dense core, leaving the gas to float around them.

This helps astronomers understand why we see dense cores in the centers of galaxies, even if those galaxies are full of gas and stars. It suggests that the dark matter core is the "boss" of the galaxy's center, and it's very hard to knock it off its throne.

Technical Summary

Here is a detailed technical summary of the paper "Effects of Self-Interaction and of an Ideal Gas in Binary Mergers of Bosonic Dark Matter Cores."

1. Problem Statement

The paper investigates the dynamics of binary mergers of solitonic cores in Bose-Einstein Condensate Dark Matter (BECDM), also known as Fuzzy Dark Matter (FDM). While previous studies (e.g., Schwabe et al., 2016) established that binary mergers of non-interacting solitons result in a "nearly universal" final core mass fraction ($M_{c,final} \approx 0.6-0.7$ of the initial total mass) due to gravitational cooling, two critical physical extensions remained unexplored in controlled 3D simulations:

1. Self-Interaction: How does the inclusion of a quartic self-interaction term (repulsive or attractive) in the Gross-Pitaevskii-Poisson (GPP) system alter the merger efficiency and final mass retention?
2. Coupling to Baryonic Matter: How does the presence of a gravitationally coupled ideal gas (modeling Fermion-Boson Stars) affect the merger process? Specifically, does the gas component disrupt the formation of the solitonic core or alter the mass loss mechanisms?

2. Methodology

The authors performed high-resolution 3D numerical simulations using the CAFE-FDM code to model two distinct physical scenarios:

A. Self-Interaction Scenario (GPP System)

• Governing Equations: The Gross-Pitaevskii-Poisson (GPP) system, where the Schrödinger equation includes a self-interaction term $g|\Psi|^2$.
  • $i\partial_t\Psi = (-\frac{1}{2}\nabla^2 + V + g|\Psi|^2)\Psi$
  • $\nabla^2 V = |\Psi|^2$
• Parameters: The self-interaction constant $g$ was varied across values: $g = -0.1$ (attractive), $0$(non-interacting),$0.1, 0.2, 0.3$ (repulsive).
• Initial Conditions: Binary mergers of stationary ground-state solutions (Newtonian Boson Stars). Asymmetry was introduced via a scaling factor $\lambda$ to control mass ratios, impact parameters, and initial velocities.
• Numerical Setup:
  • Domain: $[-20, 20]^3$ with $N=128$ grid points ($h=0.3125$).
  • Time integration: 3rd-order Runge-Kutta.
  • Boundary Conditions: Periodic for the wavefunction (to simulate an isolated system with re-entering modes) and Dirichlet ($V=0$) for the gravitational potential.

B. Ideal Gas Coupling Scenario (GPPE System)

• Governing Equations: The Schrödinger-Poisson-Euler (GPPE) system, coupling the bosonic field to a compressible ideal fluid.
  • Includes continuity, Euler momentum, and energy equations for the gas, coupled via a shared gravitational potential $V$.
  • Equation of State: Polytropic ($p = K\rho^{1+1/n}$) initially, evolving to an ideal gas ($p = (\gamma-1)\rho e$).
• Initial Conditions: Mergers of Fermion-Boson Stars (FBS), where the bosonic core is surrounded by a gas envelope.
• Parameters: The mass ratio between gas and bosons ($M_R = M_{gas}/M_{BEC}$) was varied: $0.1, 1.0, 10.0$ (representing gas-dominated regimes).
• Numerical Scheme: The Euler equations were solved using High-Resolution Shock-Capturing (HRSC) methods (Minmod limiter, HLLE Riemann solver).

3. Key Contributions

• Extension of Merger Universality: The study generalizes the "magic fraction" of mass retention found in non-interacting FDM to regimes with self-interaction and gas coupling.
• Energetic Scaling Theory: The authors provide a unified theoretical explanation based on energy scaling relations ($E \propto -M^a$). They demonstrate that the final mass fraction is determined by the exponent $a$ of the equilibrium energy scaling, which shifts depending on the dominant pressure support (quantum pressure vs. self-interaction pressure).
• Robustness of Solitonic Cores: They establish that solitonic cores are robust attractors even in environments dominated by baryonic (gas) matter, provided the coupling is purely gravitational.

4. Key Results

A. Effects of Self-Interaction

• Mass Retention Trends:
  • Repulsive Interaction ($g > 0$): Increases the final core mass fraction. As $g$ increases, the average retention ratio $R_{Mc}$ rises from $\sim 0.60$ (for $g=0$) to $\sim 0.68$ (for $g=0.3$). Repulsion acts as an additional pressure term, stabilizing the core against tidal disruption and reducing mass ejection.
  • Attractive Interaction ($g < 0$): Decreases the final core mass fraction ($R_{Mc} \approx 0.56$). Attraction enhances gravitational binding, leading to more violent relaxation and greater mass ejection into the halo.
• Theoretical Explanation:
  • For $g=0$, the system is quantum-pressure dominated, yielding an energy scaling $E \propto -M^3$. This leads to a theoretical mass fraction of $2^{-2/3} \approx 0.63$.
  • For $g > 0$ (Thomas-Fermi regime), self-interaction dominates, shifting the scaling toward $E \propto -M^2$, which theoretically yields a higher fraction ($2^{-1/2} \approx 0.71$).
  • For $g < 0$, the scaling becomes steeper ($a > 3$), resulting in lower retention fractions.

B. Effects of Ideal Gas Coupling

• Core Formation: In all cases, including when the gas mass dominates ($M_R = 10.0$), a stable solitonic bosonic core always forms. The gas component does not form a compact core; instead, it remains diffuse and settles into the potential well created by the bosons.
• Mass Retention:
  • The bosonic component retains a nearly universal mass fraction ($R_{Mc} \approx 0.62$) regardless of the gas dominance ($M_R$).
  • The gas component shows no universal retention fraction; its $R_{Mc}$ decreases as $M_R$ increases due to hydrodynamic pressure preventing central concentration.
• Mechanism: The gas modifies the global gravitational background but does not alter the intrinsic cubic energy scaling ($E \propto -M^3$) of the bosonic sector. Since the gas is spatially more extended than the core, its potential acts as a slowly varying background, leaving the bosonic condensation physics unchanged.

5. Significance

• Astrophysical Implications: The results suggest that the central density profiles and core-halo mass relations of dark matter halos are sensitive to the fundamental particle physics (self-interaction strength) but remain robust against the presence of baryonic gas.
• Model Calibration: The findings provide critical constraints for semi-analytical models of galaxy formation. Specifically, the "magic fraction" is not a single universal constant but depends on the self-interaction parameter $g$.
• Dark Matter Nature: The study supports the idea that solitonic cores are stable attractors in BECDM models. If observations of galactic cores show mass fractions significantly deviating from the $0.6-0.7$ range, it could imply specific constraints on the self-interaction strength of dark matter particles.
• Methodological Advance: The successful simulation of coupled Fermion-Boson Stars in 3D with hydrodynamic shocks provides a new benchmark for studying the interplay between dark matter and baryonic matter in the early universe and galactic centers.