A search for nonresonant Higgs boson pair production in the $b\bar{b}b\bar{b}$ final state is presented. The analysis uses 126 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}={13}$ TeV collected with the ATLAS detector at the Large Hadron Collider, and targets both the gluon-gluon fusion and vector-boson fusion production modes. No evidence of the signal is found and the observed (expected) upper limit on the cross-section for nonresonant Higgs boson pair production is determined to be 5.4 (8.1) times the Standard Model predicted cross-section at 95% confidence level. Constraints are placed on modifiers to the $HHH$ and $HHVV$ couplings. The observed (expected) $2\sigma$ constraints on the $HHH$ coupling modifier, $\kappa_\lambda$, are determined to be $[-3.5, 11.3]$ ($[-5.4, 11.4]$), while the corresponding constraints for the $HHVV$ coupling modifier, $\kappa_{2V}$, are $[-0.0, 2.1]$ ($[-0.1, 2.1]$). In addition, constraints on relevant coefficients are derived in the context of the Standard Model effective field theory and Higgs effective field theory, and upper limits on the $HH$ production cross-section are placed in seven Higgs effective field theory benchmark scenarios.

The correlations between flow harmonics $v_n$ for $n=2$, 3 and 4 and mean transverse momentum $[p_\mathrm{T}]$ in $^{129}$Xe+$^{129}$Xe and $^{208}$Pb+$^{208}$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.44$ TeV and 5.02 TeV, respectively, are measured using charged particles with the ATLAS detector. The correlations are sensitive to the shape and size of the initial geometry, nuclear deformation, and initial momentum anisotropy. The effects from non-flow and centrality fluctuations are minimized, respectively, via a subevent cumulant method and event activity selection based on particle production in the very forward rapidity. The results show strong dependences on centrality, harmonic number $n$, $p_{\mathrm{T}}$ and pseudorapidity range. Current models describe qualitatively the overall centrality- and system-dependent trends but fail to quantitatively reproduce all the data. In the central collisions, where models generally show good agreement, the $v_2$-$[p_\mathrm{T}]$ correlations are sensitive to the triaxiality of the quadruple deformation. The comparison of model to the Pb+Pb and Xe+Xe data suggests that the $^{129}$Xe nucleus is a highly deformed triaxial ellipsoid that is neither a prolate nor an oblate shape. This provides strong evidence for a triaxial deformation of $^{129}$Xe nucleus using high-energy heavy-ion collision.

Three-body nuclear forces play an important role in the structure of nuclei and hypernuclei and are also incorporated in models to describe the dynamics of dense baryonic matter, such as in neutron stars. So far, only indirect measurements anchored to the binding energies of nuclei can be used to constrain the three-nucleon force, and if hyperons are considered, the scarce data on hypernuclei impose only weak constraints on the three-body forces. In this work, we present the first direct measurement of the p$-$p$-$p and p$-$p$-\Lambda$ systems in terms of three-particle correlation functions carried out for pp collisions at $\sqrt{s} = 13$ TeV. Three-particle cumulants are extracted from the correlation functions by applying the Kubo formalism, where the three-particle interaction contribution to these correlations can be isolated after subtracting the known two-body interaction terms. A negative cumulant is found for the p$-$p$-$p system, hinting to the presence of a residual three-body effect while for p$-$p$-\Lambda$ the cumulant is consistent with zero. This measurement demonstrates the accessibility of three-baryon correlations at the LHC.

This article presents the first measurement of the interaction between charm hadrons and nucleons. The two-particle momentum correlations of $\mathrm{pD^-}$ and $\mathrm{\overline{p}D}^+$ pairs are measured by the ALICE Collaboration in high-multiplicity pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$. The data are compatible with the Coulomb-only interaction hypothesis within (1.1-1.5)$\sigma$. The level of agreement slightly improves if an attractive nucleon(N)$\overline{\mathrm{D}}$ strong interaction is considered, in contrast to most model predictions which suggest an overall repulsive interaction. This measurement allows for the first time an estimation of the 68% confidence level interval for the isospin $\mathrm{I}=0$ inverse scattering length of the $\mathrm{N\overline{D}}$ state ${f_{0,~\mathrm{I}=0}^{-1} \in [-0.4,0.9]~\mathrm{fm^{-1}}}$, assuming negligible interaction for the isospin $\mathrm{I}=1$ channel.

