Quarkonium in non-zero isospin chemical potential environment at $T \simeq 0$

Using lattice Non-Relativistic QCD with $N_f=2+1$ dynamical staggered quarks, this study investigates the impact of isospin asymmetry on bottomonium states at near-zero temperature, revealing that the Upsilon mass increases at a high isospin chemical potential of $\mu_I a = 0.106$ while exhibiting non-monotonic behavior at lower values.

The Gist

Imagine the universe is filled with a cosmic soup called Quantum Chromodynamics (QCD). Inside this soup, tiny particles called quarks stick together to form larger particles like protons and neutrons. Usually, these quarks are very picky about their partners, but under extreme conditions—like inside a neutron star or right after the Big Bang—they get squeezed so tightly that they behave differently.

This paper is a report from a team of physicists trying to understand what happens to a very special, heavy "molecule" made of quarks, called Quarkonium (specifically, a pair of bottom quarks), when the soup gets a specific kind of "flavor imbalance."

Here is the breakdown using simple analogies:

1. The Setting: The Cosmic Soup

Think of the universe's matter as a giant dance floor.

• Normal conditions: The dancers (quarks) are paired up evenly.
• High Density: Imagine the dance floor gets so crowded that people are pressed shoulder-to-shoulder. This happens in neutron stars.
• The Problem: Scientists want to study what happens to the dancers when the room is super crowded. However, running computer simulations for this is like trying to solve a math problem where the numbers keep turning into imaginary ghosts (a "complex action" problem). The computers crash or get confused.

2. The Workaround: The "Flavor" Chemical Potential

Since they can't easily simulate the crowded neutron star directly, the scientists decided to simulate a slightly different version of the dance floor.

• Instead of crowding everyone equally, they introduced an Isospin Chemical Potential.
• The Analogy: Imagine the dance floor has two types of dancers: Red Shirts and Blue Shirts. Usually, there's an equal mix. The scientists decided to artificially add more Red Shirts than Blue Shirts. This creates a "flavor imbalance."
• Why do this? This imbalance is mathematically easier for computers to handle than the real neutron star density, but it still teaches us how the "dance" changes when the crowd gets weird.

3. The Subject: The Heavy Couple (Quarkonium)

The scientists are watching a specific couple on the dance floor: the Bottomonium (or Upsilon).

• The Analogy: Imagine a very heavy, slow-moving couple (the bottom quark and anti-quark) holding hands and spinning in the center of the room. Because they are so heavy, they don't get pushed around easily by the crowd. They act like a thermometer or a pressure gauge.
• If the crowd pushes them, they might spin faster, slower, or get heavier. By watching this heavy couple, the scientists can tell what the rest of the crowd is doing.

4. The Experiment: Squeezing the Couple

The team used a supercomputer to simulate this "Red Shirt vs. Blue Shirt" imbalance at a temperature close to absolute zero (very cold). They watched how the heavy couple behaved as they increased the number of "Red Shirts" (the isospin chemical potential).

What they found:

• The "Non-Monotonic" Surprise: Usually, you expect that if you squeeze something harder, it gets heavier or lighter in a straight line. But here, the behavior was wobbly.
  • At first, as they added more "Red Shirts," the heavy couple didn't change much.
  • Then, at a certain point, they got slightly lighter.
  • But then, when the imbalance got really strong (the highest level they tested), the couple suddenly got heavier than they were in an empty room.
• The "Ghost" Effect: To make the simulation work, they had to use a tiny "current" (a nudge) to keep the math stable. They found that the results were very sensitive to how hard they nudged the system. It's like trying to balance a pencil on its tip; a tiny breeze changes the result.

5. The Big Picture: Why Does This Matter?

• Neutron Stars: This helps us understand the "Equation of State" (the rulebook) for neutron stars. If we know how heavy quark couples behave in this weird imbalance, we can guess how matter behaves inside a collapsing star.
• Different Theories: The scientists compared their results to a simpler theory (called SU(2)). In that simpler theory, the heavy couple gets lighter when crowded. But in their more realistic simulation (QCD), the couple gets heavier at high imbalance. This proves that the "real" universe behaves differently than the simplified models, and we need to be careful which model we trust.

Summary

The paper is a preliminary report saying: "We tried to see how a heavy quark couple reacts to a crowded, imbalanced universe. The results are messy and wobbly, but at the highest levels of imbalance, the couple gets heavier. We need more data and better computers to be sure, but this is a crucial step toward understanding the guts of neutron stars."

It's like trying to predict how a giant anchor would behave if you dropped it into a pool of water that was slowly turning into jelly. The anchor gets heavier, but the process isn't smooth—it's a bumpy ride that requires careful observation.

Technical Summary

1. Problem Statement

The study addresses the behavior of heavy quarkonium states (specifically bottomonium) in a QCD medium with non-zero isospin chemical potential ($\mu_I$) at near-zero temperature ($T \simeq 0$).

• Context: Understanding the QCD phase diagram at high baryon density is crucial for interpreting data from the Compressed Baryonic Matter (CBM) experiment at FAIR and for modeling neutron star structures and mergers. However, standard Lattice QCD simulations face the "sign problem" (complex action) at non-zero baryon chemical potential ($\mu_B$), making direct simulations difficult.
• Motivation: QCD with an isospin chemical potential ($\mu_I$) does not suffer from the sign problem and serves as a viable proxy to study high-density QCD dynamics. While previous studies on $SU(2)$ gauge theory showed quarkonium masses decreasing significantly at high "baryon" density, results for $SU(3)$ QCD with $\mu_I$ have been less clear.
• Objective: To investigate how isospin asymmetry affects the mass and correlators of $S$-wave and $P$-wave bottomonium states (specifically $\Upsilon$ and $\eta_b$) in the vacuum and in a pion-condensed medium.

