Quark-meson coupling model and heavy-ion collision

This paper implements the quark-meson coupling model within the Daejeon Boltzmann-Uehling-Uhlenbeck transport framework to simulate Au+Au collisions at intermediate energies, comparing the resulting maximum densities with those from a conventional quantum hadrodynamics model to highlight differences in predicted nuclear matter properties.

The Gist

Imagine you are trying to understand how a crowd of people behaves when they are squeezed into a tiny, overcrowded room. Do they push back hard? Do they squish together easily? Do they bounce off each other?

This is essentially what physicists are trying to figure out with nuclear matter—the stuff inside the center of an atom (the nucleus). Under normal conditions, this stuff is spread out. But in a heavy-ion collision, scientists smash two heavy atoms (like Gold) together at nearly the speed of light. For a split second, they create a "super-dense" ball of nuclear matter, similar to what exists inside neutron stars.

The problem is, we can't see this dense matter directly. It's like trying to understand the inside of a black box just by listening to the noise it makes when you shake it. To solve this, scientists use computer simulations called transport models. Think of these models as video game engines that predict how the particles will move and interact.

The Two "Rulebooks" for the Simulation

In this paper, the researchers are testing two different "rulebooks" (theories) to see which one better predicts what happens during these atomic crashes.

1. The "Hard Shell" Rulebook (QHD Model):
This is the traditional way of thinking. It treats the proton and neutron (nucleons) as solid, indivisible marbles. When they interact, they just bounce off each other or push against invisible fields, like billiard balls on a pool table. This is the Quantum Hadrodynamics (QHD) model.

2. The "Soft Interior" Rulebook (QMC Model):
This is the new, more sophisticated approach. It realizes that protons and neutrons aren't actually solid marbles; they are made of smaller particles called quarks.

• The Analogy: Imagine a proton isn't a solid marble, but a water balloon filled with three tiny marbles (quarks).
• When you squeeze this balloon (the nuclear matter), the water inside moves, and the shape of the balloon changes. The "Quark-Meson Coupling" (QMC) model accounts for this squishiness. It says the internal structure of the particle changes when it's under pressure, which changes how the whole system behaves.

The Experiment: The Atomic Smash-Up

The researchers took their computer simulation engine (called DJBUU, named after the city of Daejeon where a major particle facility is being built) and ran a simulation of two Gold atoms smashing into each other.

They ran the simulation twice:

1. Once using the "Hard Shell" rulebook (QHD).
2. Once using the "Soft Interior" rulebook (QMC).

What Did They Find?

Here is the surprising twist they discovered:

• The "Stiff" vs. "Soft" Surprise: Usually, if a material is "stiffer" (harder to compress), you expect it to resist being squeezed more, resulting in a lower maximum density. The QMC model predicted a "stiffer" nuclear matter (based on some of its numbers), so the scientists expected it to be harder to squeeze.
• The Reality: However, the QMC simulation actually allowed the atoms to get denser (squeezed tighter) than the traditional model.
• Why? It comes down to the "water balloon" effect. Because the QMC model accounts for the fact that the internal quarks change the "effective mass" (how heavy the particle feels) when squeezed, the particles actually become "heavier" and slower to react. This change in behavior allowed the system to pack itself tighter than the rigid "marble" model predicted.

The Takeaway

Think of it like this:
If you try to pack a suitcase with solid bricks (the old model), you can only fit so many before it's full.
But if you pack it with soft, squishy pillows that change shape when you press them (the new QMC model), you might be able to fit more stuff in, even though the pillows feel "stiffer" in a different way.

Why does this matter?
This helps scientists understand the Equation of State (the rulebook for how matter behaves under extreme pressure). This is crucial for:

1. Neutron Stars: These are giant balls of nuclear matter in space. Knowing how dense they can get helps us understand their size and structure.
2. Future Experiments: As new facilities like RAON in Korea come online, they will smash even stranger atoms together. This research helps physicists know what to look for and which theories are most likely to be correct.

In short, the paper shows that by treating atomic nuclei as flexible bags of quarks rather than rigid marbles, we get a different, and perhaps more accurate, picture of how the universe's densest matter behaves.

Technical Summary

1. Problem Statement

Heavy-ion collision (HIC) experiments are the primary terrestrial method for investigating dense nuclear matter. However, since dense matter cannot be observed directly, theoretical transport models are essential to interpret experimental observables and extract information about the Equation of State (EoS) of nuclear matter.

While the Quantum Hadrodynamics (QHD) model, based on hadronic degrees of freedom, is the standard framework for relativistic transport models, it treats nucleons as point-like particles interacting with meson fields. In contrast, the Quark-Meson Coupling (QMC) model treats nucleons as bags containing three quarks that interact directly with meson mean fields. Despite the QMC model's success in describing nuclear matter, finite nuclei, and neutron stars, its implementation in dynamical transport models for heavy-ion collisions has been limited. This study aims to bridge that gap by implementing the QMC model into a transport framework to compare its predictions against conventional QHD models.

2. Methodology

