ALICE is a large experiment at the CERN Large Hadron Collider. Located 52 meters underground, its detectors are suitable to measure muons produced by cosmic-ray interactions in the atmosphere. In this paper, the studies of the cosmic muons registered by ALICE during Run 2 (2015--2018) are described. The analysis is limited to multimuon events defined as events with more than four detected muons ($N_\mu>4$) and in the zenith angle range $0^{\circ}<\theta<50^{\circ}$. The results are compared with Monte Carlo simulations using three of the main hadronic interaction models describing the air shower development in the atmosphere: QGSJET-II-04, EPOS-LHC, and SIBYLL 2.3d. The interval of the primary cosmic-ray energy involved in the measured muon multiplicity distribution is about $ 4 \times 10^{15}<E_\mathrm{prim}< 6 \times 10^{16}$~eV. In this interval none of the three models is able to describe precisely the trend of the composition of cosmic rays as the energy increases. However, QGSJET-II-04 is found to be the only model capable of reproducing reasonably well the muon multiplicity distribution, assuming a heavy composition of the primary cosmic rays over the whole energy range, while SIBYLL 2.3d and EPOS-LHC underpredict the number of muons in a large interval of multiplicity by more than $20\%$ and $30\%$, respectively. The rate of high muon multiplicity events ($N_\mu>100$) obtained with QGSJET-II-04 and SIBYLL 2.3d is compatible with the data, while EPOS-LHC produces a significantly lower rate ($55\%$ of the measured rate). For both QGSJET-II-04 and SIBYLL 2.3d, the rate is close to the data when the composition is assumed to be dominated by heavy elements, an outcome compatible with the average energy $E_\mathrm{prim} \sim 10^{17}$ eV of these events. This result places significant constraints on more exotic production mechanisms.

We report on the measurement of inclusive, non-prompt, and prompt J/$ψ$-hadron correlations by the ALICE Collaboration at the CERN Large Hadron Collider in pp collisions at a center-of-mass energy of $13~\text{TeV}$. The correlations are studied at midrapidity ($|y| < 0.9$) in the transverse momentum ranges $p_{\rm T} < 40~\text{GeV}/c$ for the J/$ψ$ and $0.15 < p_{\rm T} < 10~\text{GeV}/c$ and $|η|<0.9$ for the associated hadrons. The measurement is based on minimum bias and high multiplicity data samples corresponding to integrated luminosities of $L_{\text{int}} = 34~\text{nb}^{-1}$ and $L_{\text{int}} = 6.9~\text{pb}^{-1}$, respectively. In addition, two more data samples are employed, requiring, on top of the minimum bias condition, a threshold on the tower energy of $E = 4$ and $9~\text{GeV}$ in the ALICE electromagnetic calorimeters, which correspond to integrated luminosities of $L_{\text{int}} = 0.9~\text{pb}^{-1}$ and $L_{\text{int}} = 8.4~\text{pb}^{-1}$, respectively. The azimuthally integrated near and away side yields of associated charged hadrons per J/$ψ$ trigger are presented as a function of the J/$ψ$ and associated hadron transverse momentum. The measurements are discussed in comparison to PYTHIA calculations.

Differential two-particle normalized cumulants ($R_2$) and transverse momentum correlations ($P_2$) are measured as a function of the relative pseudorapidity and azimuthal angle difference $( Δη, Δφ)$ of charged particle pairs in minimum bias pp collisions at $\sqrt{\textit{s}}$ = 13 TeV. The measurements use charged hadrons in the pseudorapidity region of $|η| < 0.8$ and the transverse momentum range \mbox{0.2 $< \textit{p}_{\mathrm T} < $ 2.0 $\mathrm{GeV}/\textit{c}$} in order to focus on soft multiparticle interactions and to complement prior measurements of these correlation functions in p-Pb and Pb-Pb collisions. The correlation functions are reported for both unlike-sign and like-sign pairs and their charge-independent and charge-dependent combinations. Both the $R_2$ and $P_2$ measured in pp collisions exhibit features qualitatively similar to those observed in p--Pb and Pb--Pb collisions. The $Δη$ and $Δφ$ root mean square widths of the near-side peak of the correlation functions are evaluated and compared with those observed in p-Pb and Pb-Pb collisions and show smooth evolution with the multiplicity of charged particles produced in the collision. The comparison of the measured correlation functions with predictions from PYTHIA8 shows that this model qualitatively captures their basic structure and characteristics but feature important differences. In addition, the $R_2^{\rm CD}$ is used to determine the charge balance function of hadrons produced within the detector acceptance of the measurements. The integral of the balance function is found to be compatible with those reported by a previous measurement in Pb--Pb collisions.

