We present the first measurements of femtoscopic correlations between the K$^0_{\rm S}$ and K$^{\rm \pm}$ particles in pp collisions at $\sqrt{s}=7$ TeV measured by the ALICE experiment. The observed femtoscopic correlations are consistent with final-state interactions proceeding solely via the $a_0(980)$ resonance. The extracted kaon source radius and correlation strength parameters for K$^0_{\rm S}$K$^{\rm -}$ are found to be equal within the experimental uncertainties to those for K$^0_{\rm S}$K$^{\rm +}$. Results of the present study are compared with those from identical-kaon femtoscopic studies also performed with pp collisions at $\sqrt{s}=7$ TeV by ALICE and with a K$^0_{\rm S}$K$^{\rm \pm}$ measurement in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV. Combined with the Pb-Pb results, our pp analysis is found to be compatible with the interpretation of the $a_0(980)$ having a tetraquark structure instead of that of a diquark.

The production yield of prompt D mesons and their elliptic flow coefficient $v_2$ were measured with the Event-Shape Engineering (ESE) technique applied to mid-central (10-30% and 30-50% centrality classes) Pb-Pb collisions at the centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}} =5.02$ TeV, with the ALICE detector at the LHC. The ESE technique allows the classification of events, belonging to the same centrality, according to the azimuthal anisotropy of soft particle production in the collision. The reported measurements give the opportunity to investigate the dynamics of charm quarks in the Quark-Gluon Plasma and provide information on their participation in the collective expansion of the medium. D mesons were reconstructed via their hadronic decays at mid-rapidity, $|\eta|<0.8$, in the transverse momentum interval $1<p_{\rm T}<24$ GeV/$c$. The $v_2$ coefficient is found to be sensitive to the event-shape selection confirming a correlation between the D-meson azimuthal anisotropy and the collective expansion of the bulk matter, while the per-event D-meson yields do not show any significant modification within the current uncertainties.

Charged-particle spectra at midrapidity are measured in Pb-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$ = 5.02 TeV and presented in centrality classes ranging from most central (0-5%) to most peripheral (95-100%) collisions. Possible medium effects are quantified using the nuclear modification factor ($R_{\rm AA}$) by comparing the measured spectra with those from proton-proton collisions, scaled by the number of independent nucleon-nucleon collisions obtained from a Glauber model. At large transverse momenta ($8<p_{\rm T}<20$ GeV/$c$), the average $R_{\rm AA}$ is found to increase from about $0.15$ in 0-5% central to a maximum value of about $0.8$ in 75-85% peripheral collisions, beyond which it falls off strongly to below $0.2$ for the most peripheral collisions. Furthermore, $R_{\rm AA}$ initially exhibits a positive slope as a function of $p_{\rm T}$ in the $8$-$20$ GeV/$c$ interval, while for collisions beyond the 80% class the slope is negative. To reduce uncertainties related to event selection and normalization, we also provide the ratio of $R_{\rm AA}$ in adjacent centrality intervals. Our results in peripheral collisions are consistent with a PYTHIA-based model without nuclear modification, demonstrating that biases caused by the event selection and collision geometry can lead to the apparent suppression in peripheral collisions. This explains the unintuitive observation that $R_{\rm AA}$ is below unity in peripheral Pb-Pb, but equal to unity in minimum-bias p-Pb collisions despite similar charged-particle multiplicities.

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.

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.

Azimuthally-differential femtoscopic measurements, being sensitive to spatio-temporal characteristics of the source as well as to the collective velocity fields at freeze out, provide very important information on the nature and dynamics of the system evolution. While the HBT radii oscillations relative to the second harmonic event plane measured recently reflect mostly the spatial geometry of the source, model studies have shown that the HBT radii oscillations relative to the third harmonic event plane are predominantly defined by the velocity fields. In this Letter, we present the first results on azimuthally-differential pion femtoscopy relative to the third harmonic event plane as a function of the pion pair transverse momentum $k_{\rm T}$ for different collision centralities in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV. We find that the $R_{\rm side}$ and $R_{\rm out}$ radii, which characterize the pion source size in the directions perpendicular and parallel to the pion transverse momentum, oscillate in phase relative to the third harmonic event plane, similar to the results from 3+1D hydrodynamical calculations. The observed radii oscillations unambiguously signal a collective expansion and anisotropy in the velocity fields. A comparison of the measured radii oscillations with the Blast-Wave model calculations indicate that the initial state triangularity is washed-out at freeze out.

We present the first ever measurements of femtoscopic correlations between the K$^0_{\rm S}$ and K$^{\rm \pm}$ particles. The analysis was performed on the data from Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV measured by the ALICE experiment. The observed femtoscopic correlations are consistent with final-state interactions proceeding via the $a_0(980)$ resonance. The extracted kaon source radius and correlation strength parameters for K$^0_{\rm S}$K$^{\rm -}$ are found to be equal within the experimental uncertainties to those for K$^0_{\rm S}$K$^{\rm +}$. Comparing the results of the present study with those from published identical-kaon femtoscopic studies by ALICE, mass and coupling parameters for the $a_0$ resonance are constrained. Our results are also compatible with the interpretation of the $a_0$ having a tetraquark structure over that of a diquark.