The authors integrated the QMC model into the Daejeon Boltzmann-Uehling-Uhlenbeck (DJBUU) transport model. The methodology involved three main stages:

• Theoretical Framework Implementation:
  • QHD Model: Used as a baseline, employing the standard relativistic mean-field (RMF) Lagrangian where nucleons interact with $\sigma$, $\omega$, and $\rho$ mesons. Two parameter sets were used: Liu$\rho$ (commonly used in transport models) and NL3 (known for a stiff EoS).
  • QMC Model: Implemented using a hadron-level Lagrangian derived from the quark-level interaction. The key distinction is the effective nucleon mass ($m^*_N$), which includes a quadratic term dependent on the scalar field ($\sigma$) and a parameter $a_N$:
    $$m^*_N = m_N - (g_\sigma \sigma - a_N \frac{g_\sigma^2}{2} \sigma^2)$$
    This term arises from the quark structure of the nucleon, eliminating the need for nonlinear self-interaction terms ($U(\sigma)$) often required in QHD to fit incompressibility.
• Transport Simulation:
  • The DJBUU model solves the relativistic Boltzmann-Uehling-Uhlenbeck equation, incorporating the Vlasov term (mean-field evolution) and the collision term (nucleon-nucleon collisions with Pauli blocking).
  • Simulations were performed for $^{197}\text{Au} + ^{197}\text{Au}$ collisions at a beam energy of $E_{\text{beam}} = 400$ A MeV and an impact parameter $b = 4.7$ fm.
• Benchmarking:
  • Before dynamic simulations, the authors calculated static properties of symmetric nuclear matter and pure neutron matter (pressure, binding energy, symmetry energy, effective mass) to ensure the parameter sets reproduced known nuclear constraints.

3. Key Contributions

• First Implementation of QMC in DJBUU: The paper successfully adapts the QMC model, traditionally used for static nuclear systems, into a dynamical transport model capable of simulating time-evolving heavy-ion collisions.
• Comparative Analysis of Degrees of Freedom: It provides a direct comparison between hadronic degrees of freedom (QHD) and quark degrees of freedom (QMC) within the same transport framework, isolating the effects of the underlying nuclear structure on collision dynamics.
• Parameter Set Verification: The study validates the QMC parameter set (from Ref. [12]) against the widely used Liu$\rho$ and NL3 sets, confirming that QMC satisfies constraints from heavy-ion collision experiments regarding pressure and incompressibility.

4. Results

The simulations yielded several critical findings regarding the time evolution of central baryon density and nuclear matter properties:

• Nuclear Matter Properties (Static):

  • Pressure: Both Liu$\rho$ and QMC parameter sets fall within the experimental constraints derived from heavy-ion collisions. The NL3 set predicts a "stiff" equation of state (EOS) that exceeds these constraints.
  • Effective Mass: The QMC model predicts the largest effective nucleon mass ($m^*/m \approx 0.8$) at saturation density, followed by Liu$\rho$ ($0.75$), while NL3 yields the smallest ($0.6$).
  • Incompressibility ($K_0$) and Symmetry Energy ($S$): The QMC model predicts higher values for $K_0$ (280 MeV) and $S$ (35 MeV) compared to Liu$\rho$ (240 MeV and 30.5 MeV, respectively).

• Heavy-Ion Collision Dynamics (Au+Au):

  • Density Evolution: The central baryon density evolution differed significantly among the models.
    • NL3: Due to its stiff EOS, the density rose rapidly but reached the lowest maximum density and peaked earliest.
    • Liu$\rho$ vs. QMC: These two models showed similar density evolution up to $\approx 8$ fm/c. However, QMC reached a slightly higher maximum density than Liu$\rho$.
  • Interpretation of Maximum Density: The authors noted a counter-intuitive result: despite QMC having a stiffer EOS (higher $K_0$ and $S$), which usually suppresses maximum density, it produced a higher density than Liu$\rho$.
    • Explanation: The authors attribute this to the larger effective mass in the QMC model. A larger effective mass tends to "soften" the EOS dynamically during the collision, allowing for greater compression. This effect of the effective mass ratio ($m^*/m$) dominated over the stiffening effects of $K_0$ and $S$.
  • Directed Flow: The directed flow results for Au+Au collisions using both Liu$\rho$ and QMC were similar and consistent with experimental data.

5. Significance

This work demonstrates that the internal quark structure of nucleons (modeled via QMC) significantly influences the dynamical evolution of heavy-ion collisions, particularly the maximum density achieved.

• Theoretical Insight: It highlights that the effective nucleon mass is a critical parameter in transport models, potentially more influential than bulk incompressibility in determining peak densities in intermediate-energy collisions.
• Future Applications: By validating the QMC model within DJBUU, the authors enable future studies on isospin asymmetry using rare-isotope beams (e.g., at RAON in Korea). This allows for a more fundamental investigation of dense matter properties, bridging the gap between quark-level physics and macroscopic nuclear phenomena.
• Model Robustness: The study confirms that QMC-based simulations are viable alternatives to standard QHD models, offering a different perspective on the Equation of State that is consistent with existing experimental constraints.