Inclusive J/$\psi$ yields and average transverse momenta in p-Pb collisions at a center-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}$ = 8.16 TeV are measured as a function of the charged-particle pseudorapidity density with ALICE. The J/$\psi$ mesons are reconstructed at forward $(2.03 < y_{\rm cms} < 3.53)$ and backward ($-4.46 < y_{\rm cms} < -2.96$) center-of-mass rapidity in their dimuon decay channel while the charged-particle pseudorapidity density is measured around midrapidity. The J/$\psi$ yields at forward and backward rapidity normalized to their respective average values increase with the normalized charged-particle pseudorapidity density, the former showing a weaker increase than the latter. The normalized average transverse momenta at forward and backward rapidity manifest a steady increase from low to high charged-particle pseudorapidity density with a saturation beyond the average value.

The first measurement of the $\Upsilon(1{\rm S})$ elliptic flow coefficient ($v_2$) is performed at forward rapidity (2.5 $<$ $y$ $<$ 4) in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV with the ALICE detector at the LHC. The results are obtained with the scalar product method and are reported as a function of transverse momentum ($p_{\rm{T}}$) up to 15 GeV/$c$ in the 5-60% centrality interval. The measured $\Upsilon(1{\rm S})$ $v_2$ is consistent with zero and with the small positive values predicted by transport models within uncertainties. The $v_2$ coefficient in 2 $<$ $p_{\rm T}$ $<$ 15 GeV/$c$ is lower than that of inclusive J/$\psi$ mesons in the same $p_{\rm{T}}$ interval by 2.6 standard deviations. These results, combined with earlier suppression measurements, are in agreement with a scenario in which the $\Upsilon$(1S) production in Pb-Pb collisions at LHC energies is dominated by dissociation limited to the early stage of the collision whereas in the J/$\psi$ case there is substantial experimental evidence of an additional regeneration component.

Measurements of anisotropic flow coefficients ($v_n$) and their cross-correlations using two- and multi-particle cumulant methods are reported in collisions of pp at $\sqrt{s} = 13$ TeV, p-Pb at $\sqrt{s_{_{\rm NN}}} = 5.02$ TeV, Xe-Xe at $\sqrt{s_{_{\rm NN}}} = 5.44$ TeV, and Pb-Pb at $\sqrt{s_{_{\rm NN}}} = 5.02$ TeV recorded with the ALICE detector. The multiplicity dependence of $v_n$ is studied in a very wide range from 20 to 3000 particles produced in the mid-rapidity region $|\eta|<0.8$ for the transverse momentum range $0.2 < p_{\rm T} < 3.0$ GeV/$c$. An ordering of the coefficients $v_2 > v_3 > v_4$ is found in pp and p-Pb collisions, similar to that seen in large collision systems, while a weak $v_2$ multiplicity dependence is observed relative to nucleus-nucleus collisions in the same multiplicity range. Using a novel subevent method, $v_{2}$ measured with four-particle cumulants is found to be compatible with that from six-particle cumulants in pp and p-Pb collisions. The magnitude of the correlation between $v_n^2$ and $v_m^2$, evaluated with the symmetric cumulants SC$(m,n)$ is observed to be positive at all multiplicities for $v_2$ and $v_4$, while for $v_2$ and $v_3$ it is negative and changes sign for multiplicities below 100, which may indicate a different $v_n$ fluctuation pattern in this multiplicity range. The observed long-range multi-particle azimuthal correlations in high multiplicity pp and p-Pb collisions can neither be described by PYTHIA 8 nor by IP-Glasma+MUSIC+UrQMD model calculations, and hence provide new insights into the understanding of collective effects in small collision systems.

We present the measurement of a new set of jet shape observables for track-based jets in central Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 2.76$ TeV. The set of jet shapes includes the first radial moment or angularity, $g$; the momentum dispersion, $p_{\rm T}D$; and the difference between the leading and sub-leading constituent track transverse momentum, $LeSub$. These observables provide complementary information on the jet fragmentation and can constrain different aspects of the theoretical description of jet-medium interactions. The jet shapes were measured for a small resolution parameter $R = 0.2$ and were fully corrected to particle level. The observed jet shape modifications indicate that in-medium fragmentation is harder and more collimated than vacuum fragmentation as obtained by PYTHIA calculations, which were validated with the measurements of the jet shapes in proton-proton collisions at $\sqrt{s} = 7$ TeV. The comparison of the measured distributions to templates for quark and gluon-initiated jets indicates that in-medium fragmentation resembles that of quark jets in vacuum. We further argue that the observed modifications are not consistent with a totally coherent energy loss picture where the jet loses energy as a single colour charge, suggesting that the medium resolves the jet structure at the angular scales probed by our measurements ($R=0.2$). Furthermore, we observe that small-$R$ jets can help to isolate purely energy loss effects from other effects that contribute to the modifications of the jet shower in medium such as the correlated background or medium response.

