Study of the p–p–K + and p–p–K − dynamics using the femtoscopy technique

The interactions of kaons (K) and antikaons (K) with few nucleons (N) were studied so far using kaonic atom data and measurements of kaon production and interaction yields in nuclei. Some details of the three-body KNN and KNN dynamics are still not well understood, mainly due to the overlap with multi-nucleon interactions in nuclei. An alternative method to probe the dynamics of three-body systems with kaons is to study the ﬁnal state interaction within triplet of particles emitted in pp collisions at the Large Hadron Collider, which are free from effects due to the presence of bound nucleons. This Letter reports the ﬁrst femtoscopic study of p–p–K + and p–p–K − correlations measured in high-multiplicity pp collisions at √ s = 13 TeV by the ALICE Collaboration. The analysis shows that the measured p–p–K + and p–p–K − correlation functions can be interpreted in terms of pairwise interactions in the triplets, indicating that the dynamics of such systems is dominated by the two-body interactions without signiﬁcant contributions from three-body effects or bound states.

An alternative method to explore the three-body dynamics of K + and K − with nucleons is to employ the femtoscopy technique at high-energy collider facilities.In small colliding systems, such as pp and p-Pb at the Large Hadron Collider (LHC), the inter-hadron distances at the time of the particle emission range from a few femtometers down to scales compatible with the nucleon size.This leads to an enhancement of the strength of the signal due to the short-range strong interaction in the measured correlation function [45].The femtoscopy method was proven to be able to test and constrain the hadron-hadron interaction for various two-particle systems [46][47][48][49][50][51][52][53][54][55][56][57], providing data with unprecedented precision on the hadronic interactions with strangeness at low relative momenta, down to the energy threshold of the produced pairs.Furthermore, the measured femtoscopic correlation functions are sensitive to the presence of bound states in the energy region below the threshold.Thus, they were also used to constrain the parameters of bound hadron-hadron systems [49,58].Recently, the ALICE Collaboration has extended the method to explore the correlation among three baryons, such as p-p-p and p-p-Λ, to study the dynamics of three-body systems [59].The analysis exploited Kubo's cumulant expansion method [60] and the projector method [61] to isolate the genuine three-particle correlations from the measured correlation functions.As a result, a first experimental hint of genuine three-body effects in the unbound p-p-p system was found.The study in this Letter applies the same analysis procedure adopted in [59] to the case of p-p-K + and p-p-K − particle triplets to explore possible genuine three-body effects in the correlation functions induced by the strong interaction and bound state formation.The main advantage of the femtoscopy method with respect to the previous experimental techniques is the possibility to investigate for the first time the K + and K − three-body dynamics with nucleons, free from additional effects induced by the presence of the surrounding nucleons in nuclei.
The analysed data sample of pp collisions at √ s = 13 TeV was recorded by ALICE [62][63][64] during the LHC Run 2 (2015-2018) data-taking period.In the following, the information from the V0 detector system [65], the Inner Tracking System (ITS) [66], the Time Projection Chamber (TPC) [67] and the Time-Of-Flight (TOF) detector [68] is used.The V0 detector is employed to trigger on high-multiplicity (HM) events.This trigger selects events within the 0.17% largest charged-track multiplicity of the INEL > 0 class, which is defined as inelastic collisions with at least one measured charged particle in the pseudorapidity range |η| < 1 [65].This condition results in an average of about 30 charged particles in the range |η| < 0.5 [51], hence increasing the probability of finding triplets of the desired particle species with respect to the minimum-bias sample.The primary vertex position is reconstructed with the combined information of the ITS and the TPC, and, independently, with track segments in the two innermost layers of the ITS.Only events with a reconstructed primary vertex position along the beam axis within 10 cm from the centre of the ALICE detector are selected.A total of 10 9 HM events are used in this analysis.The p-p-K + and p-p-K − data samples are built by combining three charged-particle tracks (triplets) reconstructed with the TPC.Assuming the same interactions in the particle and antiparticle systems, triplets of particles and the corresponding antiparticles are combined: p-p- The agreement of the corresponding correlation functions confirmed the validity of this assumption.
Particles and antiparticles are identified using the same kinematic and topological selections.Protons and kaons are selected in the range |η| < 0.8 and in the transverse momentum intervals p T ∈ (0.5-4.05) GeV/c and p T ∈ (0.2-2.5) GeV/c, respectively.A minimum of 80 space points (hits) inside the TPC, out of the total 159, is required to guarantee track quality and good momentum resolution.Particle identification (PID) is conducted by requiring that the measured energy loss (dE/dx) in the TPC gas is compatible with the expected one from protons and kaons within three standard deviations (σ ).For high momentum particles, the dE/dx information is combined with the time-of-flight measurement provided by the TOF, using a 3σ selection on the expected value for a given particle hypothesis at a given momentum.The PID selection for protons is described in detail in [46].The selection of kaons is based on the procedure described in [54] with several changes as explained in the following.TPC reconstructed tracks are identified as kaons either by using TPC PID information for momenta lower than 0.85 GeV/c, or, if a signal in the TOF is matched to the TPC track, by using the combined PID information from TPC and TOF up to a momentum of 2.5 GeV/c.To improve the purity of the kaon selection, the TPC and TOF information are also used to reject candidates that are compatible with the pion or electron hypothesis.
To reject particles that are non-primary or come from pile-up collisions, the distance of closest approach (DCA) to the primary vertex 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 the proton and kaon candidates is estimated employing Monte Carlo simulations based on PYTHIA 8 [69] (Monash 2013 Tune) event generator with a dedicated high-multiplicity selection to mimic the V0 high-multiplicity trigger and the GEANT4 package [70,71].The purity averaged over the p T ranges of the identified protons and kaons is about 98% for both particle species.
The number of selected and analysed triplets amounts to 4530 for p-p-K + , 3161 for p-p-K − , 6200 for p-p-K − and 4937 for p-p-K + in the femtoscopic region Q 3 < 0.4 GeV/c.The kinematic variable Q 3 , which is used in three-body analyses [59,72], is a Lorentz-invariant scalar defined as , where q i j is the norm of the relative four-momentum of the pair i j in the triplet q The systematic uncertainties of the correlation functions are evaluated by performing simultaneous variations of the selection criteria for particles and antiparticles.For protons and antiprotons, the same variations for the track selection and PID criteria as in [59] are performed.Similar variations are used for the kaons and antikaons.In order to account for the correlations between the systematic uncertainties, the variations are randomly combined in sets in which at least one selection criterion is changed.Large statistical fluctuations introduced by the variations of the selection criteria are avoided by requiring that the yield of the triplets changes by less than 10% with respect to the standard selection in the femtoscopic region Q 3 < 0.4 GeV/c.The main contribution at low Q 3 (≈ 0.15 GeV/c) is given by the variation of the DCA selection and it is found to be smaller than 5%.
The final state interactions (FSIs) among the hadrons emitted in the collisions can be explored using correlation functions [73,74] in momentum space, defined as where p i is the momentum vector of the i-th particle and P 3 (p 1 , p 2 , p 3 ) and P 1 (p i ) are the probabilities of finding three particles and one particle with the corresponding momentum, respectively.If the particles are not correlated P 3 (p 1 , p 2 , p 3 ) = P 1 (p 1 )P 1 (p 2 )P 1 (p 3 ) and then the correlation function becomes identical to unity.The presence of FSIs induces a correlation signal imprinted in the correlation function which causes deviations from unity depending on the nature of the interactions, attractive or repulsive, as well as on the properties of the emitting source (see [45,75] for the details).In the case of three-particle systems, the correlation function can be evaluated in terms of Q 3 .Hence, three-particle correlation functions are experimentally obtained as where N is a normalisation parameter, N same (Q 3 ) is the Q 3 distribution of the particle triplets emitted in the same collision and N mixed (Q 3 ) is the distribution of uncorrelated triplets.The latter are obtained by taking the three particles needed to form a triplet from three different collisions, thus called mixed event sample.Only events with similar multiplicity and position of the primary vertex along the beam axis are mixed.The number of events used for the mixing is set to 30.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 the same charge tracks on the azimuthal-polar angle plane ∆η-∆ϕ is applied to both the same and mixed event samples.The default selection is ∆η 2 + ∆ϕ 2 ≥ (0.017) 2 for p-p tracks and ∆η 2 /(0.012) 2 + ∆ϕ 2 /(0.04) 2 ≥ 1 for p-K + tracks.A systematic variation of about +10% for the values of the minimum distance is applied in the analysis.The correlation function is normalised to unity in the range 1.0 < Q 3 < 1.2 GeV/c, where no signal from FSIs is expected.Variations of the normalisation range are performed and included in the systematic uncertainties.In the case of a particle triplet X-Y-Z, pairwise correlations with one uncorrelated particle (X-Y)-Z, X-(Y-Z) and (Z-X)-Y also occur in the system.Genuine three-particle correlations were isolated in [59] by applying the cumulant expansion method to the femtoscopic correlation functions.The three-particle cumulant c 3 (Q 3 ), which incorporates information about the three-body effects in the system, is computed from the measured correlation function C 3 (Q 3 ) as follows where pairwise correlations are evaluated following the cumulant decomposition as Each component C (i−j)−k (Q 3 ) of the lower-order contributions in Eq. ( 4) is computed using two methods [59,61]: i) a data-driven approach based on the event mixing technique, building triplets in which the correlated (i-j) pairs are emitted in the same collision and the uncorrelated particle -k from another collision; ii) the projector method which uses as input the measured two-particle correlation function C(k * ) for the correlated pair (i-j), evaluated in terms of the relative momentum k * in the pair rest frame (PRF), and calculates all the possible k * configurations in the phase space for each Q 3 value of the i-j-k triplet.The input correlation functions C(k * ) used in the projector method are obtained by selecting p-p, p-K + and p-K − pairs using the same dataset as for the three-body analysis.Unlike the p-p and p-K + pairs, the p-K − correlation is significantly affected by the contribution from jet-like events [76].To use particles emitted in collisions with the same event shape as in the case of the analysed triplets, p-K − pairs are selected from p-p-K − triplets with Q 3 < 1 GeV/c, while for p-p and p-K + pairs such additional requirement is not used as the correlation functions with and without the Q 3 selection are in agreement.A ±0.1 GeV/c variation of the Q 3 limit is included in the systematic uncertainties of the projector method and it is combined with the uncertainties propagated from the data points of the p-p, p-K + and p-K − correlation functions.The latter include the variations in the selection criteria for particles, variations in the track splitting/merging rejection criteria, and normalisation range.
Figure 1: Correlation functions for (p-p)-K + (panel a 1 ), (p-p)-K − (panel a 2 ), p-(p-K + ) (panel b 1 ), p-(p-K − ) (panel b 2 ) as well as the total lower-order contributions for p-p-K + (panel c 1 ) and p-p-K − (panel c 2 ) as a function of Q 3 .The points represent the results obtained using the data-driven approach, with the statistical and systematic uncertainties represented by the error bars and the green boxes, respectively.The grey bands are the expectations of the projector method, with the band width representing the combined statistical and systematic uncertainties.
For the p-p-K + and p-p-K − systems, the lower-order contributions to the three-particle correlation functions are shown in Fig. 1, using the data-driven approach (data points) and the projector method (grey bands).The three-particle correlation functions of two correlated protons and an uncorrelated kaon, denoted as (p-p)-K + and (p-p)-K − , are shown in panels a 1 ) and a 2 ), respectively.In both cases, the correlation functions are larger than unity in the low Q 3 region as a result of the attractive p-p strong interaction, combined with the repulsive Coulomb and quantum statistics effects (see [46,48] for more details).In the case of correlated p-K + pairs with uncorrelated protons, p-(p-K + ), the correlation function shown in panel b 1 ) is lower than unity as a result of the repulsive p-K + strong and Coulomb interactions, consistent with the measurement reported in [47].In the p-(p-K − ) correlation function, shown in panel b 2 ), the main features of the p-K − interaction are visible [47]: the cusp structure due to the opening of the K 0 n channel as well as the bump due to the Λ(1520) → pK − decay appear at Q 3 ≈ 0.15 GeV/c and Q 3 ≈ 0.7 GeV/c, respectively.The total lower-order contributions to the p-p-K + , panel c 1 ), and p-p-K − , panel c 2 ), correlation functions evaluated with the data-driven approach and the projector method are in agreement.The reduced χ 2 between the two methods, evaluated in the Q 3 range shown in Fig. 1 (21 degrees of freedom) by adding in quadrature the statistical and systematic uncertainties of the data, is 1.6 and 1.4 for p-p-K + and p-p-K − , respectively.Since the projector method does not depend on the third particle mixing, it provides significantly smaller uncertainties than the data-driven approach and hence is used to extract the cumulants.This choice does not affect the final results of the analysis.Figure 2 shows the measured p-p-K + (left panel) and p-p-K − (right panel) correlation functions (data points) compared to the lower-order contributions evaluated using the projector method (grey bands).
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 Correlation functions (data points) for p-p-K + (left panel) and p-p-K − (right panel) compared to the lower-order contributions evaluated using the projector method (grey bands).Statistical and systematic uncertainties are represented by error bars and green boxes, respectively.The band widths represent the combined statistical and systematic uncertainties propagated from the two-particle correlation functions used as input to the projector method.
The corresponding cumulants are extracted using Eq.(3).The measured cumulants include the correlations of the primary particles as well as the feed-down from resonances and particle misidentifications, which need to be accounted for as discussed in [59].The number of correctly identified primary triplets is about 66% of the total sample for both p-p-K + and p-p-K − .The remaining contribution mainly stems from the feed-down of the Λ (17%) and Σ + (7%) hyperon decays into protons and from the φ meson (4%) decay into charged kaons.These numbers are obtained combining the primary and secondary fractions with particle purities.The primary and secondary fractions are extracted using Monte Carlo template fits to the measured distributions of the DCA to the primary vertex [46].The fraction of charged kaons from φ decays is calculated using the expectations from a thermal model [77].The correction of the measured cumulants to obtain the contribution from primary triplets is performed following the decomposition procedure adopted in [59].Due to the absence of theoretical and experimental information on the genuine p-Λ-K ± , p-Σ + -K ± and p-p-φ correlations, the corresponding feed-down contributions to the p-p-K ± cumulants are considered to be flat, without any specific dependence on Q 3 .This assumption was tested in the case of p-p-p correlations [59], where a similar contribution (about 19%) for the Λ feed-down into protons was found.The results obtained by using the flat and non-flat corrections are found to be in agreement within the uncertainties.A flat feed-down contribution in the cumulant corresponds to negligible correlations from genuine three-body effects in the mother particle systems X-Y-Z decaying into p-p-K ± .The dominant feed-down contributions coming from pairwise correlations in the triplets, such as (X-Y)-Z and X-(Y-Z), are already removed as they are included in the lower-order correlations.The resulting cumulants of the correctly identified primary p-p-K + and p-p-K − particle triplets are shown in left and right panels of Fig. 3, respectively.
In the absence of genuine three-particle correlations in the triplets, the measured three-particle correlation function would be identical to the lower-order contributions C 3 (Q 3 ) = C two-body (Q 3 ).Consequently, from Eq. ( 3), the cumulant would be compatible with zero within the experimental uncertainties.In the case of the p-p-K + and p-p-K − systems, the agreement of the measured cumulants with zero is evaluated by performing a χ 2 test in the region Q 3 < 0.4 GeV/c.In this region, the relative momentum for all the pairs in the particle triplets is lower than 0.2 GeV/c, which is the kinematic region sensitive to the strong interaction.Larger Q 3 values are dominated by kinematic configurations in which one of the pairs has a relative momentum larger than 0.2 GeV/c, meaning that only two of the three particles in the triplet can be in a kinematic configuration which is favourable for the strong interaction.Thus, at Q 3 > 0.4 GeV/c the contributions from the two-body interaction dominate.The local statistical significance at Q 3 < 0.4 GeV/c is obtained from the p-value of the χ 2 distribution, which is converted in a number n σ of Gaussian standard deviations.As the systematic uncertainties in different Q 3 bins are correlated, the χ 2 is calculated using the cumulants extracted from each systematic variation with the corresponding statistical uncertainties.Finally, the average χ 2 in the interval Q 3 < 0.4 GeV/c is extracted.The corresponding n σ values are 2.5 and 1.5 for p-p-K + and p-p-K − , respectively.As the Q 3 range of the genuine three-body effects is not known a priori, the n σ values have been evaluated as well in the region Q 3 < 0.2 GeV/c and they are found to be 2.7 and 1.4 for p-p-K + and p-p-K − , respectively.Such results indicate that the measured p-p-K + and p-p-K − correlation functions are compatible, within the quoted significance levels, with the assumption of pairwise correlations in the triplets without additional contributions from genuine three-body effects.More solid conclusions require a larger data sample to reduce the statistical uncertainties and full-fledged theoretical calculations for the three-particle correlation functions.In order to provide a comparison of the kinematic range accessible with this measurement to the study of K + and K − interactions in light nuclei, the momenta of the three particles in the triplets with Q 3 < 0.4 GeV/c are evaluated in the PRF of the two protons.The momentum of each proton in the PRF is p * p < 180 MeV/c, compatible with the typical Fermi momenta of nucleons in light nuclei (e.g. the Fermi momentum in Carbon-12 is about 220 MeV/c).In the rest frame of the two protons, the momentum of the K + and K − ranges from 30 MeV/c (in the first bin in Fig. 3) to 130 MeV/c (eighth bin in Fig. 3).This demonstrates that even for low kaon momenta, i.e. at energies close to the p-p-K + and p-p-K − thresholds, three-body effects such as kaonic bound state formation below threshold or three-body interactions do not contribute significantly to the measured correlation functions.These results provide additional experimental information for theoretical models aiming to understand the role of genuine three-body effects in p-p-K + and p-p-K − systems.The LHC Run 3 data taking will deliver a larger data set for more detailed studies of the K + and K − three-body dynamics in the low momentum region, down to the energy threshold of the triplets.

Figure 3 :
Figure 3: Cumulants for the p-p-K + (left panel) and p-p-K − (right panel) primary triplets.The error bars and the blue boxes on the data points represent the statistical and systematic uncertainties, respectively.The n σ deviations from zero in each bin are shown in the bottom panels, adding in quadrature the statistical and systematic uncertainties of the data.