Three- and four-pion Bose-Einstein correlations are presented in pp, p-Pb, and Pb-Pb collisions at the LHC. We compare our measured four-pion correlations to the expectation derived from two- and three-pion measurements. Such a comparison provides a method to search for coherent pion emission. We also present mixed-charge correlations in order to demonstrate the effectiveness of several analysis procedures such as Coulomb corrections. Same-charge four-pion correlations in pp and p-Pb appear consistent with the expectations from three-pion measurements. However, the presence of non-negligible background correlations in both systems prevent a conclusive statement. In Pb-Pb collisions, we observe a significant suppression of three- and four-pion Bose-Einstein correlations compared to expectations from two-pion measurements. There appears to be no centrality dependence of the suppression within the 0-50% centrality interval. The origin of the suppression is not clear. However, by postulating either coherent pion emission or large multibody Coulomb effects, the suppression may be explained.

The multi-strange baryon yields in Pb--Pb collisions have been shown to exhibit an enhancement relative to pp reactions. In this work, $\Xi$ and $\Omega$ production rates have been measured with the ALICE experiment as a function of transverse momentum, ${p_{\rm T}}$, in p-Pb collisions at a centre-of-mass energy of ${\sqrt{s_{\rm NN}}}$ = 5.02 TeV. The results cover the kinematic ranges 0.6 GeV/$c<{p_{\rm T}} <$7.2 GeV/$c$ and 0.8 GeV/$c<{p_{\rm T}}<$ 5 GeV/$c$, for $\Xi$ and $\Omega$ respectively, in the common rapidity interval -0.5 $<{y_{\rm CMS}}<$ 0. Multi-strange baryons have been identified by reconstructing their weak decays into charged particles. The ${p_{\rm T}}$ spectra are analysed as a function of event charged-particle multiplicity, which in p-Pb collisions ranges over one order of magnitude and lies between those observed in pp and Pb-Pb collisions. The measured ${p_{\rm T}}$ distributions are compared to the expectations from a Blast-Wave model. The parameters which describe the production of lighter hadron species also describe the hyperon spectra in high multiplicity p-Pb. The yield of hyperons relative to charged pions is studied and compared with results from pp and Pb-Pb collisions. A statistical model is employed, which describes the change in the ratios with volume using a canonical suppression mechanism, in which the small volume causes a species-dependent relative reduction of hadron production. The calculations, in which the magnitude of the effect depends on the strangeness content, show good qualitative agreement with the data.

The measurement of the azimuthal-correlation function of prompt D mesons with charged particles in pp collisions at $\sqrt{s}$ = 5.02 TeV and p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV with the ALICE detector at the LHC is reported. The D$^{\rm 0}$, D$^{\rm +}$, and D$^{\rm *+}$ mesons, together with their charge conjugates, were reconstructed at midrapidity in the transverse momentum interval 3 < $p_{\rm T}$ < 24 GeV/c and correlated with charged particles having $p_{\rm T}$ > 0.3 GeV/c and pseudorapidity $|\eta| <$ 0.8. The properties of the correlation peaks appearing in the near- and away-side regions (for $\Delta \varphi \approx$ 0 and $\Delta \varphi \approx \pi$, respectively) were extracted via a fit to the azimuthal correlation functions. The shape of the correlation functions and the near- and away-side peak features are found to be consistent in pp and p-Pb collisions, showing no modifications due to nuclear effects within uncertainties. The results are compared with predictions from Monte Carlo simulations performed with the PYTHIA, POWHEG+PYTHIA, HERWIG, and EPOS 3 event generators.

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.

In this Letter, comprehensive results on ${\rm\pi}^\pm$, K$^\pm$, K$^0_S$, p, $\rm\bar{p}$, $\rm \Lambda$ and $\rm \bar{\Lambda}$ production at mid-rapidity ($0 < y_{\rm cms} < 0.5$) in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, measured by the ALICE detector at the LHC, are reported. The transverse momentum distributions exhibit a hardening as a function of event multiplicity, which is stronger for heavier particles. This behavior is similar to what has been observed in pp and Pb-Pb collisions at the LHC. The measured $p_{\rm T}$ distributions are compared to results at lower energy and with predictions based on QCD-inspired and hydrodynamic models.

The pseudorapidity density of charged particles, $\rm{d}\it{N}_\rm{ch}/\rm{d}\it{\eta}$, in p-Pb collisions has been measured at a centre-of-mass energy per nucleon-nucleon pair of $\sqrt{s_{\rm{NN}}}$ = 8.16 TeV at mid-pseudorapidity for non-single-diffractive events. The results cover 3.6 units of pseudorapidity, $|\eta|<1.8$. The $\rm{d}\it{N}_\rm{ch}/\rm{d}\it{\eta}$ value is $19.1\pm0.7$ at $|\eta|<0.5$. This quantity divided by $\langle N_\rm{part} \rangle/2$, is $4.73\pm0.20$, which is 9.5% higher than the corresponding value for p-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV. Measurements are compared with models based on different mechanisms for particle production. All models agree within uncertainties with data in the Pb-going side, while HIJING overestimates, showing a symmetric behaviour, and EPOS underestimates the p-going side of the $\rm{d}\it{N}_\rm{ch}/\rm{d}\it{\eta}$ distribution. Saturation-based models reproduce the distributions well for $\eta>-1.3$. The $\rm{d}\it{N}_\rm{ch}/\rm{d}\it{\eta}$ is also measured for different centrality estimators, based both on the charged-particle multiplicity and on the energy deposited in the Zero-Degree Calorimeters. A study of the implications of the large multiplicity fluctuations due to the small number of participants for systems like p-Pb in the centrality calculation for multiplicity-based estimators is discussed, demonstrating the advantages of determining the centrality with energy deposited near beam rapidity.

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.

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.

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.