Cosmic Himalayas in CROCODILE : Probing the Extreme Quasar Overdensities by Count-in-Cells analysis and Nearest Neighbor Distribution

Here is an explanation of the paper, translated into everyday language with some creative analogies.

The Big Mystery: The "Cosmic Himalayas"

Imagine the universe as a giant, foggy ocean. Most of the time, the fog is evenly spread out. But every now and then, you find a massive, towering mountain rising out of the water. In astronomy, these "mountains" are huge clusters of galaxies and super-bright quasars (active black holes).

Recently, astronomers found a specific cluster called the "Cosmic Himalayas" (CH). It's a region packed with so many quasars that it looked like a statistical miracle.

The Problem:
When the scientists who found it calculated the odds of this happening using standard math (which assumes everything is a perfect bell curve, like a normal distribution of heights in a classroom), they got a number that made everyone gasp. They said the odds were 1 in 10^68.

To put that in perspective:

There are only about $10^{80}$ atoms in the entire observable universe.
The odds of winning the lottery every day for a year are roughly $1 $in$ 10^{30}$.
The "Cosmic Himalayas" was claimed to be so rare that it seemed impossible to exist in our current understanding of the universe (the $\Lambda$ CDM model). It was like finding a snowball made of pure gold in a blizzard.

The Investigation: The "CROCODILE" Simulation

The authors of this paper decided to play detective. They asked: "Is this thing actually impossible, or did we just use the wrong ruler to measure it?"

To find out, they used a super-computer simulation called CROCODILE. Think of this simulation as a massive, virtual universe where they can grow galaxies and black holes from scratch, following the laws of physics.

They didn't just look at one spot; they ran the simulation and looked for "Cosmic Himalayas" analogs—regions that looked just as crowded as the real one.

The Two Tools: Counting and Neighbors

To analyze the data, they used two different statistical tools, which they explain using two metaphors:

1. The "Count-in-Cells" (CIC) Method: The Cookie Jar
Imagine you have a giant jar of cookies (the universe). You take a small, fixed-size cup (a "cell") and scoop out cookies.

The Old Way (Gaussian): The old scientists assumed the cookies were distributed perfectly evenly, like sand. If you scooped up a handful with way more cookies than average, they calculated the odds based on a "perfect bell curve." This made the "Cosmic Himalayas" look like a one-in-a-trillion freak accident.
The New Way (AGND): The authors realized the universe isn't like sand; it's like a jar of cookies where some are clumped together in big piles. They used a new, more flexible math model (called the Asymmetric Generalized Normal Distribution or AGND).
The Result: When they used this new model, the "Cosmic Himalayas" wasn't a miracle anymore. It turned out to be a rare event, but not impossible. It's like finding a pile of cookies in the jar. It's unusual, but it happens naturally when you have clumps. The odds went from "impossible" ($10^{-68} $) to "very rare but possible" ($ 10^{-4}$).

2. The "Nearest Neighbor" (NND) Method: The Party
Imagine you are at a crowded party.

The Old View: You look at the crowd and say, "Wow, everyone is standing right next to each other! This is impossible!"
The New View: The authors realized that if you pick a specific group of people (a sample), they might naturally stand closer together just by chance or because of how you chose them.
The Result: They found that in their simulation, the "clumping" of quasars (how close they stand to each other) matched the "Cosmic Himalayas" perfectly. This proved that the tight grouping wasn't a sign of new physics; it was just a natural result of how these objects form and how we happen to look at them.

The "Goldilocks" Selection Bias

The paper also points out a subtle trick called Selection Bias.

Imagine you are looking for the tallest people in a city.

Case A (Strict): You only measure people wearing hats.
Case B (Loose): You measure everyone.

The authors found that depending on how you choose which quasars to count (like choosing people with hats vs. everyone), the "rarity" of the Cosmic Himalayas changes.

If you pick the "perfect" quasars (Case A), the cluster looks incredibly dense and rare.
If you pick a slightly broader group (Case B), the cluster looks less extreme.