The ALICE Collaboration has measured the energy dependence of exclusive photoproduction of $J/\psi$ vector mesons off proton targets in ultra-peripheral p-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}} = 5.02$ TeV. The e$^+$e$^-$ and $\mu^+\mu^-$ decay channels are used to measure the cross section as a function of the rapidity of the $J/\psi$ in the range $-2.5 < y < 2.7$, corresponding to an energy in the $\gamma$p centre-of-mass in the interval $40 < W_{\gamma\mathrm{p}}<550$ GeV. The measurements, which are consistent with a power law dependence of the exclusive $J/\psi$ photoproduction cross section, are compared to previous results from HERA and the LHC and to several theoretical models. They are found to be compatible with previous measurements.

The measurements of the inclusive J/$\psi$ yield at midrapidity ($\left | y \right | < 0.9$) and forward rapidity (2.5 $< y <$ 4) in Pb$-$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV with the ALICE detector at the LHC are reported. The inclusive J/$\psi$ production yields and nuclear modification factors, $R_{\rm AA}$, are measured as a function of the collision centrality, J/$\psi$ transverse momentum ($p_{\rm T}$), and rapidity. The J/$\psi$ average transverse momentum and squared transverse momentum ($\langle p_{\mathrm{T}}\rangle$ and $\langle p_{\mathrm{T}}^{\mathrm{2}}\rangle$) are evaluated as a function of the centrality at midrapidity. Compared to the previous ALICE publications, here the entire Pb$-$Pb collisions dataset collected during the LHC Run 2 is used, which improves the precision of the measurements and extends the $p_{\rm T}$ coverage. The $p_{\rm T}$-integrated $R_{\rm AA}$ shows a hint of an increasing trend towards unity from semicentral to central collisions at midrapidity, while it is flat at forward rapidity. The $p_{\rm T}$-differential $R_{\rm AA}$ shows a strong suppression at high $p_{\rm T}$ with less suppression at low $p_{\rm T}$ where it reaches a larger value at midrapidity compared to forward rapidity. The ratio of the $p_{\rm T}$-integrated yields of J/$\psi$ to those of D$^{0}$ mesons is reported for the first time for the central and semicentral event classes at midrapidity. Model calculations implementing charmonium production via the coalescence of charm quarks and antiquarks during the fireball evolution (transport models) or in a statistical approach with thermal weights are in good agreement with the data at low $p_{\rm T}$. At higher $p_{\rm T}$, the data are well described by transport models and a model based on energy loss in the strongly-interacting medium produced in nuclear collisions at the LHC.

The Quark Gluon Plasma (QGP) produced in ultra relativistic heavy-ion collisions at the Large Hadron Collider (LHC) can be studied by measuring the modifications of jets formed by hard scattered partons which interact with the medium. We studied these modifications via angular correlations of jets with charged hadrons for jets with momenta 20 < $p_{\rm{T}}^{\rm{jet}}$ < 40 GeV/$c$ as a function of the associated particle momentum. The reaction plane fit (RPF) method is used in this analysis to remove the flow modulated background. The analysis of angular correlations for different orientations of the jet relative to the second order event plane allows for the study of the path length dependence of medium modifications to jets. We present the dependence of azimuthal angular correlations of charged hadrons with respect to the angle of the axis of a reconstructed jet relative to the event plane in Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 2.76 TeV. The dependence of particle yields associated with jets on the angle of the jet with respect to the event plane is presented. Correlations at different angles relative to the event plane are compared through ratios and differences of the yield. No dependence of the results on the angle of the jet with respect to the event plane is observed within uncertainties, which is consistent with no significant path length dependence of the medium modifications for this observable.

