The azimuthal anisotropy coefficient $v_2$ of prompt D$^0$, D$^+$, D$^{*+}$ and D$_s^+$ mesons was measured in mid-central (30-50% centrality class) Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}} = 5.02$ TeV, with the ALICE detector at the LHC. The D mesons were reconstructed via their hadronic decays at mid-rapidity, $|y|<0.8$, in the transverse momentum interval $1<p_{\rm T}<24$ GeV/$c$. The measured D-meson $v_2$ has similar values as that of charged pions. The D$_s^+$ $v_2$, measured for the first time, is found to be compatible with that of non-strange D mesons. The measurements are compared with theoretical calculations of charm-quark transport in a hydrodynamically expanding medium and have the potential to constrain medium parameters.

The invariant differential cross sections for inclusive $\pi^{0}$ and $\eta$ mesons at midrapidity were measured in pp collisions at $\sqrt{s}=2.76$ TeV for transverse momenta $0.4<p_{\rm T}<40$ GeV/$c$ and $0.6<p_{\rm T}<20$ GeV/$c$, respectively, using the ALICE detector. This large range in $p_{\rm T}$ was achieved by combining various analysis techniques and different triggers involving the electromagnetic calorimeter (EMCal). In particular, a new single-cluster, shower-shape based method was developed for the identification of high-$p_{\rm T}$ neutral pions, which exploits that the showers originating from their decay photons overlap in the EMCal. The measured cross sections are found to exhibit a similar power-law behavior with an exponent of about $6.3$. Next-to-leading-order perturbative QCD calculations differ from the measured cross sections by about $30$% for the $\pi^0$, and between $30$-$50$% for the $\eta$ meson, while generator-level simulations with PYTHIA 8.2 describe the data to better than $10$-$30$%, except at $p_{\rm T}<1$ GeV/$c$. The new data can therefore be used to further improve the theoretical description of $\pi^{0}$ and $\eta$ meson production.

The measurement of azimuthal correlations of charged particles is presented for Pb-Pb collisions at $\sqrt{s_{\rm NN}}=$ 2.76 TeV and p-Pb collisions at $\sqrt{s_{\rm NN}}=$ 5.02 TeV with the ALICE detector at the CERN Large Hadron Collider. These correlations are measured for the second, third and fourth order flow vector in the pseudorapidity region $|\eta|<0.8$ as a function of centrality and transverse momentum $p_{\rm T}$ using two observables, to search for evidence of $p_{\rm T}$-dependent flow vector fluctuations. For Pb-Pb collisions at 2.76 TeV, the measurements indicate that $p_{\rm T}$-dependent fluctuations are only present for the second order flow vector. Similar results have been found for p-Pb collisions at 5.02 TeV. These measurements are compared to hydrodynamic model calculations with event-by-event geometry fluctuations in the initial state to constrain the initial conditions and transport properties of the matter created in Pb-Pb and p-Pb collisions.

The centrality dependence of the mean charged-particle multiplicity as a function of pseudorapidity is measured in approximately 1 $\mu$b$^{-1}$ of proton--lead collisions at a nucleon--nucleon centre-of-mass energy of $\sqrt{s_{_{\rm{NN}}}} = 5.02$ TeV using the ATLAS detector at the Large Hadron Collider. Charged particles with absolute pseudorapidity less than 2.7 are reconstructed using the ATLAS pixel detector. The $p$+Pb collision centrality is characterised by the total transverse energy measured in the Pb-going direction of the forward calorimeter. The charged-particle pseudorapidity distributions are found to vary strongly with centrality, with an increasing asymmetry between the proton-going and Pb-going directions as the collisions become more central. Three different estimations of the number of nucleons participating in the $p$+Pb collision have been carried out using the Glauber model as well as two Glauber--Gribov inspired extensions to the Glauber model. Charged-particle multiplicities per participant pair are found to vary differently for these three models, highlighting the importance of including colour fluctuations in nucleon--nucleon collisions in the modelling of the initial state of $p$+Pb collisions.

