Correlations of Feed-down Hadrons in a Thermal Model

This paper demonstrates that strong decays of hadrons significantly alter the yields and correlations of stable particles in a thermal hadron gas model, thereby substantially impacting measurements of net quantum number cumulants and balance functions that are crucial for identifying the QCD critical point and determining chemical susceptibilities.

The Gist

The Big Picture: The Cosmic Particle Explosion

Imagine a high-energy collision between two heavy atoms (like gold or lead) as a massive, microscopic explosion. When these atoms smash together, they create a tiny, super-hot "soup" of fundamental particles called a Quark-Gluon Plasma. As this soup cools down, it freezes into a gas of particles (hadrons) like protons, neutrons, pions, and kaons.

Scientists study this explosion to understand the fundamental rules of the universe, specifically looking for a "Critical Point"—a special state of matter where the rules of physics change dramatically. To find this, they measure fluctuations (how much the number of particles varies from one explosion to the next) and correlations (how likely it is that two specific particles appear together).

The Problem: The "Family Tree" Confusion

Here is the catch: The particles scientists can actually see in their detectors are not always the ones that were originally created in the hot soup.

Think of the initial explosion as a family having children.

• The Parents: Heavy, unstable particles (resonances) created in the soup.
• The Children: Lighter, stable particles (like pions and protons) that the parents decay into.

The paper argues that scientists have been looking at the "children" (the stable particles) and assuming they represent the "parents" perfectly. But this is misleading. A single heavy parent particle can decay into multiple children. Sometimes, a parent splits into two children of the same type; sometimes, it splits into different types.

This process is called Feed-Down (FD). It's like a genealogy tree where you are trying to count the total number of people in a family, but you are only counting the grandchildren. If you don't account for the fact that one grandparent can have five grandchildren, your math will be way off.

The Main Findings: Why "Feed-Down" Matters

The authors used a computer model (a "thermal model") to simulate this particle soup and calculate exactly how much the "children" (decays) mess up the measurements. Here is what they found:

1. The "Ghost" Protons

Scientists often use protons as a stand-in to measure the total "baryon number" (a type of matter count) of the explosion. They assume that if they count the protons, they know the total amount of matter.

• The Analogy: Imagine you are counting the number of "Red Cars" in a parking lot to estimate the total number of cars. But, you discover that many "Blue Trucks" (heavy particles) break apart and turn into "Red Cars" right before you count them.
• The Result: The paper shows that the number of protons you see is heavily inflated by these "Blue Trucks" turning into protons. In fact, the protons you see are mostly "orphans" from decays, not the original "parents." This means using proton counts to measure the total matter is like trying to guess the size of a crowd by counting only the people wearing red hats, when half the crowd got those hats from a hat factory explosion.

2. The "Chaos" of Pions

Pions are the most common particles produced. The study found that the number of pions is dominated by decays.

• The Analogy: Imagine a fireworks show. You see a lot of sparks (pions). You might think the fireworks were small, but actually, a few giant shells exploded, creating thousands of sparks. The "sparks" (pions) are mostly the result of the big explosions (heavy resonances), not the initial small sparks.
• The Result: The fluctuations (random ups and downs) in the number of pions are mostly caused by these giant explosions, not by the underlying physics of the soup itself. If you ignore this, you might think the soup is more chaotic than it really is.

3. The "Balance Sheet" Error

In physics, there are "balance functions." If you create a positive particle, you usually expect to see a negative particle nearby to balance the charge (like a bank account: +$10 and -$10).

• The Analogy: Imagine a party where people arrive in pairs (a man and a woman). If you look at the dance floor, you expect to see equal numbers of men and women. But, if a "couple" (a heavy particle) splits up and the man runs off to dance with two other women, your count gets weird.
• The Result: The study shows that decays create "unbalanced" pairs. A heavy particle might decay into two positive pions. This breaks the simple "one positive, one negative" rule scientists use to test their theories. The "balance" is hidden in the complex family tree of decays.

The Temperature Twist

The paper also looked at how temperature changes things.

• The Analogy: Think of the particle soup as a melting pot. At lower temperatures, you mostly have simple ingredients (light particles). As it gets hotter, you start melting down complex, heavy ingredients (heavy resonances).
• The Result: As the temperature rises, the number of heavy particles explodes. This means the "Feed-Down" effect gets much stronger. The "children" (decays) become the dominant feature of the data. This is crucial because if scientists are trying to find the "Critical Point" by changing the temperature of the collision, they might mistake the increase in decays for a phase transition in the matter itself.

The Bottom Line

This paper is a warning label for physicists. It says: "Don't just count the particles you see; ask who their parents were."

If you want to understand the fundamental properties of the universe (like the Critical Point or how matter behaves near a phase transition), you cannot ignore the "feed-down" from heavy particles decaying into lighter ones. If you do, your measurements of fluctuations and correlations will be distorted, potentially leading you to the wrong conclusions about the nature of the universe.

In short: The "stable" particles we measure are mostly the grandchildren of the real action. To understand the family history, we have to trace the whole tree, not just count the leaves.

Technical Summary

1. Problem Statement

The paper addresses a critical challenge in high-energy nuclear physics: distinguishing between intrinsic properties of the Quark-Gluon Plasma (QGP) and background effects caused by hadronic decays.

