Insight into thermal modifications of quarkonia from a comparison of continuum-extrapolated lattice results to perturbative QCD$^{\dagger}$

TL;DR

This paper compares continuum-extrapolated lattice QCD results with perturbative calculations to understand thermal effects on quarkonium states and heavy quark diffusion, highlighting qualitative agreement and systematic uncertainties.

Contribution

It provides a detailed comparison between lattice and perturbative spectral functions for quarkonia at finite temperature, improving understanding of their thermal modifications.

Findings

Qualitative agreement between lattice and perturbative correlators.

Systematic uncertainties account for quantitative differences.

Analysis extends from pseudoscalar to vector channels for transport insights.

Abstract

In this work, we strive to gain insight into thermal modifications of charmonium and bottomonium bound states as well as the heavy quark diffusion coefficient. The desired information is contained in the spectral function which can not be calculated on the lattice directly. Instead, the correlator given by an integration over the spectral function times an integration kernel is obtained. Extracting the spectral function is an ill-posed inversion problem and various different solutions have been proposed. We focus on a comparison to a spectral function obtained from combining perturbative and pNRQCD calculations. In order to get precise results, continuum extrapolated correlators originating from large and fine lattices are used. We first analyze the pseudoscalar channel since the absence of a transport peak simplifies the analysis. The knowledge gained from this is then used to extend the analysis to the vector channel, where information on heavy quark transport is encoded in the low frequency regime of the spectral function. The comparison shows a qualitatively good agreement between perturbative and lattice correlators. Quantitative differences can be explained by systematic uncertainties.

Figures

Tables

Table 1. Table 1: Results from the fit of the model spectral function to the lattice data in the pseudoscalar channel from Author1 ( year ) .

	Charmonium			Bottomonium
$T/T_c$	$A$	$B/T$	${\chi }^2/d.o.f.$	$A$	$B/T$	${\chi }^2/d.o.f.$
1.1	1.04	0.52	0.01	0.85	-0.11	0.02
1.3	1.04	0.37	0.01	0.87	-0.13	0.04
1.5	1.02	0.33	0.02	0.87	-0.11	0.10
2.25	1.06	0.16	0.08	0.93	-0.04	0.28

Equations

$$G_{PS}(\tau )=M_B^2{\int }_{\vec{x}}{⟨(\bar{\psi }i{\gamma }_5\psi )(\tau ,\vec{x})(\bar{\psi }i{\gamma }_5\psi )(0,\vec{0})⟩}_c,$$

$$G_{ii}(\tau )={\int }_{\vec{x}}{⟨(\bar{\psi }{\gamma }_i\psi )(\tau ,\vec{x})(\bar{\psi }{\gamma }_i\psi )(0,\vec{0})⟩}_c.$$

$$G_{PS,ii}(\tau )=\underset{0}{\overset{\infty }{\int }}\frac{\mathrm{d}\omega }{\pi }{\rho }_{PS,ii}(\omega )K(\omega ,\tau )\text{with }K(\omega ,\tau )=\frac{\cosh (\omega (\tau -\frac{1}{2T}))}{\sinh (\frac{\omega }{2T})},$$

Full text

\pubvolume

xx \issuenum1 \articlenumber5

\historyReceived: date; Accepted: date; Published: date \TitleInsight into thermal modifications of quarkonia from a comparison of continuum-extrapolated lattice results to perturbative QCD*†* \AuthorAnna-Lena Kruse 1* with H.-T. Ding2, O. Kaczmarek1,2, H. Ohno3 and H. Sandmeyer1 \AuthorNamesFirstname Lastname, Firstname Lastname and Firstname Lastname

\corresCorrespondence: alkruse@physik.uni-bielefeld.de \firstnotePresented at Hot Quarks 2018 - Workshop for young scientists on the physics of ultrarelativistic nucleus-nucleus collisions, Texel, The Netherlands, September 7-14 2018

1 Introduction

Information on the in-medium properties like quarkonium bound states and heavy quark transport coefficients is contained in the spectral function. In perturbation theory, a spectral function can be constructed by combining different energy regimes. On the lattice, the spectral function can not be calculated directly. Instead, we compare lattice and perturbative correlators to gain insight into thermal modifications of quarkonia. Since the pseudoscalar channel does not contain a transport peak, it serves as an ideal probe for this comparison method. The knowledge obtained from pseudoscalar correlators can then be used to extend the analysis to the vector channel.

The connected (i.e. flavour non-singlet) imaginary time correlator in the pseudoscalar channel is given by

[TABLE]

where $M_{B}$ is a bare quark mass and the vector correlator by

[TABLE]

These correlators are related to the corresponding spectral functions through an integral equation,

[TABLE]

