$D^0$--$D_s^+$ elliptic-flow splitting from sequential hadronization in O--O collisions at $\sqrt{s_{NN}} = 5.36$ TeV

This paper predicts that sequential hadronization, where $D_s^+$ mesons form later than $D^0$ mesons, reproduces the observed elliptic-flow splitting in O--O collisions at $\sqrt{s_{NN}} = 5.36$ TeV and establishes this splitting as a universal chronometer for the quark-gluon plasma's hadronization timeline.

The Gist

The Big Picture: A Race in a Hot Soup

Imagine two heavy-ion collisions (like smashing two heavy atoms together) as creating a tiny, incredibly hot drop of "soup" called the Quark-Gluon Plasma (QGP). This soup exists for a split second before cooling down and turning back into normal particles (hadrons).

Inside this soup, there are heavy "race cars" called charm quarks. As the soup expands and cools, these race cars eventually stop and combine with other particles to form new vehicles:

1. $D^0$ mesons (made of a charm quark and a light quark).
2. $D^+_s$ mesons (made of a charm quark and a strange quark).

The scientists in this paper are trying to figure out when these two types of vehicles are built. Do they get built at the exact same time, or does one get built before the other?

The Mystery: The "Flow" Split

When the soup expands, it doesn't just get bigger; it stretches in a specific oval shape. The particles inside start to flow along this oval. Physicists measure this flow as elliptic flow ($v_2$).

• The Observation: Recent data from the ALICE experiment showed something strange. In the middle of the race, the $D^0$ mesons were flowing more strongly than the $D^+_s$ mesons.
• The Problem: Most standard theories said they should be built at the same time. If they are built at the same time, the physics of how they combine suggests the $D^+_s$ should actually flow more than the $D^0$. This was a contradiction.

The Solution: A "Staggered" Construction Site

The authors propose a new idea: Sequential Hadronization. Think of it like a construction site with two different deadlines.

1. The Early Finishers ($D^+_s$): Because the $D^+_s$ meson is very tightly bound (like a strong magnet), it can form while the soup is still very hot (at a temperature of $1.2 T_c$). It gets built early and leaves the construction site immediately.
2. The Late Finishers ($D^0$): The $D^0$ meson is less tightly bound. It has to wait until the soup cools down a bit more (to temperature $T_c$) before it can be built.

The Analogy:
Imagine a group of runners (the charm quarks) running on a track that is slowly shrinking.

• The $D^+_s$ runners are told to stop and get on a bus at 10:00 AM. They stop running and get on the bus while the track is still wide.
• The $D^0$ runners are told to keep running until 10:15 AM. They stay on the track for those extra 15 minutes.
• Because the track is shrinking and twisting, the runners who stay longer (the $D^0$) get pushed around more by the crowd and end up with a more "twisted" path (higher flow) by the time they finally get on their bus.

This explains why the $D^0$ has more flow than the $D^+_s$: the $D^0$ had more time to get swept up in the chaos of the expanding soup.

Testing the Theory: Small vs. Big Collisions

The authors tested this idea in two different scenarios:

1. Pb-Pb Collisions (Big System): Smashing two Lead nuclei. This creates a large, long-lasting soup.
2. O-O Collisions (Small System): Smashing two Oxygen nuclei. This creates a tiny, short-lived soup (like a spark that dies out quickly).

The Findings:

• In the Big System (Lead): The "time gap" between the two construction deadlines is long (about 2–3 femtoseconds). The $D^0$ runners have plenty of time to get swept up. The difference in flow is large.
• In the Small System (Oxygen): The soup disappears so fast that the "time gap" is crushed. The $D^0$ runners barely have time to run before the soup vanishes.
• The Result: Even in the tiny Oxygen collision, the $D^0$ still flows more than the $D^+_s$, but the difference is much smaller. This matches the new preliminary data from the ALICE experiment perfectly.

