Investigating the role of strangeness in baryon$-$antibaryon annihilation at the LHC

Annihilation dynamics plays a fundamental role in the baryon$-$antibaryon interaction (B$-\rm{\overline{B}}$) at low-energy and its strength and range are crucial in the assessment of possible baryonic bound states. Experimental data on annihilation cross sections are available for the p$-\rm{\overline{p}}$ system but not in the low relative momentum region. Data regarding the B$-\rm{\overline{B}}$ interaction with strange degrees of freedom are extremely scarce, hence the modeling of the annihilation contributions is mainly based on nucleon$-$antinucleon (N$-\rm{\overline{N}}$) results, when available. In this letter we present a measurement of the p$-\rm{\overline{p}}$, p$-\overline{\Lambda}\oplus\rm{\overline{p}}-\Lambda$ and $\Lambda-\overline{\Lambda}$ interaction using correlation functions in the relative momentum space in high-multiplicity triggered pp collisions at $\sqrt{s} = 13$ TeV recorded by ALICE at the LHC. In the p$-\rm{\overline{p}}$ system the couplings to the mesonic channels in different partial waves are extracted by adopting a coupled-channel approach with recent $\chi\rm{EFT}$ potentials. The inclusion of these inelastic channels provides good agreement with the data, showing a significant presence of the annihilation term down to zero momentum. Predictions obtained using the Lednick\'y$-$Lyuboshits formula and scattering parameters obtained from heavy-ion collisions, hence mainly sensitive to elastic processes, are compared with the experimental p$-\overline{\Lambda}\oplus\rm{\overline{p}}-\Lambda$ and $\Lambda-\overline{\Lambda}$ correlations. The model describes the $\Lambda-\overline{\Lambda}$ data and underestimates the p$-\overline{\Lambda}\oplus\rm{\overline{p}}-\Lambda$ data in the region of momenta below 200 MeV/$c$. The observed deviation indicates a different contribution of annihilation channels to the two systems containing strange hadrons.

of measurements in pp collisions can probe inter-particle distances of around 1 fm and are sensitive to the presence of inelastic channels, below and above threshold [27,32,33]. In this letter we present the measurements of the correlation functions of p-p, p-Λ ⊕ p-Λ and Λ-Λ pairs in pp collisions at √ s = 13 TeV with the ALICE detector [34,35]. To better constrain the interaction, a differential analysis in pair-transverse-mass (m T ) intervals has been performed for the p-Λ ⊕ p-Λ and Λ-Λ pairs. The work presented in this Letter delivers the most precise data at low momenta for p-p, p-Λ and Λ-Λ systems and provides additional experimental constraints for the modeling of the B-B interaction.

Data analysis
The main ALICE subdetectors [34,35] used in this analysis are: the V0 detectors [36] used as trigger detectors, the Inner Tracking System (ITS) [37], the Time Projection Chamber (TPC) [38] and the Time-of-Flight (TOF) detector [39]. The last three are used to track and identify charged particles. The high-multiplicity (HM) sample used in this analysis corresponds to 0.17% of all inelastic pp collisions with at least one measured charged particle within |η| < 1 (referred to as INEL>0) [29,30]. The corresponding HM trigger is defined by coincident hits in both V0 detectors synchronous with the collider bunch crossing and by requiring as well that the sum of the measured signal amplitudes in the V0 exceeds a multiple of the average value in minimum bias collisions. The rejection of pile-up events have been applied by evaluating the presence of additional event vertices as done in [29,31] and a total of 1.0 × 10 9 HM events are selected. Protons and antiprotons are reconstructed using the procedure described in Refs. [29,31]. Primary protons and antiprotons are selected in the transverse momentum range 0.5 < p T < 4.05 GeV/c and pseudorapidity |η| < 0.8. A minimum of 80 out of the 159 available spatial points (hits) inside the TPC are required to obtain high-quality tracks. The TPC and TOF detectors select p (p) candidates by the deviation n σ between the signal hypothesis for the considered particle and the experimental measurement, normalized by the detector resolution σ [29,31]. For candidates with p < 0.75 GeV/c, the particle identification (PID) is performed with the TPC only. For larger momenta, the PID information of TPC and TOF are combined. The candidates are accepted if their |n σ | < 3. To reject non-primary protons (antiprotons), the distance of closest approach (DCA) of the candidates tracks to the primary vertex is required to be less than 0.1 cm in the xy-plane and less than 0.2 cm along the beam axis. Contributions of secondary (anti)protons stemming from weak decays and misidentified candidates are extracted using Monte Carlo (MC) template fits to the measured distance of closest approach (DCA) distributions of the to the primary vertex [25]. The resulting p (p) purity is 99.4% (98.9%). The corresponding fraction of primary particles is 82.2% (82.3%). The reconstruction of the Λ (Λ) candidates, via their weak decay Λ → pπ − (Λ → pπ + ) [40], is performed following the procedures described in Refs. [29,31]. A final selection is applied based on the reconstructed invariant mass [29,31]. The obtained Λ (Λ) purity is 95.2% (96.1%). Primary and secondary contributions for Λ and Λ are extracted in a similar way as for protons, via fits to the cosine of the pointing angle distributions using MC templates. The fraction of primary Λ (Λ) hyperons is about 57%. Secondary contributions from weak decays of neutral and charged Ξ baryons amount to 22%. The remaining fractions are attributed to Σ 0 (Σ 0 ) particles. Systematic uncertainties on the data are evaluated by varying the kinematic and topological selection criteria following [29,31].

