Comparing effective temperatures in standard and Tsallis distributions from transverse momentum spectra in small collision systems

This study investigates transverse momentum spectra of light charged hadrons in RHIC d+Au and p+p collisions at $\sqrt{s_{NN}}=200$ GeV, revealing that effective temperatures derived from standard and Tsallis distributions exhibit systematic decreases across distribution types and collision centralities, while maintaining perfect linear relationships between them.

The Gist

Imagine you are a detective trying to figure out how hot a pot of soup is, but you can't stick a thermometer in it. Instead, you have to guess the temperature by watching how fast the ingredients (carrots, potatoes, and peas) are flying around inside the pot.

This paper is about doing exactly that, but with subatomic particles instead of vegetables, and particle colliders instead of soup pots.

Here is the story of what the scientists did, explained simply:

1. The Big Experiment: The "Smash"

The scientists looked at data from two types of high-speed crashes:

• Proton-Proton (p+p): Like smashing two tiny marbles together.
• Deuteron-Gold (d+Au): Like smashing a small marble into a heavy bowling ball.

They did this at the Relativistic Heavy Ion Collider (RHIC), which is basically a giant racetrack where particles zoom around at nearly the speed of light before crashing. When they crash, they create a tiny, super-hot fireball of energy that instantly cools down and sprays out particles (like pions, kaons, and protons).

2. The Mystery: How Hot Was It?

The scientists wanted to know the "Effective Temperature" of this fireball. But here's the catch: in the quantum world, there isn't just one way to measure "temperature." It depends on the math tool you use.

Think of it like measuring the "spiciness" of a curry.

• Tool A (Boltzmann): A standard, old-school ruler. It assumes everyone behaves like a calm, average person.
• Tool B (Bose-Einstein/Fermi-Dirac): A high-tech laser scanner. It knows that some particles are "bosons" (like social butterflies who love to clump together) and others are "fermions" (like introverts who hate sharing space).
• Tool C (Tsallis): A futuristic AI that accounts for chaos. It knows that in a high-energy crash, things aren't perfectly calm; there are wild fluctuations and "hot spots."

3. The Investigation

The team took the same data (the flying particles) and ran it through all three math tools to see what temperature they each calculated.

What they found:

• The "Hotness" Ranking: The tools didn't agree on the exact number.

  • The Bose-Einstein/Fermi-Dirac tools (the high-tech scanners) gave the "truest" baseline temperature.
  • The Boltzmann tool (the old ruler) was a bit too simple. It underestimated the heat for the "social butterfly" particles and overestimated it for the "introvert" particles.
  • The Tsallis tool (the AI) consistently said the system was cooler than the others. Why? Because it accounts for the chaos and disorder, which smooths out the peaks of extreme heat.

• The "Crowdedness" Factor: They looked at how "central" the crash was.

  • Central collisions (hitting the bullseye) are like a packed mosh pit. Everyone is bumping into everyone else, creating a very hot, chaotic, but surprisingly equilibrated (balanced) mess.
  • Peripheral collisions (glancing blows) are like a sparse dance floor. There's less heat and less interaction.
  • Result: The hotter the collision (more central), the higher the temperature. As the collisions got "grazing" (peripheral), the temperature dropped.

4. The Big Discovery: The "Magic Line"

The most exciting part of the paper is what happened when they compared the results.

Even though the three math tools gave different numbers, they weren't random. They were perfectly linked.

• If you plotted the "Boltzmann temperature" against the "Bose-Einstein temperature," the dots formed a perfect straight line.
• If you plotted the "Tsallis temperature" against the "Bose-Einstein temperature," it was also a straight line.

The Analogy: Imagine you have three different currency converters (Dollars, Euros, Yen). They all give you different numbers for the price of a coffee. But if you know the exchange rate, you can perfectly convert one to the other. The scientists found the "exchange rates" between these different physics models.

5. Why Does This Matter?

In the past, scientists mostly studied huge crashes (like smashing two bowling balls together). This paper is special because it looks at small crashes (marbles and bowling balls).

• The Surprise: Even in these tiny, small-system crashes, the particles act like they are in a hot, boiling soup, not just a random scatter.
• The Conclusion: The scientists propose that we should use the "Bose-Einstein/Fermi-Dirac" numbers as the standard ruler for everyone. If you use the "Tsallis" ruler, you just need to apply a simple math formula to convert it to the standard.

Summary

This paper is like a guidebook for physicists. It says: "Hey, we have different ways to measure the heat of a particle crash. They all give different numbers, but they are all related by simple lines. If you want to compare your results with mine, just use our conversion chart. Also, even in tiny crashes, the universe gets surprisingly hot and organized."

It helps scientists stop arguing about which math tool is "right" and start using them together to understand the fundamental building blocks of our universe.

Technical Summary

1. Problem Statement

In high-energy nuclear physics, temperature is a fundamental parameter used to understand collision mechanisms, particle production, and phase transitions (specifically from hadronic matter to Quark-Gluon Plasma, QGP). However, the concept of "temperature" in these systems is model-dependent. Different statistical distributions (Bose-Einstein, Fermi-Dirac, Boltzmann, Tsallis) applied to the same experimental transverse momentum ($p_T$) spectra yield different values for the effective temperature ($T_{eff}$).

Currently, there is no standardized baseline to compare these temperatures across different models. This lack of a unified reference makes it difficult to interpret physical phenomena, such as kinetic freeze-out conditions or the nature of phase transitions, consistently. The authors aim to address this by systematically comparing $T_{eff}$ values derived from standard quantum statistics and the non-extensive Tsallis distribution within small collision systems (deuteron-gold and proton-proton collisions), where non-equilibrium effects are significant.

