Première modélisation numérique de l’évolution des quarks lourds dans le QGP formé dans les collisions d’ions lourds au RHIC et au LHC, couplé à la production de quarkonia