Dilepton Correlations from Heavy Flavor Decays

✨

This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: Hunting for Ghosts in a Storm

Imagine you are trying to study the weather inside a massive, chaotic thunderstorm (a heavy-ion collision). You want to measure the "thermal radiation" (the heat and light) coming from the very center of the storm to understand how it formed. This is the goal of physicists studying the Quark-Gluon Plasma (QGP), a super-hot soup of particles that existed right after the Big Bang.

However, there's a problem. The storm is so loud and bright that the signal you are looking for is drowned out by other, louder noises. In this case, the "noise" comes from heavy particles called Charm and Bottom quarks. When these heavy particles decay, they shoot out pairs of electrons (or muons), creating a "dilepton" signal that looks very similar to the thermal signal you want to study.

The Goal of this Paper:
The authors want to figure out exactly how these "noise" particles (the heavy flavor decays) behave. If we understand their behavior perfectly, we can subtract them from the data and finally see the "ghost" (the thermal signal) we are looking for.

The Analogy: The "Back-to-Back" Dance vs. The "Spin"

To understand the physics, let's use a dance analogy.

1. The Heavy Quark Dance (The Parents)
When two heavy quarks (Charm or Bottom) are created in a collision, they are like dance partners who are forced to spin away from each other.

At low energy: They usually spin directly away from each other, back-to-back (180 degrees apart).
At high energy: Sometimes, they get hit by a third dancer (a gluon) and end up spinning in the same direction or at a different angle.

2. The Decay (The Children)
These heavy quarks don't last long. They quickly "decay" (break apart) into lighter particles, including electrons.

The Question: If the parents were dancing back-to-back, do the children (the electrons) also fly off back-to-back? Or does the act of breaking apart scramble their directions so much that we can't tell how the parents were dancing?

3. The "kT Broadening" (The Wind)
In the real world, there is also "wind" (called kT broadening). This is a random jostling that happens to particles as they move through the dense medium.

The Fear: Physicists worried that this "wind" might blow the heavy quarks around so much that they lose their original dance formation. If the wind is strong, the back-to-back pattern disappears, and the particles look like they are flying in random directions.

What the Authors Did

The authors used a sophisticated computer program (called HVQMNR) to simulate these collisions at two different energy levels:

RHIC (USA): Like a medium-sized explosion (200 GeV).
LHC (Europe): Like a massive, high-speed collision (13 TeV).

They simulated the entire process:

Creating the heavy quarks.
Letting them "dance" (move and interact).
Watching them decay into electron pairs.
Checking the angle between the two electrons.

The Surprising Findings

1. The "Memory" is Still There
The authors found that even though the heavy quarks decay into lighter particles, the children (electrons) still remember how the parents were dancing.

If the parents were created back-to-back, the electrons tend to be found back-to-back.
If the parents were created at high speeds, the electrons can end up closer together.
The Metaphor: Even if a parent spins a child around before letting go, the child still flies off in a direction that hints at the parent's original spin. The "dance memory" survives the decay.

2. The "Wind" (kT Broadening) Matters Less Than We Thought
Previously, scientists thought the random "wind" (kT broadening) would completely scramble the angles, making it impossible to tell the difference between different types of collisions.

The Discovery: The authors found that the decay process itself acts like a filter. Because the decay happens in a specific way, it actually reduces the sensitivity to the wind. The "memory" of the original dance is surprisingly robust. The electrons don't get as scrambled as the heavy quarks themselves might.

3. Energy Changes the Story

At lower energies (RHIC): The signal is dominated by Charm quarks. The electrons are mostly found back-to-back.
At higher energies (LHC): As the energy increases, Bottom quarks become more important. Interestingly, at very high energies, the electrons from Bottom quarks start to cluster together (angle near 0) rather than flying apart. This is because high-energy collisions often produce a heavy quark pair shooting off alongside a very fast "third dancer" (a light parton), pushing them in the same direction.

Why Does This Matter?

Think of the heavy-ion collision data as a noisy radio station.

The Thermal Dileptons are the faint, beautiful music you want to hear (the signature of the early universe).
The Heavy Flavor Decays are the static and commercials drowning out the music.

This paper is like a manual on how to tune the radio. By understanding exactly how the "static" (the heavy flavor electrons) behaves—knowing that they keep their "dance memory" and aren't totally scrambled by the wind—physicists can build a better filter.