where the temperature $T$ is given by the inverse of the lattice extent. Due to the structure of the integration kernel $K(\omega,\tau)$, the low frequency regime of the spectral function influences the shape of the correlator at larger $\tau T$ while the high frequency part dominates the small distance part of the correlator. Fig. 1 shows the contribution of different frequency regimes to the correlator. In the vector channel, a transport peak is expected at $\omega\approx 0$. Its contribution mainly influences the correlator at large $\tau T$, but has negligible effect on the low $\tau T$ regime. Still, the transport peak complexifies the analysis and it is thus easier to test the method in the pseudoscalar channel, where no transport contribution is present, and use the gained knowledge for the vector channel analysis.

2 Continuum Extrapolation

In order to realize the large lattices needed for this work, we employed the standard plaquette gauge action and used the quenched approximation with Clover-improved Wilson valence quarks. Five temperatures in the range from $0.75~{}T_{c}$ to $2.25~{}T_{c}$ were generated on four different lattice spacings and the correlators were continuum extrapolated according to the method described in Author1 (year). The continuum extrapolation requires a renormalization of the correlators. For the pseudoscalar channel one- and two-loop perturbative renormalization constants are available. The continuum extrapolation is carried out using the two-loop expression and the difference is used to estimate a systematic error. In addition to the perturbative renormalization constants, non-perturbatively determined renormalization constants are available in the vector channel Author1 (year). The preferred option is to get rid of renormalization constants by building the renormalization independent ratio with the quark number susceptibility $\chi_{q}$. Since $\chi_{q}$ is given by the zeroth component of the vector correlator, it has the same renormalization that thus drops out in the ratio. We chose $\chi_{q}$ at the temperature $2.25T_{c}$ as normalization since the quark number susceptibility is more reliable at higher temperatures.

To ensure that the different lattices have the same vector meson mass, the correlators are mass-interpolated by fitting a quadratic exponential ansatz to the data points from different masses. The values of the fit function at the masses of $J/\psi$ and $\Upsilon$ at each point in $\tau T$ are taken as the new correlator. After the mass interpolation, the normalized correlators are continuum extrapolated with an ansatz quadratic in the lattice spacing. The continuum extrapolation is done for every point in $\tau T$ (see fig. 2 for examples) and leads to the continuum results shown in section 4.

3 The Perturbative Spectral Function

The spectral function to which the lattice data should be compared is derived from an interpolation of different calculations for different energy regimes as decribed in Author1 (year, year). Well above the threshold, the ultraviolet asymptotics are known up to high loop order in both channels at zero temperature. For frequencies around the threshold, thermal contributions arise and the physics for non-relativistic quark masses can be described using pNRQCD with a perturbative real time static potential. The threshold behaviour and the vacuum asymptotics are interpolated in such a way that the transition as well as its first derivative is continuous. For frequencies way below the threshold, this description overestimates the spectral function and the spectral function is modified to contain the expected exponential suppression.

4 Comparison of Lattice and Perturbative Correlators

A first look reveals a good qualitative agreement in the pseudoscalar channel. For a quantitative comparison, systematic uncertainties on both sides need to be taken into account. On the lattice side, the renormalization constants might be off. To account for this uncertainty, an overall factor $A$ is introduced. On the perturbative side, the relation between the pole mass and the $\overline{\text{MS}}$ mass is not exactly known which might result in a slightly smaller or larger threshold location. A mass shift $B$ is introduced, which can also account for any non-perturbative thermal mass shifts.

With these corrections, we construct a model spectral function $\rho^{model}(\omega)=A\rho^{pert}(\omega-B)$, where $A$ and $B$ are variables in a fit to the lattice data. The fit results match the lattice data almost perfectly as shown in fig. 3. In addition, $A$ is close to 1 and $B$ is small (see tab. 1). Overall, a good agreement between lattice and perturbative correlators is observed in the pseudoscalar channel.

In the vector channel we apply the same method as in the pseudoscalar channel, i.e. building a model spectral function that is fitted to the lattice correlator. But here the transport peak complexifies the analysis. The shape of the transport peak can not be described perturbatively, only the constant contribution to the correlator can be estimated. Since the transport peak dominates the behaviour of the high $\tau T$ region of the correlator but has only little influence on the small $\tau T$ regime, we only fit the correlator for distances smaller than $\tau T=~{}$0.25. In the vector channel, the results at $2.25~{}T_{c}$ in fig. 4 also show good agreement between the lattice and perturbative correlators. The small deviation at large distances in the charmonium results may indicate a small contribution of a transport peak.

5 Conclusions

Overall, a good agreement between perturbation theory and lattice data is observed when taking systematic errors into account. It should be noted that this might partly be due to the smaller effective strong coupling in the quenched approximation. In the future, a full QCD study is required. The next steps are to crosscheck our results for the spectral function with the result of Bayesian reconstruction methods Author1 (year) and to fit suitable ansätze for the transport peak in the vector channel correlators.

\funding

Funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project number 315477589 – TRR 211

\reftitle

References