Calculation of Particle Pair Correlation Functions with Classical Trajectory Approximation

This paper presents a novel Monte Carlo model using a classical trajectory approximation to calculate particle pair correlation functions in heavy-ion collisions, demonstrating that the method successfully fits experimental data and is highly sensitive to the source's spatio-temporal extent while being largely independent of the temperature parameter, thus advancing femtoscopic interferometry in the Fermi energy domain.

The Gist

Imagine you are trying to figure out the shape and size of a mysterious, invisible balloon that just popped in a dark room. You can't see the balloon itself, but you can see the confetti (particles) flying out in all directions. By looking at how the pieces of confetti land relative to each other, you can guess how big the balloon was and how fast it was expanding when it burst.

This is essentially what physicists do in heavy-ion collisions. They smash atomic nuclei together at high speeds, creating a tiny, super-hot "fireball" of nuclear matter. As this fireball cools down, it spits out particles. The goal is to understand the spatio-temporal evolution of that fireball—basically, its size, shape, and how long it existed before disappearing.

Here is a breakdown of the paper's story, using simple analogies:

1. The Problem: The "Ghost" in the Machine

Scientists have a powerful tool called Femtoscopic Interferometry (or "Femtoscopy"). It's like a super-microscope that uses the wave nature of particles to measure things smaller than an atom.

However, there's a catch. When two particles fly out of the fireball, they don't just travel in straight lines. They are influenced by two things:

