Next Article in Journal
Measurements of the Υ Meson Production in Au + Au Collisions by the STAR Experiment
Previous Article in Journal
Bulk Observables within a Hybrid Approach for Heavy Ion Collisions with SMASH Afterburner
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Insight into Thermal Modifications of Quarkonia From a Comparison of Continuum-Extrapolated Lattice Results to Perturbative QCD †

1
Fakultät für Physik, Universität Bielefeld, 33615 Bielefeld, Germany
2
Key Laboratory of Quark and Lepton Physics (MOE) and Institute of Particle Physics, Central China Normal University, Wuhan 430079, China
3
Center for Computational Sciences, University of Tsukuba, Tsukuba, Ibaraki 305-8577, Japan
*
Author to whom correspondence should be addressed.
Presented at Hot Quarks 2018—Workshop for Young Scientists on the Physics of Ultrarelativistic Nucleus-Nucleus Collisions, Texel, The Netherlands, 7–14 September 2018.
Published: 5 May 2019

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.

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
G P S ( τ ) = M B 2 x ( ψ ¯ i γ 5 ψ ) ( τ , x ) ( ψ ¯ i γ 5 ψ ) ( 0 , 0 ) c ,
where M B is a bare quark mass and the vector correlator by
G i i ( τ ) = x ( ψ ¯ γ i ψ ) ( τ , x ) ( ψ ¯ γ i ψ ) ( 0 , 0 ) c .
These correlators are related to the corresponding spectral functions through an integral equation,
G P S , i i ( τ ) = 0 d ω π ρ P S , i i ( ω ) K ( ω , τ ) with K ( ω , τ ) = cosh ( ω ( τ 1 2 T ) ) sinh ( ω 2 T ) ,
where the temperature T is given by the inverse of the lattice extent. Due to the structure of the integration kernel K ( ω , τ ) , the low frequency regime of the spectral function influences the shape of the correlator at larger τ T while the high frequency part dominates the small distance part of the correlator. Figure 1 shows the contribution of different frequency regimes to the correlator. In the vector channel, a transport peak is expected at ω 0 . Its contribution mainly influences the correlator at large τ T , but has negligible effect on the low τ 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 [1]. 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 [2]. The preferred option is to get rid of renormalization constants by building the renormalization independent ratio with the quark number susceptibility χ q . Since χ 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 χ q at the temperature 2.25 T 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 / ψ and Υ at each point in τ 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 τ T (see Figure 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 [3,4]. 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 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 ρ m o d e l ( ω ) = A ρ p e r t ( ω 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 Figure 3. In addition, A is close to 1 and B is small (see Table 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 τ T region of the correlator but has only little influence on the small τ T regime, we only fit the correlator for distances smaller than τ T = 0.25. In the vector channel, the results at 2.25 T c in Figure 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 [5] 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.

References

  1. Ding, H.-T.; Kaczmarek, O.; Kruse, A.-L.; Ohno, H.; Sandmeyer, H. Continuum extrapolation of quarkonium correlators at non-zero temperature. EPJ Web. Conf. 2018, 175, 07010. [Google Scholar] [CrossRef]
  2. Lüscher, M.; Sint, S.; Sommer, R.; Wittig, H. Non-perturbative determination of the axial current normalization constant in O(a) improved lattice QCD. Nucl. Phys. 1997, B491, 344. [Google Scholar] [CrossRef]
  3. Burnier, Y.; Ding, H.-T.; Kaczmarek, O.; Kruse, A.-L.; Laine, M.; Ohno, H.; Sandmeyer, H. Thermal quarkonium physics in the pseudoscalar channel. JHEP 2017, 11, 206. [Google Scholar] [CrossRef]
  4. Burnier, Y.; Laine, M.; Vepsalainen, M. Heavy quarkonium in any channel in resummed hot QCD. JHEP 2008, 1, 043. [Google Scholar] [CrossRef]
  5. Ding, H.-T.; Kaczmarek, O.; Mukherjee, S.; Ohno, H.; Shu, H.-T. Stochastic reconstructions of spectral functions: Application to lattice QCD. Phys. Rev. 2018, D97, 094503. [Google Scholar] [CrossRef]
Figure 1. The influence of different parts of the spectral function to the different regimes of the correlator.
Figure 1. The influence of different parts of the spectral function to the different regimes of the correlator.
Proceedings 10 00045 g001
Figure 2. Example of the continuum extrapolation in the pseudoscalar channel for three points in τ T .
Figure 2. Example of the continuum extrapolation in the pseudoscalar channel for three points in τ T .
Proceedings 10 00045 g002
Figure 3. Fit of the model spectral function (blue) to the lattice data in the charmonium pseudoscalar channel. The yellow line corresponds to the perturbative correlator calculated with a 10% uncertainty in the mass.
Figure 3. Fit of the model spectral function (blue) to the lattice data in the charmonium pseudoscalar channel. The yellow line corresponds to the perturbative correlator calculated with a 10% uncertainty in the mass.
Proceedings 10 00045 g003
Figure 4. Fit of the model spectral function (blue) to the lattice data up to τ T = 0.25 in the vector channel and the perturbative correlator (yellow).
Figure 4. Fit of the model spectral function (blue) to the lattice data up to τ T = 0.25 in the vector channel and the perturbative correlator (yellow).
Proceedings 10 00045 g004
Table 1. Results from the fit of the model spectral function to the lattice data in the pseudoscalar channel from [3].
Table 1. Results from the fit of the model spectral function to the lattice data in the pseudoscalar channel from [3].
CharmoniumBottomonium
T / T c A B / T χ 2 / d . o . f . A B / T χ 2 / d . o . f .
1.11.040.520.010.85−0.110.02
1.31.040.370.010.87−0.130.04
1.51.020.330.020.87−0.110.10
2.251.060.160.080.93−0.040.28
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kruse, A.-L.; Ding, H.-T.; Kaczmarek, O.; Ohno, H.; Sandmeyer, H. Insight into Thermal Modifications of Quarkonia From a Comparison of Continuum-Extrapolated Lattice Results to Perturbative QCD. Proceedings 2019, 10, 45. https://0-doi-org.brum.beds.ac.uk/10.3390/proceedings2019010045

AMA Style

Kruse A-L, Ding H-T, Kaczmarek O, Ohno H, Sandmeyer H. Insight into Thermal Modifications of Quarkonia From a Comparison of Continuum-Extrapolated Lattice Results to Perturbative QCD. Proceedings. 2019; 10(1):45. https://0-doi-org.brum.beds.ac.uk/10.3390/proceedings2019010045

Chicago/Turabian Style

Kruse, Anna-Lena, H.-T. Ding, O. Kaczmarek, H. Ohno, and H. Sandmeyer. 2019. "Insight into Thermal Modifications of Quarkonia From a Comparison of Continuum-Extrapolated Lattice Results to Perturbative QCD" Proceedings 10, no. 1: 45. https://0-doi-org.brum.beds.ac.uk/10.3390/proceedings2019010045

Article Metrics

Back to TopTop