PaperPanorama
FAQ

Nuclear Theory·nucl-th

Wednesday·July 14, 2021

5 papers—3 primary·2 cross-listed

All5Primary3Cross-lists2
  1. 04

    QTRAJ 1.0: A Lindblad equation solver for heavy-quarkonium dynamics

    Hisham Ba Omar🇺🇸 · Miguel Ángel Escobedo🇪🇸 · Ajaharul Islam🇺🇸 · Michael Strickland🇺🇸 · Sabin Thapa🇺🇸 · Peter Vander Griend🇩🇪 · Johannes Heinrich Weber🇩🇪

    We introduce an open-source package called QTraj that solves the Lindblad equation for heavy-quarkonium dynamics using the quantum trajectories algorithm. The package allows users to simulate the suppression of heavy-quarkonium states using externally-supplied input from 3+1D hydrodynamics simulations. The code uses a split-step pseudo-spectral method for updating the wave-function between jumps, which is implemented using the open-source multi-threaded FFTW3 package. This allows one to have manifestly unitary evolution when using real-valued potentials. In this paper, we provide detailed documentation of QTraj 1.0, installation instructions, and present various tests and benchmarks of the code.

    ↳ physics.comp-phhep-phnucl-thquant-phComput.Phys.Commun.(2022)·39 citations
  2. 05

    Bottomonium production in heavy-ion collisions using quantum trajectories: Differential observables and momentum anisotropy

    Nora Brambilla🇩🇪 · Miguel Ángel Escobedo🇪🇸 · Michael Strickland🇺🇸 · Antonio Vairo🇩🇪 · Peter Vander Griend🇩🇪 · Johannes Heinrich Weber🇺🇸

    We report predictions for the suppression and elliptic flow of the , , and 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 TeV Pb-Pb collisions by the ALICE, ATLAS, and CMS collaborations.

    ↳ hep-phnucl-exnucl-thphysics.comp-phPRD(2021)·71 citations

Affiliations

first authorsco-authorsvia INSPIRE