If the "Simultaneous" theory (everyone builds at the same time) were true, the Oxygen data would look completely different, and the $D^+_s$ would flow more. Since the data matches the "Staggered" theory, the staggered theory is likely correct.

The "Chronometer" Discovery

The most exciting part of the paper is a discovery about timekeeping.

The authors found a universal rule: The difference in flow between the two particles is directly linked to how long the soup exists between the two construction deadlines.

• The Analogy: Think of the flow difference as a clock.
  • If the soup lasts a long time, the clock shows a big number (big flow difference).
  • If the soup lasts a short time, the clock shows a small number (small flow difference).

They tested this across nine different collision setups (from small Oxygen to large Lead). No matter the size of the collision or the shape of the initial smash, all the data points fell onto a single straight line.

Conclusion:
The difference in how $D^0$ and $D^+_s$ particles flow acts as a "Hadronization Chronometer" (a clock for particle formation). It allows scientists to measure exactly how long the "late stage" of the quark-gluon plasma lasts, simply by looking at the difference in flow between these two specific particles.

Summary

1. The Problem: Experiments showed $D^0$ particles flow more than $D^+_s$ particles, which old theories couldn't explain.
2. The Fix: The authors suggest $D^+_s$ forms early (hot soup) and $D^0$ forms late (cool soup). The $D^0$ gets more flow because it stays in the soup longer.
3. The Proof: This theory works perfectly for both big (Lead) and small (Oxygen) collisions, matching new experimental data.
4. The Takeaway: The difference in flow between these particles is a universal "clock" that tells us how long the hot soup lasts before turning into normal matter.

Technical Summary

Technical Summary: D0–D+s elliptic-flow splitting from sequential hadronization in O–O collisions at √sNN = 5.36 TeV

Problem Statement
Recent preliminary data from the ALICE Collaboration in Pb–Pb collisions at √sNN = 5.02 TeV revealed a counter-intuitive ordering in the elliptic flow ($v_2$) of open-charm mesons: $v_2(D^0) > v_2(D_s^+)$ at intermediate transverse momentum ($p_T$). This observation contradicts standard phenomenological models based on simultaneous hadronization, where all charmed hadrons form on a single freeze-out hypersurface at the critical temperature ($T_c$). In such simultaneous scenarios, coalescence kinematics—driven by the distinct thermal momentum distributions of strange versus light quarks—predict the opposite ordering ($v_2(D_s^+) > v_2(D^0)$).

The central problem addressed in this work is whether the "sequential hadronization" paradigm, previously proposed to explain Pb–Pb data, remains viable in smaller, shorter-lived collision systems like Oxygen–Oxygen (O–O). Specifically, the authors investigate if the time separation between the formation of $D_s^+$ (predicted at $T \simeq 1.2 T_c$) and $D^0$ (at $T_c$) survives the compressed space-time evolution of O–O collisions at $\sqrt{s_{NN}} = 5.36$ TeV, and whether this mechanism can reproduce the observed $v_2$ ordering and yield ratios in this new system.

Methodology
The authors employ a multi-stage theoretical framework integrating initial state generation, hydrodynamic expansion, heavy-quark transport, and a hybrid hadronization scheme:

1. Initial Conditions and Hydrodynamics: The bulk medium evolution is simulated using the TRENTo model for initial geometry (incorporating $\alpha$-clustering structures for O–O via PGCM configurations) coupled with (3+1)D viscous hydrodynamics (CLVisc). The equation of state is based on HotQCD2014 lattice results, with temperature-dependent shear and bulk viscosities.
2. Heavy-Quark Transport: Charm quark evolution is modeled via a Langevin approach that includes both quasi-elastic scattering and medium-induced radiative energy loss. The framework incorporates cold nuclear matter effects (EPPS21 nPDFs) and tunes transport coefficients ($2\pi T D_s$ and $\hat{q}_0$) based on previous Pb–Pb constraints, with specific uncertainty bands applied for the O–O system.
3. Sequential Hadronization: The core of the study is a hybrid coalescence-plus-fragmentation model.
   • Sequential Scenario: $D_s^+$ mesons form on an earlier hypersurface at $T \simeq 1.2 T_c$. If a charm quark fails to coalesce at this stage, it continues to propagate through the QGP until $T_c$, where it may form $D^0$ or other hadrons. Remaining quarks hadronize via fragmentation.
   • Simultaneous Baseline: Both $D_s^+$ and $D^0$ are forced to form at the single $T_c$ hypersurface.
