Bottomonium suppression in an open quantum system using the quantum trajectories method

Abstract

We solve the Lindblad equation describing the Brownian motion of a Coulombic heavy quark-antiquark pair in a strongly coupled quark-gluon plasma using the highly efficient Monte Carlo wave-function method. The Lindblad equation has been derived in the framework of pNRQCD and fully accounts for the quantum and non-Abelian nature of the system. The hydrodynamics of the plasma is realistically implemented through a 3+1D dissipative hydrodynamics code. We compute the bottomonium nuclear modification factor and compare with the most recent LHC data. The computation does not rely on any free parameter, as it depends on two transport coefficients that have been evaluated independently in lattice QCD. Our final results, which include late-time feed down of excited states, agree well with the available data from LHC 5.02 TeV PbPb collisions.

Non-Abelian electric field correlator at NLO for dark matter relic abundance and quarkonium transport

“…In this work, we consider mediuminduced bound state formation and dissociation in non-Abelian gauge theories. In order to compute these processes precisely in a proper thermal field theoretical approach, we apply the open quantum system framework, which has been intensively used to investigate quarkonium transport [56][57][58][59][60][61][62][63][64][65][66][67][68][69][70] and jet quenching [71][72][73] in the QGP. In particular, we consider a subsystem consisting of a pair of heavy particles, interacting weakly with a thermal environment, i.e., the hot plasma.…”

The imaginary part of the heavy-quark potential from real-time Yang-Mills dynamics

“…Fundamentally, the computation of the survival probability of a given bottomonium state can be cast into the framework of open quantum systems (OQS) in which there is JHEP10(2021)083 a probe (bottomonium states) and medium (light quarks and gluons). Within the OQS framework, in order to describe the in-medium evolution of bottomonium states one must trace over the medium degrees of freedom and obtain evolution equations for the reduced density matrix of the system [6][7][8][9][10][11][12][13][14][15][16][17][18][19]. In the limit that the medium relaxation time scale and the intrinsic time scale of the probe are much smaller than the probe relaxation time scale, the resulting dynamical equation for the reduced density matrix can be cast into a so-called Lindblad form [20,21].…”

In medium Langevin dynamics of heavy quarkonium

“…We expand to linear order in / the master equation derived in [12] which at order 0 in / takes the form of a Lindblad equation which was solved in Refs. [12][13][14][15]. We follow the procedure of Blaizot and Escobedo of Ref.…”

“…[11,12], is the dispersive counterpart of . In this limit, the master equation can be written in Lindblad form; previous works solved this Lindblad equation and extracted experimental observables including the nuclear modification factor and the elliptic flow 2 [13][14][15].…”