• The Ghost (The Source): The remaining "residual" nucleus (the leftover debris of the explosion) still has a gravitational-like pull (actually, it's an electric and nuclear force) that tugs on the particles.
• The Dance (Final State Interactions): The two particles themselves might push or pull on each other (like two magnets) as they fly away.

Previous models often tried to ignore the "Ghost" or treated the "Dance" too simply. It's like trying to predict where two dancers will end up in a room, but ignoring the fact that the floor is sticky and they are holding hands. The results were often messy or inaccurate.

2. The Solution: The "CTA-I" Model

The authors (Sheng Xiao and colleagues) built a new computer simulation called CTA-I (Classical Trajectory Approximation). Think of this as a high-tech video game engine for nuclear physics.

Instead of using complex quantum math that gets messy, they used Classical Trajectories. Imagine you are playing a game of billiards, but the balls are charged and the table is curved.

• The Setup: They create a virtual "fireball" that is in thermal equilibrium (like a pot of boiling water where everything is mixed evenly).
• The Rules: They programmed the rules of the game:
  • The particles are born with random speeds based on the "temperature" of the pot.
  • They are born in a specific shape (a Gaussian cloud, which is like a fuzzy ball).
  • As they fly out, they feel the pull of the leftover nucleus (the "Ghost") and the push/pull of their partner particle.
• The Simulation: The computer runs millions of these "games," tracking the path of every particle pair step-by-step, just like a physics engine in a video game.

3. The Big Discovery: Size Matters, Temperature Doesn't

After running their simulation and comparing it to real data from particle detectors, they found something surprising:

• The Temperature is a Red Herring: Changing the "temperature" of the fireball (how hot the particles are) barely changed the final pattern of the particles. It's like turning up the heat in a kitchen; the smell of the food changes, but the way the smoke swirls in the air stays mostly the same.
• The Size is the Key: Changing the size of the fireball had a huge effect. If the fireball was slightly smaller, the particles interacted much more strongly, creating a very different pattern.

The Metaphor: Imagine throwing two tennis balls from a small tent versus a giant stadium.

• In the small tent, the balls are close together, and the walls (the source) push on them immediately. They interact heavily.
• In the giant stadium, the balls are far apart and the walls are far away. They barely notice each other or the walls.

The paper shows that by looking at how the "tennis balls" (particles) interact, we can measure the size of the "tent" (the source) with incredible precision, even if we don't know exactly how "hot" the tent was.

4. Why This Matters

This new model is a universal translator for nuclear physicists.

• Before: It was hard to separate the effects of the source's size from the effects of the leftover nucleus.
• Now: The CTA-I model handles both the source and the interactions "self-consistently" (meaning it doesn't cheat or ignore parts of the physics).

This allows scientists to finally extract accurate information about the equation of state of nuclear matter (how nuclear matter behaves under pressure). This is crucial for understanding:

• How stars explode (supernovae).
• What happens inside neutron stars.
• The fundamental forces that hold the universe together.

Summary

The authors built a sophisticated virtual billiard table to simulate how particles fly out of a nuclear explosion. They discovered that while the "heat" of the explosion doesn't change the outcome much, the size of the explosion changes everything. Their new tool allows them to measure the size of these microscopic fireballs with high precision, helping us understand the building blocks of the universe.

Technical Summary

Here is a detailed technical summary of the paper "Calculation of Particle Pair Correlation Functions with Classical Trajectory Approximation" by Xiao et al.

1. Problem Statement

Femtoscopic interferometry is a critical tool for probing the spatio-temporal evolution of emission sources in heavy-ion collisions (HIRs), particularly in the Fermi energy domain. However, a significant challenge exists in formulating a self-consistent description that simultaneously accounts for:

• The emission source (assumed to be in thermal equilibrium).
• Final-state interactions (FSI) between the particle pair.
• The influence of the source's own potential field on the particles.

Existing models have limitations:

• CRAB: Neglects the source's potential field influence on final-state pairs.
• Lednicky-Lyuboshits (LL): Treats two-body scattering exactly but neglects the residual nucleus (three-body effects).
• MENEKA: Uses classical trajectories and includes three-body dynamics but relies on specific assumptions (e.g., surface emission, exponential time delay) that may not fully capture the self-consistency between the source geometry and the mean field.

The authors aim to develop a unified framework that treats the source geometry and the potential field experienced by emitted particles in a self-consistent manner to extract precise spatio-temporal information.

2. Methodology: The CTA-I Model

The authors developed a novel Monte Carlo model based on the Classical Trajectory Approximation (CTA-I). The model operates through the following workflow:

A. Theoretical Foundation

• Thermal Equilibrium Source: The emission source is defined as a region of homogeneity in thermal equilibrium. The single-particle momentum distribution follows a Boltzmann form, and the spatial distribution is assumed to be a Gaussian source ($f_x(\vec{r})$).
• Self-Consistent Mean Field: By applying Liouville's theorem and requiring thermal equilibrium, the authors derived a constraint linking the source distribution to the mean field potential $V(r)$.
  • A Gaussian source implies a harmonic oscillator potential near the origin ($r \approx 0$).
  • The total mean field ($V_{MF}$) is constructed as a sum of three components:
    1. Electric Part: A smoothed Coulomb potential derived from a Gaussian charge distribution (to avoid singularities at $r=0$).
    2. Volume Part: A Wood-Saxon potential representing the attractive nuclear force.
    3. Surface Part: The derivative of the Wood-Saxon function, representing surface effects.
  • The parameters of this potential (depth, diffuseness, charge radius) are constrained by the source temperature ($T$) and source size ($\sigma_R$) to ensure self-consistency.

B. Simulation Dynamics

• Initialization: Particles are sampled from the Gaussian source with momenta drawn from the thermal distribution.
• Trajectory Evolution: The motion of correlated particle pairs is simulated using classical mechanics. The equations of motion include:
  • The interaction between the two particles (Coulomb repulsion).
  • The interaction of each particle with the residual source (the self-consistent mean field).
  • A time-stepping algorithm (similar to a predictor-corrector method) is used to integrate the trajectories until particles exit the interaction range.
• Correlation Function Calculation: The correlation function $C(q)$ is calculated as the ratio of the coincidence yield of pairs in the same event to the yield in mixed events, normalized at large relative momenta.

C. Parameterization

The model is controlled by specific parameter sets defining:

• Reaction System: Beam energy, projectile/target properties, and linear momentum transfer.
• Emission Source: Temperature ($k_BT$), Gaussian source size ($\sigma_R$), and emission time distribution (Poisson rate).
• Particle Pair: Charge, mass, and mass numbers of the emitted particles.
• Experimental Filtering: Detector acceptance and resolution are applied to the simulated events to allow direct comparison with data.

3. Key Contributions

1. Self-Consistent Framework: The model uniquely couples the Gaussian source distribution with a mean field potential that satisfies thermal equilibrium constraints, ensuring the source geometry and the potential field are mutually consistent.
2. Three-Body Dynamics: Unlike analytical models that often neglect the residual nucleus, CTA-I explicitly simulates the three-body dynamics (Source + Particle 1 + Particle 2) using classical trajectories.
3. Refined Potential Construction: The derivation of the mean field parameters ($V_0, S_0, d, \gamma_c, \gamma_d$) directly from the source temperature and size provides a rigorous physical basis for the potential shape, transitioning from harmonic near the center to Coulomb at large distances.

4. Results

The model was applied to experimental data from Fermi energy heavy-ion collisions:

• IMF-IMF Correlations (Boron pairs in $^{36}$Ar + $^{197}$Au):
  • The model successfully reproduced the experimental correlation functions.
  • Sensitivity Analysis: The correlation function showed negligible dependence on the temperature parameter ($k_BT$) but was highly sensitive to the source size parameter ($\sigma_R$). A change of only 1 fm in source size resulted in a noticeable change in correlation strength.
• LCP-LCP Correlations (Triton-Triton and $^3$He-$^3$He in $^{86}$Kr + $^{nat}$Pb):
  • The model reproduced the trends for both triton and $^3$He pairs.
  • Similar to IMFs, the correlation function was insensitive to temperature variations but strongly dependent on the source size.
  • The model captured the subtle differences in correlation strength between t-t and $^3$He-$^3$He pairs, attributing this to the stronger Coulomb anti-correlation in the $^3$He case.

5. Significance and Conclusion

• Extraction of Spatio-Temporal Information: The study demonstrates that femtoscopic correlation functions in the Fermi energy domain are primarily a probe of the spatial extent (source size) of the emission source, rather than the thermodynamic temperature (which is better determined via energy spectra).
• Validation of CTA-I: The CTA-I model provides a robust theoretical tool for interpreting experimental data without relying on the simplifying assumptions of previous models (like neglecting the source potential).
• Future Applications: The framework is flexible and can be extended to:
  • Non-spherical source geometries (e.g., disks, rings, rotating systems).
  • Non-equilibrium emission scenarios (Gaussian processes) if the emission rate is high.
  • Isospin effect studies by comparing different particle pairs (e.g., t-t vs. $^3$He-$^3$He).

In summary, the paper presents a significant advancement in femtoscopy by providing a self-consistent, classical trajectory-based model that effectively decouples source geometry from thermal parameters, offering a precise method to constrain the size and dynamics of emission sources in nuclear reactions.