Quark-Meson Coupling Model in Heavy-Ion Collision Simulations

This study integrates the quark-meson coupling (QMC) model into the DJBUU transport code to demonstrate its effectiveness in simulating intermediate-energy heavy-ion collisions, showing comparable performance to traditional QHD for bulk observables while requiring specific adjustments to in-medium $\Delta$ production cross-sections to accurately reproduce pion yields in neutron-rich systems.

The Gist

Imagine you are trying to understand how a giant, chaotic crowd behaves when two massive groups of people crash into each other. In the world of physics, these "people" are atomic nuclei, and the "crash" is a heavy-ion collision. Scientists smash these nuclei together at incredible speeds to recreate the conditions that existed just after the Big Bang or inside the cores of neutron stars.

For decades, physicists have used a set of rules called Quantum Hadrodynamics (QHD) to predict what happens in these crashes. Think of QHD as a map that treats every particle (like a proton or neutron) as a simple, solid billiard ball. It works pretty well, but it ignores the fact that these "billiard balls" are actually made of smaller, squishier things called quarks.

This paper introduces a new, more detailed map called the Quark-Meson Coupling (QMC) model. Instead of treating particles as solid balls, QMC treats them as bags of quarks that can squish and change shape when they get squeezed by other particles.

Here is a breakdown of what the researchers did and what they found, using some everyday analogies:

1. The New Map vs. The Old Map

The researchers took a sophisticated computer simulation code (called DJBUU) that acts like a traffic controller for these nuclear crashes. They swapped out the old "billiard ball" rules (QHD) for the new "squishy bag" rules (QMC).

• The Goal: They wanted to see if this new, more complex map could predict the results of nuclear crashes just as well as the old one.
• The Test: They simulated two types of crashes:
  1. Gold on Gold (Au+Au): A heavy, messy crash at high speed.
  2. Tin on Tin (Sn+Sn): A crash involving different types of "heavy" atoms to see how the balance of neutrons and protons affects the outcome.

2. The "Squish" Factor (Density)

When two nuclei crash, they get squeezed incredibly tight, like a sponge being crushed in a vice.

• The Finding: The new QMC model predicted that the nuclei got squeezed slightly tighter (higher density) than the old QHD model did.
• The Analogy: Imagine two teams of wrestlers pushing against each other. The old model says they push until they hit a wall. The new model says, "Wait, their muscles (quarks) are actually compressing under the pressure, so they can get even closer together."
• Why it matters: Even though the QMC model predicts a "stiffer" material (it resists squishing more), the internal compression of the quarks allowed the nuclei to reach a higher density. It's a bit counter-intuitive, like a spring that is harder to compress but actually gets shorter when you push it because of how its coils interact.

3. The Flow of Traffic (Directed and Transverse Flow)

After the crash, the debris flies out in specific patterns. Physicists measure this "flow" to understand the pressure inside the crash.

• The Result: Both the old model (QHD) and the new model (QMC) predicted the flow patterns almost identically. Both matched the real-world data collected by experiments like FOPI.
• The Takeaway: The new "squishy bag" model works just as well as the old "billiard ball" model for predicting how the debris flies out. This proves the new model is a valid tool for the job.

4. The "Pion" Puzzle (The Tricky Part)

The most interesting part of the paper involves pions (tiny particles created during the crash). The amount of pions produced depends heavily on a temporary particle called the Delta baryon (a heavy, unstable cousin of the proton).

• The Problem: When the researchers used the new QMC model with the exact same rules for creating pions as the old model, it produced too few pions.
• The Analogy: Imagine you have a recipe for baking cookies (pions). The old recipe works perfectly. You switch to a new oven (QMC) that bakes hotter and faster. If you use the exact same recipe, you might burn the cookies or end up with fewer of them because the oven reacts differently to the ingredients.
• The Fix: The researchers realized that because the QMC model treats the particles as "squishy bags," the rules for how they interact inside the dense crowd need a slight adjustment. They tweaked a "suppression factor" (a dial that controls how much the dense environment slows down the creation of these particles).
• The Result: After turning that dial slightly (changing a number from 2.5 to 2.2), the new model predicted the number of pions and their ratios perfectly, matching the real-world data from the SπRIT experiment.