The Conclusion:
Now that we know how to predict the "noise" (the heavy flavor decays) so accurately, we can subtract it from the data. This leaves us with a much clearer view of the "music" (the thermal radiation), allowing us to finally study the earliest moments of the universe's creation with much greater precision.

Summary in One Sentence

The paper proves that even after heavy particles break apart, their children (electrons) still remember their parents' dance moves, allowing scientists to filter out background noise and finally hear the faint signal of the early universe.

1. Problem Statement and Motivation

The study addresses the challenge of isolating thermal dilepton production (a signature of the Quark-Gluon Plasma's early-time dynamics) from the overwhelming background of non-prompt dileptons produced by the semileptonic decays of heavy-flavor hadrons (charm and bottom) in relativistic heavy-ion collisions.

The Background: In the low-to-intermediate invariant mass region ( $1 < m_{ll} < 3$ GeV), heavy-flavor decays dominate the dilepton continuum. While like-sign subtraction can remove uncorrelated pairs, correlated decays from $c\bar{c}$ and $b\bar{b}$ pairs remain an irreducible background.
The Gap: Previous studies established the azimuthal correlations of heavy-flavor hadrons in $p+p$ collisions. However, it was unclear whether the kinematic correlations (specifically azimuthal separation, $\Delta\phi$ ) of the parent heavy quarks survive the semileptonic decay process to the final-state leptons.
The Goal: To extend previous hadron-level correlation studies to the lepton level, establishing a baseline in $p+p$ collisions to distinguish heavy-flavor decay signals from thermal radiation in future heavy-ion ( $A+A$ ) studies.

2. Methodology

The authors employed a Next-to-Leading Order (NLO) perturbative QCD framework using the exclusive HVQMNR code to calculate heavy quark pair production and subsequent decay.

Production Model:
- Calculations were performed at NLO, including real and virtual corrections.
- Fragmentation: Heavy quarks were hadronized using the Peterson fragmentation function.
- $k_T$ Broadening: Intrinsic transverse momentum broadening ( $k_T$ ) was incorporated via a Gaussian distribution applied to the final state. The width $\langle k_T^2 \rangle$ was parameterized as energy-dependent, scaling with the heavy quark mass ( $n=12$ for charm, $n=3$ for bottom).
Decay Simulation:
- The code simulated semileptonic decays ( $D \to lX$ and $B \to lX$ ).
- Like-Sign Subtraction: A critical component of the methodology was the explicit calculation and subtraction of like-sign lepton pairs arising from $b\bar{b}$ decays (e.g., cascade decays like $B^- \to D^0 e^- \nu_e \to K^- e^+ \nu_e e^-$ ). This isolates the correlated opposite-sign signal.
- Branching ratios were set based on experimental data (approx. 9.6% for $D \to eX$ , and a weighted average of ~10% for $B \to eX$ including primary and secondary decays).
Experimental Acceptance:
- RHIC ( $\sqrt{s} = 200$ GeV): Calculations were restricted to the PHENIX detector acceptance (central electrons $|\eta| < 0.35$ , forward/backward muons $1.2 < |\eta| < 2.4$ ).
- LHC ( $\sqrt{s} = 13$ TeV): Calculations were restricted to the ALICE central electron acceptance ( $|\eta| < 0.8$ ).

3. Key Contributions

Extension to Lepton Level: This is the first comprehensive NLO study extending heavy-flavor azimuthal correlations from the hadron level to the dilepton level, accounting for the specific kinematics of semileptonic decays.
Decay-Induced Decorrelation: The paper quantifies how the decay process acts as a "smearing" mechanism. While heavy quarks produced back-to-back ( $\Delta\phi \sim \pi$ ) at low $p_T$ lose some correlation due to decay kinematics, the study demonstrates that the parent correlations are not entirely erased.
Flavor Separation via Kinematics: The study provides a method to distinguish charm and bottom contributions based on the transverse momentum ( $p_T$ ) of the lepton pair and the resulting $\Delta\phi$ distribution shape.
Sensitivity Analysis: The work rigorously tests the sensitivity of the dilepton correlations to the $k_T$ broadening parameter ( $\Delta$ ), a crucial variable for modeling nuclear effects in $p+A$ and $A+A$ collisions.

4. Key Results

A. Azimuthal Correlations ( $\Delta\phi$ )

Low $p_T$ Regime: At low lepton pair $p_T$ , the distribution is dominated by charm decays and peaks at $\Delta\phi \sim \pi$ (back-to-back), though the peak is broader than at the quark level due to decay kinematics.
High $p_T$ Regime: As lepton pair $p_T$ $p_{T}$ increases, the distribution undergoes a significant transformation:
- The peak shifts from $\Delta\phi \sim \pi$ toward $\Delta\phi \sim 0$ .
- Bottom dominance: The contribution from bottom decays becomes dominant, particularly at small $\Delta\phi$ . This is attributed to the production of $b\bar{b}$ pairs against a hard parton (NLO $2 \to 3$ process), which is more probable at high energies.
Detector Effects: In PHENIX, the acceptance geometry (separate arms) forces specific angular constraints. For example, pairs detected in the same arm are restricted to $\theta < \pi/2$ , enhancing the apparent back-to-back nature of charm decays in specific configurations.

B. Mass Distributions

Charm vs. Bottom: Charm decays dominate the low-mass region ( $1.1 < m_{ee} < 2.4$ GeV), while bottom decays dominate the higher mass region ( $> 3$ GeV).
$p_T$ Dependence: At higher $p_T$ , the bottom contribution increases significantly even in the intermediate mass region, altering the shape of the mass spectrum.

C. Sensitivity to $k_T$ Broadening

Reduced Sensitivity: A major finding is that the decay process significantly reduces the sensitivity of the dilepton signal to $k_T$ broadening compared to the parent hadron signal.
Mechanism: The isotropic nature of the decay in the heavy quark's rest frame, combined with the boost to the lab frame, washes out the initial state $k_T$ effects. Consequently, the dependence on the broadening parameter $\Delta$ is much weaker for dileptons than for reconstructed $D$ or $B$ mesons.
Implication: This suggests that using dilepton correlations to probe cold nuclear matter effects (like $k_T$ broadening in $p+A$ ) is less effective than using hadronic correlations, but the remaining memory of the parent correlation is still valuable.

D. Comparison with Data

PHENIX: The NLO calculations reproduce the trends of PHENIX muon and electron-muon correlation data well. The calculated cross sections for $c\bar{c}$ and $b\bar{b}$ are consistent with PHENIX extractions within uncertainties, though the HVQMNR central values are slightly higher for charm and lower for bottom compared to some PYTHIA/POWHEG tunes used by PHENIX.
ALICE: Predictions for ALICE at 13 TeV show a clear transition in the $\Delta\phi$ distribution as $p_T$ increases, shifting from a charm-dominated back-to-back peak to a bottom-dominated peak near zero.

5. Significance and Conclusions

Baseline for Heavy-Ion Physics: The paper establishes a robust $p+p$ baseline for heavy-flavor dilepton correlations. This is essential for interpreting data from $A+A$ collisions, where the goal is to extract the thermal dilepton signal.
Thermal Signal Isolation: The authors conclude that to study thermal dileptons in heavy-ion collisions, one must separate them from heavy-flavor decays. This can be achieved by:
1. Utilizing impact parameter or decay length information.
2. Focusing on $e\mu$ pairs, which cannot be produced via electromagnetic processes (unlike $e^+e^-$ or $\mu^+\mu^-$ ), thereby reducing the background from light-flavor resonances and Drell-Yan processes.
Memory of Production: Despite the decorrelation introduced by the decay process, the lepton pairs retain a "memory" of the parent heavy quark correlations. This allows for the study of heavy-flavor production mechanisms even in the presence of decay kinematics.
Limitations: The study highlights that while $k_T$ broadening effects are reduced in dileptons, they are not eliminated. However, the dominant uncertainty in the signal shape comes from the interplay between $p_T$ selection and the changing dominance of charm vs. bottom contributions.

In summary, this work provides a critical theoretical tool for the heavy-ion community, demonstrating that while heavy-flavor decays complicate the dilepton spectrum, their specific kinematic signatures (especially the $p_T$ -dependent shift in $\Delta\phi$ ) can be modeled and subtracted to reveal the underlying thermal physics of the QGP.