
==== Front
Eur Phys J C Part Fields
Eur Phys J C Part Fields
The European Physical Journal. C, Particles and Fields
1434-6044
1434-6052
Springer Berlin Heidelberg Berlin/Heidelberg

12319
10.1140/epjc/s10052-023-12319-8
Letter
Searching for a dark matter particle with anti-protonic atoms
Doser Michael michael.doser@cern.ch

1
Farrar Glennys 2
Kornakov Georgy 3
1 grid.9132.9 0000 0001 2156 142X CERN, Esplanade des Particules 1, 1211 Geneva, Switzerland
2 https://ror.org/0190ak572 grid.137628.9 0000 0004 1936 8753 Center for Cosmology and Particle Physics, New York University, New York, NY 10003 USA
3 grid.1035.7 0000000099214842 Warsaw University of Technology, ul. Koszykowa 75, 00-662 Warsaw, Poland
18 12 2023
18 12 2023
2023
83 12 114916 5 2023
3 12 2023
© The Author(s) 2023
https://creativecommons.org/licenses/by/4.0/ 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/.
Funded by SCOAP3. SCOAP3 supports the goals of the International Year of Basic Sciences for Sustainable Development.
A wide range of dark matter candidates have been proposed and are actively being searched for in a large number of experiments, both at high (TeV) and low (sub meV) energies. One dark matter candidate, a deeply bound uuddss sexaquark, S, with mass ∼2 GeV (having the same quark content as the hypothesized H-dibaryon, but long lived) is particularly difficult to explore experimentally. In this paper, we propose a scheme in which such a state could be produced at rest through the formation of p¯–3He antiprotonic atoms and their annihilation into S + K+K+π-, identified both through the unique tag of a S=+2,Q=+1 final state, as well as through full kinematic reconstruction of the final state recoiling against it.

Polish Ministry of Education and Science2022/WK/06 http://dx.doi.org/10.13039/100000001 National Science Foundation PHY-2013199 Polish National Science Centre2022/45/B/ST2/02029 http://dx.doi.org/10.13039/100000893 Simons Foundation issue-copyright-statement© EDP Sciences, Societa Italiana di Fisica (SIF) and Springer-Verlag GmbH, DE, part of Springer Nature 2023
==== Body
pmcIntroduction

The identification of the nature of dark matter is one of the most pressing concerns of particle physics. Candidates for dark matter particles cover an exceedingly wide range, with masses predicted to lie between 10-22 eV and 1010 eV, depending on the specific type of candidate and not including primordial black holes. Searches for candidates with masses above the GeV scale are generally accelerator-based: recent reviews of these can be found in [1, 2]. Putative low-mass dark matter candidates (e.g. axions, ALP’s, chameleons,...), whose detection requires highly sensitive devices, among them quantum sensors which can detect discrete state changes from one quantum state to another that interactions with dark matter can trigger, are equally actively being searched for [3, 4]. One dark matter candidate, the sexaquark (S) [5, 6], would however not have been detected by either approach, while prior searches for similar states in the GeV region (notably, the H-dibaryon [7]) would not have been sensitive to the S [8].

The sexaquark is a hypothesized deeply bound, long-lived or stable 0+ state of uuddss quarks with B=+2,S=-2 and Q=0. The quark content and quantum numbers are the same as the H-dibaryon, however the properties are distinctly different since the H-dibaryon was envisaged to be relatively weakly bound with a lifetime O(10-10s). Accelerator experiments have excluded the existence of the H-dibaryon, unless it is a weakly bound “molecule” like the deuteron, but an S would not have been discovered by experiments to date (see [8] for a survey and suggestions for strategies). Three attributes of a stable S make it very difficult to detect:The S is neutral and a flavor singlet, so it does not couple to photons, pions and most other mesons, nor does it leave a track in a detector.

The S has no pion cloud and is expected to be more compact than ordinary baryons. This means the amplitude for interconversion between S and baryons is small.

The mass of the S makes it difficult to distinguish from the much more copious neutron.