5. Why This Matters

This paper is a "proof of concept." It shows that we can successfully upgrade our understanding of nuclear physics from "solid balls" to "bags of quarks" without breaking our computer simulations.

• The Big Picture: By successfully integrating the QMC model, scientists now have a better tool to study the Equation of State (EOS) of nuclear matter. This is essentially the "rulebook" for how matter behaves under extreme pressure.
• Future Applications: This is crucial for understanding neutron stars (which are basically giant atomic nuclei) and supernovas. If we can model how quarks behave inside these stars, we can better predict how they collapse, how they spin, and what happens when they crash into each other (creating gravitational waves).

Summary

The researchers took a complex new theory (QMC) that looks at the tiny quarks inside atoms, plugged it into a heavy-ion collision simulator, and found that:

1. It predicts nuclei get squeezed slightly tighter than we thought.
2. It predicts the "traffic flow" of the crash just as well as the old theory.
3. It needed a tiny tweak to its "recipe" for creating pions to match reality.

The Bottom Line: The "squishy bag" model works! It's a more fundamental way of looking at the universe, and now we have the software tools to use it for simulating the most violent events in the cosmos.

Technical Summary

1. Problem Statement

Understanding the Equation of State (EOS) of dense nuclear matter is critical for nuclear physics and astrophysics (e.g., neutron stars, supernovae). While the Relativistic Mean Field (RMF) theory is a standard framework for describing nuclear interactions, most transport models used in heavy-ion collision simulations rely on Quantum Hadrodynamics (QHD). In QHD, nucleons are treated as point-like particles interacting via meson fields.

However, the Quark-Meson Coupling (QMC) model offers a more fundamental description by treating baryons as clusters of confined quarks (based on the MIT bag model) that interact directly with external mean fields. While QMC has been successfully applied to static nuclear systems (infinite matter, neutron stars, finite nuclei), its applicability to non-equilibrium, intermediate-energy heavy-ion collisions had remained largely unexplored. The primary challenge was to integrate the QMC model into a dynamic transport code to determine if it could reproduce experimental observables comparable to traditional QHD models.

2. Methodology

The authors implemented the QMC model into the DaeJeon Boltzmann-Uehling-Uhlenbeck (DJBUU) transport code, a covariant semi-classical transport model.

• Model Implementation:

  • DJBUU+QHD: The baseline model using the standard QHD framework where nucleons are point-like with scalar ($\sigma$) and vector ($\omega, \rho$) meson couplings.
  • DJBUU+QMC: The modified model where the mean field is derived from quark degrees of freedom. In this framework, the effective mass of a baryon ($m^*_B$) includes a quadratic term in the scalar field ($g_\sigma \sigma$) derived from quark-level calculations:
    $$m^*_B = m_B - g_\sigma \sigma + a_B (g_\sigma \sigma)^2$$
    Here, $a_B$ is a baryon-species-dependent coefficient (absent in QHD) derived from the bag model stability condition.
  • Two Schemes for QMC:
    1. QMC$_{iter.}$: Solves the meson field equations self-consistently using iterative calculations including the $a_B$ coefficients.
    2. QMC$_{param.}$: Uses density-dependent parameterizations for the scalar and vector potentials fitted to baryon density to reduce computational cost.