Ultrarelativistic heavy-ion collisions create a quark-gluon plasma (QGP), a hot and dense state of strongly interacting QCD matter. In ultracentral collisions, the QGP volume remains nearly constant event-by-event, while its total entropy can fluctuate due to quantum effects, leading to temperature variations. These features allow the correlation between the mean transverse momentum $(\langle p_{\mathrm{T}} \rangle)$ of charged hadrons and their multiplicity to serve as a probe of the QGP's speed of sound, $c_{s}$. This study extracts $c_{s}$ by analyzing the relative increase in $\langle p_{\mathrm{T}} \rangle$ with respect to the charged-particle density $(\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle)$ at midrapidity in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$, using data from the ALICE detector. Centrality is determined with estimators based on multiplicity and transverse energy, applying a pseudorapidity gap to reduce selection biases. The extracted value of $c_{s}^{2}$ is found to strongly depend on the employed centrality estimator and ranges between $0.1146 \pm 0.0028 \,\mathrm{(stat.)} \pm 0.0065 \,\mathrm{(syst.)}$ and $0.4374 \pm 0.0006 \mathrm{(stat.)} \pm 0.0184 \mathrm{(syst.)}$ in natural units. Additionally, the event-by-event $[p_{\mathrm{T}}]$ distribution is studied through its variance, skewness, and kurtosis. A pronounced decrease in the self-normalized variance and a peak followed by a drop in skewness suggest the suppression of impact-parameter fluctuations in ultracentral collisions. These observations provide new insights into the thermodynamic properties and initial-state fluctuations of the QGP.

The transverse momentum spectra and integrated yields of $\overlineΣ^{\pm}$ have been measured in pp and p-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV with the ALICE experiment. Measurements are performed via the newly accessed decay channel $\overlineΣ^{\pm} \rightarrow {\rm\overline{n}}π^{\pm}$. A new method of antineutron reconstruction with the PHOS electromagnetic spectrometer is developed and applied to this analysis. The $p_{\rm T}$ spectra of $\overlineΣ^{\pm}$ are measured in the range $0.5 < p_{\rm T} < 3$ GeV/$c$ and compared to predictions of the PYTHIA 8, DPMJET, PHOJET, EPOS LHC and EPOS4 models. The EPOS LHC and EPOS4 models provide the best descriptions of the measured spectra both in pp and p-Pb collisions, while models which do not account for multiparton interactions provide a considerably worse description at high $p_{\rm T}$. The total yields of $\overlineΣ^{\pm}$ in both pp and p-Pb collisions are compared to predictions of the Thermal-FIST model and dynamical models PYTHIA 8, DPMJET, PHOJET, EPOS LHC and EPOS4. All models reproduce the total yields in both colliding systems within uncertainties. The nuclear modification factors $R_{\rm pPb}$ for both $\overlineΣ^{+}$ and $\overlineΣ^{-}$ are evaluated and compared to those of protons, $Λ$ and $Ξ$ hyperons, and predictions of EPOS LHC and EPOS4 models. No deviations of $R_{\rm pPb}$ for $\overlineΣ^{\pm}$ from the model predictions or measurements for other hadrons are found within uncertainties.