4. Hadronic Rescattering: Post-hadronization, formed mesons undergo rescattering in the hadronic phase until kinetic freeze-out at 137 MeV.

Key Contributions and Results

• Reproduction of $v_2$ Ordering in O–O: The study presents predictions for $p_T$-differential $v_2$ in 0–20% central O–O collisions. The sequential scenario successfully reproduces the experimentally observed ordering $v_2(D^0) > v_2(D_s^+)$ at intermediate $p_T$. In contrast, the simultaneous scenario yields the inverted ordering ($v_2(D_s^+) > v_2(D^0)$), failing to match preliminary ALICE data. This confirms that the time separation between formation hypersurfaces is the essential dynamical ingredient for describing heavy-flavor flow in small systems.
• System-Size Dependence of Flow Splitting: The magnitude of the flow splitting $\Delta v_2 \equiv v_2(D^0) - v_2(D_s^+)$ is found to be substantially smaller in O–O (0–20%) than in Pb–Pb (30–50%). The authors attribute this to the competition between the initial geometric eccentricity and the duration of the $1.2 T_c \to T_c$ evolution window. In O–O, the fireball lifetime is significantly shorter, compressing the time available for $D^0$-parent charm quarks to accumulate additional momentum anisotropy before hadronization.
• Yield Ratio Predictions: The sequential scenario predicts an enhanced $D_s^+/D^0$ yield ratio at low $p_T$ compared to the simultaneous case, driven by the higher strange-quark phase-space density at the earlier $1.2 T_c$ hypersurface. The ratio is systematically higher in Pb–Pb than in O–O across the full $p_T$ range, consistent with greater strangeness abundance and longer fireball lifetimes in larger systems.
• Partonic Origin and Universal Scaling: By decomposing the flow splitting, the authors identify that the reversal of the $v_2$ ordering is kinematic, driven by the late-stage partonic flow accumulated by $D^0$ parents during the $1.2 T_c \to T_c$ interval. A systematic scan across nine collision configurations (varying centrality in O–O and Pb–Pb) reveals a universal linear scaling between the final hadronic splitting $\Delta v_2$ and the partonic flow increment $\Delta v_{2,c}$ accumulated during this specific temperature window.
  • This scaling holds regardless of system size, initial eccentricity, or collision energy (including a projection for FCC energies at 39 TeV).
  • The data points for all configurations collapse onto a single line, indicating that $\Delta v_2$ is dictated entirely by the partonic anisotropy accumulated during the sequential window, rather than the global collision history.

Significance
The paper establishes the $D^0$–$D_s^+$ flow splitting as a robust "hadronization chronometer" of the Quark-Gluon Plasma (QGP). The study demonstrates that:

1. The sequential hadronization mechanism is not limited to large systems (Pb–Pb) but persists in small systems (O–O), providing a consistent explanation for the $v_2$ ordering across different collision sizes.
2. The magnitude of the splitting is a direct probe of the late-stage QGP evolution time ($1.2 T_c \to T_c$). The smaller splitting in O–O is a natural consequence of the compressed evolution window, not a failure of the mechanism.
3. The observed universal linear scaling between hadronic splitting and partonic flow increment offers a predictive framework. Future measurements at varying energies or with different ion species are expected to follow this universal trend, allowing for a quantitative probe of the space-time structure of heavy-flavor hadronization and the dynamics of the late-stage QGP.