A measurement of $J/\psi$ and $\psi(2\mathrm{S})$ production is presented. It is based on a data sample from Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV and $pp$ collisions at $\sqrt{s}$ = 5.02 TeV recorded by the ATLAS detector at the LHC in 2015, corresponding to an integrated luminosity of $0.42\mathrm{nb}^{-1}$ and $25\mathrm{pb}^{-1}$ in Pb+Pb and $pp$, respectively. The measurements of per-event yields, nuclear modification factors, and non-prompt fractions are performed in the dimuon decay channel for $9 < p_{T}^{\mu\mu} < 40$ GeV in dimuon transverse momentum, and $-2.0 < y_{\mu\mu} < 2.0$ in rapidity. Strong suppression is found in Pb+Pb collisions for both prompt and non-prompt $J/\psi$, as well as for prompt and non-prompt $\psi(2\mathrm{S})$, increasing with event centrality. The suppression of prompt $\psi(2\mathrm{S})$ is observed to be stronger than that of $J/\psi$, while the suppression of non-prompt $\psi(2\mathrm{S})$ is equal to that of the non-prompt $J/\psi$ within uncertainties, consistent with the expectation that both arise from \textit{b}-quarks propagating through the medium. Despite prompt and non-prompt $J/\psi$ arising from different mechanisms, the dependence of their nuclear modification factors on centrality is found to be quite similar.

To assess the properties of the quark-gluon plasma formed in heavy-ion collisions, the ATLAS experiment at the LHC measures a correlation between the mean transverse momentum and the magnitudes of the flow harmonics. The analysis uses data samples of lead-lead and proton-lead collisions obtained at the centre-of-mass energy per nucleon pair of 5.02 TeV, corresponding to total integrated luminosities of $22 ~\mu b^{-1}$ and $28~nb^{-1}$, respectively. The measurement is performed using a modified Pearson correlation coefficient with the charged-particle tracks on an event-by-event basis. The modified Pearson correlation coefficients for the $2^{nd}$-, 3$^{rd}$-, and 4$^{th}$-order harmonics are measured as a function of event centrality quantified as the number of charged particles or the number of nucleons participating in the collision. The measurements are performed for several intervals of the charged-particle transverse momentum. The correlation coefficients for all studied harmonics exhibit a strong centrality evolution in the lead-lead collisions, which only weakly depends on the charged-particle momentum range. In the proton-lead collisions, the modified Pearson correlation coefficient measured for the second harmonics shows only weak centrality dependence. The data is qualitatively described by the predictions based on the hydrodynamical model.

We report the triton ($t$) production in mid-rapidity ($|y| <$ 0.5) Au+Au collisions at $\sqrt{s_\mathrm{NN}}$= 7.7--200 GeV measured by the STAR experiment from the first phase of the beam energy scan at the Relativistic Heavy Ion Collider (RHIC). The nuclear compound yield ratio ($\mathrm{N}_t \times \mathrm{N}_p/\mathrm{N}_d^2$), which is predicted to be sensitive to the fluctuation of local neutron density, is observed to decrease monotonically with increasing charged-particle multiplicity ($dN_{ch}/d\eta$) and follows a scaling behavior. The $dN_{ch}/d\eta$ dependence of the yield ratio is compared to calculations from coalescence and thermal models. Enhancements in the yield ratios relative to the coalescence baseline are observed in the 0%-10% most central collisions at 19.6 and 27 GeV, with a significance of 2.3$\sigma$ and 3.4$\sigma$, respectively, giving a combined significance of 4.1$\sigma$. The enhancements are not observed in peripheral collisions or model calculations without critical fluctuation, and decreases with a smaller $p_{T}$ acceptance. The physics implications of these results on the QCD phase structure and the production mechanism of light nuclei in heavy-ion collisions are discussed.

