La production de quarkonia lourd est l'une des sondes les plus prometteuses pour l'exploration du QGP. Des observables tels que la forte perte d'énergie du quarkonium et sa suppression peuvent offrir une dynamique directe de boule de feu. Pour cette raison, la construction de modèles capables d'expliquer et de traiter la dynamique des quarkoniums lourds et l'évolution temporelle est cruciale pour étudier la matière dans ces conditions extrêmes. Cependant, certaines approches telles que les modèles statistiques d'équation du taux d'hadronisation ont grandement contribué à la communauté des collisions d'ions lourds. Ils ont également souligné que pour obtenir une description complètement satisfaisante du problème de la production de quarkonia lourds, une modélisation dynamique de la recombinaison tout en préservant certaines caractéristiques quantiques, en particulier lorsque nous abordons la formation d'états liés relativistes, est nécessaire. Suivant cette ligne, nos recherches ont porté sur la construction d'un modèle de coalescence dans l'espace des phases basé sur le formalisme de formation de particules composées de Remler pour explorer la formation de quarkonium lourd dans un environnement de collision d'ions lourds. Notre modèle est un formalisme multi-particulaire capable de traiter la probabilité de recombinaison de toutes les paires de quarks lourds possibles (c ou b et ses quarks antimatière respectifs) voyageant dans le QGP à un instant donné. Par le calcul d'un taux différentiel effectif (effets de dissociation et de recombinaison mis ensemble), nous pouvons suivre la probabilité tout au long du temps d'évolution ainsi que d'étudier plusieurs observables qui dépendent de la probabilité comme le spectre pT, RAA (le facteur de modification nucléaire) et les coefficients de débit (v2). De plus, ces modèles offrent une fenêtre intéressante pour aborder la dynamique interne des paires, y compris le potentiel d'interaction entre les quarks et l'influence du milieu sur la dynamique des paires.
