Transport coefficients from in-medium quarkonium dynamics

This is a self-archived version of an original article. This version may differ from the original in pagination and typographic details. Author(s): Title: Year: Version: Copyright: Rights: Rights url: Please cite the original version: CC BY 4.0 https://creativecommons.org/licenses/by/4.0/ Transport coefficients from in-medium quarkonium dynamics © Authors, 2019 Published version Brambilla, Nora; Escobedo, Miguel A.; Vairo, Antonio; Vander Griend, Peter Brambilla, N., Escobedo, M. A., Vairo, A., & Vander Griend, P. (2019). Transport coefficients from in-medium quarkonium dynamics. Physical Review D, 100(5), Article 054025. https://doi.org/10.1103/PhysRevD.100.054025 2019 Transport coefficients from in-medium quarkonium dynamics Nora Brambilla,1,2,* Miguel A. Escobedo,3,4,†Antonio Vairo,1,‡and Peter Vander Griend 1,§ 1Physik Department, Technische Universität München, 85748 Garching, Germany 2Institute for Advanced Study, Technische Universität München, Lichtenbergstrasse 2 a, 85748 Garching, Germany 3Instituto Galego de Física de Altas Enerxías (IGFAE), Universidade de Santiago de Compostela, 15782 Galicia, Spain 4Department of Physics, University of Jyväskylä, P.O. Box 35, FI-40014 Jyväskylä, Finland (Received 22 July 2019; published 18 September 2019) The in-medium dynamics of heavy particles are governed by transport coefficients. The heavy quark momentum diffusion coefficient, κ, is an object of special interest in the literature, but one which has proven notoriously difficult to estimate, despite the fact that it has been computed by weak-coupling methods at next-to-leading order accuracy, and by lattice simulations of the pure SU(3) gauge theory. Another coefficient, γ, has been recently identified. It can be understood as the dispersive counterpart of κ. Little is known about γ. Both κand γare, however, of foremost importance in heavy quarkonium physics as they entirely determine the in and out of equilibrium dynamics of quarkonium in a medium, if the evolution of the density matrix is Markovian, and the motion, quantum Brownian; the medium could be a strongly or weakly coupled plasma. In this paper, using the relation between κ,γand the quarkonium in-medium width and mass shift respectively, we evaluate the two coefficients from existing 2þ1flavor lattice QCD data. The resulting range for κis consistent with earlier determinations, the one for γis the first nonperturbative determination of this quantity. DOI: 10.1103/PhysRevD.100.054025 I. INTRODUCTION Heavy quarkonium has long been theorized to serve as a probe of the medium formed in heavy-ion collisions with the purpose to detect a new state of matter, the quark gluon plasma (QGP) [1]. In turn, the study of the QGP offers a unique window on the Universe at an early time. In this paper, we focus on the out of equilibrium dynamics of heavy quarkonium in the medium. Our aims are twofold. Making use of recent results [2,3], we further elaborate on the out of equilibrium dynamics under the assumptions that the evolution of the quarkonium density is Markovian, and the motion, quantum Brownian. We emphasize the relation that exists, under these conditions, between the quarkonium dynamics and the transport coefficient κ, describing the momentum diffusion of a heavy quark in a medium, and γ, the dispersive counterpart of κ. Finally, exploiting this relation, and under the conditions of its validity, we provide a method for extracting the coefficients κand γfrom existing 2þ1 flavor lattice QCD data. A recent work on the extraction of heavy-quark transport coefficients is Ref. [4]. The remainder of the paper is structured as follows. In Sec. II, we summarize some recent progress in the study of the out of equilibrium dynamics of heavy particles and, in particular, heavy quarkonium in a medium, supplying the relevant background for the results to follow. Section III contains the evolution equations for the heavy quarkonium density matrix in a medium under the condition of Markovianity and quantum Brownian motion. In this section, we also relate the transport coefficients κand γ with the quarkonium in-medium width and mass shift respectively. In Sec. IV, we use existing lattice data to assign numerical values to both coefficients. We conclude in Sec. V. Technical details can be found in the Appendix. II. THEORETICAL BACKGROUND A. Heavy quarkonia via Lindblad equation: κand γ A recent body of work has sought to model heavy quarkonium evolution in the medium formed in heavy-ion collisions using the formalism of open quantum systems [2,3,5–11]. For a general review of open quantum systems, we direct the reader to [12]. Quarkonium serves as the *[email protected] †[email protected] ‡[email protected] §v[email protected] 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 SCOAP3. PHYSICAL REVIEW D 100, 054025 (2019) 2470-0010=2019=100(5)=054025(10) 054025-1 Published by the American Physical Society system, and the medium, which can be a QGP, as the environment. The system, quarkonium, is characterized by at least three energy scales: the mass Mof the heavy quark, the inverse of the Bohr radius, a0, and the binding energy E. These energy scales, quarkonium being a nonrelativistic bound state, are hierarchically ordered: M≫1=a0≫E. We identify the inverse of Ewith the intrinsic timescale of the system:τS∼1=E. The environment, the medium, may be characterized by several energy scales. We will assume just one single energy scale, πT. We identify the inverse of πTwith the correlation time of the environment:τE∼1=ðπTÞ. If the medium is in thermal equilibrium, or locally in thermal equilibrium, we may understand Tas the temperature. Since we do not exploit any further, possible, hierarchy among the energy scales of the medium, we are indeed considering that the medium may be strongly coupled. For instance, if the medium is a strongly coupled QGP in local thermal equilibrium, then πTis of the same order as the Debye mass, mD∼gT, and of the same order as the magnetic screening mass, mM∼g2T. The evolution of the system in the environment is characterized by a relaxation time τR. We assume that the quarkonium is Coulombic, which applies to the charmonium and bottomonium ground states. This requires that 1 a0 ≫πT; ΛQCD;ð1Þ and that in-medium and nonperturbative corrections to the Coulomb potential are subleading.1Furthermore, we also assume that πT≫E: ð2Þ Under these assumptions the relaxation time is given by the inverse of the self-energy diagram shown in Fig. 1, τR∼ 1 Σs ∼ 1 a2 0ðπTÞ3:ð3Þ In the case of a weakly coupled medium, the relaxation time may be enhanced by a factor 1=g2ðTÞ,gbeing the QCD gauge coupling. Under the Coulombic assumption (1), it follows that τR≫τE;ð4Þ which is a necessary condition for the system to be insensitive to the initial condition of the environment and, therefore, to show a Markovian evolution. Moreover, from (2) it follows that τS≫τE:ð5Þ This qualifies the regime of the quarkonium in the medium as quantum Brownian motion [6]. In [2,3] it has been shown that under the Coulombic, Eq. (1) or Eq. (4), and the Brownian motion assumption, Eq. (2) or Eq. (5), the evolution equation for the density matrix, ρ, of the heavy quark-antiquark system can be written in the Lindblad form [13,14]: dρ dt¼−i½H;ρþX nCnρC† n− 1 2fC† nCn;ρg;ð6Þ where His a Hermitian operator, and Cnare known as collapse operators. These operators were computed in [2,3]; we give the explicit expressions of ρ,H, and the Cnin the following: Eqs. (11) and (23) to (25) in Sec. III. The operators Hand Cnturn out to depend on only two transport coefficients, κand γ, that encode the entire inmedium dynamics. They are related to the real and imaginary parts of the heavy quarkonium self-energy, Σs, with a rigorous exposition and derivation of κand γ in the context of heavy quarkonium presented in Sec. III.In [2,3] it was further recognized that κis in fact the heavyquark momentum diffusion coefficient, while γcould be understood as its dispersive counterpart. We elaborate more on κand its role in the in-medium dynamics of heavy quarks in the following section. B. Heavy quarks via Langevin equation: κ The heavy-quark momentum diffusion coefficient κis an object of great interest in the literature [15–21] as it affects the momentum distribution of heavy flavor mesons measured in several experimental facilities [22–25]. It is a key FIG. 1. A diagrammatic representation of the leading order color-singlet self-energy diagram, Σs, in pNRQCD. Single lines represent quark-antiquark color-singlet propagators, double lines quark-antiquark color-octet propagators, curly lines gluons, and crossed circles chromoelectric dipole vertices. 1The temperature of the medium formed in a heavy-ion collision at the LHC ranges from approximately 475 MeV down to the freeze-out temperature giving a maximum of πTof approximately 1.5 GeV in the initial stages of the collision. In Sec. IV, for the ϒð1SÞstate, we calculate 1=a0≈1.5GeV. As the medium expands rapidly in the initial stages of the thermal evolution, it quickly cools to lower temperatures, and for the lowest lying bottomonium states, we expect Eq. (1) to hold for all but the earliest times. As the J=ψhas a significantly larger radius than the ϒð1SÞ, this relation will hold at lower temperatures. In Sec. IV for the J=ψ, we calculate 1=a0≈0.84 GeV giving a range of validity up to temperatures of approximately 250 MeV. NORA BRAMBILLA et al. PHYS. REV. D 100, 054025 (2019) 054025-2 component in understanding the heavy-quark diffusion in a thermal medium in the framework of the Langevin equations. Specifically, for a heavy quark of mass Min a thermal medium at a temperature T, with M≫T, the momentum of the heavy quark changes little over the characteristic timescale of the plasma due to random interactions with the medium constituents [15]. This slow evolution due to uncorrelated interactions with the medium is described by the Langevin equations: dpi dt¼−ηDpiþξiðtÞ;hξiðtÞξjðt0Þi ¼ κδijδðt−t0Þ; ð7Þ where piis the momentum of the heavy quark, ηDis the drag coefficient, and ξiencodes the random, uncorrelated interactions of the quark with the medium. Demanding that the system approaches thermal equilibrium entails an Einstein relation between the drag coefficient and the heavy-quark momentum diffusion coefficient, i.e., ηD¼ κ=ð2MTÞ. We thus see that the dynamics of the heavy quark in the thermal medium are governed by a single transport coefficient, namely the heavy-quark momentum diffusion coefficient κ. C. Determining κand γ In spite of the relevance of κin the theoretical description of heavy-quark diffusion in a thermal medium, its calculation has proven arduous and a precise determination elusive; for recent reviews see [26–29]. Calculations of κ require a number of assumptions on the dynamics and initial conditions of the medium along with its evolution and interaction with the heavy quark. Comparison with data would then allow for a discrimination among different assumptions and models. Large experimental uncertainties combined with subtle interactions among different assumptions have complicated the attempts to fix κreliably. The heavy-quark momentum diffusion coefficient may be also determined by means of lattice simulations. While the extraction from the spectral function of current-current correlators has turned out to be very difficult [30], more recently, κhas been related to the spectral function of the chromoelectric field correlator hgEa;iðt; 0ÞgEa;ið0;0Þi,ρel, in thermal QCD [16,18]. For definitions and details, see the Appendix. The relation reads κ¼T 6Nc lim ω→0 ρelðωÞ ω;ð8Þ which constrains κto be positive. Nc¼3is the number of colors. Equation (8) has allowed for determinations of κ, so far on quenched lattices in thermal QCD for temperatures between Tcand 2Tc,Tcbeing the crossover temperature to the QGP [19,20]. Finally, an analytic, perturbative estimate of κup to next-to-leading order in the hard-thermal-loop effective theory appears to suffer from poor convergence [17]. In this paper, we determine κfrom the thermal decay width of a heavy quarkonium in a strongly coupled medium [2,3]. This determination uses a different observable, the quarkonium thermal width, a different set of assumptions, Eqs. (1) (Coulombic bound state) and (2) (quantum Brownian motion), and a different source of data, 2 flavor lattice QCD data from [31] and 2þ1 flavor lattice QCD data from [32]. Therefore, it is an independent determination with different systematic uncertainties, potentially competitive with other determinations. We estimate our main sources of systematic uncertainties in this determination of κto be higherorder corrections inherent in our effective field theory approach and the systematic uncertainties inherited from the specific lattice data used in our calculations. As discussed in Sec. II A, for bottomonium, we expect our hierarchy of scales in Eq. (1) to be fulfilled and these higher-order corrections to be small. In contrast to the theoretical understanding of the role of κin the dynamics of heavy quarks in a thermal medium and the progress towards its calculation, comparatively little is known about γ. Since γmay be understood as a correction to the heavy quark-antiquark potential, no similar object arises in the description of the in-medium heavy-quark dynamics. A proper definition relates γto the chromoelectric field correlator in such a way that it may be considered the dispersive counterpart of κ[2,3]. Like κ, γcan be written in thermal QCD in terms of the chromoelectric spectral function ρel: γ¼− 1 3NcZ∞ 0 dω 2π ρelðωÞ ω:ð9Þ A derivation of Eq. (9) is in the Appendix. Differently from κ, however, the coefficient γis a function of ρelðωÞ=ωover the whole spectrum of frequencies. Since ρelðωÞ∼ω3for large frequencies, the above integral is ultraviolet divergent and needs to be regularized and renormalized. The large frequencies behavior of ρelðωÞis entirely given by the in vacuum (T¼0) contributions. These are known up to nextto-leading order [33]. From Ref. [33] it also follows that the thermal part of γis finite. Just as κis the parameter of central importance in the study of in-medium heavy quarks, κand γappear to be the parameters of central importance in the study of the quantum Brownian motion of Coulombic quarkonia in a strongly coupled medium. In this paper, taking advantage of the relation between γand the quarkonium thermal mass shift in a strongly coupled medium [2,3],we determine the thermal part of γfrom the 2þ1flavor lattice QCD data of [32]. The procedure will be similar to the one used to extract κ, as well as the underlying assumptions (1) and (2). TRANSPORT COEFFICIENTS FROM IN-MEDIUM …PHYS. REV. D 100, 054025 (2019) 054025-3 III. QUARKONIUM IN THE QUANTUM BROWNIAN REGIME In [3], a set of master equations governing the time evolution of heavy quarkonium in a medium were derived. The equations follow from assuming the inverse Bohr radius of the quarkonium to be greater than the energy scale of the medium, Eq. (1), and model the quarkonium as evolving in the vacuum up to a time t¼t0, at which point interactions with the medium begin. The equations express the time evolution of the density matrices of the heavy quark-antiquark color-singlet, ρs, and octet states, ρo,in terms of the color-singlet and octet Hamiltonians, hs¼ p2=M −CFαs=r þ and ho¼p2=Mþαs=ð2NcrÞþ, and interaction terms with the medium, which, at order r2in the multipole expansion, are encoded in the selfenergy diagram shown in Fig. 1. These interactions account for the mass shift of the heavy quark-antiquark pair induced by the medium, its decay width induced by the medium, the generation of quark-antiquark color-singlet states from quark-antiquark color-octet states interacting with the medium and the generation of quark-antiquark color-octet states from quark-antiquark (color-singlet or octet) states interacting with the medium. The color-singlet and octet Hamiltonians, hsand ho, describe particles of mass Mand momentum pinteracting at a distance rthrough a Coulomb potential; CF¼ðN2 c−1Þ=ð2NcÞ¼4=3is the Casimir of the fundamental representation, and αs¼g2=ð4πÞis the strong coupling. The dots in our expressions of hsand ho stand for higher-order terms that are irrelevant for the present analysis. The effective field theory framework in which the nonrelativistic heavy quark-antiquark dynamics can be systematically described in terms of quark-antiquark color-singlet and octet fields, whose interactions with the medium are organized in powers of 1=M and r, is potential nonrelativistic QCD (pNRQCD) [34–36]. The leading order interaction between a heavy quark-antiquark field and the medium is encoded in pNRQCD in a chromoelectric dipole interaction, which appears at order r=M0in the effective field theory Lagrangian. Further assuming that any energy scale in the medium is larger than the heavy quark-antiquark binding energy,2 Eq. (2), leads to the following evolution equations [3]: dρ dt¼−i½H;ρþX nm hnmLn iρLm i †− 1 2fLm i †Ln i;ρg; ð10Þ where ρ¼ρs0 0ρo;ð11Þ H¼hs0 0hoþrirj 2 ˜γijðtÞ10 0N2 c−2 2ðN2 c−1Þ;ð12Þ L0 i¼ri00 01 ;ð13Þ L1 i¼rj 2ð˜ κijðtÞþi˜γijðtÞÞ00 0N2 c−4 2ðN2 c−1Þ;ð14Þ L2 i¼ri01 10 ;ð15Þ L3 i¼rj 2ð˜ κijðtÞþi˜γijðtÞÞ01 N2 c−1 10 :ð16Þ Despite the fact that hnm are the elements of a matrix, h¼0 B B B @ 0100 1000 0001 0010 1 C C C A ;ð17Þ which is not positive definite, it is straightforward to show that after a redefinition of Ln ithe eigenvectors of h associated to negative eigenvalues, which are proportional to L1 i−L0 iand L3 i−L2 i, can be set to zero. Hence, according to Refs. [13,14], we can map the above evolution equations into the Lindblad form (6). We emphasize that the evolution equations in this form hold for a Coulombic bound state in quantum Brownian motion in a medium. We do not require the medium to be weakly coupled. The tensors ˜ κijðtÞand ˜γijðtÞhave a field theoretical definition: ˜ κijðtÞ¼ 1 2NcZt t0 dt0hfgEa;iðt; 0Þ;gE a;jðt0;0Þgi;ð18Þ ˜γijðtÞ¼−i 2NcZt t0 dt0h½gEa;iðt; 0Þ;gE a;jðt0;0Þi;ð19Þ where h  i stands for the in-medium average, the curly brackets signify anticommutator, the square ones commutator, and Eis the chromoelectric field. In the above expressions, the chromoelectric field has to be understood as Ω†×ðusual EÞ×Ω, where Ωis a Wilson line going from −∞ to t:Ω¼exp ½−ig Rt −∞ dsA0ðs; 0Þ. The Wilson lines guarantee that the definitions of ˜ κijðtÞand ˜γijðtÞare gauge invariant. The tensors ˜ κijðtÞand ˜γijðtÞmay be related 2Note that, in pNRQCD, quark-antiquark pairs in a color-octet configuration have energy p2=M, where pis a relative momentum of the same order as 1=a0. Hence the typical energy of the color-octet pair is of order E, which is the binding energy of the color-singlet pair. NORA BRAMBILLA et al. PHYS. REV. D 100, 054025 (2019) 054025-4 to the real and imaginary parts of the quark-antiquark colorsinglet self-energy diagram shown in Fig. 1: ΣsðtÞ¼rirj1 2NcZt t0 dt0hgEa;iðt; 0ÞgEa;jðt0;0Þi ¼rirj 2½˜ κijðtÞþi˜γijðtÞ:ð20Þ If the medium is isotropic, then ˜ κijðtÞ¼δij ˜ κðtÞand ˜γijðtÞ¼δij ˜γðtÞ, where ˜ κðtÞ¼ 1 6NcZt t0 dt0hfgEa;iðt; 0Þ;gE a;iðt0;0Þgi;ð21Þ ˜γðtÞ¼−i 6NcZt t0 dt0h½gEa;iðt; 0Þ;gE a;iðt0;0Þi:ð22Þ The Lindblad equation describing the time evolution of the density matrix has then a particularly simple form [2,3],as the Hermitian operator His given by H¼hs0 0hoþr2 2 ˜γðtÞ10 0N2 c−2 2ðN2 c−1Þ;ð23Þ and the collapse operators Ciby: C0 i¼ffiffiffiffiffiffiffiffiffiffiffiffiffiffi ˜ κðtÞ N2 c−1 sri01 ffiffiffiffiffiffiffiffiffiffiffiffiffiffi N2 c−1 p0;ð24Þ C1 i¼ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ðN2 c−4Þ˜ κðtÞ 2ðN2 c−1Þ sri00 01 :ð25Þ At short times after the formation of the medium, t≳t0, ˜ κðtÞscales like T4ðt−t0Þand the thermal part of ˜γðtÞlike T5ðt−t0Þ2. In the opposite, large time limit, t−t0is the largest timescale in the problem. At this point, it is convenient to assume that the chromoelectric correlators appearing in (21) and (22) are, at least approximately, time translation invariant:hEa;iðt;0ÞEa;iðt0;0Þi¼hEa;iðt−t0;0Þ× Ea;ið0;0Þi. This is the case, for instance, at thermal equilibrium or close to it, if the variation in time of the temperature is slow. In the large time limit and under the assumption of (approximate) time translation invariance ˜ κðtÞand ˜γðtÞapproach the asymptotic values κ¼˜ κð∞Þand γ¼˜γð∞Þrespectively: κ¼1 6NcZ∞ 0 dthfgEa;iðt; 0Þ;gE a;ið0;0Þgi;ð26Þ γ¼−i 6NcZ∞ 0 dth½gEa;iðt; 0Þ;gE a;ið0;0Þi:ð27Þ The quantity κis the heavy-quark momentum diffusion coefficient. The above form of the heavy-quark momentum diffusion coefficient was first derived in [16] (see also [18]) in the context of the diffusion of a heavy quark in a thermal medium according to the Langevin equations (7). The coefficient κcan also be written as the real part of the time ordered correlator 1=ð6NcÞRþ∞ −∞ dthTgEa;iðt;0ÞgEa;ið0;0Þi; γis then its imaginary part. In the large time limit, for an isotropic medium and under approximate time translation invariance, the color-singlet self-energy (20) becomes Σs¼r2 2ðκþiγÞ:ð28Þ This allows us to write r2κ¼ΣsþΣ† s¼−2Imð−iΣsÞ;ð29Þ r2γ¼−iΣsþiΣ† s¼2Reð−iΣsÞ;ð30Þ and eventually relate κand γto the quarkonium in-medium width, Γ, and the in-medium mass shift, δM.For1S Coulombic quarkonium states, these relations read [2,3] Γð1SÞ¼3a2 0κ;ð31Þ δMð1SÞ¼3 2a2 0γ;ð32Þ where 3a2 0is the expectation value of r2on a 1SCoulombic bound state. The Bohr radius is a0¼2=ðMCFαsÞ. Γand δMdo not contain, by definition, in vacuum contributions. Also κ, as defined in (26), does not contain in vacuum contributions, reflecting the fact that energy conservation prohibits the decay of a heavy quark-antiquark color singlet into a heavy quark-antiquark color octet in vacuum. In contrast, γ, as defined in (27), does contain in vacuum contributions. Hence, Eq. (32) relates δMð1SÞto γ subtracted of its vacuum (T¼0) part. The coefficient γ should be understood in this subtraction scheme in Eq. (32) and in the next section. Now that we have explicit relations for κand γin terms of a0,Γð1SÞ, and δMð1SÞ, we can proceed to extract κand γfrom available lattice estimates of Γð1SÞand δMð1SÞ. IV. RESULTS AND COMPARISONS Equations (31) and (32) fix the ratio γ=κto be γ κ¼2δMð1SÞ Γð1SÞ;ð33Þ which may turn out to be useful once both δMð1SÞand Γð1SÞare reliably determined, for the ratio does not depend on the Bohr radius. The quantities δMð1SÞand Γð1SÞ cannot be accessed by experiments, as in heavy-ion collisions the quarkonium decays, at a late time, in the TRANSPORT COEFFICIENTS FROM IN-MEDIUM …PHYS. REV. D 100, 054025 (2019) 054025-5 vacuum. These quantities can instead by computed by lattice QCD, with the thermal mass shift, δMð1SÞ, clearly in a more reliable fashion than the thermal width, Γð1SÞ. In order to determine γand κfrom δMð1SÞand Γð1SÞ using Eqs. (32) and (31), we need to calculate a0. Since the system is Coulombic, we can do it by solving the selfconsistency equation a0¼2 MCFαsð1=a0Þ;ð34Þ where αsð1=a0Þis the strong coupling evaluated at the scale 1=a0. For the bottom and charm masses, we take M¼ Mb¼4.78 GeV and M¼Mc¼1.67 GeV, respectively. These are the central values for the pole masses quoted by the Particle Data Group [37].3We solve for a0using the one-loop, 3-flavor running of αswith ΛMS ¼332 MeV, also from Ref. [37]. We take the running of αsat one loop for consistency with the fact that the radius of a 1S Coulombic bound state is given by Eq. (34) only at leading order. For the bottomonium ground state we obtain a0¼0.67 GeV−1¼0.13 fm, while for the charmonium ground state we obtain a0¼1.19 GeV−1¼0.23 fm. For the above choice of heavy-quark masses, the ϒð1SÞbinding energy is E¼Mðϒð1SÞÞ −2Mb¼−0.1GeV, and the J=ψbinding energy is E¼MðJ=ψÞ−2Mc¼−0.24 GeV.4 We evaluate γ=T3using the thermal mass shifts computed for the J=ψat T¼251 MeV and for the ϒð1SÞat T¼407 MeV and T¼251 MeV in the 2þ1flavor lattice simulation of Ref. [32]. After rescaling for the mass, see Footnote 3, the mass shifts are ð−85  29ÞMeV for the J=ψat T¼251 MeV, and ð−48  16ÞMeV and ð−30 12ÞMeV for the ϒð1SÞat T¼ 407 MeV and T¼251 MeV, respectively.5The results for γare shown by the three first entries (black bars) of Fig. 2. For the J=ψ, we have 1=a0¼0.84 GeV, which fulfills (somewhat marginally) the hierarchy of Eq. (1) at T¼251 MeV, since π×ð251 MeVÞ¼0.79 GeV. It also fulfills the condition (2). In the case of the ϒð1SÞ,1=a0¼ 1.5GeV and both conditions (1) and (2) are fulfilled at both temperatures T¼251 MeV and T¼407 MeV although at the lower temperature more clearly than at the higher one. The two conditions (1) and (2) guarantee that the J=ψ and ϒð1SÞremain Coulombic also in the medium and that their motion through the medium is a quantum Brownian one. Because of this, we consider all three extractions of γ, from the J=ψ(one temperature) and the ϒð1SÞ(two temperatures), reliable, as they are consistent with our assumptions. We take their range as an estimate of γ=T3for 251 ≲T≲407 MeV: −3.8≲γ T3≲−0.7:ð35Þ It is particularly significant to see that the extraction of γ=T3from the J=ψat T¼251 MeV overlaps perfectly with the extraction of γ=T3from the ϒð1SÞat the same temperature. This shows that, as expected, γ=T3depends only on the temperature, while it does not depend on the quarkonium state. Concerning the temperature dependence, the extraction of γ=T3from the ϒð1SÞat T¼407 MeV could suggest that −γ=T3tends towards smaller values at higher temperatures. The last two entries (blue bars) of Fig. 2refer to γdetermined from the leading order expression of the thermal mass shift computed in Ref. [38] and reported in Refs. [2,3]. We see that perturbation theory also gives a negative value for γ, as our nonperturbative estimate above. Moreover, there is a partial overlap between the perturbative result at the highest temperature, where a weak-coupling treatment is expected FIG. 2. The first three entries (black bars) show γ=T3as obtained from Eq. (32) using lattice data of Ref. [32] for the thermal mass shift of the J=ψand of the ϒð1SÞat two different temperatures. The error bars account for the lattice uncertainties only. The last two entries (blue bars) provide γ=T3from the perturbative, leading order, expression of the thermal mass shift with the strong coupling computed at πtimes the two different temperatures 251 and 407 MeV. We assign a 50% uncertainty to these results. The gray band gives our final range for γ=T3, see text. 3Since in Ref. [32] the masses Mb¼4.65 GeV and Mc¼ 1.275 GeV were used, we have rescaled their values for δMð1SÞ and Γð1SÞby ð4.65=4.78Þ2in the bottomonium case and by ð1.275=1.67Þ2in the charmonium case. Similarly, since in Ref. [31] the mass Mb¼5GeV was used, we have rescaled their value for Γð1SÞby ð5=4.78Þ2. The coefficients γand κare mass independent, so the choice of the mass should not affect them. There is, however, a residual dependence due to having truncated the expressions in the right-hand sides of Eqs. (31) and (32) at leading order in the various expansions underlying the effective field theory. We have checked that this residual dependence is, indeed, well accounted for by the quoted errors. 4Defining the binding energy as the Coulombic Bohr level, E¼−1=ðMa2 0Þ, changes the numerical value (to E¼ −0.46 GeV for the ϒð1SÞand E¼−0.42 GeV for the J=ψ), but not the hierarchy of energy scales. Hence, the following arguments and extractions of γand κ, which only depend on that hierarchy, would remain unchanged. 5Consistently with the Coulombic assumption, the J=ψand ϒð1SÞin vacuum binding energies are negative and at least a factor 2 larger than the corresponding thermal mass shifts. NORA BRAMBILLA et al. PHYS. REV. D 100, 054025 (2019) 054025-6 to work better, and the range given in Eq. (35).Itis possible, however, that higher-order corrections will spoil the leading order result at the temperatures considered here, as it is the case for the weak-coupling expression of κ. Finally, we remark that a small and negative value of γ=T3 is also phenomenologically favored by the comparison of the ϒð1SÞnuclear modification factor, as computed from the Lindblad equation (6),(23),(24),(25), with the most recent CMS data [2,3,39]. Concerning κ, the available lattice data for the thermal width are less precise than those available for the mass shift. The width given in the 2þ1flavor lattice simulation of Ref. [32] is preliminary and should be understood as a lower bound rather than the width itself [40]. We combine this lower bound with the somewhat older 2 flavor lattice results from [31] that supply an upper bound to the thermal width. More specifically, we take the ϒð1SÞwidth at T¼407 MeV from [32] [Γð1SÞ≈22.3MeV, after rescaling for the mass] as a lower bound and the highest temperature estimate of Γð1SÞfrom Fig. 5 of Ref. [31] [Γð1SÞ=T ≈1.1, after rescaling for the mass, for T=Tc≈2 with Tc≈220 MeV] as an upper bound. We obtain 0.24 ≲κ T3≲4.2:ð36Þ The above range is the first entry (black bar) in Fig. 3.We note that we do not plot an estimate of κobtained from the thermal width of the J=ψas Ref. [31] performed measurements only on bottomonium states; this leaves us with only a lower bound from the measurement of the width of the J=ψat T¼251 MeV performed in Ref. [32]. We find this gives a lower limit of κ=T3≳0.235 0.208 (with purely statistical uncertainties); we note the agreement with the lower bound of κ=T3≳0.24 obtained from the ϒð1SÞat T¼407 MeV. The second entry (brown bar) in Fig. 3 reports the result, 1.8≲κ=T3≲3.4, of the lattice study done in a pure SU(3) plasma at T≈1.5Tcwith Tc≈ 313 MeV in Ref. [20]. Indirect bounds can also be placed on κfrom the experimental measure of the D-meson azimuthal anisotropy coefficient v2. In particular, measurements from the ALICE [24] and STAR [25] Collaborations can be compared with theoretical models to place bounds on the heavy-quark spatial diffusion coefficient. Relating the heavy-quark spatial diffusion coefficient Dto the heavy-quark momentum diffusion coefficient (κ=T3¼ 2=DT), we find that ALICE data give 1.8≲κ=T3≲8.4 at T¼Tcand STAR data give 1.0≲κ=T3≲6.3for Tc≲T≲2Tc. The third and fourth entries (green bars) in Fig. 3represent these two bounds, respectively. Calculations of the nuclear modification factor RAA and the elliptic flow v2from phenomenological models also place bounds on the momentum diffusion coefficient; for a collection of recent results, see [27]. The results of this reference give 3.1≲κ=T3≲6.3at T≈155–160 MeV and are the fifth entry (red bar) in Fig. 3.6The sixth entry (blue bar) in Fig. 3follows from computing κin weak-coupling perturbation theory at T¼407 MeV. In the perturbative expression we include the complete order g4contribution and the order g5term 7=ð48π2ÞCFNcg4ðmD=TÞ(this is, truncated at order g5, what is called leading order expression in [17]); we note that the order g4contribution alone would give unphysical negative values of κfor realistic couplings. The perturbative expression of κ=T3, in the above sense, gives 2≲κ=T3≲6by assigning a 50% uncertainty. This uncertainty may, however, be underestimated since the complete g5correction is known and very large [17]. In fact it may increase the leading order result even by an order of magnitude under some circumstances [17,19], which obviously questions the reliability of a naive weakcoupling expansion for κat the considered temperatures. All determinations of κ=T3are consistent with each other, and, in particular, with the range presented in Eq. (36). This is noteworthy as these determinations are very different, some of them do not even rely on full QCD, FIG. 3. The first entry (black bar) shows κ=T3as obtained from Eq. (31) using lattice data of Refs. [31,32] for the upper and lower bounds of the thermal decay width of the ϒð1SÞ. The second entry (brown bar) reports the (quenched) lattice estimate of Ref. [20]. The third and fourth entries (green bars) are the determinations based on the ALICE [24] and STAR [25] measurements of the D-meson azimuthal anisotropy coefficient v2, respectively. The fifth entry (red bar) is the determination based on a range of values of the heavy-quark spatial diffusion coefficient obtained in [27] from an analysis of phenomenological models. The sixth entry (blue bar) is a perturbative result (see text) with the strong coupling computed at the scale π×ð407 MeVÞ¼1.28 GeV. We assign a 50% uncertainty to it. 6We note that the phenomenological models used in the last three extractions may be more complicated than the Langevin equation given in Eq. (7). Furthermore, these extractions are not based on a nonrelativistic expansion and may include into κ dynamics that in our effective field theory framework occur as relativistic corrections of higher order in 1=M; these corrections scale as the square of the relative heavy-quark velocity in the quarkonium possibly contributing to a systematic effect by up to 10% for bottomonium and 30% for charmonium when comparing with the mass-independent extractions from the chromoelectric correlator. TRANSPORT COEFFICIENTS FROM IN-MEDIUM …PHYS. REV. D 100, 054025 (2019) 054025-7 and have been obtained at different temperatures. Indeed, κ=T3, as well as γ=T3, does depend on the temperature. The lattice study of [19] could suggest κ=T3assumes lower values at higher temperatures, which seems consistent with studies from phenomenological models [41]. Finally, we mention that the same relationship between the width of the state and κgiven by Eq. (31) has also been exploited in [42], but in the opposite direction, to supply an estimate of the quarkonium thermal width from κ. V. CONCLUSION The Lindblad equation describing the heavy quarkonium evolution in the hot medium created at the early stages of high-energy heavy-ion collisions requires, under some conditions, only two parameters to describe the interaction of the heavy quark-antiquark pair with the medium [2,3,43]. The conditions are that the quarkonium is a Coulombic bound state, which holds when Eq. (1) or Eq. (4) are satisfied, and that its motion in the medium is a quantum Brownian motion, which holds when Eq. (2) or Eq. (5) are satisfied. The two parameters are the heavyquark momentum diffusion coefficient, κ, which crucially enters also the Langevin equation describing the heavyquark diffusion in the medium, whose field theoretical definition is in Eq. (26), and its dissipative counterpart, the coefficient γ, whose field theoretical definition is in Eq. (27). No assumption is required on the nature of the medium. In this paper, we have estimated γ, using its relation to the in-medium quarkonium mass shift and the 2þ1flavor lattice data of Ref. [32]. The result is given in Eq. (35). This is the first nonperturbative determination of γ. Its sign is consistent with the weak-coupling, leading order thermal mass shift, which, however, is affected by large uncertainties. We have also computed κ, using its relation to the inmedium quarkonium decay width and the two different sets of lattice data [31,32]. The data are not precise enough to pin down a narrow range of κ, nevertheless they allow us to establish an upper and a lower limit for this transport coefficient. They are given in Eq. (36). The range of κis consistent with other determinations, see Fig. 3. Once more precise lattice data will become available, this method has the potential to provide a competitive determination of the heavy-quark momentum diffusion coefficient. Already now, it gives a range for κthat is based on 2þ1and 2 flavor lattice data, while current lattice determinations rely on pure SU(3) gauge theory simulations [19,20]. From the above considerations, it is clear that this work calls for several further lattice analyses. Concerning γ,it would be important to have a direct evaluation based on the spectral function of the chromoelectric field correlator and Eq. (9) in a given renormalization scheme. Concerning κ, besides direct full QCD determinations based on Eq. (8), this work aims also at motivating further lattice computations of the in-medium quarkonium decay width. Finally, we remark that the logic of the present work may be reversed and the agreement noticed between the present determination of κand other independent ones used to support our starting assumptions on the nature of the studied heavy quarkonia and of their diffusion in the medium. Indeed, our determination of κfrom the in-medium decay width and the explicit computation of the relevant energy scales provide evidence that at least the ϒð1SÞis a Coulombic bound state that propagates with a quantum Brownian motion in the medium formed, and at the temperature attained, by present day heavy-ion colliders. ACKNOWLEDGMENTS The authors thank Jacopo Ghiglieri, Viljami Leino, Peter Petreczky and Alexander Rothkopf for discussions and a reading of the manuscript. N. B. and A. V. thank Yukinao Akamatsu and Joan Soto for an inspiring discussion at the Institute for Nuclear Theory while attending the program “Multi-Scale Problems Using Effective Field Theories.” This work was funded by the Bundesministerium für Bildung und Forschung Project No. 05P2018 and by the DFG cluster of excellence “Universe”(www.universecluster.de). The work of M. A. E. was supported by the Academy of Finland Project No. 297058, by Ministerio de Ciencia e Innovacion of Spain under Project No. FPA201783814-P and Maria de Maetzu Unit of Excellence MDM20160692, by Xunta de Galicia and FEDER. We thank the Munich Institute for Astroand Particle Physics (MIAPP) of the DFG cluster of excellence “Universe”during which activity “Probing the Quark-Gluon Plasma with Collective Phenomena and Heavy Quarks”the idea of this paper was first discussed. Note added.—While finishing this paper, we became aware of the preprint of A. M. Eller, J. Ghiglieri and G. D. Moore “Thermal quarkonium mass shift from Euclidean correlators”[44] that proposes a way to determine γdirectly from a proper Euclidean correlator. We thank the authors for sharing with us their results prior to publication. APPENDIX: SPECTRAL FUNCTION AND γ In this Appendix, we define the spectral function of the chromoelectric correlator, ρel, and derive expression (9) for γ. We start considering the real time chromoelectric correlator hgEa;iðt; 0ÞgEa;ið0;0Þi;ðA1Þ where the chromoelectric fields are understood in the convention specified after Eq. (18), i.e., with Wilson lines attached. The Wilson lines make the correlator gauge invariant. The medium average h  i is normalized by NORA BRAMBILLA et al. PHYS. REV. D 100, 054025 (2019) 054025-8