This suggests that the "Cosmic Himalayas" might just be a case of looking at the right (or wrong) place at the right time. It's like taking a photo of a crowded room and zooming in on the one spot where everyone is huddled together. It looks crazy in the photo, but if you zoom out, it's just a normal crowd.

The Conclusion: No Need to Rewrite Physics

The big takeaway is simple: We don't need to throw out our current model of the universe.

The "Cosmic Himalayas" is not a glitch in the matrix. It is a natural, albeit rare, outcome of how the universe forms structures. The previous claim that it was "impossible" was simply because the scientists used a math tool (the Gaussian bell curve) that was too rigid to handle the messy, clumpy reality of the cosmos.

In short:

Old Math: "This mountain is impossible! The universe is broken!"
New Math: "This mountain is rare, but it's just a big pile of rocks. The universe is fine."

The authors successfully showed that the "Cosmic Himalayas" fits perfectly into our standard understanding of the universe, provided we use the right kind of math to count the stars.

Here is a detailed technical summary of the paper "Cosmic Himalayas in CROCODILE: Probing the Extreme Quasar Overdensities by Count-in-Cells analysis and Nearest Neighbor Distribution."

1. Problem Statement

The paper addresses a significant tension in modern cosmology regarding the "Cosmic Himalayas" (CH), an extreme overdensity of quasars at redshift $z \sim 2$ reported by Liang et al. (2025).

The Anomaly: The CH was reported to have an overdensity significance of $\delta = 16.9\sigma$ under the assumption of a Gaussian distribution, corresponding to a probability of $P \sim 10^{-68}$ .
The Challenge: Such an event is statistically impossible within the standard $\Lambda$ CDM framework, suggesting either a failure of the cosmological model, a misunderstanding of quasar formation physics, or a flaw in the statistical methodology used to quantify the overdensity.
The Hypothesis: The authors propose that the extreme significance arises from the inappropriate use of Gaussian statistics to describe the highly non-Gaussian, heavy-tailed distribution of quasar counts in the early universe.

2. Methodology

The study utilizes the CROCODILE suite of cosmological hydrodynamic simulations, which self-consistently models galaxy-black hole coevolution, alongside a cross-check using the Uchuu N-body simulation.

A. Simulation Setup (CROCODILE)

Code: gadget4-osaka (SPH code).
Volume: $200 , h^{-1} \text{cMpc} $box with$ 1024^3$ dark matter and gas particles.
Physics: Includes star formation, supernova feedback, and AGN feedback.
Key Model Updates (BF50NEL):
- No Eddington Limit (NEL): Allows for brief episodes of super-Eddington accretion, crucial for reproducing high-redshift quasar populations.
- Sound Speed (V2): Updated to use the maximum of the effective Equation of State temperature and the actual thermodynamic gas temperature, preventing unphysically high Bondi accretion rates in cold, dense gas.
- Feedback: Thermal AGN feedback coupled stochastically to gas.
Validation: The simulation reproduces observed Eddington ratio distributions, X-ray luminosity functions (XLF), and Black Hole Mass Functions (BHMF) at $z \sim 2$ .

B. Quasar Selection

Two selection criteria were applied to the simulation data to bracket observational uncertainties:

Case A: Narrower selection (Mass $8.0 \le \log(M_{BH}/M_\odot) \le 9.5 $,$ L_{bol} > 10^{45.47} $erg/s). Matches SDSS mean density ($ \bar{N} \approx 0.25 $per$ 30 , h^{-1}\text{Mpc}$ cell).
Case B: Broader selection (Mass $7.1 \le \log(M_{BH}/M_\odot) \le 9.9 $,$ L_{bol} > 10^{45} $erg/s). Higher density ($ \bar{N} \approx 1.5$).

C. Statistical Analysis

Two complementary statistical tools were employed:

Count-in-Cells (CIC):
- Divides the volume into cells of $L_{cell} = 30 \, h^{-1}\text{Mpc}$ .
- Constructs the probability distribution function (PDF) of quasar counts, $P_{CIC}(N)$ .
- Fitting Models: Compared a standard Gaussian fit against an Asymmetric Generalized Normal Distribution (AGND).
Nearest-Neighbor Distribution (NND):
- Analyzes the distance to the closest quasar pair ( $P_{NN}(<r)$ ).
- Fits the data using a power-law two-point correlation function (TPCF), $\xi(r) = (r/r_0)^{-\gamma}$ , to derive an effective correlation length ( $r_{eff}^0$ ). This isolates intrinsic clustering strength from number density effects.

3. Key Contributions

Statistical Reframing: Demonstrates that the extreme significance of the Cosmic Himalayas is an artifact of assuming Gaussianity. The paper introduces the AGND as a physically motivated, superior model for quasar counts in the high-density regime.
Simulation Validation: Confirms that CROCODILE, with updated accretion physics (NEL and V2 sound speed), successfully reproduces the observed population of luminous quasars at $z \sim 2$ .
Cross-Simulation Robustness: Validates findings using the Uchuu simulation (a $2 , h^{-1}\text{Gpc} $volume), proving that the non-Gaussian nature of the CIC distribution is a generic feature of$ \Lambda$CDM structure formation, not an artifact of the simulation box size.

4. Key Results

A. CIC Analysis: Non-Gaussianity

Gaussian Failure: The Gaussian model fails to reproduce the skewness and heavy tails of the quasar count distribution. It significantly overestimates the rarity of extreme events.
AGND Success: The AGND model provides an excellent fit to the simulation data.
Re-evaluation of Significance:
- Under Gaussian assumptions, the most extreme regions in the simulation appear as $12\sigma $outliers ($ P \sim 10^{-33}$).
- Under the AGND model, these same regions have a probability of $P \sim 10^{-4}$ .
- Conclusion: The "16.9 $\sigma$ " claim for the CH is a statistical misinterpretation. When modeled correctly, such overdensities are rare but not impossible within $\Lambda$ CDM.

B. NND Analysis: Clustering Properties

Correlation Lengths: The NND analysis reveals that the strong clustering observed in the CH ( $r_{eff}^0 \approx 30 \, h^{-1}\text{Mpc}$ ) is naturally reproduced in the simulation, specifically in Case A (the more biased selection).
Selection Bias: The apparent agreement between the CH and Case B in raw NND plots was driven by higher number density, not intrinsic clustering. When corrected for density, Case A matches the CH's clustering strength.
Implication: The strong clustering is a natural outcome of sample selection biases combined with cosmic variance, rather than a unique physical anomaly.

C. Probability Estimates

Using the AGND model to calculate the probability of finding such a structure in a Gpc-scale survey (like SDSS):

For Case A (similar to SDSS density), the probability of finding at least one cell with $N \ge 8$ quasars in a $1 , h^{-3}\text{Gpc}^3 $volume is **$ P \approx 0.71$**.
This confirms that the existence of the Cosmic Himalayas is highly probable in a standard $\Lambda$ CDM universe.

5. Significance

Resolution of the Tension: The paper resolves the apparent conflict between the existence of the Cosmic Himalayas and the $\Lambda$ CDM model. The "anomaly" is a statistical illusion caused by applying Gaussian statistics to a non-Gaussian distribution.
Methodological Shift: It establishes that extreme overdensities in the high-redshift universe must be quantified using non-Gaussian statistics (like AGND) rather than standard $\sigma$ -based Gaussian metrics.
Cosmological Consistency: The findings reinforce the robustness of the $\Lambda$ CDM paradigm, showing that extreme structures are natural, albeit rare, outcomes of hierarchical structure formation and biased tracer selection.
Future Directions: The study highlights the need for careful treatment of selection functions and statistical models when interpreting future high-redshift quasar surveys (e.g., from LSST or Euclid).