2. Methodology

The authors employed Lattice Non-Relativistic QCD (NRQCD) to simulate bottom quarks on gauge field ensembles generated with dynamical light quarks.

• Lattice Setup:

  • Gauge Ensembles: Generated with $N_f = 2+1$ flavors of dynamical staggered quarks.
  • Parameters: Lattice size $32^3 \times 48$, lattice spacing $a \approx 0.1535$ fm, and pion mass $m_\pi \approx 135$ MeV.
  • Chemical Potential: Simulations were performed at various isospin chemical potentials ($\mu_I a$) ranging from $0.000$ to $0.106$. This range covers the region below and above the pion condensation threshold ($\mu_I a \approx m_\pi a / 2 \approx 0.053$).
  • Current Sources: To handle the Goldstone mode associated with pion condensation, simulations were run with non-zero $(u, d)$ current sources ($\lambda_a$) of magnitudes $0.0010, 0.0018,$ and $0.0036$. Physical observables require extrapolation to the $\lambda_a \to 0$ limit.

• NRQCD Formalism:

  • The bottom quark dynamics were modeled using the $O(v^4)$ NRQCD Lagrangian, including relativistic corrections up to order $1/M_b^3$.
  • The Lagrangian includes terms for kinetic energy, chromo-electric/magnetic interactions, and spin-dependent forces.
  • Propagator: The lattice propagator was formulated as an initial value problem, utilizing a Lepage parameter ($n=2$) to ensure stability for high momentum behavior.
  • Tuning: The bottom quark mass ($M_b$) was tuned using the dispersion relation $aE(\hat{P}^2) = aM + a\hat{P}^2/2M^2$ to match the experimental spin-averaged $1S$ bottomonium mass.

• Analysis Strategy:

  • Constructed $S$-wave ($\Upsilon, \eta_b$) and $P$-wave correlators.
  • Calculated the ratio of correlators at non-zero $\mu_I$ to those at $\mu_I = 0$. This ratio serves as a direct probe of the medium effect, minimizing systematic errors.
  • Analyzed the behavior of these ratios as a function of Euclidean time ($\tau$) and $\mu_I$.

3. Key Contributions

• First Lattice NRQCD Study at $T \simeq 0$: This work provides preliminary lattice results for bottomonium in an isospin medium at zero temperature, a regime previously unexplored with this specific formalism.
• Systematic Variation of Parameters: The study systematically varies both the isospin chemical potential ($\mu_I$) and the symmetry-breaking current source ($\lambda$) to disentangle physical effects from artifacts of the simulation setup.
• Comparison with $SU(2)$ Theory: The paper explicitly contrasts the behavior of quarkonium in $SU(3)$ QCD (isospin medium) with $SU(2)$ gauge theory (baryon medium), highlighting fundamental differences in how heavy quark bound states respond to density.

4. Results

The preliminary results yield several critical observations:

• Non-Monotonic Behavior: The effect of isospin asymmetry on the $\Upsilon$ mass is non-monotonic with respect to $\mu_I$.
  • Below Condensation ($\mu_I a < 0.053$): The correlator ratios are generally close to 1. In some cases (e.g., $\mu_I a = 0.048$ with small $\lambda$), the ratio is slightly $>1$, suggesting a slight mass reduction (lighter state), though statistical errors are large.
  • Above Condensation ($\mu_I a \geq 0.106$): At the highest chemical potential tested ($\mu_I a = 0.106$), the ratio becomes consistently less than 1. Since the correlator $G(\tau) \sim e^{-M\tau}$, a ratio $<1$ implies the exponent (mass) in the medium is larger than in the vacuum. Thus, the $\Upsilon$ mass increases (becomes heavier) in the high-density isospin medium.
• Sensitivity to Current Source ($\lambda$): The correlator ratios are highly sensitive to the strength of the $(u, d)$ current source $\lambda$, particularly at lower $\mu_I$. This indicates that the extrapolation to $\lambda \to 0$ is crucial for accurate physical predictions.
• Magnitude of Effect: The mass shift is small (less than 3% at $\tau/a = 48$) but statistically distinguishable at high $\mu_I$.
• Contrast with $SU(2)$: Unlike $SU(2)$ gauge theory where quarkonium masses drop significantly at high density, $SU(3)$ QCD with isospin shows a mass increase at high $\mu_I$.

5. Significance and Conclusion

• Probe of QCD Dynamics: The results suggest that heavy quarkonium states can serve as a "densimeter" for high-density QCD matter. The distinct behavior (mass increase vs. mass decrease) compared to $SU(2)$ theories implies that the underlying dynamics of $SU(3)$ QCD at high density are qualitatively different, potentially related to the nature of the baryon (three-quark vs. di-quark).
• Implications for Astrophysics: Understanding the Equation of State (EoS) and the behavior of hadronic states at high density is vital for neutron star physics. These findings provide constraints on how heavy quark states behave in the extreme environments found in neutron star cores.
• Future Work: The authors conclude that while the trend of mass increase at high $\mu_I$ is clear, the statistical precision needs improvement. Future work will involve increasing the number of Monte Carlo configurations to reduce errors and performing simulations at even higher $\mu_I$ to confirm the non-monotonic trend and the $\lambda \to 0$ extrapolation.

In summary, this paper establishes that in an isospin-asymmetric QCD medium at zero temperature, the bottomonium mass initially remains stable or slightly decreases but eventually increases significantly as the isospin chemical potential exceeds the pion condensation threshold, offering a unique window into the non-perturbative dynamics of dense QCD matter.