2. Methodology

The study utilizes experimental data from the STAR Collaboration at the Relativistic Heavy Ion Collider (RHIC) for collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_{NN}} = 200$ GeV.

• Data Sources:

  • Systems: Deuteron-gold (d+Au) and proton-proton (p+p) collisions.
  • Centrality: d+Au collisions are categorized into central (0–20%), semi-central (20–40%), and peripheral (40–100%) classes.
  • Particles: Identified light charged hadrons: Bosons ($\pi^\pm$, $K^\pm$) and Fermions ($p, \bar{p}$).
  • Kinematic Range: Mid-rapidity ($|y| < 0.5$).

• Distribution Models:
  The authors fit the invariant yield spectra $(1/2\pi p_T) d^2N/dydp_T$ using three distinct approaches:

  1. Standard Distributions:
     • Bose-Einstein (BE) for bosons and Fermi-Dirac (FD) for fermions.
     • Boltzmann distribution (the classical approximation).
     • Implementation: To achieve satisfactory fits for the complex spectra, a three-component multi-source thermal model was employed for BE, FD, and Boltzmann distributions. This accounts for temperature fluctuations from different emission sources.
  2. Tsallis Distribution:
     • A single-component Tsallis distribution incorporating the entropy index $q$.
     • Rationale: The Tsallis distribution effectively smooths the temperature fluctuations inherent in the multi-component standard models, allowing for a good fit with fewer parameters.

• Analysis Strategy:

  • Extract $T_{eff}$ values from all models for each centrality and particle type.
  • Analyze the trends of $T_{eff}$ and the entropy index $q$ as a function of collision centrality.
  • Establish mathematical relationships (linear correlations) between the $T_{eff}$ values derived from different distributions.

3. Key Contributions

• Establishment of a Standard Baseline: The paper proposes using temperatures derived from Bose-Einstein (for bosons) and Fermi-Dirac (for fermions) distributions as the standard baseline for comparing effective temperatures across different statistical models.
• Small System Analysis: Unlike previous studies focusing on large heavy-ion systems (Au+Au, Pb+Pb), this work specifically targets small systems (d+Au, p+p), demonstrating that multi-component thermal models and Tsallis statistics are equally applicable and reveal similar trends in these environments.
• Quantification of Linear Relationships: The study provides precise linear equations linking $T_{Boltzmann}$ and $T_{Tsallis}$ to the standard $T_{BE/FD}$ baselines, allowing for the conversion of temperature values between different statistical frameworks.

4. Key Results

• Systematic Trends in $T_{eff}$:

  • Magnitude Order: For a given particle and centrality, the effective temperatures follow a systematic decreasing trend:
    $$T_{Bose-Einstein} > T_{Boltzmann} > T_{Fermi-Dirac} > T_{Tsallis}$$
    (Note: $T_{Boltzmann}$ underestimates $T_{BE}$ but overestimates $T_{FD}$; $T_{Tsallis}$ is consistently the lowest due to the influence of the entropy index $q$.)
  • Centrality Dependence: $T_{eff}$ values decrease as collision centrality decreases (moving from central to peripheral collisions).
  • Particle Dependence: $T_{eff}$ shows mass dependence and isospin dependence (e.g., differences between $K^+$ and $K^-$ due to absorption rates, and $p$ vs. $\bar{p}$ due to pre-existing protons in nuclei).

• Entropy Index ($q$) Behavior:

  • The entropy index $q$ approaches unity ($q \to 1$) in central collisions, indicating a system closer to thermal equilibrium due to increased multi-scattering.
  • $q$ increases in peripheral collisions and for lighter particles (pions), reflecting stronger non-equilibrium features. Pions, being the lightest hadrons, freeze out earlier in a high-energy density environment, retaining more non-extensive characteristics ($q > 1$).

• Linear Correlations:

  • A perfect linear relationship was observed between $T_{Boltzmann}$ and $T_{BE/FD}$, and between $T_{Tsallis}$ and $T_{BE/FD}$.
  • Example linear fits (for $\pi^+$ in central d+Au):
    • $T_{BE} \approx 0.819 T_{Boltzmann} + 0.009$
    • $T_{BE} \approx 0.628 T_{Tsallis} + 0.007$
  • These correlations hold across different centrality classes and particle types, validating the use of a unified baseline.

• Small vs. Large Systems:

  • Parameters from p+p collisions are found to be very similar to those from peripheral d+Au collisions, suggesting comparable energy deposition and initial-state geometries with minimized hadronic re-scattering.

5. Significance and Implications

• Unified Thermodynamic Framework: By establishing linear relationships between different statistical models, the paper enables a unified interpretation of temperature fluctuations across statistical paradigms (quantum vs. non-extensive statistics).
• Insight into QCD Evolution: The study highlights the coexistence of pronounced non-equilibrium effects (captured by Tsallis $q$) and local thermal equilibrium (captured by multi-component standard distributions) in small systems. This challenges the universality of pure equilibrium assumptions in high-energy collisions.
• Phase Transition Probes: The extracted effective temperatures serve as indicators for the kinetic freeze-out stage. The non-monotonic dependence of $T_{eff}$ on energy/centrality provides indirect evidence for phase transitions and potential critical endpoints in the QCD phase diagram.
• Future Experimental Design: The findings guide the design of precision momentum spectrum measurements at RHIC, the Large Hadron Collider (LHC), and the future Electron-Ion Collider (EIC), suggesting that hybrid statistical frameworks are necessary to decode the evolution of the QCD medium.

In conclusion, this work bridges the gap between traditional thermal models and non-extensive statistics, providing a robust, standardized method for comparing effective temperatures in high-energy collisions, particularly in the complex environment of small collision systems.