The coherent J/$\psi$ photoproduction cross section is measured for the first time at midrapidity in peripheral to semicentral Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV. The centrality differential cross section ${\rm d} \sigma/ {\rm d}y$ is reported for the centrality range 40-90%, together with the doubly-differential cross section ${\rm d}^2 \sigma /{\rm d}y {\rm d} p_{\rm T}$, extracted in two peripheral centrality classes. The J/$\psi$ mesons are reconstructed in the dielectron channel, in the rapidity interval $|y| <$ 0.9 using the ALICE central barrel detectors. The J/$\psi$ cross section at midrapidity is statistically compatible to the earlier ALICE measurement at forward rapidity and at the same centre-of-mass energy, and shows only a mild centrality dependence over the covered range. Several sets of theoretical calculations taking into account the hadronic overlap in the collisions but ignoring possible final-state effects from a hot expanding medium are found to give a fairly good description of the current measurements within uncertainties.

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.

A search for heavy Higgs bosons produced in association with a vector boson and decaying into a pair of vector bosons is performed in final states with two leptons (electrons or muons) of the same electric charge, missing transverse momentum and jets. A data sample of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded with the ATLAS detector at the Large Hadron Collider between 2015 and 2018 is used. The data correspond to a total integrated luminosity of 139 fb$^{-1}$. The observed data are in agreement with Standard Model background expectations. The results are interpreted using higher-dimensional operators in an effective field theory. Upper limits on the production cross-section are calculated at 95% confidence level as a function of the heavy Higgs boson's mass and coupling strengths to vector bosons. Limits are set in the Higgs boson mass range from 300 to 1500 GeV, and depend on the assumed couplings. The highest excluded mass for a heavy Higgs boson with the coupling combinations explored is 900 GeV. Limits on coupling strengths are also provided.

This paper presents a search for a new Z' vector gauge boson with the ATLAS experiment at the Large Hadron Collider using pp collision data collected at $\sqrt{s} = 13$ TeV, corresponding to an integrated luminosity of 139 fb$^{-1}$. The new gauge boson Z' is predicted by $L_{\mu}-L_{\tau}$ models to address observed phenomena that can not be explained by the Standard Model. The search examines the four-muon (4$\mu$) final state, using a deep learning neural network classifier to separate the Z' signal from the Standard Model background events. The di-muon invariant masses in the $4\mu$ events are used to extract the Z' resonance signature. No significant excess of events is observed over the predicted background. Upper limits at a 95% confidence level on the Z' production cross-section times the decay branching fraction of $pp \rightarrow Z'\mu\mu \rightarrow 4\mu$ are set from 0.31 to 4.3 fb for the Z' mass ranging from 5 to 81 GeV. The corresponding common coupling strengths, $g_{Z'}$, of the Z' boson to the second and third generation leptons above 0.003 - 0.2 have been excluded.

Inclusive and differential measurements of the top-antitop ($t\bar{t}$) charge asymmetry $A_\text{C}^{t\bar{t}}$ and the leptonic asymmetry $A_\text{C}^{\ell\bar{\ell}}$ are presented in proton-proton collisions at $\sqrt{s} = 13$ TeV recorded by the ATLAS experiment at the CERN Large Hadron Collider. The measurement uses the complete Run 2 dataset, corresponding to an integrated luminosity of 139 fb$^{-1}$, combines data in the single-lepton and dilepton channels, and employs reconstruction techniques adapted to both the resolved and boosted topologies. A Bayesian unfolding procedure is performed to correct for detector resolution and acceptance effects. The combined inclusive $t\bar{t}$ charge asymmetry is measured to be $A_\text{C}^{t\bar{t}} = 0.0068 \pm 0.0015$, which differs from zero by 4.7 standard deviations. Differential measurements are performed as a function of the invariant mass, transverse momentum and longitudinal boost of the $t\bar{t}$ system. Both the inclusive and differential measurements are found to be compatible with the Standard Model predictions, at next-to-next-to-leading order in quantum chromodynamics perturbation theory with next-to-leading-order electroweak corrections. The measurements are interpreted in the framework of the Standard Model effective field theory, placing competitive bounds on several Wilson coefficients.