The production of non-prompt ${\rm D}^{0}$ mesons from beauty-hadron decays was measured at midrapidity ($\left| y \right| < 0.5$) in Pb-Pb collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_{\rm NN}}=5.02~\mathrm{TeV}$ with the ALICE experiment at the LHC. Their nuclear modification factor ($R_{\rm AA}$), measured for the first time down to $p_{\rm T}=1~\mathrm{GeV}/c$ in the $0-10$% and $30-50$% centrality classes, indicates a significant suppression, up to a factor of about three, for $p_{\rm T} > 5~\mathrm{GeV}/c$ in the $0-10$% central Pb-Pb collisions. The data are described by models that include both collisional and radiative processes in the calculation of beauty-quark energy loss in the quark-gluon plasma, and quark recombination in addition to fragmentation as a hadronisation mechanism. The ratio of the non-prompt to prompt ${\rm D}^{0}$-meson $R_{\rm AA}$ is larger than unity for $p_{\rm T} > 4~\mathrm{GeV}/c$ in the $0-10$% central Pb-Pb collisions, as predicted by models in which beauty quarks lose less energy than charm quarks in the quark-gluon plasma because of their larger mass.

The production yields of non-prompt $\mathrm{D_s}^{+}$ mesons, namely $\mathrm{D_s}^{+}$ mesons from beauty-hadron decays, were measured for the first time as a function of the transverse momentum ($p_{\rm T}$) at midrapidity ($|y|<0.5$) in central and semi-central Pb$-$Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}=5.02$ TeV with the ALICE experiment at the LHC. The $\mathrm{D_s}^{+}$ mesons and their charge conjugates were reconstructed from the hadronic decay channel $\mathrm{D_s}^{+} \to \phi\pi^+$, with $\phi \to {\rm K}^-{\rm K}^+$, in the $4 < p_{\rm T}<36$ GeV/$c$ and $2 < p_{\rm T}<24$ GeV/$c$ intervals for the 0$-$10% and 30$-$50% centrality classes, respectively. The measured yields of non-prompt $\mathrm{D_s}^{+}$ mesons are compared to those of prompt $\mathrm{D_s}^{+}$ and non-prompt $\mathrm{D^0}$ mesons by calculating the ratios of the production yields in Pb$-$Pb collisions and the nuclear modification factor $R_\mathrm{AA}$. The ratio between the $R_\mathrm{AA}$ of non-prompt $\mathrm{D_s}^{+}$ and prompt $\mathrm{D_s}^{+}$ mesons, and that between the $R_\mathrm{AA}$ of non-prompt $\mathrm{D_s}^{+}$ and non-prompt $\mathrm{D^0}$ mesons in central Pb$-$Pb collisions are found to be on average higher than unity in the $4< p_{\rm T}<12$ GeV/$c$ interval with a statistical significance of about $1.6\,\sigma$ and $1.7\,\sigma$, respectively. The measured $R_\mathrm{AA}$ ratios are compared with the predictions of theoretical models of heavy-quark transport in a hydrodynamically expanding QGP that incorporate hadronisation via quark recombination.

The production of prompt $D^{0}$, $D^{+}_{\rm s}$, and $\Lambda_{\rm c}^{+}$ hadrons, and their ratios, $D^{+}_{\rm s}$/$D^{0}$ and $\Lambda_{\rm c}^{+}$/$D^{0}$, are measured in proton-proton collisions at $\sqrt{s}$ = 13 TeV at midrapidity ($|y| <0.5$) with the ALICE detector at the LHC. The measurements are performed as a function of the charm-hadron transverse momentum ($p_{\rm T}$) in intervals of charged-particle multiplicity, measured with two multiplicity estimators covering different pseudorapidity regions. While the strange to non-strange $D^{+}_{\rm s}$/$D^{0}$ ratio indicates no significant multiplicity dependence, the baryon-to-meson $p_{\rm T}$-differential $\Lambda_{\rm c}^{+}$/$D^{0}$ ratio shows a multiplicity-dependent enhancement, with a significance of 5.3$\sigma$ for $1< p_{\rm T} < 12$ GeV/$c$, comparing the highest multiplicity interval with respect to the lowest one. The measurements are compared with a theoretical model that explains the multiplicity dependence by a canonical treatment of quantum charges in the statistical hadronisation approach, and with predictions from event generators that implement colour reconnection mechanisms beyond the leading colour approximation to model the hadronisation process. The $\Lambda_{\rm c}^{+}$/$D^{0}$ ratios as a function of $p_{\rm T}$ present a similar shape and magnitude as the $\Lambda/K^{0}_{s}$ ratios in comparable multiplicity intervals, suggesting a potential common mechanism for light- and charm-hadron formation, with analogous multiplicity dependence. The $p_{\rm T}$-integrated ratios, extrapolated down to $p_{\rm T}$=0, do not show a significant dependence on multiplicity within the uncertainties.

