Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBottomonium suppression and elliptic flow using Heavy Quarkonium Quantum Dynamics
139
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We introduce a framework called Heavy Quarkonium Quantum Dynamics (HQQD) which can be used to compute the dynamical suppression of heavy quarkonia propagating in the quark-gluon plasma using real-time in-medium quantum evolution. Using HQQD we compute large sets of real-time solutions to the Schr{o}dinger equation using a realistic in-medium complex-valued potential. We sample 2 million quarkonia wave packet trajectories and evolve them through the QGP using HQQD to obtain their survival probabilities. The computation is performed using three different HQQD model parameter sets in order to estimate our systematic uncertainty. After taking into account final state feed down we compare our results to existing experimental data for the suppression and elliptic flow of bottomonium states and find that HQQD predictions are good agreement with available data for $R_{AA}$ as a function of $N_{rm part}$ and $p_T$ collected at $sqrt{s_{rm NN}} =$ 5.02 TeV. In the case of $v_2$ for the various states, we find that the path-length dependence of $Upsilon(1s)$ suppression results in quite small $v_2$ for $Upsilon(1s)$. Our prediction for the integrated elliptic flow for $Upsilon(1s)$ in the $10{-}90$% centrality class, which now includes an estimate of the systematic error, is $v_2[Upsilon(1s)]$ = 0.003 $pm$ 0.0007 $pm,^{0.0006}_{0.0013}$. We also find that, due to their increased suppression, excited bottomonium states have a larger elliptic flow. Based on this observation we make predictions for $v_2[Upsilon(2s)]$ and $v_2[Upsilon(3s)]$ as a function of centrality and transverse momentum.
rate research
Read More
We compute the suppression and elliptic flow of bottomonium using real-time solutions to the Schr{o}dinger equation with a realistic in-medium complex-valued potential. To model the initial production, we assume that, in the limit of heavy quark masses, the wave-function can be described by a lattice-smeared (Gaussian) Dirac delta wave-function. The resulting final-state quantum-mechanical overlaps provide the survival probability of all bottomonium eigenstates. Our results are in good agreement with available data for $R_{AA}$ as a function of $N_{rm part}$ and $p_T$ collected at $sqrt{s_{rm NN}} =$ 5.02 TeV. In the case of $v_2$ for the various states, we find that the path-length dependence of $Upsilon(1s)$ suppression results in quite small $v_2$ for $Upsilon(1s)$. Our prediction for the integrated elliptic flow for $Upsilon(1s)$ in the $10{-}90$% centrality class is $v_2[Upsilon(1s)] = 0.0026 pm 0.0007$. We additionally find that, due to their increased suppression, excited bottomonium states have a larger elliptic flow and we make predictions for $v_2[Upsilon(2s)]$ and $v_2[Upsilon(3s)]$ as a function of centrality and transverse momentum. Similar to prior studies, we find that it is possible for bottomonium states to have negative $v_2$ at low transverse momentum.
Many prior studies of in-medium quarkonium suppression have implicitly made use of an adiabatic approximation in which it was assumed that the heavy quark potential is a slowly varying function of time. In the adiabatic limit, one can separately determine the in-medium breakup rate and the medium time evolution, folding these together only at the end of the calculation. In this paper, we relax this assumption by solving the 3d Schrodinger equation in real-time in order to compute quarkonium suppression dynamically. We compare results obtained using the adiabatic approximation with real-time calculations for both harmonic oscillator and realistic complex heavy quark potentials. Using the latter, we find that, for the Upsilon(1s), the difference between the adiabatic approximation and full real-time evolution is at the few percent level, however, for the Upsilon(2s), we find that the correction can be as large as 18% in low temperature regions. For the J/Psi, we find a larger difference between the dynamical evolution and the adiabatic approximation, with the error reaching approximately 36%.
The strong suppression of bottomonia production in ultra-relativistic heavy-ion collisions is a smoking gun for the creation of a deconfined quark-gluon plasma (QGP). In this proceedings contribution, I review recent work that aims to provide a more comprehensive and systematic understanding of bottomonium dynamics in the QGP through the use of pNRQCD and an open quantum systems approach. This approach allows one to evolve the heavy-quarkonium reduced density matrix, taking into account non-unitary effective Hamiltonian evolution of the wave-function and quantum jumps between different angular momentum and color states. In the case of a strong coupled QGP in which E << T,m_D << 1/a_0, the corresponding evolution equation is Markovian and can therefore be mapped to a Lindblad evolution equation. To solve the resulting Lindblad equation, we make use of a stochastic unraveling called the quantum trajectories algorithm and couple the non-abelian quantum evolution to a realistic 3+1D viscous hydrodynamical background. Using a large number of Monte-Carlo sampled bottomonium trajectories, we make predictions for bottomonium R_AA and elliptic flow as a function of centrality and transverse momentum and compare to data collected by the ALICE, ATLAS, and CMS collaborations.
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.
We report predictions for the suppression and elliptic flow of the $Upsilon(1S)$, $Upsilon(2S)$, and $Upsilon(3S)$ as a function of centrality and transverse momentum in ultra-relativistic heavy-ion collisions. We obtain our predictions by numerically solving a Lindblad equation for the evolution of the heavy-quarkonium reduced density matrix derived using potential nonrelativistic QCD and the formalism of open quantum systems. To numerically solve the Lindblad equation, we make use of a stochastic unraveling called the quantum trajectories algorithm. This unraveling allows us to solve the Lindblad evolution equation efficiently on large lattices with no angular momentum cutoff. The resulting evolution describes the full 3D quantum and non-abelian evolution of the reduced density matrix for bottomonium states. We expand upon our previous work by treating differential observables and elliptic flow; this is made possible by a newly implemented Monte-Carlo sampling of physical trajectories. Our final results are compared to experimental data collected in $sqrt{s_{NN}} = 5.02$ TeV Pb-Pb collisions by the ALICE, ATLAS, and CMS collaborations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