• Simulation Setup:

  • System 1 (Flow Observables): $^{197}\text{Au} + ^{197}\text{Au}$ collisions at $E_{beam} = 400$ A MeV.
    • Goal: Compare transverse and directed flow ($v_1$) with FOPI experimental data.
  • System 2 (Pion Production): $^{132}\text{Sn} + ^{124}\text{Sn}$ (neutron-rich, $N/Z=1.56$) and $^{108}\text{Sn} + ^{112}\text{Sn}$ (less neutron-rich, $N/Z=1.2$) at $E_{beam} = 270$ A MeV.
    • Goal: Investigate pion multiplicities ($\pi^-, \pi^+$) and their ratios (Single Ratio and Double Ratio) inspired by S$\pi$RIT experiments.
  • In-Medium Effects: The production of $\Delta$ resonances (the primary source of pions) was modeled using a density-dependent suppression factor in the cross-section: $\sigma^* = \sigma \times \exp(-C \rho_B/\rho_0)$.

3. Key Contributions

• First Integration: This work presents the first implementation of the QMC model within a covariant transport code (DJBUU) for intermediate-energy heavy-ion collisions.
• Benchmarking: It provides a rigorous benchmark comparing QMC against the standard QHD model and experimental data for both collective flow and particle production.
• Parameter Sensitivity: The study highlights the sensitivity of pion production to the in-medium modification of the $\Delta$-production cross-section when switching from QHD to QMC, specifically regarding the effective mass scaling.

4. Key Results

A. Au+Au Collisions (400 A MeV)

• Central Density: The QMC models (both iterative and parametric) predicted a higher maximum central baryon density (~3–6% higher) than QHD, despite QMC having a stiffer incompressibility ($K_0 = 280$ MeV vs. 240 MeV). The authors attribute this counter-intuitive result to the larger Dirac effective mass in QMC ($m^*/m \approx 0.8$ vs. 0.75), which induces a stronger softening effect on the EOS.
• Collective Flow:
  • Transverse Flow: QMC results matched QHD and experimental data well. The parametric version showed a slightly larger slope due to a weaker attractive force.
  • Directed Flow ($v_1$): Both QHD and QMC successfully reproduced the experimental directed flow of free protons, validating the QMC model's ability to describe collective dynamics.

B. Sn+Sn Collisions (270 A MeV) & Pion Production

• Initial Discrepancy: When using the same in-medium suppression parameters (specifically the coefficient $C=2.5$) tuned for QHD, the QMC model produced fewer pions and $\Delta$ baryons than QHD, despite the higher central density. This was caused by the strong density suppression factor ($e^{-C\rho/\rho_0}$) acting on the QMC system, which has a different effective mass profile.
• Refitting the Suppression: The authors found that the in-medium modification tuned for QHD is not directly transferable to QMC. By reducing the density suppression coefficient from $C=2.5$ to $C=2.2$, the QMC model successfully reproduced the experimental pion yields and ratios.
  • This adjustment aligns with the effective-mass scaling approach ($\sigma^*/\sigma_{free} \approx (m^*/m)^2$), where the larger effective mass in QMC ($0.8$) suggests a weaker suppression of the cross-section compared to QHD ($0.75$).
• Double Ratio (DR): The QMC model (with $C=2.2$) reproduced the double pion ratio (DR) as well as QHD, demonstrating its capability to handle isospin-dependent observables.
• Iterative vs. Parametric: The self-consistent iterative approach (QMC$_{iter.}$) showed slightly better agreement with experimental data than the parametric approach, emphasizing the importance of self-consistent meson field treatment.

5. Significance and Outlook

• Validation: The study confirms that the QMC model, which incorporates quark degrees of freedom, is a viable and effective framework for simulating intermediate-energy heavy-ion collisions.
• Quark Effects: It demonstrates that quark-level effects (via the bag model and effective mass modifications) significantly influence the EOS and particle production in non-equilibrium environments, challenging the assumption that such effects are negligible at these energies.
• Future Applications: The successful integration of QMC into DJBUU opens the door for future investigations into:
  • In-medium baryon magnetic moments.
  • Short-range quark-quark correlations.
  • The role of quark degrees of freedom in the formation of hypernuclei and $\Delta$-matter in dense environments.

In conclusion, the paper establishes a robust bridge between quark-level nuclear models and macroscopic transport simulations, providing a new tool to explore the structure of dense nuclear matter under extreme conditions.