• Context: Measurements of fluctuations in net conserved quantum numbers (baryon number, charge, strangeness) and balance functions are primary tools for searching for the QCD critical point and determining chemical susceptibilities near the phase transition (e.g., at RHIC's Beam Energy Scan and LHC).
• The Issue: These measurements are heavily contaminated by Feed-Down (FD) from resonance decays. While strong and weak decays conserve global quantum numbers, they alter the yields and correlations of observable "stable" particles (e.g., protons, pions, kaons).
• The Gap: Experimental identification of specific parent resonances in heavy-ion collisions is often impossible due to high multiplicity and combinatorial backgrounds. Consequently, there is a lack of baseline theoretical understanding regarding how much resonance decays distort fluctuation cumulants and balance functions, potentially masking or mimicking critical point signals.

2. Methodology

The authors employ a Grand Canonical Thermal Hadron Gas Model (GCM) to quantify the impact of feed-down.

• Model Setup:
  • System: A thermalized hadron gas at vanishing baryo-chemical potential ($\mu_B = 0$).
  • Temperature Range: $T = 140$ to $200$ MeV, covering typical chemical freeze-out temperatures.
  • Particle List: A comprehensive database of 353 known hadrons with masses up to $2.5$ GeV/$c^2$ (excluding charm and bottom due to low thermal yield).
  • Decay Chains: The model includes all known two, three, and four-prong decay channels with experimental branching fractions.
• Calculation Procedure:
  1. Thermal Density: Calculate the initial thermal density ($\rho^{TH}$) of all species using the standard thermal distribution (involving modified Bessel functions).
  2. Iterative Feed-Down: Compute the probability ($P^\alpha_a$) of a parent species $\alpha$ decaying into a stable daughter $a$. This is done iteratively from the lightest to the heaviest species, summing all possible decay paths.
  3. Total Yield: The final observable density is the sum of the primary thermal yield and the feed-down contribution: $\rho^{TH+FD} = \rho^{TH} + \rho^{FD}$.
  4. Correlations: The authors calculate two-particle cumulants ($C_2$) and balance functions. They distinguish between "primary daughters" (directly from a parent) and "cousins" (produced via sequential decays) to determine correlated pair densities ($\rho^{FD}_{ab}$).
• Validation: The thermal model results are compared against Pythia 8.3 Monte Carlo simulations for $pp$ collisions to check the robustness of balance function fractions.

3. Key Contributions

• Quantification of FD Impact: The paper provides a rigorous, quantitative baseline for how resonance decays modify the yields and correlations of stable hadrons, moving beyond qualitative assumptions.
• Decoupling of Net-Quantum Numbers: The study explicitly demonstrates that while decays conserve net baryon number, they do not conserve the number of specific species (e.g., net protons). A baryon decay can yield a neutron or a $\Lambda$ instead of a proton, causing large fluctuations in the observable proton count.
• Balance Function Analysis: The authors introduce a method to normalize correlated pair densities by their integrals to analyze charge/baryon balancing, revealing that non-balancing correlations are highly sensitive to temperature.

4. Key Results

• Yield Enhancement:
  • Feed-down drastically increases the yields of stable particles. For pions, the yield increases by a factor of 2.64 at $T=140$ MeV and 10.6 at $T=200$ MeV.
  • Proton yields increase by a factor of ~3.1 at $T=160$ MeV.
  • The density of observable protons is only ~1/3 of the total thermal baryon density, making protons a poor proxy for total baryon number fluctuations.
• Fluctuation Distortion:
  • The standard deviation of the feed-down contribution to the proton yield is nearly as large as the thermal yield itself. This implies that net-proton fluctuation cumulants are dominated by decay stochasticity rather than the underlying baryon number physics.
  • Correlated pair densities from decays ($C_2$) often exceed uncorrelated thermal expectations by factors of 10 to 20 (e.g., $C_2(p|\pi^+)$ is ~23 times larger than the product of thermal densities).
• Temperature Dependence:
  • Charge Balancing: Correlations that balance charge (e.g., $\pi^+$ with $\pi^-$) show weak sensitivity to temperature.
  • Non-Balancing Correlations: Correlations that do not balance charge (e.g., $\pi^+$ with $\pi^+$, or $K^+$ with $K^+$) increase significantly with temperature (up to 108% increase for $p+\pi^+$ from 140 to 200 MeV). This is driven by the population of heavier resonances at higher temperatures.
• Balance Functions:
  • The fractional integral of balance functions (e.g., the fraction of $\pi^+$ balanced by $\pi^-$) remains relatively constant (~89% for $\pi^-$) across temperatures and matches Pythia predictions for $pp$ collisions. This suggests that balance functions are robust against thermalization assumptions but may not be sensitive enough to distinguish the degree of thermalization in A-A collisions.

5. Significance and Implications

• Critical Point Search: The results warn that observed fluctuations in net-proton numbers at RHIC BES may be driven largely by the reduction of feed-down at lower temperatures rather than a phase transition or critical point. Correcting for FD is essential to extract true chemical susceptibilities.
• Experimental Strategy: The paper suggests that "non-balancing" correlations (e.g., same-sign pairs) serve as a better "thermometer" for the chemical freeze-out temperature than traditional charge-balancing pairs, due to their strong dependence on the population of heavy resonances.
• Modeling Limitations: The study highlights that current thermal models (which assume no initial correlations) and experimental measurements (limited by acceptance and $p_T$ cuts) need more sophisticated modeling to fully disentangle decay effects from medium dynamics.
• General Applicability: The findings apply to both elementary ($pp$) and heavy-ion ($AA$) collisions, indicating that feed-down is a fundamental background process that must be accounted for in all fluctuation and correlation analyses.

In conclusion, the paper establishes that feed-down is not a minor correction but a dominant factor in the measurement of net quantum number fluctuations and balance functions. Ignoring these effects risks misinterpreting experimental data regarding the QCD phase diagram and the search for the critical point.