Probing the photon emissivity of the quark-gluon plasma without an inverse problem in lattice QCD

Marco Cè,Tim Harris,Ardit Krasniqi,Harvey B. Meyer,and Csaba Török

DOI: https://doi.org/10.1103/physrevd.109.014507

IF: 5.407

2024-01-19

Physical Review D

Abstract:The thermal photon emissivity of the quark-gluon plasma is determined by the in-medium spectral function of the electromagnetic current at lightlike kinematics, σ(ω) . In this work, we present the first lattice QCD results on moments of σ(ω)/ω , defined by the weight function 1/(ω2+(2πTn)2) , n∈Z and computed without encountering an inverse problem. We employ two dynamical flavors of O(a) -improved Wilson fermions at a temperature T≈250 MeV and perform the continuum limit. We compare our results for the first two moments to those obtained dispersively by integrating over the spectral function computed at weak coupling by Arnold, Moore and Yaffe. https://doi.org/10.1103/PhysRevD.109.014507 Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI. Funded by SCOAP 3 . Published by the American Physical Society

astronomy & astrophysics,physics, particles & fields

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to directly determine the moments of the in - medium spectral function of the photon emission rate in Quark - Gluon Plasma (QGP) by the Lattice QCD (Quantum Chromodynamics) method without encountering the inverse problem. Specifically, the authors calculated the first two moments of the spectral function \(\sigma(\omega)\), which are defined as the integrals under the weight function \(\frac{1}{\omega^2+(2\pi Tn)^2}\), where \(n\in\{1, 2\}\).

### Background and Motivation

1. **Importance of Photon Emission Rate**
- As the quantum of electromagnetic radiation, photons can escape from the strongly interacting medium due to their relatively large free path compared to strongly interacting matter at high temperatures, thus providing information about the space - time evolution of the medium.
- Experimentally, by measuring the photon spectrum, physical properties such as the temperature of the medium can be inferred.

2. **Limitations of Existing Models**
- Existing models usually rely on the predictions of QCD perturbation theory, but these predictions have difficulties in dealing with non - equilibrium quark - gluon plasma.
- The difference between experimental data and theoretical models, especially in the low transverse momentum region, is known as the "direct photon puzzle".

### Main Contributions of the Paper

1. **Avoiding the Inverse Problem**
- The traditional Lattice QCD method needs to obtain the spectral density from the correlation function in Euclidean space through analytic continuation, which is a numerically unstable problem.
- This paper proposes a new method to obtain the integral of the spectral function by calculating a correlation function that can be directly accessed in Lattice QCD without facing the inverse problem.

2. **Calculation Method**
- The authors used two dynamical \(O(a)\)-improved Wilson fermions and carried out the continuous - limit calculation at a temperature of about 250 MeV.
- By calculating the transverse Euclidean correlation function in the imaginary space - momentum, the first two moments of the spectral function were obtained and compared with the predictions of the weak - coupling theory.

### Conclusion

- This research successfully calculated the moments of the spectral function of the photon emission rate of quark - gluon plasma by the Lattice QCD method for the first time without involving the inverse problem.
- These results provide an important non - perturbative test for verifying the predictions of the weak - coupling QCD theory and help to better understand the properties of quark - gluon plasma.

### Mathematical Expressions

- **Definition of Spectral Function**
\[
\rho_{\mu\nu}(\omega, k)=\int d^{4}x e^{i(\omega t - k\cdot x)}\langle[J_{\text{em}}^{\mu}(x), J_{\text{em}}^{\nu}(0)]\rangle
\]
where \(J_{\text{em}}^{\mu}(x)=\sum_{f}Q_{f}\bar{\psi}_{f}(x)\gamma^{\mu}\psi_{f}(x)\) is the electromagnetic current, \(Q_{f}\) is the quark charge, and \(\gamma^{\mu}\) is the Euclidean Dirac matrix.

- **Formula for Photon Emission Rate**
\[
\frac{d\Gamma_{\gamma}(\omega)}{d\omega}=\frac{\alpha_{\text{em}}}{\pi^{2}\omega}\sigma(\omega)\frac{1}{e^{\omega/T}- 1}+O(\alpha_{\text{em}}^{2})
\]
where \(\sigma(\omega)\) is the spectral function of the two - point correlation function of the electromagnetic current with zero imaginary part.

- **Definition of Moments**
\[
H_E(\omega_n)=-\frac{\o