If mS<mΛ+mp+me=2.054 GeV, its decay must be doubly-weak and its lifetime would be greater than the age of the Universe [9], making it an interesting dark matter candidate. Direct detection experiments have not yet probed the relevant mass range with sufficient sensitivity to have excluded this possibility [8]. It has also been shown that the stability of neutron stars is not affected by the existence of sexaquarks, due to quark deconfinement [10]. In addition to its significance as a Standard Model dark matter candidate, the existence of the S and production of SS¯ pairs in e+e- collisions could make an undetected contribution to Rhad and increase the Hadronic Vacuum Polarization (HVP) [11]. This would reduce at some level the tension between the dispersive determination of the HVP based on experimental measurements of Rhad [12, 13] and the predictions of lattice QCD [14, 15], which are much closer to the value needed to explain the measured muon anomalous magnetic moment [16] (Note however that the recent CMD-3 data has greatly reduced the discrepancy between the dispersive value and g-2 [17] and lattice QCD predictions [15]; moreover the change in g-2 from SS¯ production in e+e- collisions is expected to be small). These considerations, combined with the importance of fully exploring the hadron spectrum, makes the existence of the S an important question even if the S does not explain Dark Matter.

Formation of multiquark states (e.g. tetraquarks or pentaquarks [18]) is well established in high energy particle interactions [19], with their dense environment of numerous quark-antiquark pairs, and such a production channel for sexaquarks has also been proposed and investigated [20, 21]. While in such high energy collisions, production rates of multiquark states may well be high, detection of all involved participants and complete kinematical reconstruction of the potentially produced sexaquark is very difficult, given the generally high multiplicity of participants and the low mass of the proposed sexaquark. In contrast, production at low energy, which facilitates full reconstruction, requires establishing a situation with high numbers and high local density of quarks, and in which furthermore, at least two strange quarks are produced. While searches in e.g. Υ decays have been attempted [22, 23], annihilation of antiprotons (or antineutrons) with nuclei, with a concomitant larger number of quarks and antiquarks, may have a better potential to achieve these conditions. On one hand, annihilation between an antiproton and a nucleon results - in the naive quark model - in the interaction of three quark-antiquark pairs, leading potentially to the formation of multiple strange-antistrange pairs, as has been pointed out in proposals to investigate double-strangeness production with antiprotons [24]. On the other hand, the spectator nucleons potentially provide additional light quarks; such multi-nucleon interactions involving antiprotons have been observed [25] in antiproton-nucleus annihilations at rest and have been studied theoretically [26]. This strategy has been applied in the past to (unsuccessful) searches for similar states, in particular for the H dibaryon [7] with identical quark composition as the S(uuddss), but with different mass, physical structure, and thus decay characteristics. In particular, previous searches for the H-dibaryon (with its ΛΛ quark composition) have either relied on decays into ΛΛ, Λpπ or Λn [27–29] or on exclusive co-production with Λ¯Λ¯ in Υ decays [22] (requiring formation and binding of exactly 12 quarks and antiquarks simultaneously, with an estimated branching fraction of O(10-11) [8]).

Here, we propose a novel approach that, while having similarities to the above, would lead to a very clean experimental tag for the production of stable sexaquarks in a low energy environment characterized by a high local density of quarks and antiquarks. We will investigate in the following a specific antiproton-nucleus system, first conceptually, to extract the salient characteristics of a reaction potentially resulting in production of a uuddss sexaquark, and then investigate with a simulation the sensitivity of such a production process against potential background processes.

Formation process

Our proposed scheme for producing S(uuddss) sexaquark states at low energy starts with antiprotonic atoms of 3He. Simply considering the masses of the involved participants results in a total energy budget of mp¯+mp+mp+mn ∼ 3750 MeV for a system at rest. As we shall show below, the formation of S(uuddss) in the process p¯–3He →SK+K+π- should be detectable for relative production rates down to 10-9 of all annihilations, with a clean experimental signature that allows determining precisely the mass of the (invisible) recoiling dark matter candidate S(uuddss) system.

Annihilation of the p¯ with an n or a p of the nucleus results – in the naive quark model – in quark annihilation and/or quark rearrangement; for simplicity, we will first focus on the case of p¯p annihilation with no additional spectator nucleons. If two quark-antiquark annihilations occur, then two ss¯ pairs can be produced (albeit with low probability, and ignoring energy limitations in a first step), with either the uu¯ or the dd¯ surviving, in other words, uud + u¯u¯d¯ → ussu¯s¯s¯ or dssd¯s¯s¯. Such processes, leading to the production of ϕϕ pairs, have been observed to occur in p¯p annihilations in flight [30], once the available rest-frame energy is sufficient.