Analysis of the correlation function
The main observable in the analysis presented here is the two-particle correlation function C(k * ), which depends on the relative momentum k * evaluated in the pair rest frame [25]. In femtoscopy measurements, the final state is fixed to the measured particle pair and the corresponding correlation function is sensitive to all the available initial, elastic and inelastic, channels produced in the collision [32,33]. For the study of the B-B interaction, the single-channel Koonin-Pratt equation [41] has to be modified in order to accommodate the inelastic contributions stemming from the annihilation channels [32,33]. Assuming that the interaction of the pair in the final state i is affected by the inelastic channels j, the Koonin-Pratt formula is modified by the introduction of an additive term related to the processes j → i [32,33,42]: The first integral on the right-hand side describes the elastic contribution where initial and final state coincide, while the second integral is responsible for the remaining inelastic processes j → i. This last integral depends on two main ingredients: the wave function ψ j (k * , r * ) for channel j going to the final state i and the conversion weights ω j . These latter quantities can be written as is related to the amount of j pairs produced in the initial collision and kinematically available to be converted to the final measured state. Quantitative estimates on these production weights can be obtained combining thermal model calculations of particle yields [43] with kinematics constraints from transport models [44]. If the assumed inelastic wave function ψ j (k * , r * ) is properly accounting for the coupling strength, the corresponding ω s j weight is equal to unity. Recent femtoscopic measurements by the ALICE Collaboration performed on K − -p pairs in pp [27] and in Pb-Pb collisions [45] showed that by changing the colliding system, and hence the size of the emitting source S(r * ), the effects on the C(k * ) due to the inelastic contributions given by the last term in Eq. (1) are enhanced or suppressed. The wave functions ψ j (k * , r * ) related to the inelastic channels are localized at distances r * approximately below 1.5-2 fm and equal to zero above. Hence, performing femtoscopic measurements with a large emitting source, as it occurs in central heavy-ion collisions (r * above 5 fm), results in a correlation function mainly dominated by the elastic contribution, given by the first term of Eq. (1). For this reason, the C(k * ) measured in Pb-Pb can be modeled with the single-channel Lednický-Lyuboshits formula [42] assuming a complex scattering length f 0 , in which the imaginary part I f 0 accounts for an average inelastic contribution from all j channels. These inelastic contributions become more relevant when performing the same measurement in small colliding systems as pp, where the emitting source size is of the order of 1 fm [46] and the modeling of the C(k * ) requires the knowledge of the exact elastic ψ i (k * , r * ) and inelastic ψ j (k * , r * ) wave functions obtained from the solution of a coupled-channel approach [46]. If the theoretical modeling of the interaction properly accounts for the inelastic channels (ω s j = 1), the use of the modified Koonin-Pratt formula in Eq. (1) with a proper estimate of the production weights ω prod j will describe the data in both small and large colliding systems as shown in [45] for the K − -p system. The use of the single-channel Lednický-Lyuboshits model will only be applicable if the wave functions ψ j (k * , r * ) would be strongly suppressed, corresponding to a very weak coupling to the inelastic channels. The B-B interaction investigated in this work is less known with respect to the K − -p case in [27,45], hence two different approaches have been used to calculate the theoretical correlation for the p-p and the Λ-Λ, p-Λ pairs, respectively. For both approaches the CATS framework is used [47]. The genuine p-p correlation is modeled either by assuming a Coulomb-only interaction or by also including a strong interaction from N-N chiral effective (χEFT) potentials at next-to-next-to-next-to-leading order (N 3 LO) [11]. The p-p wave functions, available for S ( 1 S 0 , 3 S 1 ) and P ( 1 P 1 , 3 P 0 , 3 P 1 , 3 P 2 ) partialwaves (PW), have been evaluated within a coupled-channel formalism in which only the coupling to the charge-exchange n-n channel is explicitly included. The formula in Eq. (1) is used for the genuine p-p correlation function with the chiral wave functions for the elastic i = p-p and the charge-exchange channel j = n-n [11]. The wave functions ψ PW X→pp , accounting for the multi-meson annihilation channels j = X, are not currently available. The annihilation contribution is implicitly present in the χEFT potentials in [11] since the parameters of the model are constrained to the most-recent partial-wave analysis on the available p-p and N-N cross sections [4]. The Migdal-Watson approximation [48] is used as an approximate way to explicitly include the additional j = X annihilation channels. This approximation relies on the fact that these X multi-meson channels open below the p-p threshold and hence the momentum dependence of the annihilation potential V X→pp around the p-p threshold can be neglected. The wave functions ψ PW X→pp for each PW can be rewritten in terms of the elastic component as ω PW ψ PW pp→pp , with the weights ω PW to be determined from data. These latter weights are directly connected to the conversion weights ω j in Eq. (1) with the strong coupling term ω s j extended to the different PW states. A detailed estimate of the yields and kinematics of the annihilation channels (ω prod j ), necessary to isolate the strong coupling term in each PW, is not trivial since it involves contributions stemming from multi-pions channels and it should include also intermediate states of resonances strongly decaying into pions. For this reason, the extracted weights ω PW in this work contain informations not only on the coupling strength of the mesonic channels to p-p, but also on the abundances of the contributing multi-meson channels produced in the initial state. The modeled correlation function reads [33]: The first and second terms describe the elastic and n-n contributions, while the last term accounts for the annihilation channels. The degeneracy in spin and angular momentum is embedded in the statistical factors ρ PW . To reduce the number of ω PW weights to be fitted, a study on the shape of the single inelastic correlation terms C PW X→pp (k * ) is performed in each partial wave. The correlations with a different profile in k * are selected, allowing to determine three representative contributions: the 1 S 0 for S states, the 1 P 1 and 3 P 0 for P states. For the two systems containing strangeness, p-Λ and Λ-Λ, no theoretical wave functions are currently available, hence the single-channel Lednický-Lyuboshits analytical formula with a complex scattering length f 0 is used to evaluate the theoretical correlations [23,42]. As mentioned above, in this singlechannel approach, only the elastic contributions are explicitly accounted for in the R f 0 , corresponding to the first term in Eq. 1. The imaginary part I f 0 accounts for an average over all the inelastic contributions of the B-B interaction, mainly dominated by annihilation. The same approach has been used in the ALICE femtoscopic measurements in Pb-Pb collisions [23], which delivered the only available scattering parameters on both the p-Λ and Λ-Λ interaction. For the latter, theoretical predictions are available [16], providing values for the scattering parameters compatible with the ALICE Pb-Pb results [23]. The emitting source in Eq. 1 can be determined as a function of the pair-transverse-mass m T with a data-driven model based on proton-proton correlations [31]. This allows us to investigate the interaction for different particle pairs. The properties of the underlying interaction in p-Λ and Λ-Λ systems do not depend on m T and can hence be better constrained using a m T differential analysis. Considering the available sample, 6 and 3 m T intervals are used for the p-Λ and Λ-Λ measured correlations, respectively. These experimental correlations are compared, in each m T interval, to the Lednický-Lyuboshits model by assuming at first the scattering parameters obtained in the Pb-Pb analysis [23]. Secondly, a simultaneous fit for each pair in all the available m T bins is performed leaving the I f 0 to vary in order to test if a better agreement with the data is achieved. Further discussions on the two different fitting procedures can be found in the next section.
Experimentally, the correlation function is defined as Here N SE (k * ) is the distribution of pairs measured in the same event, N ME (k * ) is the reference distribution of uncorrelated pairs sampled from different (mixed) events and N is a normalization parameter determined by requiring that particle pairs with large k * are not correlated. The mixed-event sample is obtained by pairing particles stemming from events with a similar number of charged particles at midrapidity and a close-by primary vertex position along the beam direction as done in [27,29,30]. The correlation functions of baryon-antibaryon and antibaryon-baryon pairs are combined to enhance the statistical significance for the p-Λ pairs, hence in the following p-Λ denotes the sum p-Λ ⊕ p-Λ. The p-p, p-Λ and Λ-Λ data are fitted with a total correlation function where N D is a normalization constant fitted to data. The default fit range is 0 < k * < 500 MeV/c. The mod- from Eq. 2 and using the Lednický-Lyuboshits model, and the residual secondary contributions weighted by the λ i parameters [25]. The genuine contributions for p-p, p-Λ and Λ-Λ amount to λ p-p = 66.5%, λ p-Λ = 45.8% and λ Λ-Λ = 30.9%, respectively. Residual contributions involving pairs measured in this work are modeled assuming the corresponding theoretical predictions mentioned above. Contributions involving Σ ±,0 (Σ ±,0 ) and Ξ −,0 (Ξ +,0 ) are considered to be constant in k * due to the limited theoretical knowledge, and amount to 10.1%, 44.6% and 65.7% for p-p, p-Λ and Λ-Λ, respectively. A crosscheck on these residuals by assuming a strong interaction based on the scattering parameters extracted in Pb-Pb measurements [23] was performed and differences in the extracted results with respect to the constant assumption are found to be negligible. A variation of ±10% to the upper limit of the default fit range is applied for evaluating the systematic uncertainties. Additionally, the systematic uncertainties related to the λ i parameters are evaluated based on variations of the amount of secondary contributions to each measured particle species, where the largest source of uncertainty stems from the ratio Σ 0 : Λ = 0.33 ± 0.07 [31,43,[49][50][51]. In addition to the feed-down contributions, a correction for finite experimental momentum resolution has to be taken into account for a direct comparison with data [25]. The size of the emitting source employed in the calculation of C model (k * ) for the three B-B pairs is fixed from the data-driven analysis of p-p pairs, which demonstrates the existence of a common Gaussian core as a function of m T for all baryon-baryon pairs when contributions from short-lived strongly decaying resonances are properly included [31]. For the p-p pairs, the core source size at the corresponding m T = 1.45 GeV/c 2 is r core = 1.06 ± 0.04 fm and the associated effective Gaussian source size is r 0 = 1.22 fm. The core radii for the p-Λ and Λ-Λ m T bins presented in this letter are r core ( m T = 1.75 GeV/c 2 ) = 0.95 ± 0.04 fm (r 0 = 1.15 fm) and r core ( m T = 2.12 GeV/c 2 ) = 0.87 ± 0.04 fm (r 0 = 1.11 fm), respectively. The second term in Eq. 4, C background (k * ), accounts for non-femtoscopic effects due to energy-momentum conservation at large k * [25] and to minijet phenomena arising from hard processes at the parton level, largely present in the measurement of B-B correlations: A data-driven approach is employed using PYTHIA 8.2 [52] to model the mini-jet part contained in C background (k * ). The particle production in such simulations is associated to two processes: particles stemming from a common parton (common ancestors), leading to the minijet component, and particles coming from different partons (non-common ancestors), responsible for the non-jet part. The C minijet (k * ) in Eq. (5) is given by a linear combination of the common (C C (k * )) and non-common (C NC (k * )) contributions weighted by a factor w C and (1 − w C ), respectively. The ancestor weight w C is a free parameter in the fit of C tot (k * ) to the data. The common and non-common correlations obtained from PYTHIA 8.2 are fitted with a product of three Gaussian functions up to k * = 2500 MeV/c, providing good agreement with the simulated data. To account for remaining non-femtoscopic effects at large k * [25], a linear baseline C baseline (k * ) = a + bk * is added to the ancestors term in Eq. (5). The coefficients a and b are fixed by fitting C background (k * ) to the data in the region of 400 < k * < 2500 MeV/c. The results for p-p pairs are shown in Fig. 1. The band represents the 1σ uncertainty associated to the template fitting. The shape of C background (k * ) agrees within uncertainties with the data in the region above k * ≈ 200 MeV/c, where the non-flat behavior of minijet contributions is visible. A change of ±10% in this range and a quadratic polynomial are included to estimate the systematic uncertainty related to the total background. Similar results and conclusions are obtained for the p-Λ and Λ-Λ systems.

Results
The correlation functions for p-p and for two representative m T bins of p-Λ and Λ-Λ are shown in Fig. 2 and in Fig. 3, respectively. The results for the remaining m T bins are presented in Figs. A.1, A.2 of Appendix A.1. The lower panels show the statistical deviation between data and model expressed in terms of numbers of standard deviation n σ . The width of the band represents the total uncertainty of the fit. The grey boxes correspond to the systematic uncertainties of the data. They are maximal at the lowest k * bin and amount to 1%, 4% and 10% for p-p, p-Λ and Λ-Λ pairs, respectively. The C p-p (k * ) correlation is compared first to a Coulomb-only interaction and secondly to a Coulomb + strong interaction from N-N χEFT potentials with wave functions for the n-n → p-p process explicitly included [11]. Results for this latter scenario are obtained by evaluating the genuine p-p correlation in Eq. (2) with only the first two terms and shown in blue in Fig. 2. The opening of the n-n channel above threshold, expected as a cusp structure in the C(k * ) at k * ≈ 50 MeV/c, is not visible in agreement with the weak coupling already measured in scattering experiments [4]. The chiral model underestimates the data in the region below 200 MeV/c and it cannot reproduce the enhancement above unity of the C(k * ) as k * approaches zero. This increase is not described either by assuming only the Coulomb attraction (green band), showing that annihilation is largely present close to threshold as k * → 0 MeV/c. The contributions to the p-p correlation from the multi-meson annihilation channels, produced as initial states which feed into the measured p-p system, are not explicitly accounted for in the chiral potential and hence the last term ∑ PW C PW X→pp (k * ) in Eq. (2) is currently missing in the fit shown by the blue band. The red band in Fig. 2 represents the results obtained from the explicit inclusion of the annihilation channels in the third term of Eq. (2) via the Migdal-Watson approximation. The corresponding fit provides a better description of the data in the low k * region where annihilation is dominant. The extracted coupling weights ω PW from this femtoscopic fit are ω1 S 0 = 1.19 ± 0.10 (stat) ±0.19 (syst) and ω3 P 0 = 40.04±4.06 (stat) ±4.24 (syst), while ω1 P 1 is compatible with zero. The hierarchy of the coupling weights in the different PW agrees with the inelasticity parameters η obtained in the recent partial-wave analysis [4].  Figure 2: (Color online) Measured correlation function of p-p pairs. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The Coulomb only interaction is shown by the green band. The blue band represents the fit performed using N 3 LO χEFT potentials [11] with elastic and n-n coupled-channel. The inclusion of annihilation channels is shown by the red band, along with the C background (k * ), multiplied by the normalization constant N D obtained in the fit. The reported average χ 2 /NDF is evaluated in the k * interval [0, 400] MeV/c and it includes correlations between the data points. Lower panel: n σ deviation between data and model in terms of numbers of statistical standard deviations.
For the systems containing strangeness, the Migdal-Watson approach cannot be employed since only scattering parameters for the p-Λ and Λ-Λ interaction are available [23]. The values of R f 0 and I f 0 obtained in Pb-Pb measurements are employed in the Lednický-Lyuboshits analytical formula [23,42] to model the p-Λ and Λ-Λ genuine correlation functions. In Fig. 3, the results obtained modeling the p-Λ and Λ-Λ theoretical correlations with the Lednický-Lyuboshits model described in Sec. 3 are shown. The first tested scenario assumes the scattering parameters extracted from Pb-Pb results [23] and the results are denoted by light green bands. It can be expected that if the direct contributions of the annihilation channels are negligible, the values extracted in Pb-Pb will reproduce well also the pp data in this analysis. As can be seen from the right panel in Fig. 3, this first approach reproduces the measured Λ-Λ correlation function, with an average χ 2 /NDF = 2.8 evaluated in the k * interval [0, 400] MeV/c but it clearly underestimates the p-Λ correlation data in the k * region below 200 MeV/c. A similar trend is observed when performing the fit to the p-p measured correlation with the Lednický-Lyuboshits approach used in the Pb-Pb results [23] as shown in Fig. A.3 in the Appendix A.2. The discrepancy hence, as in the p-p case, has to be attributed to a larger amount of annihilation channels feeding into the p-Λ system with respect to the Λ-Λ pairs. To validate this interpretation, a simultaneous fit in all the m T bins is performed leaving free to vary the imaginary part of the scattering length I f 0 , accounting for inelastic channels, and the effective range d 0 . The negative real part of the scattering length R f 0 , indicating either a repulsive elastic interaction or a possible bound state, is kept fixed to the Pb-Pb results [23]. To reach a reasonable agreement of the model with p-Λ data, I f 0 has to be increased by approximately a factor 5.3, while the change in the extracted d 0 is negligible. Such a discrepancy can be attributed to the failure of the single-channel Lednický-Lyuboshits model to properly accommodate the direct contribution of inelastic channels (last term in Eq. (1) The different results for the p-Λ and Λ-Λ systems may also be related to a different amount of initially produced multi-meson states feeding into the two B-B pairs. To substantiate this scenario, a study of the two-meson channel contributions (ππ, πK) is performed using the EPOS transport model [44]. The fraction of two-mesons pairs f 2M→BB produced in the initial collision and kinematically available to produce B-B pairs with low k * is estimated. The latter is obtained by dividing the amount of mesonmeson pairs initially produced, having a center-of-mass energy above the B-B threshold and leading to B-B pairs at low k * , by the total number of produced two-meson pairs kinematically allowed to create the B-B pairs. Based on this study and considering these kinematics considerations, a similar amount (≈ 6.4%) is found for p-Λ and Λ-Λ pairs, indicating that the above effect is related to the properties of the p-Λ and Λ-Λ interaction. To quantify the final relative amount of annihilation channels feeding to the p-Λ and Λ-Λ systems, the fractions have to be multiplied by the corresponding coupling constant g, obtained within an SU(3) Lagrangian by evaluating the trace of the meson-baryon interaction term [53]. Within this simplified calculations, the coupling strength for the p-Λ system is found to be approximately 3.3 times larger than for the Λ-Λ pairs. The estimated contribution g × f 2M→BB , although limited to only two-meson channels, for p-Λ pairs is found to be about 6 times larger than for Λ-Λ pairs, indicating a different annihilation contributions occurring in p-Λ and Λ-Λ interaction which is confirmed by the measured correlation functions in Fig. 3. These estimations, even though based on a qualitative approach, clearly indicates that the annihilation for the Λ-Λ interaction is present but it should not be largely dominant over the elastic part. More input from theory is needed in order to claim if such a condition is ideal for the formation of bound-states in the Λ-Λ system. The results for the p-Λ system, however, clearly point to a much larger presence of the annihilation channels, which might reduce the possibility to create baryonia. The data presented in Fig. 3 represent the most precise data currently available on p-Λ and Λ-Λ pairs and can provide constraints for theoretical models on these interactions.
In conclusion, femtoscopic techniques have been adopted to study the annihilation dynamics in p-p, p-Λ and Λ-Λ systems. A quantitative determination of the effective coupling weights, connected to the annihilation channels present in p-p, has been obtained adopting a coupled-channel approach with N 3 LO χEFT potentials [11]. The largest couplings have been obtained in the spin triplet P ( 3 P 0 ) and singlet S ( 1 S 0 ) state. The inclusion of these inelastic channels leads to a better agreement between data and model in the region of k * below 50 MeV/c, indicating a wide presence of annihilation channels close to threshold. The scattering parameters obtained in Pb-Pb collisions [23] have been used to model the p-Λ and Λ-Λ data using the Lednický-Lyuboshits formula. A consistent description of the Λ-Λ correlation is achieved while an increase of the I f 0 in the p-Λ interaction is needed to improve the agreement with the p-Λ data. These results, confirmed by kinematics and SU (3)  a larger contribution in p-Λ from annihilation channels in comparison to Λ-Λ. The ALICE data shown in this work delivered the most precise measurements on p-p, p-Λ and Λ-Λ systems at low momenta and suggest that baryonia are unlikely to occur in p-p and p-Λ systems due to the large annihilation contributions present for these pairs. A modeling of the B-B interaction for systems as p-Λ, based on optical potentials, and a quantitative estimate of production of the multi-meson annihilation channels in the collisions, can provide a better understanding of the elastic and the annihilation term which can help to strengthen final conclusions on possible bound states.   [23], are shown by the violet band. The scattering parameters used as input for the Lednický-Lyuboshits calculations include only the n-n contribution as coupled-channel [54][55][56]. The model completely underestimates the k * region from 50 to 150 MeV/c. As can be seen in Fig. 2, the annihilation channels play a role in this intermediate region and a better description of the data is achieved when using the Migdal-Watson approximation to include them. This is a clear indication that the multi-meson channels are explicitly needed to model the current measured p-p correlation function. The Lednický-Lyuboshits calculation also overestimates the coupling to the n-n channel, as can be seen from the large cusp structure at k * ≈ 50 MeV/c not present in the data. Lower panel: n σ deviation between data and model in terms of numbers of statistical standard deviations.