Azimuthal anisotropies of muons from charm and bottom hadron decays are measured in Pb+Pb collisions at $\sqrt{s_\mathrm{NN}}= 5.02$ TeV. The data were collected with the ATLAS detector at the Large Hadron Collider in 2015 and 2018 with integrated luminosities of $0.5~\mathrm{nb}^{-1}$ and $1.4~\mathrm{nb^{-1}}$, respectively. The kinematic selection for heavy-flavor muons requires transverse momentum $4 < p_\mathrm{T} < 30$ GeV and pseudorapidity $|\eta|<2.0$. The dominant sources of muons in this $p_\mathrm{T}$ range are semi-leptonic decays of charm and bottom hadrons. These heavy-flavor muons are separated from light-hadron decay muons and punch-through hadrons using the momentum imbalance between the measurements in the tracking detector and in the muon spectrometers. Azimuthal anisotropies, quantified by flow coefficients, are measured via the event-plane method for inclusive heavy-flavor muons as a function of the muon $p_\mathrm{T}$ and in intervals of Pb+Pb collision centrality. Heavy-flavor muons are separated into contributions from charm and bottom hadron decays using the muon transverse impact parameter with respect to the event primary vertex. Non-zero elliptic ($v_{2}$) and triangular ($v_{3}$) flow coefficients are extracted for charm and bottom muons, with the charm muon coefficients larger than those for bottom muons for all Pb+Pb collision centralities. The results indicate substantial modification to the charm and bottom quark angular distributions through interactions in the quark-gluon plasma produced in these Pb+Pb collisions, with smaller modifications for the bottom quarks as expected theoretically due to their larger mass.

The azimuthal anisotropic flow of identified and unidentified charged particles has been systematically studied in Cu+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV for harmonics $n=$ 1-4 in the pseudorapidity range $|\eta|<1$. The directed flow in Cu+Au collisions is compared with the rapidity-odd and, for the first time, the rapidity-even components of charged particle directed flow in Au+Au collisions at $\sqrt{s_{_{NN}}}$ = 200~GeV. The slope of the directed flow pseudorapidity dependence in Cu+Au collisions is found to be similar to that in Au+Au collisions, with the intercept shifted toward positive $\eta$ values, i.e., the Cu-going direction. The mean transverse momentum projected onto the spectator plane, $\langle p_x\rangle$, in Cu+Au collision also exhibits approximately linear dependence on $\eta$ with the intercept at about $\eta\approx-0.4$, closer to the rapidity of the Cu+Au system center-of-mass. The observed dependencies find natural explanation in a picture of the directed flow originating partly due the "tilted source" and partly due to the rapidity dependent asymmetry in the initial density distribution. Charge-dependence of the $\langle p_x\rangle$ was also observed in Cu+Au collisions, indicating an effect of the initial electric field created by charge difference of the spectator protons in two colliding nuclei. The rapidity-even component of directed flow in Au+Au collisions is close to that in Pb+Pb collisions at $\sqrt{s_{_{NN}}}$ = 2.76 TeV, indicating a similar magnitude of dipole-like fluctuations in the initial-state density distribution. Higher harmonic flow in Cu+Au collisions exhibits similar trends to those observed in Au+Au and Pb+Pb collisions and is qualitatively reproduced by a viscous hydrodynamic model and a multi-phase transport model. For all harmonics with $n\ge2$ we observe an approximate scaling of $v_n$ with the number of constituent quarks.

New measurements of directed flow for charged hadrons, characterized by the Fourier coefficient \vone, are presented for transverse momenta $\mathrm{p_T}$, and centrality intervals in Au+Au collisions recorded by the STAR experiment for the center-of-mass energy range $\mathrm{\sqrt{s_{_{NN}}}} = 7.7 - 200$ GeV. The measurements underscore the importance of momentum conservation and the characteristic dependencies on $\mathrm{\sqrt{s_{_{NN}}}}$, centrality and $\mathrm{p_T}$ are consistent with the expectations of geometric fluctuations generated in the initial stages of the collision, acting in concert with a hydrodynamic-like expansion. The centrality and $\mathrm{p_T}$ dependencies of $\mathrm{v^{even}_{1}}$, as well as an observed similarity between its excitation function and that for $\mathrm{v_3}$, could serve as constraints for initial-state models. The $\mathrm{v^{even}_{1}}$ excitation function could also provide an important supplement to the flow measurements employed for precision extraction of the temperature dependence of the specific shear viscosity.