Since even the lightest final states containing s¯s¯ss quark combinations (4K, Ξ0Ξ0¯, ϕϕ and ϕKK) all exceed 2mp, energy conservation requires that for annihilation at rest, this process occur together with at least one additional nucleon, either a p or an n. For p¯–3He, this thus requires a three-nucleon interaction with a single spectator nucleon (in the case of the p¯d subsystem of p¯–3He) or a four-nucleon interaction (in the case of p¯–3He); here, we will consider the later process in the form of the quark picture involving all four nucleons (Fig. 1).Fig. 1 Quark rearrangement and annihilation graph for the formation of a uuddss sexaquark state in p¯–3He annihilations. S denotes the putative S(uuddss) sexaquark state

For multi-nucleon annihilations, the question of whether cross-nucleon quark clusters can form is of relevance. In a study of deep-inelastic electron scattering from 3He [31], a geometric probability for a six-quark cluster configuration to occur in the 3He nucleus was found to be 16% (and even that of a 9-quark cluster remains non negligible, at 1%). In deuterons, the same probability is estimated to be 10% [32]; theoretical considerations [26] indicate both the production rate of individual two-body final states at a level of O(10-5) as well as an enhanced strange particle yield.

Not surprisingly, corresponding experimental indications for multi-nucleon annihilations have been found in several systems, such as the process p¯d→pπ which requires all nine (anti)quarks of the system to interact, and has a branching ratio of 10-5 [25]. Similarly, three-nucleon processes with additional spectator nucleons have been argued on the basis of the high yield of Λ production, both in deuterium and in nitrogen [25, 33]. Similarly, while searching for deeply bound K¯K¯-nuclear states in p¯4He annihilations, the OBELIX experiment has identified final states containing two K+, accompanied by two π- and missing momentum and energy, implying necessarily the involvement in the annihilation of more than one nucleon. This final state was interpreted as 2K+Σ-Σ-ps→2K+2π-(2nps) with a yield of (0.172±0.038) × 10-4 [34]. Annihilation processes involving multiple nucleons are thus well established, with production rates in the order of 10-6∼10-4 and with indications of the presence of multi-nucleon annihilations.

In the case of p¯–3He (Fig. 1), after annihilation of two of the antiquarks of the p¯ with two of the quarks of other nucleons to produce two s¯s pairs, the resulting intermediate quark composition is then [ussu¯s¯s¯][ddu][uud] or [dssd¯s¯s¯][ddu][uud], depending on which quarks of the antiproton-proton pair annihilate (the former occurring two times more frequently than the later).

Pairing off the 3rd and 4th nucleons’ quarks results in [ussddu][u¯d][s¯u][s¯u] (in 2/3 × 1/3 of the cases), in [dssuud][d¯u][s¯d][s¯d] (in 1/3 × 1/3 of the cases), in [ussddu][u¯u][s¯d][s¯u] (in 2/3 × 2/3 of the cases) and in [dssuud][d¯d][s¯d][s¯u] (in 1/3 × 2/3 of the cases).

The outcome is then production of the (putative) S(uuddss) state in coincidence with two kaons and one pion:1 (uuddss)+π-+K++K+(in2/3×1/3of the cases)

2 (uuddss)+π++K0+K0(in1/3×1/3of the cases)

3 (uuddss)+π0+K0+K+(in2/3×2/3of the cases)

4 (uuddss)+π0+K0+K+(in1/3×2/3of the cases)

In all cases, the total strangeness of the final state mesons is +2, while the sum of charges ΣQ=+1. Final state masses sum to m(S) + 2 m(K) + m(π), allowing producing S with mS ∈[0,2600] MeV. The experimentally cleanest tag is (1), since all final state particles (except for the invisible S(uuddss) state) are easily detectable and uniquely identifiable. Furthermore, kinematic reconstruction of the missing mass is unambiguously possible, since the initial energy and momentum of the p¯–3He system are known. Note that also p¯d annihilation can lead to the production of the sexaquark (p¯d→Ξ¯0S, Fig. 2, left), but results in a much reduced mass window (mS< 1500 GeV) than is the case for p¯–3He and is experimentally much more difficult to detect due to the dominant Ξ¯0→Λ¯π0,Λ¯→p¯π+ decay.Fig. 2 Three-nucleon rearrangement and annihilation graphs for p¯d leading to the formation of a final state containing the K+K+π- topology or of the S(uuddss). In both cases, the spectator p would also form part of the final state. Top: graph for p¯d→Ξ0¯S. Bottom: p¯d→Ξ-K+K+π-. In p¯–3He annihilations, this channel, in which the Ξ- can be fully reconstructed via its Ξ-→Λπ-,Λ→pπ- decay, and which is accompanied by a spectator proton ps, is a useful cross check for the presence of multinucleon annihilation, but also for triggering, tracking and particle identification capabilities of the detector used to search for S(uuddss)