Understanding the production mechanism of light (anti)nuclei is one of the key challenges of nuclear physics and has important consequences for astrophysics, since it provides an input for indirect dark-matter searches in space. In this paper, the latest results about the production of light (anti)nuclei in pp collisions at $\sqrt{s} = 13$ TeV are presented, focusing on the comparison with the predictions of coalescence and thermal models. For the first time, the coalescence parameters $B_2$ for deuterons and $B_3$ for helions are compared with parameter-free theoretical predictions that are directly constrained by the femtoscopic measurement of the source radius in the same event class. A fair description of the data with a Gaussian wave function is observed for both deuteron and helion, supporting the coalescence mechanism for the production of light (anti)nuclei in pp collisions. This method paves the way for future investigations of the internal structure of more complex nuclear clusters, including the hypertriton.

The first measurements of elliptic flow of $\pi^\pm$, ${\rm K}^\pm$, p+$\overline{\rm p}$, ${\rm K_{S}^0}$, $\Lambda$+$\overline{\Lambda}$, $\phi$, $\Xi^-$+$\Xi^+$, and $\Omega^-$+$\Omega^+$ using multiparticle cumulants in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV are presented. Results obtained with two- ($v_2\{2\}$) and four-particle cumulants ($v_2\{4\}$) are shown as a function of transverse momentum, $p_{\rm T}$, for various collision centrality intervals. Combining the data for both $v_2\{2\}$ and $v_2\{4\}$ also allows us to report the first measurements of the mean elliptic flow, elliptic flow fluctuations, and relative elliptic flow fluctuations for various hadron species. These observables probe the event-by-event eccentricity fluctuations in the initial state and the contributions from the dynamic evolution of the expanding quark-gluon plasma. The characteristic features observed in previous $p_{\rm T}$-differential anisotropic flow measurements for identified hadrons with two-particle correlations, namely the mass ordering at low $p_{\rm T}$ and the approximate scaling with the number of constituent quarks at intermediate $p_{\rm T}$, are similarly present in the four-particle correlations and the combinations of $v_2\{2\}$ and $v_2\{4\}$. In addition, a particle species dependence of flow fluctuations is observed that could indicate a significant contribution from final state hadronic interactions. The comparison between experimental measurements and CoLBT model calculations, which combine the various physics processes of hydrodynamics, quark coalescence, and jet fragmentation, illustrates their importance over a wide $p_{\rm T}$ range.

The study of nuclei and antinuclei production has proven to be a powerful tool to investigate the formation mechanism of loosely bound states in high-energy hadronic collisions. The first measurement of the production of ${\rm ^{3}_{\Lambda}\rm H}$ in p-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV is presented in this Letter. Its production yield measured in the rapidity interval $-1 < y < 0$ for the 40% highest multiplicity p-Pb collisions is ${\rm d} N /{\rm d} y =[\mathrm{6.3 \pm 1.8 (stat.) \pm 1.2 (syst.) ] \times 10^{-7}}$. The measurement is compared with the expectations of statistical hadronisation and coalescence models, which describe the nucleosynthesis in hadronic collisions. These two models predict very different yields of the hypertriton in charged particle multiplicity environments relevant to small collision systems such as p-Pb and therefore the measurement of ${\rm d} N /{\rm d} y$ is crucial to distinguish between them. The precision of this measurement leads to the exclusion with a significance larger than 6.9$\sigma$ of some configurations of the statistical hadronization model, thus constraining the theory behind the production of loosely bound states at hadron colliders.