Hot QCD physics studies the nuclear strong force under extreme temperature and densities. Experimentally these conditions are achieved via high-energy collisions of heavy ions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). In the past decade, a unique and substantial suite of data was collected at RHIC and the LHC, probing hydrodynamics at the nucleon scale, the temperature dependence of the transport properties of quark-gluon plasma, the phase diagram of nuclear matter, the interaction of quarks and gluons at different scales and much more. This document, as part of the 2023 nuclear science long range planning process, was written to review the progress in hot QCD since the 2015 Long Range Plan for Nuclear Science, as well as highlight the realization of previous recommendations, and present opportunities for the next decade, building on the accomplishments and investments made in theoretical developments and the construction of new detectors. Furthermore, this document provides additional context to support the recommendations voted on at the Joint Hot and Cold QCD Town Hall Meeting, which are reported in a separate document.

Collective flow observables are known to be a sensitive tool to gain insights on the equation of state of nuclear matter from heavy-ion collision observations. Towards more quantitative constraints one has to carefully assess other influences on the collective behaviour. In this work a hadronic transport approach SMASH (Simulating Many Accelerated Strongly-interacting Hadrons) is applied to study the first four anisotropic flow coefficients in Au+Au collisions at $E_{\rm lab}=1.23A$ GeV in the context of the recently measured data by the HADES collaboration. In particular, the formation of light nuclei is important in this energy regime. Two different approaches are contrasted to each other: A clustering algorithm inspired by coalescence as well as microscopic formation of deuterons via explicit cross-sections. The sensitivity of directed and elliptic flow observables to the strength of the Skyrme mean field is explored. In addition, it is demonstrated that the rapidity-odd $v_3$ coefficient is practically zero in this energy regime and the ratio of $v_4/v_2^2$ is close to the value of 0.5 expected from hydrodynamic behaviour. This study establishes the current understanding of collective behaviour within the SMASH approach and lays the ground for future more quantitative constraints on the equation of state of nuclear matter within improved mean field calculations.

The nature of the QCD phase transition in the chiral limit presents a challenging problem for lattice QCD. However, its study provides constraints on the phase diagram at the physical point. In this work, we investigate how the order of the chiral phase transition depends on the number of light quark flavours. To approach the lattice chiral limit, we map out and extrapolate the chiral critical surface that separates the first-order region from the crossover region in an extended parameter space, which includes the gauge coupling, the number of quark flavours, their masses, and the lattice spacing. Lattice simulations with standard staggered quarks reveal that for each N_f < 8, there exists a tricritical lattice spacing $a^\text{tric}(N_f)$, at which the chiral transition changes from first order (a>a^\text{tric}) to second order ($a

Fermions become polarized in a vortical fluid due to spin-vorticity coupling. The spin polarization density is proportional to the local fluid vorticity at the next-to-leading order of a gradient expansion in a quantum kinetic theory. Spin correlations of two $\Lambda$-hyperons can therefore reveal the vortical structure of the dense matter in high-energy heavy-ion collisions. We employ a (3+1)D viscous hydrodynamic model with event-by-event fluctuating initial conditions from A MultiPhase Transport (AMPT) model to calculate the vorticity distributions and $\Lambda$ spin correlations. The azimuthal correlation of the transverse spin is shown to have a cosine form plus an offset due to a circular structure of the transverse vorticity around the beam direction and global spin polarization. The longitudinal spin correlation shows a structure of vortex-pairing in the transverse plane due to the convective flow of hot spots in the radial direction. The dependence on colliding energy, rapidity, centrality and sensitivity to the shear viscosity are also investigated.

(Abridged) We here continue our effort to model the behaviour of matter when orbiting or accreting onto a generic black hole by developing a new numerical code employing advanced techniques geared solve the equations of in general-relativistic hydrodynamics. The new code employs a number of high-resolution shock-capturing Riemann-solvers and reconstruction algorithms, exploiting the enhanced accuracy and the reduced computational cost of AMR techniques. In addition, the code makes use of sophisticated ray-tracing libraries that, coupled with general-relativistic radiation-transfer calculations, allow us to compute accurately the electromagnetic emissions from such accretion flows. We validate the new code by presenting an extensive series of stationary accretion flows either in spherical or axial symmetry and performed either in 2D or 3D. In addition, we consider the highly nonlinear scenario of a recoiling black hole produced in the merger of a supermassive black hole binary interacting with the surrounding circumbinary disc. In this way we can present, for the first time, ray-traced images of the shocked fluid and the light-curve resulting from consistent general-relativistic radiation-transport calculations from this process. The work presented here lays the ground for the development of a generic computational infrastructure employing AMR techniques to deal accurately and self-consistently with accretion flows onto compact objects. In addition to the accurate handling of the matter, we provide a self-consistent electromagnetic emission from these scenarios by solving the associated radiative-transfer problem. While magnetic fields are presently excluded from our analysis, the tools presented here can have a number of applications to study accretion flows onto black holes or neutron stars.