We now discuss a number of standard model processes that result either in similar final state topologies, or that produce final states containing the characteristic K+K+π- topology (with additional mesons or baryons, but which might escape detection in a real detector) to determine whether these could result in potential backgrounds. No standard model processes result in the signal topology with a single missing particle; in all cases, particle misidentification plays a role in potential backgrounds:The most challenging background channel stems from p¯p(d)→K+K0¯π+π-(d), with the π+ being misidentified as a K+ and both the K0¯ and the spectator d escaping detection. This channel with a production rate of O(10-3) will be studied in more detail below (section 3).

Additional potential backgrounds stem from incomplete detector solid angle coverage. Annihilation of antiprotons on a deuteron-like two-nucleon system inside 3He with the fourth nucleon (a proton ps) remaining as spectator, p¯d(ps)→Ξ-K+K+π-(ps) is kinematically allowed (Fig. 2, right). Here, the decay of the Ξ-→Λπ-;Λ→nπ0 produces a recoiling π-, a neutron and two photons. Misidentification of this topology as signal requires missing the π- of the Ξ- decay and the spectator proton, but also the missing mass recoiling against the K+K+π- system will not result in a unique recoil mass.

A related background with a spectator neutron (ns) stems from: p¯d(ns)→Λ0K+K+K-(ns), which is is kinematically allowed, as is the related process p¯d(n)→Σ0K+K+K-n. While mis-identification of the K- as π- can lead to the signal topology, the missing mass recoiling against the K+K+“π-” system will not result in a unique recoil mass.

A more challenging, also related, background stems from p¯pp(ns)→Ξ0K+K+π-(ns). Fully neutral decays of the Ξ0→Λπ0;Λ→nπ0 produce a recoiling neutron and four photons, in addition to the spectator neutron. Here, only the missing mass recoiling against the K+K+π- system, which will not result in a unique recoil mass, can be used to discriminate this potential background.

p¯n+p+p, with p¯n→K-K+π-. This background requires both misidentification of the K- as K+ as well as both protons escaping detection.

To summarize: the potential backgrounds to the signal, either via rescattering or multi-nucleon annihilation processes, cannot result in a K+K+π- final state recoiling against a single state; only particle misidentification (or particles escaping detection) can result in such a final state. In the later case, the resulting topology requires a study of the missing mass distribution to establish down to which level sensitivity to the searched for S(uuddss) is feasible. The fake and misidentification rates are evaluated in the following with a simple Monte Carlo simulation and using known branching ratios for p¯p and p¯n annihilations into final states with at least one K+.

Simulations

The p¯–3He annihilation process has been studied using the Geant4 simulation framework [35–37], obtaining an estimate of the production rate of the topologies of interest with final states containing K+K+π- and possible backgrounds to the process. In addition, for rare multi-nucleon processes we have studied the final states using Monte-Carlo phase space generators implemented in the ROOT framework [38, 39]. Fig. 3 provides a sketch of the searched-for topology within a generic tracking and particle identification detector, such as a compact TPC.Fig. 3 Sketch of the searched-for π-K+K+mm topology (where mm is the recoiling missing mass) with a generic compact cylindrical TPC-like tracking detector surrounding Penning trap structures within which the formation and annihilation process would occur. The overall device would be embedded in a multi-T solenoidal magnetic field

The annihilation at rest of antiprotons on nuclei is simulated in several steps, following the procedure as in [40]. First, 1 keV antiprotons are shot at a nm-thin 3He target where via electromagnetic interaction they are slowed down until captured by an atom. The thin target prevents new interactions of any produced states with the bulk material, providing similar conditions to those of ultra-high vacuum. Then, the annihilation of the antiproton with one of the nucleons is described by the Geant4 FTFP_BERT_HP physics list [41], which assumes a quark gluon string model for high energy interactions.

