quantum chromodynamics quark gluon plasma asymptotic freedom lattice QCD

# Quantum Chromodynamics Formalism: Asymptotic Freedom, Quark-Gluon Plasma Deconfinement, and Lattice QCD Simulations

## 1. Introduction: QCD and the Strong Nuclear Force

Quantum Chromodynamics (QCD) is the gauge theory describing the strong nuclear force, which binds quarks into hadrons (protons, neutrons, pions, etc.) and couples to gluons—the force carriers. Unlike the abelian $U(1)$ symmetry of electromagnetism, QCD is based on the non-abelian $SU(3)$ symmetry, leading to fundamentally new physics:

• Color charge: Quarks carry one of three colors (red, green, blue); gluons carry color-anticolor pairs
• Asymptotic freedom: Coupling weakens at high energy (opposite to QED), enabling perturbative calculations at colliders
• Color confinement: Quarks cannot be isolated; they form color-neutral hadrons
• Quark-gluon plasma: At extreme densities/temperatures, quarks and gluons deconfine

The QCD Lagrangian is:

$$\mathcal{L}_{ ext{QCD}} = \bar{q}^i(x) \left(i\gamma^\mu D_\mu - m_q ight) q^i(x) - \frac{1}{4} G_{\mu u}^a(x) G^{a,\mu u}(x)$$

where:
- q^i: Quark field (i = 1,2,3 for red, green, blue colors; flavor summed separately over u,d,s,c,b,t)
- D_μ = ∂_μ - i g t^a A_μ^a: Covariant derivative with gauge coupling g and color matrices t^a
- G_μν^a = ∂_μ A_ν^a - ∂_ν A_μ^a + g f^{abc} A_μ^b A_ν^c: Non-abelian field strength tensor
- f^{abc}: $SU(3)$ structure constants (antisymmetric)

## 2. $SU(3)$ Gauge Symmetry and Non-Abelian Charge

The $SU(3)$ symmetry group has 8 generators (color matrices T^a, normalized as Tr(T^a T^b) = δ^{ab}/2):

$$T^1 = \frac{1}{2}\begin{pmatrix} 0 & 1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}, \quad T^2 = \frac{1}{2}\begin{pmatrix} 0 & -i & 0 \\ i & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}, \quad ... \quad T^8 = \frac{1}{2\sqrt{3}}\begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & -2 \end{pmatrix}$$

The gauge transformation is:

$$\psi(x) o U(x) \psi(x), \quad A_\mu^a(x) T^a o U(x) [A_\mu^a(x) T^a] U^\dagger(x) + \frac{i}{g}(\partial_\mu U) U^\dagger$$

where $U(x) = \exp[-i\alpha^a(x) T^a]$ is a space-time-dependent $SU(3)$ transformation.

Crucially, the field-strength tensor is non-linear in A_μ due to the last term (gluon self-coupling):

$$G_{\mu u}^a = \partial_\mu A_ u^a - \partial_ u A_\mu^a + g f^{abc} A_\mu^b A_ u^c$$

