Abstract
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.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
One of the open challenges of nuclear physics is the understanding of many-particle dynamics. Studies of the nuclear structure have unambiguously shown that calculations based only on nucleon–nucleon (N–N) interactions fail to accurately describe many experimental observables, such as nuclear binding energies along the periodic table of elements [1], the position of the neutron drip line for neutron-rich nuclei [2] or the properties of the recently observed four-neutrons resonance [3]. A significant improvement in the modelling of nuclear bound objects has been achieved by including three-body forces in theoretical calculations. These three-body forces are implemented in chiral effective field theories [4] and in a number of ab initio many-body methods such as no-core shell model [5], coupled-cluster theory [6, 7], self-consistent Green’s function theory [8], similarity renormalisation group [9, 10], and quantum Monte Carlo [11]. Studies conducted on intermediate mass neutron-rich nuclei proved that the sensitivity to the three-body forces increases with the number of neutrons in the system [6]. Three-body forces within light and medium-mass nuclei, where the nuclear saturation density corresponds to typical inter-particle distances of 2 fm, contribute about 10–15% to the total interaction strength [12, 13]. However, at higher densities and shorter inter-particle distances their contribution might increase [2], but no data are available in such a regime and the properties of nuclear matter can be only extrapolated using the available information at saturation densities. The experimental information on the three-body forces involving \(\Lambda \) hyperons is even more scarce since the data available for hypernuclei are much less than the data for nuclei. Recent hypetriton measurements in several colliding systems at RHIC and LHC [14,15,16,17] provide important input to the understanding of N–N–\(\Lambda \) forces and future measurements will resolve the current tensions among the different estimations of the binding energy and life-time. Theoretical works assign to the hypertriton a radius of the order of 5 fm [18] and hence a N–\(\Lambda \) distance of 10 fm [13, 19] within this state. Heavier hypernuclei are more compact and their size is comparable to that of normal nuclei so that they represent an optimal test bed for the N–N–\(\Lambda \) interaction [20]. However, good fits of the theoretical models to the available hypernuclear data, from \({}_\Lambda ^{7}\)Li to \({}_\Lambda ^{208}\)Pb [21,22,23,24], require a full understanding of the shell-structure of such bound objects as well as accurate experimental constraints on the spin-dependent N–\(\Lambda \) interaction, in particular for the p-wave and higher partial waves. The lack of precise data as well as the limitations in the microscopic description of the structure of hypernuclei cause large ambiguities on the strength of the N–N–\(\Lambda \) three-body force. Further opportunities are provided by the recently observed \(^{3}_{\Lambda }\)n bound state [25] and planned experimental programs focused on neutron-rich hypernuclei [26]. Nevertheless, the contribution from three-body forces in bound objects such as nuclei and hypernuclei cannot be separated from the lower-order two-body interactions, hence, complementary experimental methods to investigate three-baryon systems could provide an important contribution to this field.
Neutron-rich and dense baryonic matter constitutes an interesting system also because of its connections to the physics of neutron stars (NS) [27]. The structure and composition of the innermost part of NS is not known. Amongst many possible scenarios, some models support the appearance of various hadronic particle species with increasing baryon density inside the star [27, 28]. The presence of hadronic degrees of freedom and their relative abundances are sensitive to the two- and three-body interaction models which are used to compute the equation of state (EoS) of NS matter. The different hypotheses can be tested by deriving the masses and radii of NS for a specific EoS and comparing them with the corresponding astrophysical observations [28]. The suggestion of strange baryons inside NS is motivated by the fact that central densities of NS might become sufficiently large (\(\rho \approx 3-4 \rho _0\), where \(\rho _0\) is the nuclear saturation density) to provide favourable conditions for the onset of strangeness production processes leading to, in particular, the formation of hyperons. The appearance of \(\Lambda \) hyperons in NS matter results in a softening of the EoS which is at variance with astrophysical observations of two solar mass stars [29, 30]. However, in Ref. [31], it was shown that by adding a strongly repulsive N–N–\(\Lambda \) interaction, tuned to reproduce the separation energies of \(\Lambda \) hyperons in several hypernuclei, a sufficiently stiff EoS can be obtained and even the massive NS observables can be reproduced. This indicates that three-body forces may have a significant contribution in models that describe the structure of NS. Hence, a direct measurement of the three-body forces involving nucleons and hyperons at small inter-particle distances is required.
The femtoscopy technique can be used as a tool to investigate the strong interaction amongst hadrons produced in particle collisions [32,33,34,35] and recently has been successfully employed to analyse experimental data. The produced hadrons may undergo final state interactions (FSIs) and the resulting correlation in the momentum space can be studied to test the underlying dynamics using correlation functions [34, 35]. The method has been applied by the STAR Collaboration to measure hadron-hadron correlations in Au–Au collisions with a centre-of-mass energy of 200 GeV per nucleon pair [36,37,38]. In such ultra-relativistic heavy-ion collisions, the average relative distances of emitted particles is about 7–8 fm [35]. In small colliding systems, such as pp and p–Pb collisions at the LHC, particles are produced at distances of the order of 1 fm, hence, the sensitivity of the correlation function to the short-range strong interaction is enhanced. Recently, the method has been employed by ALICE to study FSIs of hadrons produced in such small colliding systems. The large data samples allowed for the precise measurement of correlation functions for multiple hadronic pairs (p–p [39], p–K\(^+\) and p–K\(^-\) [40], p–\(\mathrm {\Lambda }\) [39], p–\(\mathrm {\Sigma }^0\) [41], \(\mathrm {\Lambda }\)–\(\mathrm {\Lambda }\) [42], p–\(\mathrm {\Xi }^-\) [43], p–\(\mathrm {\Omega }^-\) [44], p–\(\phi \) [45] and baryon–antibaryon [46]). By using these results, several models for the two-body strong interaction could be validated (for a complete review see Ref. [47]).
The femtoscopy technique was also employed in the analysis of three- and four-pion correlations measured in pp, p–Pb, Pb–Pb collision systems by ALICE [48, 49] to probe coherent hadron production. The Kubo’s cumulant expansion method [50] was used to isolate the genuine three-particle correlation from the two-body contributions where the latter were evaluated by combining two particles from the same event and a third particle taken from another event. Alternatively, the recently developed projector method [51], where either the theoretical or the measured two-body correlation functions are used to obtain the lower-order contributions, can be employed. This method allows a significant reduction of the statistical uncertainties.
In this article, the first femtoscopic study of three-baryon correlations is performed for the p–p–p and p–p–\(\Lambda \) systems measured in high-multiplicity (HM) pp collisions at \(\sqrt{s} = 13\) TeV. The Kubo’s formalism and the projector method are employed to isolate the genuine three-body correlation and the choice of the reaction system aims to study the interaction at small distances. The article is organised as follows: in Sect. 2.1 the data analysis procedure is presented starting from the event selection; in Sect. 2.2 the definition of the two-particle correlation function is extended to the three-particle case; in Sect. 2.3 the femtoscopic three-particle cumulant is defined; the lower-order two-particle correlation contributions in the measured correlation functions are evaluated in Sect. 2.4 and the decomposition of the cumulant to account for misidentifications and particle feed-down are presented in Section 2.5; the final results are discussed in Sect. 3 and the conclusions are given in Sect. 4.
2 Analysis
2.1 Event selection and particle identification
The data sample of pp collisions at a centre of mass energy \(\sqrt{s}~=~13\) TeV was recorded with the ALICE detector [52, 53] during the LHC Run 2 (2015–2018). The sample has been collected employing a HM trigger. The trigger is based on the measured amplitude in the V0 detector system, consisting of two arrays of plastic scintillators located at forward (\(2.8<\eta <5.1\)) and backward (\(-3.7<\eta <-1.7\)) pseudorapidities [54]. The selected HM events correspond to the highest 0.17% multiplicity interval with respect to all inelastic collisions with at least one measured charged particle within \(|\eta |<1\) (INEL\(>0\)). This condition results in an average of 30 charged particles in the range \(|\eta |<0.5\) [44]. Charged-particle tracking in the midrapidity region is conducted with the Inner Tracking System (ITS) [52] and the Time Projection Chamber (TPC) [55]. These detectors are immersed in a homogeneous 0.5 T magnetic field parallel to the beam direction. The ITS consists of six cylindrical layers of high position-resolution silicon detectors placed radially between 3.9 and 43 cm around the beam vacuum tube. The TPC consists of a 5 m long, cylindrical gaseous detector with full azimuthal coverage in the pseudorapidity range \(|\eta | < 0.9\).
Particle identification (PID) is conducted via the measurement of the specific ionisation energy loss (\(\text {d}E/\text {d}x\)) in the TPC gas with up to 159 reconstructed space points along the particle trajectory. For high momentum particles, the TPC measurement is combined with information provided by the time-of-flight (TOF) [56] detector system, which is located at a radial distance of 3.7 m from the nominal interaction point and consists of multigap resistive plate chambers covering the full azimuthal angle in \(|\eta |<0.9\).
The primary vertex (PV) of the event is reconstructed with the combined track information of the ITS and the TPC, and independently with track segments in the two innermost layers of the ITS. The reconstructed PV of the event is required to have a maximal displacement with respect to the nominal interaction point of 10 cm along the beam axis, in order to ensure a uniform acceptance. Pile-up events with multiple primary vertices are removed following the procedure described in Refs. [39, 43, 57]. This rejects the events with pile-up of collisions occurring in the same or nearby bunch crossings. However, additional clean-up has to be applied on the track selection level to reject particles produced in pile-up collisions in the long TPC readout time.
A total of \(1.0\times 10 ^9\) HM events are used for the analysis after event selection. In order to build the three-particle correlation functions of p–p–p and p–p–\(\Lambda \) systems, particle and antiparticle distributions are combined. In the following, p–p–p refers to p–p–p \(\oplus \) \(\overline{\textrm{p}}\text{-- }\overline{\textrm{p}}\text{-- }\overline{\textrm{p}}\) and p–p–\(\Lambda \) refers to p–p–\(\Lambda \) \(\oplus \) \(\overline{\textrm{p}}\text{-- }\overline{\textrm{p}}\text{-- }\overline{\Lambda } \). The proton and \(\Lambda \) candidates as well as their antiparticles need to be selected. As the particle and antiparticle selections are identical, only the particles are explicitly discussed below. Both particle species are reconstructed using the procedure described in Ref. [57], while the related systematic uncertainties are evaluated by varying the kinematic and topological selection criteria used in the reconstruction. In the following text, the systematic variations are enclosed in parentheses.
The primary protons are selected in the momentum interval \(0.5~(0.4,~0.6) {<} p_{\textrm{T}} {<} 4.05\) \(\text {Ge\hspace{-1.00006pt}V/}c\) and \(|\eta | {<} 0.8~(0.77, 0.85)\). To improve the quality of the tracks a minimum of 80 (70, 90) out of the 159 possible spatial points inside the TPC are required. The PID selections are applied by comparing the measured \(\text {d}E/\text {d}x\) and time-of-flight with the expected values for a proton candidate. The agreement is expressed in multiples (\(n^\text {PID}_\sigma \)) of the detector resolution \(\sigma \). For protons with \(p_{\textrm{T}} < 0.75\) \(\text {Ge\hspace{-1.00006pt}V/}c\) the \(n^\text {PID}_\sigma \) is evaluated only based on the specific energy loss in the TPC, while for \(p_{\textrm{T}} \ge 0.75\) \(\text {Ge\hspace{-1.00006pt}V/}c\) a combined TPC and TOF PID selection is applied \(\left( n^\text {PID}_\sigma =\sqrt{n_{\sigma ,\textrm{TPC}}^2+n_{\sigma ,\textrm{TOF}}^2}\right) \). The \(n^\text {PID}_\sigma \) of the accepted proton candidates is required to be lower than 3 (2.5, 3.5). To reject particles that are non-primary or come from pile-up collisions, the distance of closest approach (DCA) to the PV of the tracks is required to be less than 0.1 cm in the transverse plane and less than 0.2 cm along the beam axis. The purity of candidates is estimated using Monte Carlo (MC) simulations by taking the ratio of the number of reconstructed true protons produced by the generator and the number of all candidates identified as protons as a function of the reconstructed transverse momentum. The contributions of secondary protons stemming from weak decays of strange baryons and from interactions in the detector material are extracted using MC template fits to the measured distributions of the DCA to the PV [39]. The average purity of the identified protons is \(98.3\%\) and \(86.6\%\) of them are primaries.
The \(\Lambda \) candidates are reconstructed via the weak decay \(\Lambda \rightarrow \textrm{p}\pi ^-\) (the \(\overline{\Lambda }\rightarrow \overline{\textrm{p}}\pi ^+\) in case of \(\overline{\Lambda }\) reconstruction). The secondary daughter tracks are selected with similar criteria as for the primary protons regarding \(|\eta |\) and the number of hits in the TPC. However, a less strict PID requirement of \(n^\text {PID}_\sigma <5(4)\) is used. In addition, the daughter tracks are required to have a DCA to the PV of at least 0.05 (0.06) cm and the DCA between the daughter tracks at the secondary vertex must be smaller than 1.5 (1.2) cm. The cosine of the pointing angle (CPA) between the vector connecting the PV to the decay vertex and the 3-momentum of the \(\Lambda \) candidate is required to be larger than 0.99 (0.995). To reject unphysical secondary vertices, reconstructed with tracks stemming from pile-up of pp collisions occurring in different bunch crossings, the decay tracks are required to possess a hit in the two innermost or the two outermost ITS layers or a matched TOF signal [42]. Finally, a selection on the candidate invariant mass (IM) is applied by requiring it to be in a \(\pm \,4\) MeV/\(c^2\) interval around the nominal \(\Lambda \) mass [58]. The primary and secondary contributions to the yield of \(\Lambda \) are extracted employing a similar method as for protons but using the CPA as an observable for the template fits. The \(\Lambda \) hyperons produced in primary interactions contribute to about 58.5% of their total yield. About 19.5% originate from the electromagnetic decays of \(\Sigma ^0\). The number of \(\Sigma ^0\) particles is related to their ratio to the \(\Lambda \) hyperons, which is fixed to 1/3 based on predictions from the isospin symmetry and a measurement of the corresponding production ratios [59]. Further, each of the weak decays of \(\Xi ^-\) and \(\Xi ^0\) contributes about 11 % to the yield of \(\Lambda \) hyperons. The purity of \(\Lambda \) and \(\overline{\Lambda }\) has been extracted by fitting the IM spectra of candidates as a function of the three-particle kinematic variable \(Q_3\) which is defined in Eq. 5. The fits have been performed in the IM range of 1090 to 1150 MeV/\(c^2\) using a double Gaussian for the \(\Lambda \) signal and a second-order polynomial for the background. The result has been averaged for \(Q_3 <1\) \(\text {Ge\hspace{-1.00006pt}V/}c\), leading to a combined purity of \(\Lambda \) and \(\overline{\Lambda }\) of 95.6%.
The systematic uncertainties are evaluated by performing simultaneous variations of the selection criteria for protons and \(\Lambda \) candidates as well as for the corresponding antiparticles. The variations are randomly combined in 44 sets in which at least one of the selection criteria is varied. Such procedure allows to account for the correlations between the systematic uncertainties. Each random set of variations is accepted for the evaluation of the systematic uncertainties only if the yield of the triplets is varied by less than 10% with respect to the standard selection in the kinematic region \(Q_3 <0.4\) \(\text {Ge\hspace{-1.00006pt}V/}c\).
2.2 Three-particle correlation function
The observable of interest in femtoscopy is usually the two-particle momentum correlation function [35, 60], which is defined as the probability to simultaneously find two particles with momenta \({\textbf{p}}_1\) and \({\textbf{p}}_2\) divided by the product of the corresponding single particle probabilities
These probabilities are related to the inclusive Lorentz-invariant spectra \(P({\textbf{p}}_1,{\textbf{p}}_2) {\propto } E_1 E_2 \frac{\textrm{d}^6N}{\textrm{d}^3{\textbf{p}}_1 \textrm{d}^3{\textbf{p}}_2}\) and \(P({\textbf{p}}_{i}) {\propto } E_{i}\frac{\textrm{d}^3N_{i}}{\textrm{d}^3{\textbf{p}}_{i}}\). In the absence of a correlation signal, the value of \(C({\textbf{p}}_1,{\textbf{p}}_2)\) is constant and normalised to unity. A similar logic can be followed to construct the three-particle correlation functions as
Following [61, 62], Eq. 1 can also be written as
where \(S(r^*)\) is the distribution of the relative distances of particle pairs in the pair rest frame (PRF, denoted by the \(^*\)) – the so-called source function. The properties of the source in pp collisions at \(\sqrt{s}\) = 13 TeV have been evaluated in Ref. [57], including the effects of short-lived resonance decays which enlarge the effective source size. The wave function of the particle pair relative motion is denoted by \(\psi ({\textbf{r}}^*,\mathbf {k^*})\) where \(\mathbf {k^*} = ({\textbf{p}}^*_1-{\textbf{p}}^*_2)/2\) is the relative momentum. The wave function encapsulates the details of the particle interaction and drives the shape of the correlation function. In case of the three-particle correlation function, the two-particle source function and the wave function of the particle pair relative motion must be replaced by a three-particle source function and wave function. In this analysis, the measured three-particle correlation functions are not compared to theoretical predictions. The goal here is to extract the three-particle femtoscopic cumulants which provide experimental evidence of the existence, or the absence, of genuine three-particle correlations, as explained in Sect. 2.3.
The three-particle correlation function can be written as
where \(N_{\textrm{s}}(Q_3)\) and \(N_{\textrm{m}}(Q_3)\) are the same-event and mixed-event distributions of three particle combinations (triplets) as a function of \(Q_3\) and \({\mathcal {N}}\) is the normalisation parameter. The Lorentz-invariant variable \(Q_3\) is defined in [48] as
where \(q_{ij}\) is the norm of the four-vector [35]
which can be rewritten as
Here \(m_{i}\) and \(m_{j}\) are the particle i and j masses, \(p_{i}^\mu \) and \(p_{j}^\mu \) are the particle four momenta, while \(q_{ij}^{\mu }\) is the relative four-momentum of the pair ij. In the case of same mass particles, the term \(\frac{\left( p_{i}-p_{j}\right) \cdot P_{ij}}{ P_{ij}^{2}} P_{ij}^{\mu }\) becomes 0. In the non-relativistic case \(q_{ij}^2 = - 4 {k^*_{ij}}^2\), where \(k^*_{ij}\) is the relative momentum of the ij pair in the PRF.
The mixed-event sample is obtained using event-mixing techniques, in which the particle triplets of interest are generated by combining single particles stemming from three different events. To maintain the same acceptance effects as in the same event sample, the mixing procedure is conducted only for events with similar z position of the primary vertex and multiplicity [39]. Additionally, in order to correct for possible differences in terms of multiplicity distribution between same and mixed events, the yield of the latter is re-weighted in each multiplicity interval to have the same statistical weight as the distribution when particles are from the same event. To account for the two-track merging and splitting effects due to the finite two-track resolution in the same-event sample, a minimum value of the distance between two proton tracks (in case of p–\(\Lambda \) pairs, the proton from \(\Lambda \) decay is considered along with the primary proton) on the azimuthal-polar angles plane \(\Delta \eta \)–\(\Delta \varphi \) is applied to both the same- and mixed-event samples. The default selection is \(\Delta \eta ^2 + \Delta \varphi ^2 \ge 0.017^2\) and a systematic variation of \(+10\) % for the value of the minimum distance is applied in the analysis. The normalisation parameter \({\mathcal {N}}\) is chosen such that the mean value of the correlation function equals unity in a \(Q_3\) region where the effects of FSIs are negligible. The interval \(Q_3\in (1.0 - 1.2)\) GeV/\(\textit{c}\) is chosen for all triplets.
2.3 Three-particle femtoscopic cumulants
The measurable three-particle correlation function \(C({\textbf{p}}_1,{\textbf{p}}_2, {\textbf{p}}_3)\) include all interactions at work in the three-particle system: the two-body interactions among all pairs within the selected triplet and the genuine three-body interaction. To access only the genuine three-body correlations, one can use cumulants. Given random variables \(X_i\), the cumulant for a triplet is defined by Kubo [50] as
where \(\left\langle X_{i}\right\rangle \) is the expectation value of the variable \(X_i\) and \(\left\langle X_{i} X_{j}\right\rangle \), \(\left\langle X_{i} X_{j} X_{k}\right\rangle \) are the two- and three-variable joint moments. The three-particle correlation function, defined in Eq. 4, is the three-particle momentum distribution normalised to the mixed-event distribution. The cumulants method can be applied to the numerator which contains the correlated particles, and then the expression is normalised to the mixed-event distribution. The three-particle femtoscopic cumulant \({c}_{3}\) thus can be defined as
where \(N_{3}\left( {\textbf {p}}_{1}, {\textbf {p}}_{2}, {\textbf {p}}_{3}\right) \) and \(N_{2}\left( {\textbf {p}}_{i}, {\textbf {p}}_{j}\right) \) are the same-event three- and two-particle momentum distributions; \(N_{1}\left( {\textbf {p}}_{i}\right) \) is the single-particle momentum distribution; the product terms \(N_{2}\left( {\textbf {p}}_{i}, {\textbf {p}}_{j}\right) N_{1}\left( {\textbf {p}}_{k}\right) \) and \(N_{1}\left( {\textbf {p}}_{i}\right) N_{1}\left( {\textbf {p}}_{j}\right) N_{1}\left( {\textbf {p}}_{k}\right) \) indicate the mixed event distributions. Thus one can further rewrite the femtoscopic cumulant as
This method has been already successfully applied within the ALICE Collaboration to study the possibility of coherent pion production by measuring three-pion femtoscopic cumulants in Refs. [48, 49]. Theorem I from Ref. [50] enunciates that the three-particle cumulant is zero if the variables \(X_{i}, X_{j}, \ldots \) can be divided into two or more groups that are statistically independent. In case of femtoscopic cumulants, this translates into \(c_{3}\left( {\textbf {p}}_{1}, {\textbf {p}}_{2}, {\textbf {p}}_{3}\right) =0\) in the absence of genuine three-body correlations. Therefore, the measurements of non-vanishing values of \(c_3\) can be used as an experimental confirmation of the existence of genuine three-body effects.
If genuine three-body correlations are not present in the particle triplet, the three-particle correlation function can be expressed using only lower order contributions as follows
In Eq. 11, \(C([{\textbf{p}}_i,{\textbf{p}}_j],{\textbf{p}}_k)\) is built by combining particles i and j from the same event with particle k from another event to obtain the numerator \(\left( N_{2}\left( {\textbf {p}}_{i}, {\textbf {p}}_{j}\right) N_{1}\left( {\textbf {p}}_{k}\right) \right) \) of the correlation function while the denominator \(\left( N_{1}\left( {\textbf {p}}_{1}\right) N_{1}\left( {\textbf {p}}_{2}\right) N_{1}\left( {\textbf {p}}_{3}\right) \right) \) is estimated using three particles from three different events as described in Sect. 2.2.
2.4 Projector method
An alternative method to isolate the genuine three-body contribution to the measured three-particle correlation functions is the projector method [51]. This method makes use of the subtraction rule provided by the Kubo’s cumulant decomposition (Eq. 10) but, instead of evaluating them with the data-driven approach based on event mixing described above, it calculates \(C([{\textbf{p}}_i,{\textbf{p}}_j],{\textbf{p}}_k)\) using the measured or the calculated two-particle correlation function and the projection of the third non-interacting (spectator) particle. The method is described in Ref. [51]. Given the three-particle correlation function, \(C(Q_3)\), and the two-body correlation functions, \(C(k^{*}_{ij})\), the projector method provides a kinematic transformation from the relative momentum \(k^{*}_{ij}\) of the interacting pairs ij to the \(Q_3\) of the three-body system \((i-j)-k\) under study. The transformation is given by the following integral in the momentum space
where the indices ij denote the interacting pair and the projector function \(W_{ij}\) is equal to [51]
The constants \(\alpha \), \(\beta \) and \(\gamma \) depend on the particle masses.Footnote 1
The integral in Eq. 12 can be evaluated using the measured p–p and p–\(\Lambda \) correlation functions from Refs. [41, 57, 63]. The resulting correlation functions are compared to the ones obtained by employing the data-driven method (Eq. 11) and shown in Fig. 1. Panel (a) shows the (p–p)–p correlation function, the green points are obtained using the data-driven approach and the grey band is obtained with the projector method. The statistical and systematic uncertainties are shown separately for the data driven method, while the width of the grey band represents the sum in quadrature of the statistical and systematic uncertainties for the projector method. The statistical uncertainties of all the measured correlation functions have been estimated using a bootstrap [63] method by sampling same- and mixed-event counts from Poisson distributions. The statistical uncertainties shown correspond to the central 68% confidence interval and are consistent with the uncertainties obtained employing the standard error propagation method. The systematic uncertainties are estimated by varying the selection criteria of the particle candidates as described in Sect. 2.1. Panels (b) and (c) show the same comparison for the (p–p)–\(\Lambda \) and the p–(p–\(\Lambda \)) correlation functions. The number of events used for mixing to obtain (p–p)–p and (p–p)–\(\Lambda \) correlation functions is 30. The numerator of p–(p–\(\Lambda \)) correlation function requires p–\(\Lambda \) pairs in same event sample, which are less abundant than p–p pairs. Thus to obtain good statistical precision, the number of events used for mixing (to account for the third uncorrelated particle) must be increased to 100 in case of p–(p–\(\Lambda \)) triplets.
The results from the data-driven and the projector method are in good agreement between each other. The number of deviations \(n_{\sigma }\) in each bin are shown in the bottom panels of Fig. 1, where \(\sigma \) is the combined statistical and systematic uncertainty for both the experimental data and the projector. The agreement in the region \(Q_3 < 0.8\) GeV/c has been evaluated by performing a \(\chi ^2\) test. The \(\chi ^2\) is calculated combining the \(n_{\sigma }\) values of each bin. Finally, the p value from the \(\chi ^2\)-distribution is computed and the global \(n_{\sigma }\) values are extracted. The latter amount to 0.167, 0.0006 and 2.75 for (p–p)–p, (p–p)–\(\Lambda \) and p–(p–\(\Lambda \)), respectively. The data-driven method requires the usage of the third particle in the triplet from the mixed-event data sample and consequently the statistical uncertainty depends on the number of events used for mixing, while the projector method does not have this limitation. Thus, the latter significantly reduces the total uncertainty in the evaluation of the two-particle correlation effect on the three-particle correlation functions. For this reason, the projector method is used to calculate the three-particle cumulants for the p–p–p and p–p–\(\Lambda \) triplets.
The total two-particle contribution to the three-particle correlation function is obtained by substituting all terms on the right-hand side of Eq. 11 with the corresponding kinematic transformation, i.e.
where the indices refer to the label of the correlated pairs. In the case of p–p–p we have
and in the case of p–p–\(\Lambda \) we have
The resulting total lower-order contributions to the three-particle correlation functions (Eqs. 15 and 16) are shown in Fig. 2. The agreement between the data-driven approach and the projector method predictions translate into \(n_{\sigma }\) = 0.167 and \(n_{\sigma }\) = 0.0014 for the p–p–p and p–p–\(\Lambda \) lower-order contributions, respectively.
2.5 Decomposition of the three-particle cumulants
The experimental determination of the correlation function is mainly distorted by two distinct impurities in the candidate sample: misidentified particles and feed-down particles originating from weakly decaying particles. This introduces additional contributions to the correlation function of interest. These contributions are either assumed to be flat or, when the interaction is known, they are explicitly considered as discussed in Ref. [39]. The contributions to the correlation function stemming from decaying particles or impurities of the sample are weighted with the so-called \(\lambda \) parameters. By adopting this technique the residual correlations can be included in the final description of the experimental correlation function of two particles as
where the \(ij\ne 00\) denote all possible impurity and feed-down contributions and the \(ij=00\) is the correctly identified primary particle contribution. These \(\lambda \) parameters are obtained employing single particle properties such as the purity and feed-down probability. The underlying mathematical formalism is outlined in Ref. [39]. This mechanism has been extended to the three-particle case and the genuine three-particle cumulants can be obtained by subtracting the impurity and feed-down contributions from the measured cumulants. The full mathematical derivation is presented in Appendix C. The final expression of the genuine three-particle cumulants is
where X, Y and Z represent three generic particle species, the index 0 refers to correctly identified primary particle and the indexes i, j, k refer to misidentified or to secondary particles of a generic particle species. As shown in Appendix C, the specific weights \(\lambda \) depend on the purity and feed-down fraction of the single particles and are found to be equal to \(\lambda _{X_0 Y_0 Z_0}(\textrm{ppp})=\, 0.618\) and \(\lambda _{X_0 Y_0 Z_0}(\mathrm {pp\Lambda })=\,0.405\) for the p–p–p and p–p–\(\Lambda \) cumulants, respectively. Only 60% (40%) of the (p–p–\(\Lambda \)) triplets correspond to correctly identified primary particles.
In the following, the results for the p–p–p cumulants will be corrected according to the evaluated \(\lambda \) parameters assuming that all the three-particle contributions stemming from feed-down and impurities are flat in the momentum space. This assumption is supported by the observation that the measured p–p–\(\Lambda \) cumulants are consistent with zero within uncertainties (see Fig. 4 and the discussion in Sect. 3). The correction is not applied to the p–p–\(\Lambda \) cumulants because the shape of the feed-down contribution is not known and also because the statistical uncertainties are too large to provide any sensitivity to the three particle correlations.
3 Results
The measured three-particle correlation functions for p–p–p and p–p–\(\Lambda \) triplets are shown in Fig. 3 on the left and right panels, respectively. The number of events used for mixing for both cases is 30. The total number of same event triplets that are present at the range \(Q_{3}<\) 0.8 GeV/c are 17840 for p–p–p, 10980 for \(\overline{\textrm{p}}\text{-- }\overline{\textrm{p}}\text{-- }\overline{\textrm{p}}\), 9191 for p–p–\(\Lambda \) and 5886 for \(\overline{\textrm{p}}\text{-- }\overline{\textrm{p}}\text{-- }\overline{\Lambda } \). The green symbols represent the data points with their statistical and systematic uncertainties, while the grey bands correspond to the lower-order two-body interaction contributions obtained using the projector method already shown in Fig. 2. The non-femtoscopic contributions to the measured correlation functions, evaluated using Monte Carlo simulations, are found to be negligibly small (see Appendix A for a detailed discussion).
In the low \(Q_3\) region, the measured correlation functions deviate from the projected lower order contributions obtained using only two-particle correlations. The genuine three-body effects are then isolated by evaluating the cumulants
The lower-order contribution \(C^\mathrm {two\text {-}body}(Q_3)\) obtained with the projector method is used. The results for p–p–p and p–p–\(\Lambda \) triplets are shown in Fig. 4 on the left and right panels, respectively.
The p–p–p cumulant, already corrected for the feed-down contributions, is negative for \(0.16< Q_3 < 0.22\) GeV/c, while the large statistical uncertainty in the lowest \(Q_3\) interval prevents a conclusion on the sign for \(Q_3 < 0.16\) GeV/c. The agreement between the measured cumulant and the assumption that there are no genuine three-body effects is evaluated using the \(\chi ^2\) test in the region \(Q_3 < 0.4\) GeV/c, where the two-body interactions are prominent. There is no theoretical or experimental knowledge on the exact \(Q_3\) range where three-body effects become relevant, however they are expected to contribute at lower or same \(Q_3\) values as the two-body interactions. For this reason, the region of two-body forces was chosen. The obtained p-value corresponds to 6.7 standard deviations. If the cumulant is obtained using data-driven method to estimate lower order contributions, it results in 6.0 standard deviations. This result hints to the presence of an effect beyond the two-body interactions in the p–p–p system that could be either due to Pauli blocking (Fermi-Dirac quantum statistics) at the three particle level [64] or to the contribution of a three-body nuclear repulsive interactions. Long-range Coulomb interactions may also lead to significant contributions [65]. More quantitative conclusions on the interpretation of the non-zero cumulant require more sophisticated calculations for the three-body system. The present analysis demonstrates the experimental accessibility of the three-baryon cumulant in the data sample of pp collisions at \(\sqrt{s}=13\) TeV recorded by ALICE. In addition to the p–p–p system, the mixed-charge \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\) case has also been studied (see Appendix B for details). In the case of \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\), the two-body \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\) interaction contains both elastic and inelastic components, previously measured in an independent analysis [46]; the p–p interaction is well known from scattering data [66, 67] and verified by correlation measurements [39, 57, 68]; and the genuine three-body strong interaction for the \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\) triplet should be negligible. The extracted cumulant for \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\) is shown in the left panel of Fig. 4 by the red open circles. Since the number of the mixed-charge triplets is a factor four higher than the one of the same-charge triplets, it is possible to extend the measurement of the three-particle correlation to lower \(Q_3\) values. The correction for the feed-down contributions has been applied and the statistical and systematic uncertainties are shown. The \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\) cumulant evaluated using the projector method (data-driven approach) agrees with the assumption of only two-body correlations present in the system within 2.1 (2.2) standard deviations in the region \(Q_3 < 0.4\) GeV/c and within 0.9 (0.9) standard deviations in the region \(Q_3 < 0.2\) GeV/c, suggesting that genuine three-body effects are not statistically significant. This result as well demonstrates that the measured p–p–p cumulant deviation from zero is not due to detector effects.
In the case of p–p–\(\Lambda \), a positive cumulant is measured at \(Q_3 < 0.16\) GeV/c. The p-value obtained from the \(\chi ^2\) test in the region \(Q_3 < 0.4\) GeV/c corresponds to a deviation of 0.8 \(\sigma \) from the assumption that no genuine three-body correlations are present. Using data-driven method to estimate lower order contributions results in 0.8 \(\sigma \) as well. A similar value is found by repeating the significance test with the Fisher method, meaning that the measured cumulant is compatible with the assumption of no genuine three-body effects within the uncertainties. The current measurement does not allow to draw any firm conclusion yet on the three-body interaction in the p–p–\(\Lambda \) system, but since in this case only two of the particles are identical and charged, a non zero cumulant can be directly linked to the presence of a strong three-body interaction. It is estimated that employing a three-baryon event filtering during the upcoming Run 3 data taking should increase the number of triplets by a factor up to 500 for the target integrated luminosity of 200 nb\(^{-1}\) at \(\sqrt{s}=13.6\) TeV [69]. This opens up the possibility of measuring precisely the three-body correlations for both the p–p–p and p–p–\(\Lambda \) systems.
4 Conclusions
In this article, the first femtoscopic study of the p–p–p and p–p–\(\Lambda \) systems measured in high-multiplicity pp collisions at \(\sqrt{s} = 13\) TeV with the ALICE detector has been presented. In the chosen colliding system, hadrons are emitted at average relative distances of about 1 fm providing a unique environment to test three-body interactions at scales shorter than inter-particle ones in nuclei. The data collected during the LHC Run 2 enabled the measurement of the p–p–p and p–p–\(\Lambda \) correlation functions in the low \(Q_3\) region down to 0.1 GeV/c, giving access to the region where the effects of the hadronic two- and three-body interactions are more pronounced. The genuine three-particle correlations have been isolated using the Kubo’s cumulant expansion method. The lower-order two-body contributions have been estimated employing both a data-driven event mixing technique and a newly developed projector method. The two approaches have been compared and found to be in good agreement between each other, providing the first validation of the projector method using the data. The extracted p–p–p and p–p–\(\Lambda \) cumulants deviate from zero in the low \(Q_3\) region. In the case of p–p–p, a negative three-particle cumulant is measured. The p-value extracted from the \(\chi ^2\) test corresponds to a deviation of 6.7\(\sigma \) from the assumption of only two-body correlations present in the system for \(Q_3 < 0.4\) GeV/c. The obtained result provides an experimental hint to the presence of an effect beyond pairwise interactions in the p–p–p system. The observed deviation could be due to genuine three-body effects arising from: Pauli blocking, short-range strong interactions, or long-range Coulomb interactions. Refined three-body system calculations are required to give a solid interpretation of the measurement. The mixed-charge \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\) cumulant has also been measured as a benchmark and the result is consistent with the assumption that only two-body correlations are present in the system showing that the effect observed for the p–p–p system is a genuine one. In the case of p–p–\(\Lambda \), a positive cumulant is observed at low \(Q_3\). The deviation from zero at \(Q_3 < 0.4\) GeV/c is 0.8 \(\sigma \), which suggests no significant deviation from the assumption that only two-body correlations are present in the system within the current uncertainties. For this system, where one particle is uncharged and only two particles are identical, genuine three-particle correlations can be directly linked to the three-body strong interaction. The upcoming LHC Run 3 data taking will provide significantly larger samples of measured triplets, allowing more quantitative conclusions to be drawn for many-body dynamics.
The analysis presented in this article represents a first important step towards the direct measurement of the three-body interaction among baryons, demonstrating that genuine three-particle effects can be studied using three-particle correlation functions as experimental observables.
Data Availability Statement
This manuscript has associated data in a data repository. [Authors’ comment: Manuscript has associated data in a HEPData repository at https://www.hepdata.net/record/ins2092560.]
Notes
- $$\begin{aligned}\alpha&= \frac{4\ m_k^2}{(m_i+m_k)^2}+\frac{4\ m_k^2}{(m_j+m_k)^2}+ 4 \ ;\\ \beta&= \frac{4\ m_k (m_i+m_j+m_k)}{m_i+m_j} \left[ \frac{m_j}{(m_j+m_k)^2}-\frac{m_i}{(m_i+m_k)^2}\right] \ ;\\ \gamma&= \frac{4\ (m_i+m_j+m_k)^2}{(m_i + m_j)^2} \left[ \frac{ m_i^2}{(m_i + m_k )^2}+\frac{m_j^2}{( m_j + m_k)^2}\right] \ .\end{aligned}$$
References
P. Navrátil, S. Quaglioni, G. Hupin, C. Romero-Redondo, A. Calci, Unified ab initio approaches to nuclear structure and reactions. Phys. Scr. 91(5), 053002 (2016)
K. Hebeler, J. Holt, J. Menéndez, A. Schwenk, Nuclear forces and their impact on neutron-rich nuclei and neutron-rich matter. Annu. Rev. Nucl. Part. Sci. 65(1), 457–484 (2015). https://doi.org/10.1146/annurev-nucl-102313-025446
M. Duer et al., Observation of a correlated free four-neutron system. Nature 606(7915), 678–682 (2022). https://doi.org/10.1038/s41586-022-04827-6
E. Epelbaum, H.-W. Hammer, U.-G. Meißner, Modern theory of nuclear forces. Rev. Mod. Phys. 81, 1773–1825 (2009). https://doi.org/10.1103/RevModPhys.81.1773
B.R. Barrett, P. Navrátil, J.P. Vary, Ab initio no core shell model. Prog. Part. Nucl. Phys. 69, 131–181 (2013)
G. Hagen, M. Hjorth-Jensen, G. Jansen, R. Machleidt, T. Papenbrock, Evolution of shell structure in neutron-rich calcium isotopes. Phys. Rev. Lett. 109(3), 032502 (2012)
G. Hagen, T. Papenbrock, M. Hjorth-Jensen, D.J. Dean, Coupled-cluster computations of atomic nuclei. Rep. Prog. Phys. 77(9), 096302 (2014)
V. Soma, C. Barbieri, T. Duguet, Ab initio Gorkov-Green’s function calculations of open-shell nuclei. Phys. Rev. C 87(1), 011303 (2013)
R. Roth, J. Langhammer, A. Calci, S. Binder, P. Navrátil, Similarity-transformed chiral N N + 3 N interactions for the ab initio description of \(^{12}\)C and \(^{16}\)O. Phys. Rev. Lett. 107(7), 072501 (2011)
S.R. Stroberg, H. Hergert, S.K. Bogner, J.D. Holt, Nonempirical interactions for the nuclear shell model: an update. Annu. Rev. Nucl. Part. Sci. 69, 307–362 (2019)
J. Carlson, S. Gandolfi, F. Pederiva, S.C. Pieper, R. Schiavilla, K. Schmidt, R.B. Wiringa, Quantum Monte Carlo methods for nuclear physics. Rev. Mod. Phys. 87(3), 1067 (2015)
H.-W. Hammer, A. Nogga, A. Schwenk, Colloquium: three-body forces: from cold atoms to nuclei. Rev. Mod. Phys. 85, 197–217 (2013). https://doi.org/10.1103/RevModPhys.85.197
A. Kievsky, S. Rosati, M. Viviani, L.E. Marcucci, L. Girlanda, A high-precision variational approach to three- and four-nucleon bound and zero-energy scattering states. J. Phys. G Nucl. Part. Phys. 35(6), 063101 (2008). https://doi.org/10.1088/0954-3899/35/6/063101
STAR Collaboration, L. Adamczyk et al., Measurement of the \(^3_{\Lambda }\)H lifetime in Au+Au collisions at the BNL relativistic heavy ion collider. Phys. Rev. C 97(5), 054909 (2018).https://doi.org/10.1103/PhysRevC.97.054909. arXiv:1710.00436 [nucl-ex]
ALICE Collaboration, S. Acharya et al., \(^3_\Lambda H{ and ^3_{\bar{\Lambda }}\overline{H}}\) lifetime measurement in Pb-Pb collisions at \(\sqrt{s_{\rm NN}} = \) 5.02 TeV via two-body decay. Phys. Lett. B 797 134905 (2019). https://doi.org/10.1016/j.physletb.2019.134905. arXiv:1907.06906 [nucl-ex]
STAR Collaboration, J. Adam et al., Measurement of the mass difference and the binding energy of the hypertriton and antihypertriton. Nat. Phys. 16(4), 409–412 (2020). https://doi.org/10.1038/s41567-020-0799-7. arXiv:1904.10520 [hep-ex]
ALICE Collaboration, S. Acharya et al., Hypertriton production in p-Pb collisions at \(\sqrt{s_{\rm NN}}\)=5.02 TeV. Phys. Rev. Lett. 128(25), 252003 (2022). https://doi.org/10.1103/PhysRevLett.128.252003. arXiv:2107.10627 [nucl-ex]
H. Nemura, Y. Suzuki, Y. Fujiwara, C. Nakamoto, Study of light \(\Lambda \)- and \(\Lambda \Lambda \)-hypernuclei with the stochastic variational method and effective \(\Lambda \)N potentials. Prog. Theor. Phys. 103(5), 929–958 (2000). https://doi.org/10.1143/PTP.103.929
F. Hildenbrand, H.-W. Hammer, Three-body hypernuclei in pionless effective field theory. Phys. Rev. C 100, 034002 (2019). https://doi.org/10.1103/PhysRevC.100.034002
A. Feliciello, T. Nagae, Experimental review of hypernuclear physics: recent achievements and future perspectives. Rep. Prog. Phys. 78(9), 096301 (2015). https://doi.org/10.1088/0034-4885/78/9/096301
T. Hasegawa et al., Spectroscopic study of \({}_{{{\Lambda }}}^{10}{{\rm B}}\), \({}_{{{\Lambda }}}^{12}{{\rm C}}\), \({}_{{{\Lambda }}}^{28}{{\rm Si}}\), \({}_{{{\Lambda }}}^{89}{{\rm Y}}\), \({}_{{{\Lambda }}}^{139}{{\rm La}}\), and \({}_{{{\Lambda }}}^{208}{{\rm Pb}}\) by the (\({{{\pi }}}^{+}\),\({{k}}^{+}\)) reaction. Phys. Rev. C 53, 1210–1220 (1996). https://doi.org/10.1103/PhysRevC.53.1210
FINUDA Collaboration, M. Agnello et al., Hypernuclear spectroscopy with K\(^-\) at rest on \(^7\)Li, \(^9\)Be, \(^{13}\)C and \(^{16}\)O. Phys. Lett. B 698(3), 219–225 (2011)
J-PARC E13 Collaboration, S.B. Yang et al., First determination of the level structure of an \(sd\)-shell hypernucleus, \(_{\rm \Lambda \mathit{}^{19}\rm F}\). Phys. Rev. Lett. 120, 132505 (2018). https://doi.org/10.1103/PhysRevLett.120.132505
Jefferson Lab Hall A Collaboration, F. Garibaldi et al., High-resolution hypernuclear spectroscopy at jefferson lab, hall a. Phys. Rev. C 99, 054309 (2019). https://doi.org/10.1103/PhysRevC.99.054309
HypHI Collaboration, C. Rappold et al., Search for evidence of \({}_{\Lambda }^{3}n\) by observing \(d+{\pi }^{-}\) and \(t+{\pi }^{-}\) final states in the reaction of \({}^{6}\)Li+\({}^{12}\)C at \(2A\) GeV. Phys. Rev. C 88, 041001 (2013). https://doi.org/10.1103/PhysRevC.88.041001
T.R. Saito et al., New directions in hypernuclear physics. Nat. Rev. Phys. 3(12), 803–813 (2021). https://doi.org/10.1038/s42254-021-00371-w
M. Oertel, M. Hempel, T. Klähn, S. Typel, Equations of state for supernovae and compact stars. Rev. Mod. Phys. 89(1), 015007 (2017). https://doi.org/10.1103/RevModPhys.89.015007. arXiv:1610.03361 [astro-ph.HE]
L. Tolos, L. Fabbietti, Strangeness in nuclei and neutron stars. Prog. Part. Nucl. Phys. 112, 103770 (2020). https://doi.org/10.1016/j.ppnp.2020.103770. arXiv:2002.09223 [nucl-ex]
P. Demorest, T. Pennucci, S. Ransom, M. Roberts, J. Hessels, Shapiro delay measurement of a two solar mass neutron star. Nature 467, 1081–1083 (2010). https://doi.org/10.1038/nature09466. arXiv:1010.5788 [astro-ph.HE]
J. Antoniadis et al., A massive pulsar in a compact relativistic binary. Science340(6131) (2013). https://doi.org/10.1126/science.1233232. https://science.sciencemag.org/content/340/6131/1233232
D. Lonardoni, A. Lovato, S. Gandolfi, F. Pederiva, Hyperon puzzle: hints from quantum Monte Carlo calculations. Phys. Rev. Lett. 114(9), 092301 (2015). https://doi.org/10.1103/PhysRevLett.114.092301. arXiv:1407.4448 [nucl-th]
U.A. Wiedemann, U. Heinz, Particle interferometry for relativistic heavy-ion collisions. Phys. Rep. 319(4–5), 145–230 (1999)
U. Heinz, B.V. Jacak, Two-particle correlations in relativistic heavy-ion collisions. Annu. Rev. Nucl. Part. Sci. 49(1), 529–579 (1999)
R. Lednický, Correlation femtoscopy of multiparticle processes. Phys. At. Nucl. 67(1), 72–82 (2004)
M.A. Lisa, S. Pratt, R. Soltz, U. Wiedemann, Femtoscopy in relativistic heavy ion collisions. Ann. Rev. Nucl. Part. Sci. 55, 357–402 (2005). https://doi.org/10.1146/annurev.nucl.55.090704.151533. arXiv:nucl-ex/0505014
STAR Collaboration, L. Adamczyk, J. Adkins, G. Agakishiev, M. Aggarwal, Z. Ahammed, I. Alekseev, J. Alford, C. Anson, A. Aparin, D. Arkhipkin, et al., \(\Lambda \)\(\Lambda \) correlation function in Au+ Au collisions at \(\sqrt{s_{{\rm NN}}}\) = 200 GeV. Phys. Rev. Lett.114(2), 022301 (2015)
STAR Collaboration, L. Adamczyk et al., Measurement of interaction between antiprotons. Nature 527, 345–348 (2015). https://doi.org/10.1038/nature15724. arXiv:1507.07158 [nucl-ex]
STAR Collaboration, J. Adam et al., The proton-\(\Omega \) correlation function in Au+Au collisions at \(\sqrt{s_{\rm NN}}\)=200 GeV. Phys. Lett. B 790, 490–497 (2019). https://doi.org/10.1016/j.physletb.2019.01.055. arXiv:1808.02511 [hep-ex]
ALICE Collaboration, S. Acharya et al., p-p, p-\(\Lambda \) and \(\Lambda \)-\(\Lambda \) correlations studied via femtoscopy in pp reactions at \(\sqrt{s}=\)7 TeV. Phys. Rev. C 99, 024001 (2019). https://doi.org/10.1103/PhysRevC.99.024001. arXiv:1805.12455 [nucl-ex]
ALICE Collaboration, S. Acharya et al., Scattering studies with low-energy kaon-proton femtoscopy in proton-proton collisions at the LHC. Phys. Rev. Lett. 124(9), 092301 (2020). https://doi.org/10.1103/PhysRevLett.124.092301. arXiv:1905.13470 [nucl-ex]
ALICE Collaboration, S. Acharya et al., Investigation of the p–\(\Sigma \)0 interaction via femtoscopy in pp collisions. Phys. Lett. B 805, 135419 (2020). arXiv:1910.14407
ALICE Collaboration, S. Acharya et al., Study of the \(\Lambda \)-\(\Lambda \) interaction with femtoscopy correlations in pp and p-Pb collisions at the LHC. Phys. Lett. B 797, 134822 (2019). https://doi.org/10.1016/j.physletb.2019.134822. arXiv:1905.07209 [nucl-ex]
ALICE Collaboration, S. Acharya et al., First observation of an attractive interaction between a proton and a cascade baryon. Phys. Rev. Lett. 123, 112002 (2019). https://doi.org/10.1103/PhysRevLett.123.112002. arXiv:1904.12198 [nucl-ex]
ALICE Collaboration, S. Acharya et al., Unveiling the strong interaction among hadrons at the LHC. Nature 588(7837), 232–238 (2020). https://doi.org/10.1038/s41586-020-3001-6. arXiv:2005.11495
ALICE Collaboration, S. Acharya et al., Experimental evidence for an attractive p-\(\phi \) interaction. Phys. Rev. Lett. 127(17), 172301 (2021). https://doi.org/10.1103/PhysRevLett.127.172301. arXiv:2105.05578 [nucl-ex]
ALICE Collaboration, S. Acharya et al., Investigating the role of strangeness in baryon–antibaryon annihilation at the LHC. Phys. Lett. B 829, 137060 (2022). https://doi.org/10.1016/j.physletb.2022.137060. arXiv:2105.05190 [nucl-ex]
L. Fabbietti, V.M. Sarti, O.V. Doce, Study of the strong interaction among hadrons with correlations at the LHC. Ann. Rev. Nucl. Part. Sci. 71, 377–402 (2021). https://doi.org/10.1146/annurev-nucl-102419-034438. arXiv:2012.09806 [nucl-ex]
ALICE Collaboration, B.B. Abelev et al., Two- and three-pion quantum statistics correlations in Pb-Pb collisions at \(\sqrt{{s}_{\rm NN}}\) = 2.76 TeV at the CERN Large Hadron Collider. Phys. Rev. C 89(2), 024911 (2014). https://doi.org/10.1103/PhysRevC.89.024911. arXiv:1310.7808 [nucl-ex]
ALICE Collaboration, J. Adam et al., Multipion Bose-Einstein correlations in pp, p-Pb, and Pb-Pb collisions at energies available at the CERN Large Hadron Collider. Phys. Rev. C 93, 054908 (2016). https://doi.org/10.1103/PhysRevC.93.054908. arXiv:1512.08902
R. Kubo, Generalized cumulant expansion method. J. Phys. Soc. Jpn. 17(7), 1100–1120 (1962). https://doi.org/10.1143/JPSJ.17.1100
R. Del Grande, L. Šerkšnytė, L. Fabbietti, V.M. Sarti, D. Mihaylov, A method to remove lower order contributions in multi-particle femtoscopic correlation functions. Eur. Phys. J. C 82(3), 244 (2022). https://doi.org/10.1140/epjc/s10052-022-10209-z. arXiv:2107.10227 [nucl-th]
ALICE Collaboration, K. Aamodt et al., The ALICE experiment at the CERN LHC. J. Instr. 3, S08002 (2008). https://doi.org/10.1088/1748-0221/3/08/S08002
ALICE Collaboration, B. Abelev et al., Performance of the ALICE experiment at the CERN LHC. Int. J. Mod. Phys. A 29, 1430044 (2014). https://doi.org/10.1142/S0217751X14300440
ALICE Collaboration, E. Abbas et al., Performance of the ALICE VZERO system. JINST 8, P10016 (2013). https://doi.org/10.1088/1748-0221/8/10/P10016. arXiv:1306.3130 [nucl-ex]
J. Alme et al., The ALICE TPC, a large 3-dimensional tracking device with fast readout for ultra-high multiplicity events. Nucl. Instrum. Methods A 622, 316–367 (2010). https://doi.org/10.1016/j.nima.2010.04.042
A. Akindinov et al., Performance of the ALICE Time-of-Flight detector at the LHC. Eur. Phys. J. Plus 128, 44 (2013). https://doi.org/10.1140/epjp/i2013-13044-x
ALICE Collaboration, S. Acharya et al., Search for a common baryon source in high-multiplicity pp collisions at the LHC. Phys. Lett. B 811, 135849 (2020). https://doi.org/10.1016/j.physletb.2020.135849. arXiv:2004.08018 [nucl-ex]
Particle Data Group Collaboration, M. Tanabashi et al., Review of particle physics. Phys. Rev. D 98, 030001 (2018). https://doi.org/10.1103/PhysRevD.98.030001
ARGUS Collaboration, H. Albrecht et al., Observation of Octet and Decuplet Hyperons in \(e^+ e^-\) annihilation at 10-GeV center-of-mass energy. Phys. Lett. B 183, 419–424 (1987). https://doi.org/10.1016/0370-2693(87)90988-9
R. Lednický, V.L. Lyuboshits, Final state interaction effect on pairing correlations between particles with small relative momenta. Sov. J. Nucl. Phys. 35, 770 (1982)
S.E. Koonin, Proton pictures of high-energy nuclear collisions. Phys. Lett. B 70(1), 43–47 (1977)
S. Pratt, T. Csörgő, J. Zimányi, Detailed predictions for two-pion correlations in ultrarelativistic heavy-ion collisions. Phys. Rev. C 42(6), 2646 (1990)
ALICE Collaboration, S. Acharya et al., Exploring the N\(\Lambda \)-N\(\Sigma \) coupled system with high precision correlation techniques at the LHC. arXiv:2104.04427 [nucl-ex]
P. Niemann, H.W. Hammer, Pauli blocking effects and Cooper triples in three-component Fermi gases. Phys. Rev. A 86, 013628 (2012). https://doi.org/10.1103/PhysRevA.86.013628. arXiv:1203.1824 [cond-mat.quant-gas]
E. Alt, A. Mukhamedzhanov, Asymptotic solution of the Schrödinger equation for three charged particles. Phys. Rev. A 47(3), 2004 (1993)
R.A. Arndt, W.J. Briscoe, I.I. Strakovsky, R.L. Workman, Updated analysis of NN elastic scattering to 3-GeV. Phys. Rev. C 76, 025209 (2007). https://doi.org/10.1103/PhysRevC.76.025209. arXiv:0706.2195 [nucl-th]
R. Navarro Pérez, J.E. Amaro, E. Ruiz Arriola, Partial wave analysis of nucleon-nucleon scattering below pion production threshold. Phys. Rev. C 88, 024002 (2013). https://doi.org/10.1103/PhysRevC.88.024002. arXiv:1304.0895 [nucl-th]. [Erratum: Phys. Rev. C 88, 069902 (2013)]
ALICE Collaboration, S. Acharya et al., Investigation of the p-\(\Sigma ^{0}\) interaction via femtoscopy in pp collisions. Phys. Lett. B 805, 135419 (2020). https://doi.org/10.1016/j.physletb.2020.135419. arXiv:1910.14407 [nucl-ex]
ALICE Collaboration, Future high-energy pp programme with ALICE. ALICE-PUBLIC-2020-005, CERN-LHCC-2020-018, LHCC-G-179 (2020)
T. Sjöstrand, S. Ask, J.R. Christiansen, R. Corke, N. Desai, P. Ilten, S. Mrenna, S. Prestel, C.O. Rasmussen, P.Z. Skands, An introduction to PYTHIA 8.2. Comput. Phys. Commun. 191, 159–177 (2015). https://doi.org/10.1016/j.cpc.2015.01.024. arXiv:1410.3012 [hep-ph]
ALICE Collaboration, J. Adam et al., Insight into particle production mechanisms via angular correlations of identified particles in pp collisions at \(\sqrt{s}\) = 7 TeV. Eur. Phys. J. C77(8), 569 (2017). https://doi.org/10.1140/epjc/s10052-017-5129-6
Acknowledgements
The ALICE Collaboration is grateful to Prof. Alejandro Kievsky for the fruitful discussions on the theoretical aspects of three-body systems. The ALICE Collaboration would like to thank all its engineers and technicians for their invaluable contributions to the construction of the experiment and the CERN accelerator teams for the outstanding performance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources and support provided by all Grid centres and the Worldwide LHC Computing Grid (WLCG) collaboration. The ALICE Collaboration acknowledges the following funding agencies for their support in building and running the ALICE detector: A. I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation (ANSL), State Committee of Science and World Federation of Scientists (WFS), Armenia; Austrian Academy of Sciences, Austrian Science Fund (FWF): [M 2467-N36] and Nationalstiftung für Forschung, Technologie und Entwicklung, Austria; Ministry of Communications and High Technologies, National Nuclear Research Center, Azerbaijan; Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Financiadora de Estudos e Projetos (Finep), Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) and Universidade Federal do Rio Grande do Sul (UFRGS), Brazil; Bulgarian Ministry of Education and Science, within the National Roadmap for Research Infrastructures 2020¿2027 (object CERN), Bulgaria; Ministry of Education of China (MOEC) , Ministry of Science & Technology of China (MSTC) and National Natural Science Foundation of China (NSFC), China; Ministry of Science and Education and Croatian Science Foundation, Croatia; Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear (CEADEN), Cubaenergía, Cuba; Ministry of Education, Youth and Sports of the Czech Republic, Czech Republic; The Danish Council for Independent Research | Natural Sciences, the VILLUM FONDEN and Danish National Research Foundation (DNRF), Denmark; Helsinki Institute of Physics (HIP), Finland; Commissariat à l’Energie Atomique (CEA) and Institut National de Physique Nucléaire et de Physique des Particules (IN2P3) and Centre National de la Recherche Scientifique (CNRS), France; Bundesministerium für Bildung und Forschung (BMBF) and GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany; General Secretariat for Research and Technology, Ministry of Education, Research and Religions, Greece; National Research, Development and Innovation Office, Hungary; Department of Atomic Energy Government of India (DAE), Department of Science and Technology, Government of India (DST), University Grants Commission, Government of India (UGC) and Council of Scientific and Industrial Research (CSIR), India; National Research and Innovation Agency - BRIN, Indonesia; Istituto Nazionale di Fisica Nucleare (INFN), Italy; Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT) and Japan Society for the Promotion of Science (JSPS) KAKENHI, Japan; Consejo Nacional de Ciencia (CONACYT) y Tecnología, through Fondo de Cooperación Internacional en Ciencia y Tecnología (FONCICYT) and Dirección General de Asuntos del Personal Academico (DGAPA), Mexico; Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO), Netherlands; The Research Council of Norway, Norway; Commission on Science and Technology for Sustainable Development in the South (COMSATS), Pakistan; Pontificia Universidad Católica del Perú, Peru; Ministry of Education and Science, National Science Centre and WUT ID-UB, Poland; Korea Institute of Science and Technology Information and National Research Foundation of Korea (NRF), Republic of Korea; Ministry of Education and Scientific Research, Institute of Atomic Physics, Ministry of Research and Innovation and Institute of Atomic Physics and University Politehnica of Bucharest, Romania; Ministry of Education, Science, Research and Sport of the Slovak Republic, Slovakia; National Research Foundation of South Africa, South Africa; Swedish Research Council (VR) and Knut & Alice Wallenberg Foundation (KAW), Sweden; European Organization for Nuclear Research, Switzerland; Suranaree University of Technology (SUT), National Science and Technology Development Agency (NSTDA), Thailand Science Research and Innovation (TSRI) and National Science, Research and Innovation Fund (NSRF), Thailand; Turkish Energy, Nuclear and Mineral Research Agency (TENMAK), Turkey; National Academy of Sciences of Ukraine, Ukraine; Science and Technology Facilities Council (STFC), United Kingdom; National Science Foundation of the United States of America (NSF) and United States Department of Energy, Office of Nuclear Physics (DOE NP), United States of America. In addition, individual groups or members have received support from: Marie Skłodowska Curie, European Research Council, Strong 2020-Horizon 2020 (grant nos. 950692, 824093, 896850), European Union; Academy of Finland (Center of Excellence in Quark Matter) (grant nos. 346327, 346328), Finland; Programa de Apoyos para la Superación del Personal Académico, UNAM, Mexico.
Funding
Open access funding provided by CERN (European Organization for Nuclear Research)
Author information
Authors and Affiliations
Consortia
Additional information
Communicated by C. Munoz Camacho.
Appendices
Appendix A: Monte Carlo studies
One of the benchmarks of the presented analysis consists in verifying that the three-particle correlations obtained from PYTHIA 8 [70] (Monash 2013 Tune) simulations do not show any significant deviation from unity. Indeed, no FSIs – either two- or three-particle – are included in the simulated and reconstructed Monte Carlo data using the PYTHIA 8 event generator for pp collisions at \(\sqrt{s} = \, 13\) TeV, thus Monte Carlo can neither be used to estimate the lower-order contributions (two-body correlations) nor genuine three-body effects. Figure 5 shows the comparison between the measured and simulated correlation functions as a function of \(Q_3\), where the simulation includes a dedicated high-multiplicity selection to mimic the V0 high-multiplicity trigger in the real data. The green symbols represent the experimental data while the black symbols refer to the simulation. The simulated correlation functions are consistent with unity for the entire \(Q_3< 0.8 \) GeV/c range, showing that there are no effects caused by the track reconstruction in the detector, as well as no sign of mini-jets contribution (see more details in Ref. [46]) and that the energy and momentum conservation effects are eventually present at larger values of \(Q_3\) that are not relevant for the studies carried out in this work. Also the simulations of the lower-order contributions display the same behaviour.
Appendix B: Mixed-charge correlation studies
An additional benchmark for the p–p–p cumulant result and the measured deviation from zero has been considered by studying the \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\) triplets. Identical event and track selection criteria and systematic variations as those employed for the \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\) analysis have been used. The rejection in the \(\Delta \eta \)–\(\Delta \varphi \) plane (see Sect. 2.2 for the details) has been applied only for same-charge pairs in the triplet. The obtained correlation functions for \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\), \({(\textrm{p}\text{-- }\textrm{p})\text{-- }\overline{\textrm{p}}}\) and \((\textrm{p}\text{-- }{\overline{\textrm{p}}})\text{-- }\textrm{p}\) triplets with the corresponding statistical and systematic uncertainties are shown in panels (a), (b) and (c) of Fig. 6, respectively. The grey bands are obtained using the projector method following the procedure described in Sect. 2.4. The \({(\textrm{p}\text{-- }\textrm{p})\text{-- }\overline{\textrm{p}}}\) correlation functions obtained with the data-driven approach and the projector method (panel b) in Fig. 6 are in agreement with the (p–p)–p results (panel a) in Fig. 1. This is consistent with the expectation from Eqs. 12 and 13, since such correlation functions depend only on the correlation of the two particles in the same event and on the mass of the uncorrelated particle, which are identical in the \({(\textrm{p}\text{-- }\textrm{p})\text{-- }\overline{\textrm{p}}}\) and (p–p)–p systems. The \((\textrm{p}\text{-- }{\overline{\textrm{p}}})\text{-- }\textrm{p}\) correlation function reflects the interplay of FSIs and non-femtoscopic correlations measured in the study of \(\textrm{p}\text{-- }\overline{\textrm{p}}\) pairs [46]. The grey band in panel (c) of Fig. 6 is obtained by using as input of Eq. 12 the correlation function \(C(k^*)\) of \(\textrm{p}\text{-- }\overline{\textrm{p}}\) pairs emitted in \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\) triplets with \(Q_3 < 1\) GeV/c. Such selection is performed in order to use \(\textrm{p}\text{-- }{\overline{\textrm{p}}}\) pairs produced in similar shape events as those where low \(Q_3\) triplets are found. The requirement is necessary for \(\textrm{p}\text{-- }{\overline{\textrm{p}}}\) correlations due to the mini-jet contribution [46] which is instead not present in p–\(\Lambda \) and p–p (see e.g. [71]). The systematics induced by this additional selection are considered in the uncertainties.
The \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\) cumulant is obtained using the projector method and it is shown in Fig. 4 by red circular symbols. The result is compatible with zero within the statistical and systematic uncertainties, demonstrating that strong-interaction as well as Coulomb effects on three-particle level are not statistically significant in \({\textrm{p}\text{-- }\textrm{p}\text{-- }\overline{\textrm{p}}}\).
Appenndix C: Feed-down contributions
The measured femtoscopic correlations do not originate only from correctly identified primary particles but they include as well contributions from misidentified particles and feed-down particles originating from weakly decaying hadrons. In case of two-body femtoscopy, the decomposition method explained in Sect. 2.5 and the Eq. 17 can be used to account for such impurities and feed-down effects.
This method can be extended to the three-particle case. The total data sample X contains the particles which stem from feed-down \(X_F\), misidentified particles \(X_M\) and the correctly identified primary particles \(X_0\). Both feed-down and misidentified particles can originate from different channels and the contributions can be expressed as
where \(N_F\) and \(N_M\) are the number of feed-down and misidentification contributions. The purity is the fraction of correctly identified particles to the total number of particles in the data sample and can be defined as
The correctly identified particles can stem from the decays of particles and for this purpose the channel fraction \(f(X_i)\) is defined as
The fraction of specific channel in the whole data sample then can be written as
The correlation function for three particles can be written as
where N and M are the yields of XYZ triplet in same and mixed events, respectively. Using the identities introduced before one can write
The correlation function then can be rewritten as
where \(C_{i, j,k}(X YZ)\) is the correlation function of the i, j, k-th channel of origin of the particles X, Y, Z and the \( \lambda _{i, j,k}(X Y Z)\) is the weight for such correlation. This parameter depends only on the mixed event sample and can be related to previously introduced single particle quantities, channel fraction and purity, as follows
To study the lower-order contributions in the measured three-particle correlation function, one needs to define the decomposition for the measurement of two correlated particles. In such case, the origin of the third particle in the numerator does not matter as it is uncorrelated. Equation C.9 becomes
Here \(C_{i, j}(X Y,Z)\) denotes the correlation function where the two correlated particles X and Y are from origin i and j respectively and Z is from any origin. Here \(\lambda _{i, j}(X Y, Z) \) is
This notation is valid only if one wants to study the \((X-Y)-Z\) correlation. To obtain the cumulant of the primary particles which were identified correctly, one has to subtract the lower-order correlations, such as \((X-Y)-Z\), from the three-particle correlation. For this purpose, the Eq. C.11 must be rewritten. As previously explained, in case of \((X-Y)-Z\) correlation, the origin of the third uncorrelated particle is not important, which means that \(C(X Y,Z_l)=C(X Y,Z_m)\), where \(Z_l\) and \(Z_m\) have different origin. As previously shown, the fraction of particles \(Z_l\) in the whole data sample is \(P(Z_l)\), which can be as well expressed with the \(\lambda \) parameter of one particle \(\lambda _l(Z)\). Using the property \(1 = \sum _{k}\lambda _k(Z)\) of the \(\lambda \) parameters in Eq. C.12, one can write
and Eq. C.11 can be rewritten as
In the following, the above formalism is used to obtain the correctly identified primary particle cumulant. Starting with the cumulant expression
the correctly identified primary particle correlations can be isolated from the rest as follows
Written in such a way, one can already see that the cumulant of the measured correlation function can be expressed as the sum of the correctly identified primary particle cumulant and the cumulant of the rest of the contributions as follows
The terms inside the brackets are almost a cumulant expression, except the \(+2\) term is missing, but one can add and subtract 2 to obtain
This is the final result – the cumulant calculated using the measured correlation functions consists of the three correctly identified primary particle cumulant and the cumulant which consist of the rest of possible contributions. In such case, the correctly identified particle cumulant is
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Acharya, S., Adamová, D., Adler, A. et al. Towards the understanding of the genuine three-body interaction for p–p–p and p–p–\(\Lambda \). Eur. Phys. J. A 59, 145 (2023). https://doi.org/10.1140/epja/s10050-023-00998-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epja/s10050-023-00998-6