After each capture and annihilation event, information on the produced final state particles such as their type, four-momentum vector, and charge is accessible for subsequent analyses. Final states containing K+K+π- candidates are searched for in view of subsequent analysis. Zero events with this pure final state were found in 109 simulated annihilations. This is expected due to the lack of rare multi-nucleon annihilation processes in the code. However, similar final states are possible when there is a wrong tagging of the particle type (misidentification) or some of the final state particles are not measured (inefficiency) in the experiment.

For simulating these effects we have assumed a minimal momenta for particles for being reconstructed: 10 MeV/c for π, 25 MeV/c for K, 40 MeV/c for protons, and 80 MeV/c for deuterons. The inefficiency of the tracking device is of 0.5 %, momentum resolution (Δp/p) is set to 1 %, and a probability of wrong particle type assignment (confusing a pion with a kaon and vice-versa) is set to 0.5 %. With those assumptions, we have analysed the final states containing reconstructed K+K+π- triplets.

The analysis of the tagged final states produced in Geant4 shows that most of the background events have the following origin sorted from most likely to less likely: lost low energy deuteron or proton with one negative pion and two positive pions misidentified as kaons.

lost low energy proton or deuteron, kaon and two pions, one misidentified as the second kaon.

inefficient reconstruction of proton/deuteron, kaon and two pions, one of them misidentified.

The expected background distributions for 1010 antiproton annihilations on 3He and the above main sources are shown in Fig. 4. The condition of full kinematic reconstruction of the final state is able to reject most of the events with similar topologies. In the mass region around the deuteron, we expect to have 60–70 background counts per 0.052 GeV/c2, increasing towards higher masses, reaching approximately 350 counts at the maximum at 2.2 GeV/c2. In the same figure the expected S signal distribution after smearing the momentum of the reconstructed tracks is also included for an assumed production rate in p¯–3He annihilations of 10-8 and a mass of 1.89 GeV/c2. The observed width for a stable S is only due to reconstruction and tracking uncertainties and is expected to be in the order of 0.01 GeV/c2.

In addition to the Geant4 simulations we have performed checks for the rare multi-nucleon processes with branching ratios below 10-5 [24]. In order to study the expected background distribution in events with two kaons produced together with Λ or Ξ baryons, we studied the missing mass distribution assuming an inefficient reconstruction of the weakly decaying particles. The missing mass distribution starts at 2.2 GeV/c2, thus not representing a real source of background for the S searches.

We also used the phase-space generator to test the missing mass distribution for a possible channel including in addition to the S, two kaons and a charged pion further neutral pions. In this case, the sharp peak in the missing mass disappears and it is shifted towards higher values. Thus, this potential source of background will not contribute to the expected signal region.

The background distribution obtained using the Geant4 framework has been used as the basis to estimate the sensitivity to different branching ratios of S in an experiment. An example of the S with a mass of 1890 MeV/c2 and a branching ratio of 10-8 is shown in Fig. 4.

The simulated points then are used to estimate the sensitivity using the Asimov significance test [42] for 1010 annihilation events. The resulting 3 and 5-σ sensitivities to branching ratios as a function of the S mass (mS) are shown in the right panel of Fig. 4. The highest sensitivity, to a branching ratio of 10-9 is in the interval between 1200 and 1600 MeV/c2, decreasing to 10-8 for masses around 2200 MeV/c2. In the most favourable mass region between 1800 and 2000 MeV/c2 the predicted sensitivity is 3–4 ×10-9.

The simulated number of annihilations is equivalent to an experiment lasting for a year at the Antiproton Decelerator (AD) at CERN providing 100 keV antiprotons in bunches of 106∼7 particles every 100 s [43, 44]. In order to ensure the detection of any of the recoiling deuterons (or protons) which are the main source of background events the experiment should use a light target, preventing absorption of the final-state particles. For this reason, reconstruction using liquid or gaseous targets would be very difficult.

An alternative is provided by the formation of p¯–3He in Penning traps [40]. In-trap formation in ultra high vacuum (UHV) ensures that any low transverse energy nuclear remnants can be detected and identified axially (e.g. via a micro channel plate (MCP) and time-of-flight detectors), while higher energy mesons can be detected via a traditional cylindrical tracking detector embedded within the Penning trap’s solenoidal magnetic field which allows to reconstruct accurately the momentum of any produced charged particles (Fig. 3).Fig. 4 Top: Reconstructed missing mass distribution from 1010 simulated events in Geant4 of the K+K+π- system in p¯–3He →uuddss+K++K++π- together with background. The peak at 1.89 GeV/c2 corresponds to production of S(uuddss) with a branching ratio of 10-8. Bottom: Experimental sensitivity (3 and 5-σ) to the branching ratios of S(uuddss) in p¯–3He →uuddss+K++K++π- as a function of mS for 1010 events