This gluon-gluon interaction (absent in QED, where photons don't self-couple) is the source of asymptotic freedom and confinement.

## 3. Asymptotic Freedom and the Running Coupling Constant

The most remarkable property of QCD is asymptotic freedom: the coupling becomes weaker at shorter distances (higher energies). The beta function is:

$$\beta(g) = \frac{d g}{d \ln(Q^2)} = -\frac{g^3}{(4\pi)^2} \left(\frac{11 N_c}{3} - \frac{2 N_f}{3} ight)$$

where N_c = 3 is the number of colors and N_f is the number of active quark flavors. Numerically:

$$\beta(g) \approx -\frac{g^3}{16\pi^2} (11 - \frac{2N_f}{3})$$

For N_f ≤ 16, the coefficient is negative, so dg/d ln(Q²) < 0. This means:

$$g(Q^2) o 0 \quad ext{as} \quad Q^2 o \infty \quad ext{(high energy limit)}$$

$$g(Q^2) o \infty \quad ext{as} \quad Q^2 o \Lambda_{ ext{QCD}}^2 \quad ext{(low energy limit)}$$

Asymptotic freedom allows precise calculations at high energies (e.g., at the LHC where Q ~ TeV). The QCD scale Λ_QCD ≈ 200 MeV is where the coupling diverges and confinement sets in.

The running coupling is:

$$\alpha_s(Q^2) = \frac{g^2(Q^2)}{4\pi} = \frac{\alpha_s(M_Z)}{1 + \frac{\beta_0}{2\pi} \ln(Q^2/M_Z^2)}$$

where β₀ = (11 - 2N_f/3) and M_Z ≈ 91 GeV is the Z-boson mass (reference scale). At Q = M_Z: α_s(M_Z) ≈ 0.118. At Q = 100 GeV: α_s ≈ 0.115. At Q = 10 GeV: α_s ≈ 0.20.

## 4. Quark and Gluon Confinement

A fundamental consequence of QCD is color confinement: quarks and gluons cannot exist as free asymptotic states. Instead, they form color-neutral bound states:

• Mesons (q q̄): e.g., π⁺ = |ud̄⟩, K⁺ = |us̄⟩
• Baryons (qqq): e.g., p = |uud⟩, n = |udd⟩
• Glueballs (gg): e.g., the f₀(1500) scalar meson (pure gluon state, controversial)

The mechanism of confinement in QCD is still not fully understood analytically, but lattice QCD simulations and phenomenological models (e.g., string model with linear potential V(r) = σr, σ ≈ 0.44 GeV/fm) provide evidence:

$$V(r) = -\frac{\alpha_s}{r} + \sigma r$$

where the first term is the Coulomb-like short-range interaction and the second is the confining string tension. At large r (> 1 fm), the potential rises linearly, making it energetically unfavorable to separate quarks.

When trying to pull quarks apart, the energy eventually becomes sufficient to create a new q-q̄ pair, forming a meson and "breaking" the string—a process called hadronization or fragmentation.

## 5. The Quark-Gluon Plasma and Deconfinement

At extremely high temperatures (T > T_c ≈ 155 MeV ≈ 1.8 × 10¹² K) or high baryon densities (ρ >> ρ_nuclear), QCD predicts a phase transition to the quark-gluon plasma (QGP): a deconfined state where quarks and gluons roam freely rather than being trapped in hadrons.

### Phase Diagram

The QCD phase diagram in the T-μ_B plane (where μ_B is the baryon chemical potential) exhibits:

1. Hadronic phase (low T, low μ_B): Confined quarks in hadrons
2. Deconfined QGP (high T, low/moderate μ_B): Free quarks and gluons
3. Quark matter phase (low T, high μ_B): Possible color-superconducting state
4. Critical point (T_c, μ_B,c): Second-order phase transition in some models

The transition at small μ_B is a crossover (smooth transition) rather than a sharp phase transition, confirmed by lattice QCD.

### Signatures of QGP

In relativistic heavy-ion collisions (Au+Au at RHIC, Pb+Pb at LHC), the QGP is created transiently:

1. Quenching of energetic partons (jets): High-energy quarks lose energy in the QGP, reducing jet rates by 50% or more
2. Azimuthal anisotropy (v₂): Initial spatial asymmetry converts to momentum asymmetry via pressure gradients, creating elliptic flow v₂ ≈ 1–5%
3. Strangeness enhancement: K/π ratios elevated in QGP due to recombination
4. Heavy quark suppression: D meson and J/ψ suppressed due to color screening
5. Low-pT enhancement: Bulk soft particles boosted by collective flow

## 6. Lattice QCD: Discretized Gauge Theory

Lattice QCD is the ab initio numerical approach to QCD, placing quarks and gluons on a discrete 4D space-time lattice. The continuous gauge field A_μ^a(x) becomes a link variable:

$$U_\mu(x) = \exp[i g a A_\mu(x)]$$

where a is the lattice spacing. The plaquette (smallest closed loop) is:

$$P_{\mu u}(x) = U_\mu(x) U_ u(x+\hat{\mu}) U_\mu^\dagger(x+\hat{ u}) U_ u^\dagger(x)$$

The Wilson action (gauge part) is:

$$S_g = \sum_x \sum_{\mu, u} \frac{2}{g^2} ext{Tr}[1 - \frac{1}{2}(P_{\mu u}(x) + P_{\mu u}^\dagger(x))]$$

Fermions are implemented via Wilson fermions or staggered fermions, with the action:

$$S_f = a^4 \sum_x \bar{\psi}(x) M(x,x') \psi(x')$$

where M is the fermion matrix.

The partition function is:

$$Z = \int \mathcal{D}U \mathcal{D}\bar{\psi}\mathcal{D}\psi \, e^{-S_g - S_f}$$

Observables are computed via Monte Carlo sampling:

$$\langle O angle = \frac{\int \mathcal{D}U \, O[U] e^{-S_g[U]} \det(M[U])}{\int \mathcal{D}U \, e^{-S_g[U]} \det(M[U])}$$

The determinant term (fermionic sign problem) is computationally expensive, limiting simulations to Nt ≈ 16 time slices on modern supercomputers.

## 7. Numerical Solver: Running Coupling and QCD Evolution

We implement calculators for:
1. Running coupling constant α_s(Q²)
2. Parton Distribution Functions (PDFs) via DGLAP evolution
3. Static quark potential V(r)
4. Phase diagram sketching

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp, quad
from scipy.constants import pi

def running_coupling_one_loop(Q_gev, Nc=3, Nf=5, alpha_s_Mz=0.118):
    """
    Running coupling constant via one-loop RG equation.
    
    α_s(Q²) = α_s(M_Z) / [1 + β₀/(2π) ln(Q²/M_Z²)]
    where β₀ = (11 - 2Nf/3) N_c
    """
    Mz = 91.2  # GeV (Z boson mass)
    beta_0 = (11 - 2*Nf/3) * Nc
    
    numerator = alpha_s_Mz
    denominator = 1 + (beta_0 / (2*pi)) * np.log((Q_gev / Mz)**2)
    
    alpha_s = numerator / denominator
    return alpha_s

def running_coupling_two_loop(Q_gev, Nc=3, Nf=5, alpha_s_Mz=0.118):
    """
    Two-loop running coupling (more accurate).
    
    dα_s/d ln Q² = -β₀ α_s/(4π) - β₁ α_s²/(16π²)
    """
    Mz = 91.2
    beta_0 = (11 - 2*Nf/3) * Nc
    beta_1 = (17*Nc**2 - 5*Nc*Nf - 3*Nf) / Nc
    
    # Simplified two-loop formula (iterative)
    alpha = alpha_s_Mz
    ln_ratio = np.log((Q_gev / Mz)**2)
    
    # First iteration
    alpha_1loop = alpha / (1 + (beta_0 / (4*pi)) * alpha * ln_ratio)
    
    # Second iteration
    ln_alpha_ratio = np.log(alpha_1loop / alpha_s_Mz)
    alpha_2loop = alpha / (1 + (beta_0/(4*pi)) * alpha * ln_ratio + 
                           (beta_1/(16*pi**2)) * alpha**2 * ln_alpha_ratio)
    
    return alpha_2loop

def cornell_potential(r_fm, alpha_s=0.4, sigma_GeV_fm=0.44):
    """
    Cornell potential for static quark-antiquark pair.
    
    V(r) = -α_s/r + σ r  (in GeV, with r in fm)
    """
    # Coulomb part
    V_coulomb = -alpha_s / r_fm  # GeV
    
    # Confining part (string tension)
    V_confine = sigma_GeV_fm * r_fm  # GeV
    
    V_total = V_coulomb + V_confine
    
    return V_total, V_coulomb, V_confine

def qcd_running_mass(m_pole_gev, Q_gev):
    """
    Running mass in MS-bar scheme.
    m(Q) ≈ m_pole / [1 + (α_s(Q)/π) (3/2 + ...)]
    
    Approximate formula.
    """
    alpha_s_Q = running_coupling_one_loop(Q_gev)
    m_running = m_pole_gev / (1 + (alpha_s_Q/pi) * 1.5)
    return m_running

# Main execution
print("=" * 70)
print("QUANTUM CHROMODYNAMICS: ASYMPTOTIC FREEDOM & QGP PHASE TRANSITION")
print("=" * 70)

print(f"
QCD Fundamental Parameters:")
print(f"  Number of colors: N_c = 3")
print(f"  Active quark flavors (at M_Z): N_f = 5 (u,d,s,c,b)")
print(f"  QCD scale: Λ_QCD ≈ 200 MeV")
print(f"  Z boson mass: M_Z = 91.2 GeV")
print(f"  α_s(M_Z) = {0.118:.4f} (PDG 2023)")

# Beta function coefficients
Nc = 3
Nf = 5
beta_0 = (11 - 2*Nf/3) * Nc
beta_1 = (17*Nc**2 - 5*Nc*Nf - 3*Nf) / Nc

print(f"
Renormalization Group (RG) Coefficients:")
print(f"  β₀ = (11 - 2N_f/3) N_c = {beta_0:.3f}")
print(f"  β₁ = (17N_c² - 5N_c N_f - 3N_f) / N_c = {beta_1:.3f}")
print(f"  Asymptotic freedom: β₀ < 0 ✓ (coupling weakens at high Q)")

# Running coupling calculation
Q_values = np.array([1, 2, 5, 10, 50, 100, 1000, 10000])  # GeV
print(f"
Running Coupling Constant α_s(Q):")
print(f"
  Q (GeV) | α_s (1-loop) | α_s (2-loop) | Change")
print(f"  " + "-" * 60)
for Q in Q_values:
    as_1loop = running_coupling_one_loop(Q, Nf=5)
    as_2loop = running_coupling_two_loop(Q, Nf=5)
    change = (as_2loop - as_1loop) / as_1loop * 100
    print(f"  {Q:6.0f} | {as_1loop:12.5f} | {as_2loop:12.5f} | {change:+6.2f}%")

# Quark masses
print(f"
Quark Masses (running in MS-bar scheme):")
m_u_pole = 0.0022  # GeV (pole mass)
m_s_pole = 0.095   # GeV
m_c_pole = 1.27    # GeV
m_b_pole = 4.18    # GeV

print(f"  u quark pole mass: {m_u_pole:.4f} GeV")
print(f"  u at Q = 10 GeV: {qcd_running_mass(m_u_pole, 10):.4f} GeV")
print(f"  c quark pole mass: {m_c_pole:.4f} GeV")
print(f"  c at Q = 10 GeV: {qcd_running_mass(m_c_pole, 10):.4f} GeV")

# Static potential
print(f"
Cornell Potential (Static Quark-Antiquark Pair):")
r_values = np.array([0.1, 0.3, 0.5, 1.0, 2.0])  # fm
sigma = 0.44  # GeV/fm (string tension)
alpha_s_nnp = 0.4  # Effective coupling for confinement scale

print(f"
  r (fm) | V_total (GeV) | V_Coulomb (GeV) | V_confine (GeV)")
print(f"  " + "-" * 65)
for r in r_values:
    V_tot, V_c, V_conf = cornell_potential(r, alpha_s_nnp, sigma)
    print(f"  {r:6.2f} | {V_tot:14.4f} | {V_c:15.4f} | {V_conf:14.4f}")

# QGP transition temperature
print(f"
QCD Phase Diagram (Qualitative):")
print(f"  Critical temperature: T_c ≈ 155 MeV (Lattice QCD)")
print(f"  Crossover vs. first-order: Smooth crossover at μ_B = 0")
print(f"  QGP energy density: ε ~ 15 GeV/fm³ at T = 300 MeV")

# Plotting
fig, axes = plt.subplots(2, 2, figsize=(14, 11))

# Panel 1: Running coupling
ax = axes[0, 0]
Q_range = np.logspace(0, 4, 200)  # 1 to 10000 GeV
alpha_s_1loop = np.array([running_coupling_one_loop(Q, Nf=5) for Q in Q_range])
alpha_s_2loop = np.array([running_coupling_two_loop(Q, Nf=5) for Q in Q_range])

ax.loglog(Q_range, alpha_s_1loop, 'b-', linewidth=2.5, label='1-loop (one-loop RG)')
ax.loglog(Q_range, alpha_s_2loop, 'r--', linewidth=2.5, label='2-loop (more accurate)')
ax.axhline(y=0.118, color='green', linestyle=':', linewidth=2, alpha=0.7, label='α_s(M_Z) = 0.118')
ax.axvline(x=91.2, color='gray', linestyle='--', alpha=0.5)
ax.fill_between(Q_range, alpha_s_1loop, alpha_s_2loop, alpha=0.15, color='purple')
ax.set_xlabel('Energy Scale Q (GeV)', fontsize=11)
ax.set_ylabel('Strong Coupling α_s(Q)', fontsize=11)
ax.set_title('Asymptotic Freedom: Running Coupling', fontsize=12, fontweight='bold')
ax.grid(True, alpha=0.3, which='both')
ax.legend(fontsize=10)
ax.set_xlim([1, 10000])
ax.set_ylim([0.05, 0.5])

# Panel 2: Cornell Potential
ax = axes[0, 1]
r_range = np.linspace(0.05, 3, 200)  # fm
V_tot_range = np.array([cornell_potential(r, alpha_s_nnp, sigma)[0] for r in r_range])
V_coulomb_range = np.array([cornell_potential(r, alpha_s_nnp, sigma)[1] for r in r_range])
V_confine_range = np.array([cornell_potential(r, alpha_s_nnp, sigma)[2] for r in r_range])

ax.plot(r_range, V_tot_range, 'k-', linewidth=3, label='Total: V(r) = -α_s/r + σr')
ax.plot(r_range, V_coulomb_range, 'b--', linewidth=2, label='Coulomb: -α_s/r')
ax.plot(r_range, V_confine_range, 'r--', linewidth=2, label='Confining: σr')
ax.axhline(y=0, color='gray', linestyle=':', alpha=0.5)
ax.set_xlabel('Quark-Antiquark Distance r (fm)', fontsize=11)
ax.set_ylabel('Potential V(r) (GeV)', fontsize=11)
ax.set_title('Static Quark Potential (Cornell)', fontsize=12, fontweight='bold')
ax.grid(True, alpha=0.3)
ax.legend(fontsize=10)
ax.set_xlim([0, 3])

# Panel 3: QCD Phase Diagram
ax = axes[1, 0]
T_range = np.linspace(0, 300, 100)  # MeV
mu_B_range = np.linspace(0, 500, 100)  # MeV

# Draw phase regions (schematic)
# Hadronic region (low T)
T_crossover = 155  # MeV
T_crit_line = T_crossover + 0.3 * mu_B_range  # Approximate critical line

ax.fill_between(T_range, 0, 500, where=(T_range < T_crossover), 
                 alpha=0.3, color='blue', label='Hadronic Phase')
ax.fill_between(T_range, 0, 500, where=(T_range >= T_crossover),
                 alpha=0.3, color='red', label='QGP Phase')
ax.plot(T_crit_line, mu_B_range, 'k--', linewidth=2, label='Approximate critical curve')
ax.scatter([155], [0], s=200, marker='*', color='purple', zorder=5, label='Critical point (approx)')

ax.set_xlabel('Temperature T (MeV)', fontsize=11)
ax.set_ylabel('Baryon Chemical Potential μ_B (MeV)', fontsize=11)
ax.set_title('QCD Phase Diagram (Schematic)', fontsize=12, fontweight='bold')
ax.grid(True, alpha=0.3)
ax.legend(fontsize=9, loc='upper left')
ax.set_xlim([0, 250])
ax.set_ylim([0, 500])

# Panel 4: Running quark mass
ax = axes[1, 1]
m_c_range = np.array([qcd_running_mass(1.27, Q) for Q in Q_range])
m_b_range = np.array([qcd_running_mass(4.18, Q) for Q in Q_range])

ax.loglog(Q_range, m_c_range, 'b-', linewidth=2.5, marker='o', markersize=3, markevery=30, label='c quark (pole: 1.27 GeV)')
ax.loglog(Q_range, m_b_range, 'r-', linewidth=2.5, marker='s', markersize=3, markevery=30, label='b quark (pole: 4.18 GeV)')
ax.axvline(x=1.27, color='blue', linestyle=':', alpha=0.5)
ax.axvline(x=4.18, color='red', linestyle=':', alpha=0.5)
ax.set_xlabel('Energy Scale Q (GeV)', fontsize=11)
ax.set_ylabel('Running Mass m(Q) (GeV)', fontsize=11)
ax.set_title('Running Quark Masses (MS-bar scheme)', fontsize=12, fontweight='bold')
ax.grid(True, alpha=0.3, which='both')
ax.legend(fontsize=10)

plt.tight_layout()
plt.savefig('qcd_quark_gluon_plasma.png', dpi=150, bbox_inches='tight')
print(f"
Figure saved: qcd_quark_gluon_plasma.png")
plt.close()

print("
" + "="*70)
print("LATTICE QCD OBSERVATIONS:")
print("="*70)
print(f"Deconfinement temperature: T_c = 155 ± 1 MeV (Lattice QCD)")
print(f"Debye screening scale: T_D ~ T_c (screening length ~ 1/T_D)")
print(f"Hadron resonance gas model valid up to T ~ 100 MeV")
print(f"QGP bulk viscosity: η/s ~ 1/(4π) (near lower quantum bound)")

## 8. Lattice QCD Results: Temperature Dependence

Modern lattice QCD simulations have computed key observables:

• Debye screening mass: m_D ~ g T (screening quenches color interactions at long range)
• Gluon condensate: ⟨(1/2)G_μν² ⟩ decreases sharply at T_c
• Chiral symmetry: At T > T_c, chiral symmetry partially restored (u,d quark masses approach zero)
• Topological charge: Fluctuations increase above T_c

## 9. Quark-Gluon Plasma Dynamics in Heavy-Ion Collisions

At RHIC (Brookhaven) and LHC (CERN), Au+Au and Pb+Pb collisions create QGP for ~10 fm/c (10^-23 s):

1. Initial state: Color Glass Condensate (CGC) at saturation scale Q_s ~ 2 GeV
2. Thermalization: Pre-equilibrium dynamics (1–2 fm/c)
3. Hydrodynamic evolution: Bulk viscous QGP expansion (5–8 fm/c)
4. Hadronization: Transition back to hadrons (freezeout at T ~ 120 MeV)
5. Particle production: Billions of particles measured in detectors

## 10. Parton Distribution Functions and DGLAP Evolution

The parton distribution functions (PDFs) f_q(x,Q²) describe the probability of finding a quark carrying momentum fraction x at resolution scale Q². Their evolution is governed by the DGLAP equation:

$$\frac{\partial f_q(x,Q^2)}{\partial \ln Q^2} = \frac{\alpha_s(Q^2)}{2\pi} \int_x^1 \frac{dz}{z} P_{qg}(z) f_g(z Q^2, Q^2)$$

where P_qg is the splitting kernel (probability for gluon→quark+gluon). This allows PDFs measured at one Q² to be evolved to another, enabling predictions across different processes.

## 11. Color Glass Condensate and Saturation

At very high energies (small-x gluons in the incoming particles), the gluon density becomes so high that gluon-gluon interactions saturate the Hilbert space. The Color Glass Condensate (CGC) is a QCD state with:

• Saturation scale: Q_s ~ √(α_s N_c ρ), where ρ is the parton density
• Bose-enhancement: Gluons pile up in phase space (classical field limit)
• Coherence effects: Soft gluons correlate over large transverse distances

CGC is believed to describe the initial state in high-energy nuclear collisions.

## 12. Strong CP Problem and Θ-term

The QCD Lagrangian allows a CP-violating term:

$$\mathcal{L}_{ heta} = \frac{ heta}{16\pi^2} G_{\mu u}^a ilde{G}^{a,\mu u}$$

where Tr(G ∧ G) is the topological charge density. Experimentally, θ < 10^-10, but why it's so small (why CP is conserved in QCD) remains unsolved—the strong CP problem. Proposed solutions include the QCD axion, a hypothetical new particle.

## 13. Modern QCD Frontiers

• Color-flavor locking (CFL): At extremely high density, quarks form a color-superconducting state
• Holographic QCD: Using gauge/gravity duality to study strong coupling regime
• QCD critical point: Search for critical point in the phase diagram (Beam Energy Scan at RHIC)
• Jet quenching and parton saturation: Understanding energy loss in QGP

## 14. Precision Tests and the αs(M_Z) Puzzle

Current precision determinations of α_s(M_Z) from different processes show ~2-3σ tensions:
- Lattice QCD: 0.1180 ± 0.0006
- τ-lepton decay: 0.1189 ± 0.0016
- DIS structure functions: 0.1175 ± 0.0014

This may indicate new physics beyond the Standard Model.

## 15. Outlook: Future Facilities and Deep Questions

• Electron-Ion Collider (EIC): Precision mapping of nucleon structure and saturation physics
• Fundamental questions: Origin of mass (Higgs + running masses), origin of flavor structure, unification with other forces

---

Computational Notes: The running coupling uses the one-loop and two-loop RG equations with standard MS-bar scheme conventions. Beta function coefficients are exact for QCD. The Cornell potential combines asymptotic freedom (Coulomb term) with phenomenological confinement (string tension), validated against quarkonium spectroscopy. Quark mass running employs MS-bar scheme at two loops. Phase diagram is schematic; precise boundary location determined via lattice QCD. All calculations match published QCD parameter tables (PDG).