Islands in Simulated Cosmos: Probing the Hubble Flow around Groups and Clusters

Using IllustrisTNG simulations, this study demonstrates that while local Hubble flow kinematics can constrain halo masses and the Hubble constant with reasonable accuracy, intrinsic environmental variance fundamentally limits their ability to distinguish between different dark energy models or specific dynamical extensions.

The Gist

Imagine the universe as a giant, expanding balloon. As it inflates, everything on its surface moves away from everything else. This is the Hubble Flow: the general rule that the universe is stretching, and galaxies are drifting apart like dots on that balloon.

But sometimes, things get messy. If you have a heavy weight (like a massive galaxy cluster) sitting on the balloon, it creates a dip. Nearby dots might stop drifting away and start rolling toward the weight because of gravity. This creates a tug-of-war: the universe wants to pull them apart, but gravity wants to pull them together.

This paper, titled "Islands in Simulated Cosmos," is like a scientific detective story. The authors, David Benisty and Antonino Del Popolo, wanted to see if we can use this "tug-of-war" zone to measure two very important things:

1. How heavy these galaxy clusters are (their mass).
2. How fast the universe is actually expanding (the Hubble Constant, $H_0$).
3. Can we detect Dark Energy? (The mysterious force pushing the universe apart).

Here is the breakdown of their investigation using simple analogies:

1. The Cosmic Laboratory (The Simulation)