Discussion

The above discussion demonstrates that S(uuddss) sexaquarks, a specific dark matter candidate, can potentially be formed and be detected with high sensitivity and over a very low background in antiproton-nuclear annihilations at rest, with the most appropriate system for detection being p¯–3He annihilations. This process assumes that re-arrangements of quarks can happen across nucleon boundaries. Experimentally, such multi-nucleon interactions have been shown to occasionally occur in antiproton annihilation [25, 45]. The same process could also occur in heavier nuclei than 3He, resulting in final states of the type (uuddss)+K+K+π-+n×[uud]+m×[udd]. Kinematically, such a state with supernumerary protons or neutrons might decay into ΛΛN′ or would result in multiple particles recoiling against the K+K+π- system, and thus entail a less clean final state topology than for p¯–3He. Finally, while the formation of the sexaquark system in antiproton–tritium interactions should proceed as in 3He, the resulting final state does not have the same characteristic S=+2 signature.

Formation of the p¯–3He antiprotonic atom through co-trapped ionic precursor systems, e.g. p¯ + 3He+, can rely either on pulsed formation via charge exchange between co-trapped antiprotons and anionic atomic ions that are first photo-ionized and then immediately excited into Rydberg states [46], similar to the pulsed formation of antihydrogen [47], requiring production of metastable 3He- [48], or on three-body interactions between antiprotons and a trapped ensemble of cationic 3He+ atomic ions, similarly to the manner in which antihydrogen atoms can be formed [49, 50].

Alternatively, replacing 3He with tritium would result in an (Sπ-K+K0) final state in a process similar to that of Fig. 1. Contrary to antiproton annihilation on 3He to (Sπ-K+K+) however, which has no standard model background with the same final state topology, for antiproton annihilation on tritium the standard model process p¯p(nn)→π-K+K0¯(nn) with two spectator neutrons occurs at a branching ratio level of O(10-3).

In summary, we have shown that formation of S(uuddss) in the process p¯–3He →SK+K+π- should be detectable for relative production rates down to 10-9 of all annihilations, with a clean experimental signature that allows determining precisely the mass of the (invisible) recoiling dark matter candidate S(uuddss). It should be noted, however – since the novel scheme here is based on multi-nucleon interactions – that further theoretical work will be needed to assess the constraining power of any upper limits in the event that a signal is not observed.

Acknowledgements

This research of GRF has been supported by National Science Foundation Grant no. PHY-2013199 and the Simons Foundation. This work of GK was supported by Warsaw University of Technology within the Excellence Initiative: Research University (ID-UB) programme and the IDUB-POB-FWEiTE-1 project grant as well as the Polish National Science Centre under agreement no. 2022/45/B/ST2/02029, and by the Polish Ministry of Education and Science under agreement no. 2022/WK/06.

Data Availability Statement

This manuscript has no associated data or the data will not be deposited. [Authors’ comment: The simulated data of this manuscript can be obtained from the authors.]

G. Farrar and G. Kornakov contributed equally to this work.
==== Refs
References