Measurements of anisotropic flow coefficients with two- and multi-particle cumulants for inclusive charged particles in Pb-Pb collisions at $\sqrt{{\textit s}_\text{NN}} = 5.02$ and 2.76 TeV are reported in the pseudorapidity range $|\eta| < 0.8$ and transverse momentum $0.2 < p_\text{T} < 50$ GeV/$c$. The full data sample collected by the ALICE detector in 2015 (2010), corresponding to an integrated luminosity of 12.7 (2.0) $\mu$b$^{-1}$ in the centrality range 0-80%, is analysed. Flow coefficients up to the sixth flow harmonic ($v_6$) are reported and a detailed comparison among results at the two energies is carried out. The $p_\text{T}$ dependence of anisotropic flow coefficients and its evolution with respect to centrality and harmonic number $n$ are investigated. An approximate power-law scaling of the form $v_n(p_\text{T}) \sim p_\text{T}^{n/3}$ is observed for all flow harmonics at low $p_\text{T}$ ($0.2 < p_\text{T} < 3$ GeV/$c$). At the same time, the ratios $v_n/v_m^{n/m}$ are observed to be essentially independent of $p_\text{T}$ for most centralities up to about $p_\text{T} = 10$ GeV/$c$. Analysing the differences among higher-order cumulants of elliptic flow ($v_2$), which have different sensitivities to flow fluctuations, a measurement of the standardised skewness of the event-by-event $v_2$ distribution $P(v_2)$ is reported and constraints on its higher moments are provided. The Elliptic Power distribution is used to parametrise $P(v_2)$, extracting its parameters from fits to cumulants. The measurements are compared to different model predictions in order to discriminate among initial-state models and to constrain the temperature dependence of the shear viscosity to entropy-density ratio.

The elliptic ($v_2$), triangular ($v_3$), and quadrangular ($v_4$) flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, ${\rm p+\overline{p}}$, ${\Lambda+\overline{\Lambda}}$, ${\rm K}^{\rm 0}_{\rm S}$, and the $\phi$-meson are measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV. Results obtained with the scalar product method are reported for the rapidity range $\vert y \vert <$ 0.5 as a function of transverse momentum, $p_\text{T}$, at different collision centrality intervals between 0-70%, including ultra-central (0-1%) collisions for $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$. For $p_\text{T} < 3$ GeV$\kern-0.05em/\kern-0.02em c$, the flow coefficients exhibit a particle mass dependence. At intermediate transverse momenta ($3<p_\text{T}<$~8-10 GeV$\kern-0.05em/\kern-0.02em c$), particles show an approximate grouping according to their type (i.e., mesons and baryons). The $\phi$-meson $v_2$, which tests both particle mass dependence and type scaling, follows ${\rm p+\overline{p}}$ $v_2$ at low $p_\text{T}$ and $\pi^{\pm}$ $v_2$ at intermediate $p_\text{T}$. The evolution of the shape of $v_{\rm n}(p_{\mathrm{T}})$ as a function of centrality and harmonic number $n$ is studied for the various particle species. Flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$ for $p_\text{T}<3$ GeV$\kern-0.05em/\kern-0.02em c$ are compared to iEBE-VISHNU and MUSIC hydrodynamical calculations coupled to a hadronic cascade model (UrQMD). The iEBE-VISHNU calculations describe the results fairly well for $p_\text{T} < 2.5$ GeV$\kern-0.05em/\kern-0.02em c$, while MUSIC calculations reproduce the measurements for $p_\text{T} < 1$ GeV$\kern-0.05em/\kern-0.02em c$. A comparison to $v_{\rm n}$ coefficients measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV is also provided.

The transverse structure of jets was studied via jet fragmentation transverse momentum ($j_{\rm{T}}$) distributions, obtained using two-particle correlations in proton-proton and proton-lead collisions, measured with the ALICE experiment at the LHC. The highest transverse momentum particle in each event is used as the trigger particle and the region $3 < p_{\rm{Tt}} < 15$ GeV/$c$ is explored in this study. The measured distributions show a clear narrow Gaussian component and a wide non-Gaussian one. Based on Pythia simulations, the narrow component can be related to non-perturbative hadronization and the wide component to quantum chromodynamical splitting. The width of the narrow component shows a weak dependence on the transverse momentum of the trigger particle, in agreement with the expectation of universality of the hadronization process. On the other hand, the width of the wide component shows a rising trend suggesting increased branching for higher transverse momentum. The results obtained in pp collisions at $\sqrt{s}$ = 7 TeV and in p-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV are compatible within uncertainties and hence no significant cold nuclear matter effects are observed. The results are compared to previous measurements from CCOR and PHENIX as well as to Pythia 8 and Herwig 7 simulations.