Mid-rapidity production of $\pi^{\pm}$, $\rm{K}^{\pm}$ and ($\bar{\rm{p}}$)p measured by the ALICE experiment at the LHC, in Pb-Pb and inelastic pp collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV, is presented. The invariant yields are measured over a wide transverse momentum ($p_{\rm{T}}$) range from hundreds of MeV/$c$ up to 20 GeV/$c$. The results in Pb-Pb collisions are presented as a function of the collision centrality, in the range 0$-$90%. The comparison of the $p_{\rm{T}}$-integrated particle ratios, i.e. proton-to-pion (p/$\pi$) and kaon-to-pion (K/$\pi$) ratios, with similar measurements in Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 2.76 TeV show no significant energy dependence. Blast-wave fits of the $p_{\rm{T}}$ spectra indicate that in the most central collisions radial flow is slightly larger at 5.02 TeV with respect to 2.76 TeV. Particle ratios (p/$\pi$, K/$\pi$) as a function of $p_{\rm{T}}$ show pronounced maxima at $p_{\rm{T}}$ $\approx$ 3 GeV/$c$ in central Pb-Pb collisions. At high $p_{\rm{T}}$, particle ratios at 5.02 TeV are similar to those measured in pp collisions at the same energy and in Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 2.76 TeV. Using the pp reference spectra measured at the same collision energy of 5.02 TeV, the nuclear modification factors for the different particle species are derived. Within uncertainties, the nuclear modification factor is particle species independent for high $p_{\rm{T}}$ and compatible with measurements at $\sqrt{s_{\rm{NN}}}$ = 2.76 TeV. The results are compared to state-of-the-art model calculations, which are found to describe the observed trends satisfactorily.

In ultrarelativistic heavy-ion collisions, the event-by-event variation of the elliptic flow $v_2$ reflects fluctuations in the shape of the initial state of the system. This allows to select events with the same centrality but different initial geometry. This selection technique, Event Shape Engineering, has been used in the analysis of charge-dependent two- and three-particle correlations in Pb-Pb collisions at $\sqrt{s_{_{\rm NN}}} =2.76$ TeV. The two-particle correlator $\langle \cos(\varphi_\alpha - \varphi_\beta) \rangle$, calculated for different combinations of charges $\alpha$ and $\beta$, is almost independent of $v_2$ (for a given centrality), while the three-particle correlator $\langle \cos(\varphi_\alpha + \varphi_\beta - 2\Psi_2) \rangle$ scales almost linearly both with the event $v_2$ and charged-particle pseudorapidity density. The charge dependence of the three-particle correlator is often interpreted as evidence for the Chiral Magnetic Effect (CME), a parity violating effect of the strong interaction. However, its measured dependence on $v_2$ points to a large non-CME contribution to the correlator. Comparing the results with Monte Carlo calculations including a magnetic field due to the spectators, the upper limit of the CME signal contribution to the three-particle correlator in the 10-50% centrality interval is found to be 26-33% at 95% confidence level.

We report measurements of the production of prompt D$^0$, D$^+$, D$^{*+}$ and D$^+_{\rm s}$ mesons in Pb-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}=5.02$ TeV, in the centrality classes 0-10%, 30-50% and 60-80%. The D-meson production yields are measured at mid-rapidity ($|y|<0.5$) as a function of transverse momentum ($p_{\rm T}$). The $p_{\rm T}$ intervals covered in central collisions are: $1<p_{\rm T}<50$ Gev/$c$ for D$^0$, $2<p_{\rm T}<50$ GeV/$c$ for D$^+$, $3<p_{\rm T}<50$ GeV/$c$ for D$^{*+}$, and $4<p_{\rm T}<16$ GeV/$c$ for D$^+_{\rm s}$ mesons. The nuclear modification factors ($R_{\rm AA}$) for non-strange D mesons (D$^0$, D$^+$, D$^{*+}$) show minimum values of about 0.2 for $p_{\rm T}$ = 6-10 GeV/$c$ in the most central collisions and are compatible within uncertainties with those measured at $\sqrt{s_{\rm NN}}=2.76$ TeV. For D$^+_{\rm s}$ mesons, the values of $R_{\rm AA}$ are larger than those of non-strange D mesons, but compatible within uncertainties. In central collisions the average $R_{\rm AA}$ of non-strange D mesons is compatible with that of charged particles for $p_{\rm T} > 8$ GeV/$c$, while it is larger at lower $p_{\rm T}$. The nuclear modification factors for strange and non-strange D mesons are also compared to theoretical models with different implementations of in-medium energy loss.