1. Kahlhoefer F Review of lhc dark matter searches Int. J. Mod. Phys. A 2017 32 13 1730006 10.1142/S0217751X1730006X
2. Boveia A Doglioni C Dark matter searches at colliders Ann. Rev. Nucl. Part. Sci. 2018 68 1 429 459 10.1146/annurev-nucl-101917-021008
3. Ferreira EGM Ultra-light dark matter Astron. Astrophys. Rev. 2021 29 7 10.1007/s00159-021-00135-6
4. Chadha-Day F Ellis J Marsh DJ Axion dark matter: what is it and why now? Sci. Adv. 2022 8 8 3618 10.1126/sciadv.abj3618
5. Farrar, G.R.: Stable sexaquark (2017). arXiv:1708.08951 [hep-ph]
6. Farrar, G.R.: A precision test of the nature of dark matter and a probe of the qcd phase transition (2018). arXiv:1805.03723 [hep-ph]
7. Jaffe RL Perhaps a stable dihyperon Phys. Rev. Lett. 1977 38 195 10.1103/PhysRevLett.38.195
8. Farrar, G.R.: A stable sexaquark: overview and discovery strategies (2022). arXiv:2201.01334 [hep-ph]
9. Farrar GR Zaharijas G Nuclear and nucleon transitions of the h dibaryon Phys. Rev. D 2004 70 014008 10.1103/PhysRevD.70.014008
10. Shahrbaf M Blaschke D Typel S Farrar GR Alvarez-Castillo DE Sexaquark dilemma in neutron stars and its solution by quark deconfinement Phys. Rev. D 2022 105 10 103005 10.1103/PhysRevD.105.103005
11. Farrar, G.R.: The muon g-2 and lattice qcd hadronic vacuum polarization may point to new, long-lived neutral hadrons (2022). arXiv:2206.13460 [hep-ph]
12. Aoyama T The anomalous magnetic moment of the muon in the Standard Model Phys. Rep. 2020 887 1 166 10.1016/j.physrep.2020.07.006
13. Colangelo G El-Khadra AX Hoferichter M Keshavarzi A Lehner C Stoffer P Teubner T Data-driven evaluations of Euclidean windows to scrutinize hadronic vacuum polarization Phys. Lett. B 2022 833 137313 10.1016/j.physletb.2022.137313
14. Borsanyi S Fodor Z Guenther JN Leading hadronic contribution to the muon magnetic moment from lattice qcd Nature 2021 10.1038/s41586-021-03418-1
15. Blum T Update of Euclidean windows of the hadronic vacuum polarization Phys. Rev. D 2023 108 5 054507 10.1103/PhysRevD.108.054507
16. Aguillard DP Measurement of the positive muon anomalous magnetic moment to 0.20 ppm Phys. Rev. Lett. 2023 131 16 161802 10.1103/PhysRevLett.131.161802 37925710
17. F.V. Ignatov et al., Measurement of the pion formfactor with CMD-3 detector and its implication to the hadronic contribution to muon (g-2) (2023). arXiv:2309.12910 [hep-ex]
18. Ali A Lange JS Stone S Exotics: Heavy pentaquarks and tetraquarks Prog. Part. Nucl. Phys. 2017 97 123 198 10.1016/j.ppnp.2017.08.003
19. Lebed RF Mitchell RE Swanson ES Heavy-quark qcd exotica Prog. Part. Nucl. Phys. 2017 93 143 194 10.1016/j.ppnp.2016.11.003
20. Blaschke D Bravina L Bugaev L Farrar G Thermal production of sexaquarks in heavy-ion collisions Int. J. Mod. Phys. A 2021 36 25 2141005 10.1142/S0217751X21410050
21. Bruce R d’Enterria D Roeck A Drewes M Farrar GR New physics searches with heavy-ion collisions at the cern large hadron collider J. Phys. G Nucl. Part. Phys. 2020 47 6 060501 10.1088/1361-6471/ab7ff7
22. Kim BH Search for an h-dibaryon with mass near 2mλ in υ(1s) and υ(2s) decays Phys. Rev. Lett. 2013 110 22 222002 10.1103/PhysRevLett.110.222002 23767713
23. Lees JP Search for a stable six-quark state at BABAR Phys. Rev. Lett. 2019 122 7 072002 10.1103/PhysRevLett.122.072002 30848619
24. Zmeskal J Double-strangeness production with antiprotons Hyperfine Interact. 2009 194 249 10.1007/s10751-009-0082-x
25. Riedlberger J Antiproton annihilation at rest in nitrogen and deuterium gas Phys. Rev. C 1989 40 2717 10.1103/PhysRevC.40.2717
26. Cugnon J Vandermeulen J Antiproton and antilambda annihilations on several nucleons Phys. Rev. C 1989 39 181 191 10.1103/PhysRevC.39.181
27. Acharya S Study of the Λ-Λ interaction with femtoscopy correlations in pp and p-Pb collisions at the LHC Phys. Lett. B 2019 797 134822 10.1016/j.physletb.2019.134822
28. Takahashi H Observation of a (Lambda Lambda)He-6 double hypernucleus Phys. Rev. Lett. 2001 87 212502 10.1103/PhysRevLett.87.212502 11736336
29. Adam J Search for weakly decaying λn¯ and λλ exotic bound states in central pb-pb collisions at sNN=2.76 tev Phys. Lett. B 2016 752 267 277 10.1016/j.physletb.2015.11.048
30. Bertolotto L Buzzo A Debevec PT Observation of ϕϕ production in the reaction p¯p→4k± at 1.4 gev / c incident p¯ momentum Phys. Lett. B 1995 345 3 325 334 10.1016/0370-2693(94)01670-8
31. Pirner HJ Vary JP Deep Inelastic electron Scattering and the Quark Structure of He-3 Phys. Rev. Lett. 1981 46 1376 1379 10.1103/PhysRevLett.46.1376
32. Veselov A Kalashnikova Y Narodetskii I Realistic nucleon–nucleon potentials with allowance for the quark degrees of freedom Sov. J. Nucl. Phys. 1985 42 5 774 780
33. Bizzari R Three-body interactions (N¯-nn) in antiproton-deuteron annihilation at rest Lett. Nuovo Cimento 1969 2 431 10.1007/BF02755626
34. Bendiscioli G Bressani T Lavezzi L Panzarasa A Salvini P S =-2 strangeness production in p¯4he annihilations at rest Nucl. Phys. A 2007 797 109 130 10.1016/j.nuclphysa.2007.10.003
35. Agostinelli S GEANT4—a simulation toolkit Nucl. Instrum. Methods A 2003 506 250 303 10.1016/S0168-9002(03)01368-8
36. Allison J Geant4 developments and applications IEEE Trans. Nucl. Sci. 2006 53 270 10.1109/TNS.2006.869826
37. Allison J Recent developments in Geant4 Nucl. Instrum. Methods A 2016 835 186 225 10.1016/j.nima.2016.06.125
38. Brun R Rademakers F Root-an object oriented data analysis framework Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 1997 389 1 81 86 10.1016/S0168-9002(97)00048-X
39. James, F.E.: Monte Carlo Phase Space. CERN Academic Training Lecture. CERN, Geneva (1968). 10.5170/CERN-1968-015. CERN, Geneva, 1 May 1968. https://cds.cern.ch/record/275743
40. Kornakov, G., Cerchiari, G., Zieliński, J., Lappo, L., Sadowski, G., Doser, M.: Synthesis of cold and trappable fully stripped HCI’s via antiproton-induced nuclear fragmentation in traps (2022). arXiv:2209.02596 [physics.atom-ph]
41. Galoyan A Uzhinsky V Ribon A Simulation of antimatter-matter interactions in geant4 EPJ Web Conf. 2018 173 06005 10.1051/epjconf/201817306005
42. G. Cowan, K. Cranmer, E. Gross, O. Vitells, Asymptotic formulae for likelihood-based tests of new physics. Eur. Phys. J. C 71, 1554 (2011). 10.1140/epjc/s10052-011-1554-0 (Erratum: Eur.Phys.J.C 73, 2501 (2013)). arXiv:1007.1727 [physics.data-an]
43. B. Autin et al., The anti-proton decelerator (ad), a simplified anti-proton source (feasibility study) (1995)
44. Bartmann W Belochitskii P Breuker H Butin F Carli C Eriksson T Oelert W Ostojic R Pasinelli S Tranquille G The ELENA facility Philos. Trans. R. Soc. Lond. A 2018 376 2116 20170266 10.1098/rsta.2017.0266
45. Montagna P Single and multinucleon antiproton-4he annihilation at rest Nucl. Phys. A 2002 700 159 192 10.1016/S0375-9474(01)01302-1
46. Doser M Antiprotonic bound systems Prog. Part. Nucl. Phys. 2022 125 103964 10.1016/j.ppnp.2022.103964
47. C. Amsler, M. Antonello, A. Belov et al., Pulsed production of antihydrogen. Commun. Phys. 4, 19 (2021)
48. Andersen T Andersen LH Balling P Haugen HK Hvelplund P Smith WW Taulbjerg K Metastable-ion lifetime studies utilizing a heavy-ion storage ring: measurements on he- Phys. Rev. A 1993 47 890 896 10.1103/PhysRevA.47.890 9909007
49. Amoretti M Amsler C Bonomi G Production and detection of cold antihydrogen atoms Nature 2002 419 456 459 10.1038/nature01096 12368849
50. Gabrielse G Bowden NS Oxley P Speck A Storry CH Tan JN Wessels M Grzonka D Oelert W Schepers G Sefzick T Walz J Pittner H Hänsch TW Hessels EA Background-free observation of cold antihydrogen with field-ionization analysis of its states Phys. Rev. Lett. 2002 89 213401 10.1103/PhysRevLett.89.213401 12443407