Instead of looking at the real, messy universe (where we can't know the "true" answers), the authors used a super-computer simulation called IllustrisTNG.

• The Analogy: Imagine a video game like SimCity or Minecraft, but for the entire universe. They built a digital universe with perfect rules. They know exactly how heavy every digital galaxy is and exactly how fast the digital universe is expanding.
• The Goal: They wanted to see if their mathematical formulas could look at this digital universe, guess the answers, and get it right. If the formula works on the "perfect" simulation, maybe it works on the real universe too.

2. The "Turnaround" Zone

They focused on a specific boundary called the Turnaround Radius.

• The Analogy: Think of a river flowing downstream (the Hubble Flow). If you drop a leaf near a giant whirlpool (a galaxy cluster), the leaf might get pulled in. But if you drop the leaf far enough away, the river current is too strong, and it keeps flowing downstream.
• The Turnaround Radius is the exact distance where the leaf stops flowing downstream and starts getting sucked into the whirlpool. The authors measured the speed of "leaves" (galaxies) at different distances to see where this switch happens.

3. The Detective Work (The Math)

They used a set of equations (based on the Lemaître–Tolman framework) to model this motion. They tried different versions of the math:

• Version A: Simple gravity only.
• Version B: Gravity + spinning (angular momentum).
• Version C: Gravity + spinning + friction (like air resistance).
• Version D: Gravity + Dark Energy models.

They fed the data from their simulation into these equations to see which version matched the "real" digital universe best.

4. The Results: What They Found

✅ The Good News: Mass and Expansion Speed
The method worked surprisingly well for the basics!

• Mass: When they guessed the weight of the galaxy clusters, they were off by less than 1% on average. It's like weighing a truck and guessing within a few pounds.
• Expansion Rate ($H_0$): They guessed the expansion speed of the universe very accurately (within about 14%).
• Takeaway: If you look at the "tug-of-war" zone around a galaxy, you can figure out how heavy it is and how fast the universe is expanding.

❌ The Bad News: Dark Energy
This is the big twist. They hoped to use this method to figure out the nature of Dark Energy (the force pushing the universe apart).

• The Problem: The universe is messy. Even in their perfect simulation, every galaxy cluster has a slightly different history. Some are spinning, some are being squished by neighbors, some are lopsided.
• The Analogy: Imagine trying to hear a whisper (Dark Energy) in a crowded, noisy stadium (the messy local universe). The "noise" of the galaxy's own gravity and its weird shape is so loud that it drowns out the whisper.
• The Result: The different models for Dark Energy all looked almost identical once they added in the "noise" of the real environment. The math couldn't tell them which Dark Energy model was the right one because the natural differences between galaxy clusters were too big.

5. The Conclusion

The paper concludes that while the "Hubble Flow" method is a great tool for weighing galaxy clusters and measuring the general expansion of the universe, it is not precise enough to solve the mystery of Dark Energy.

• Why? Because the local environment (neighboring galaxies, weird shapes, past collisions) creates too much "static" in the data.
• The Future: To solve the Dark Energy mystery, we need better tools. The authors suggest that future telescopes (like JWST or Euclid) that can measure distances more precisely, combined with better math that accounts for the "lumpiness" of the universe, might eventually crack the code.

In a nutshell:
The authors built a digital universe to test a measuring stick. They found the stick is great for measuring weight and speed, but it's too wobbly to measure the mysterious force (Dark Energy) pushing the universe apart, because the local neighborhood is just too chaotic.

Technical Summary

1. Problem Statement

The local Hubble flow—the expansion of the universe observed in the vicinity of galaxy groups and clusters—represents a critical transition zone between the inward pull of nonlinear gravitational collapse and the outward push of cosmic expansion. While the Hubble law holds on large scales, local environments exhibit significant deviations due to gravitational interactions.

A key open question in cosmology is whether the local Hubble flow can serve as a precision diagnostic for Dark Energy (DE). Previous studies have attempted to constrain DE models or the Hubble constant ($H_0$) using velocity–distance relations around nearby structures. However, it remains unclear whether the intrinsic variance between different halos (galaxy groups/clusters) in a realistic universe is small enough to distinguish between different cosmological models (e.g., $\Lambda$CDM vs. dynamical DE) or if environmental noise fundamentally limits this method.

2. Methodology

The authors utilize the IllustrisTNG300 cosmological simulation as a controlled "cosmic laboratory" to test the robustness of the Lemaître–Tolman (LT) framework.

• Simulation Data: They selected isolated halos from the TNG300 simulation (volume $\sim 300$ Mpc, high resolution) using strict criteria to ensure isolation (no neighboring halos $>10\%$ of host mass within the turnaround radius).
• Theoretical Framework: They modeled the kinematics using an extended LT formalism. The equation of motion for a test particle includes:
  • Newtonian gravity ($-GM/r^2$).
  • Angular momentum ($l^2/r^3$).
  • Cosmological acceleration ($\ddot{a}/a$).
  • A phenomenological dynamical friction term ($-\eta \dot{r}$).
• Velocity-Radius Law: They derived an empirical velocity-radius relation:
  $$v(r) = -AH_0 \left(\frac{GM}{H_0^2}\right)^{\frac{n+1}{3}} r^{-n} + bH_0r$$
  Where $A$, $b$, and $n$ are dimensionless coefficients dependent on physical assumptions (e.g., presence of angular momentum, friction, or specific DE models).
• Model Variants: They tested four specific model configurations (Table I in the paper):
  1. $M_{LT}$: No angular momentum, no friction.
  2. $J_{LT}$: Includes specific angular momentum.
  3. $J\eta_{LT}$: Includes both angular momentum and dynamical friction.
  4. $TA$ (Timing Argument): A simplified case with fixed $n=0.5$.
• Analysis Technique:
  • Constructed median radial velocity profiles for isolated halos using dark matter and star particles.
  • Employed Bayesian inference using the PolyChord nested sampling algorithm to fit the model parameters ($A, b, n$) and recover halo masses ($M$) and $H_0$.
  • Used $R^2$ (coefficient of determination) to weight the fits and assess quality.

3. Key Contributions

• Systematic Uncertainty Quantification: Unlike previous works that relied on observational statistical errors, this study isolates systematic uncertainties inherent to the physical environment (e.g., tidal fields, triaxiality, merger history) by using a simulation where the "true" values are known.
• Model Discrimination Test: The paper rigorously tests whether the intrinsic halo-to-halo variance allows for the discrimination of different physical models (e.g., distinguishing between models with/without friction or different DE equations of state).
• Validation of the Timing Argument: It evaluates the validity of fixing the scaling index $n=0.5$ (as per the classical Timing Argument) against data-driven free-parameter fits.

4. Results

A. Parameter Recovery and Variance

• Coefficients ($A, b, n$):
  • $A$ (Infall amplitude): Shows large scatter ($\sim 80\%$), reflecting diverse accretion histories.
  • $b$ (Expansion normalization): More tightly constrained ($\sim 25\%$ scatter).
  • $n$ (Radial scaling): Exhibits substantial variation ($\sim 67\%$ scatter). Isolated, relaxed halos cluster around $n \approx 0.96$ (closer to Newtonian infall), while the Timing Argument value ($n=0.5$) is significantly lower.
• Fixed vs. Free $n$: Fixing $n=0.5$ (Timing Argument) yields consistent $A$ and $b$ values but results in systematically lower fit quality ($R^2$), particularly for massive clusters and dense environments. This suggests a universal $n=0.5$ does not capture the diversity of real halo velocity profiles.

B. Mass and $H_0$ Recovery

• Halo Mass: The recovered masses ($M_{fit}$) are unbiased on average but show significant scatter.
  • Median ratio: $\langle M_{fit}/M_{true} \rangle = 0.991 \pm 0.148$.
• Hubble Constant: The inferred $H_0$ shows a slight systematic overestimation.
  • $\langle H_{0,fit} \rangle = 68.25 \pm 9.82$ km s$^{-1}$ Mpc$^{-1}$ (True input: $67.74$).
  • This represents a $\sim 14\%$ overestimate, driven by systems with complex velocity fields where the simplified LT model fails.

C. Model Discrimination (The Core Finding)

• Indistinguishability: While different physical assumptions (e.g., adding friction or angular momentum) theoretically predict different values for $A, b,$ and $n$, the intrinsic environmental scatter in the simulation is so large that these differences are statistically indistinguishable.
• Conclusion on Dark Energy: The method cannot currently constrain the Dark Energy equation of state ($w$) because the "noise" from local environmental dynamics (tidal shear, triaxiality, assembly history) overwhelms the subtle signatures of different DE models.

5. Significance and Implications

• Fundamental Limitation: The study demonstrates that the local Hubble flow method has a fundamental floor of uncertainty due to nonlinear environmental effects. Even with perfect sampling (infinite tracer galaxies), the variance in fitted parameters remains a signature of the environment's deviation from the idealized spherical model.
• Practical Utility: Despite the inability to constrain DE, the method is robust for estimating halo masses and the Hubble constant if quality cuts are applied. Systems with high fit quality ($R^2 > 0.8$) yield mass estimates within $\sim 12\%$ and $H_0$ within $\sim 14\%$ of true values.
• Future Directions: To improve precision, future models must relax the assumption of spherical symmetry, incorporating ellipsoidal collapse models and accounting for local shear and triaxiality. The authors suggest that upcoming high-precision distance measurements (e.g., from JWST, Euclid, DESI) combined with these improved dynamical models will be necessary to push beyond current limitations.

Summary Verdict: The paper concludes that while the local Hubble flow is a valuable probe for $H_0$ and mass, it is not currently a precision tool for distinguishing Dark Energy models due to the overwhelming intrinsic variance of the cosmic web's local dynamics. Previous claims of constraining DE via local flows likely underestimated these systematic physical uncertainties.