## Searches for ultra-high energy gamma-ray at the Pierre Auger Observatory and implications on super-heavy dark matter

## Olivier Deligny , for the Pierre Auger Collaboration 1 2

1

Laboratoire de Physique des 2 Infinis Ir` ene Joliot-Curie (IJCLab) CNRS/IN2P3, Universit´ e Paris-Saclay, Orsay, France 2

Full author list: http://www.auger.org/archive/authors 2024 06.html

E-mail: deligny@ijclab.in2p3.fr

Abstract. The first interactions of photon-induced showers are of electromagnetic nature, and the transfer of energy to the hadron/muon channel is reduced with respect to the bulk of hadron-induced showers. This results in a lower number of secondary muons. Additionally, as the development of photon showers is delayed by the typically small multiplicity of electromagnetic interactions, their maximum of shower development is deeper in the atmosphere than for showers initiated by hadrons. These salient features have enabled searches for photon showers at the Pierre Auger Observatory. They have led to stringent upper limits on ultra-high-energy gamma-ray fluxes over four orders in magnitude in energy. These limits are not only of considerable astrophysical interest, but they also allow us to constrain beyond-standard-physics scenarios. For instance, dark matter particles could be superheavy, provided their lifetime is much longer than the age of the universe. Constraints on specific extensions of the Standard Model of particle physics that meet the lifetime requirement for a superheavy particle will be presented. They include limits on instanton strength as well as on mixing angle between active and sterile neutrinos.

§ 1 Searches for ultra-high energy gamma-rays from dark-matter particles. Compelling evidence for the observation of the decay of super-heavy dark matter particles would be the detection of a flux of gamma-rays with energies in excess of 10 8 GeV, in particular from regions of higher dark-matter density such as the center of our Galaxy. The identification of gamma-ray primaries relies on distinct salient features between the extensive air showers initiated by gamma-rays and those by the overwhelming background of nuclei: a lower number of secondary muons and deeper 𝑋 max values. Both the ground signal and 𝑋 max can be measured at the Pierre Auger Observatory [1] where a hybrid detection technique is employed for the observation of extensive air showers by combining fluorescence detectors with ground arrays of particle detectors. The combination of the various instruments allows ultra-high-energy gamma-rays to be detected in the energy range above 10 7 7 . GeV with a directional exposure obtained through the time and area integration of the detection efficiency 𝜖 𝛾 and selection efficiency 𝜅 𝛾 projected onto the direction perpendicular to the arrival direction,

$$\mathcal { E } _ { \gamma } ( E, \delta, \alpha ) & = \int \mathrm d { \mathbf x } \int \mathrm d { t } \cos \theta \epsilon _ { \gamma } ( E, \theta, { \mathbf x }, t ) \kappa _ { \gamma } ( E, \theta ).$$

In total, four different analyses, differing in the detector used, have been developed to cover the wide energy range probed at the Observatory [2, 3, 4, 5]. In particular, two of these analyses benefit from a direct estimate of the depth of shower maximum by the fluorescence detector as one of the discriminating variables. The use of the fluorescence-detector datasets introduces for these analyses an explicit dependence in time for 𝜖 𝛾 due to the limited

Figure 1: Directional exposure of the Pierre Auger Observatory in Galactic coordinates to ultra-high energy gamma-rays at 10 8 GeV ( left ) and 10 10 GeV ( right ). From [6].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00795v1-figure-7f16013f-70f3-47c4-9997-754f97bd17c7.jpg)

duty cycle on moonless nights that propagates into a small dependence in right ascension 𝛼 for E 𝛾 . In addition to the detection efficiency, the 𝜅 𝛾 factor accounts for the dependencies of the selection process aimed at separating hadronic-induced showers from gamma-ray-induced ones. For illustration, we show in Fig. 1 an example of directional exposure to photons at 10 8 GeV (left panel) and at 10 10 GeV (right panel). The former is built on the 1.14 km 2 array with a separation of detectors of 433 m optimized to study the range of energies around 10 8 GeV, while the latter is the 3,000 km 2 one, optimized for higher energies. This explains the large increase of exposure observed at high energies. Overall, no gamma-ray could be firmly identified up to now. The non-observation of candidates (beyond the expected background) allowed the derivation of upper bounds that can constrain various models very effectively, as discussed below.

Due to their attenuation over intergalactic distances, only ultra-high energy gamma-rays emitted in the Milky Way can survive on their way to Earth. The emission rate per unit volume and unit energy 𝑞 𝛾 from decaying dark-matter particles is shaped by the energy density of dark matter 𝜌 DM,

$$q _ { \gamma } ( E, \mathbf x _ { \odot } + s \mathbf u ) = \frac { 1 } { M _ { X } \tau _ { X } } \frac { d N _ { \gamma } } { d E } \rho _ { \text{DM} } ( \mathbf x _ { \odot } + s \mathbf u ),$$

where 𝑀 𝑋 is the mass of the dark-matter particle, 𝜏 𝑋 its lifetime, 𝑑𝑁 𝛾 / 𝑑𝐸 is the energy spectrum of the gamma-ray decay byproducts, x ⊙ is the position of the Solar system in the Galaxy, and u is a unit vector on the sphere. The energy density is normalized to 𝜌 ⊙ = 0 45 GeV cm . -3 [7]. The flux (per steradian) of ultra-high energy gamma-rays produced by the decay of dark-matter particles, 𝐽 DM ,𝛾 ( 𝐸, u ) , is then obtained by integrating the position-dependent emission rate 𝑞 𝛾 along the path of the photons in the direction u ,

$$J _ { \gamma } ( E, \mathbf u ) = \frac { 1 } { 4 \pi } \int _ { 0 } ^ { \infty } d s \ q _ { \gamma } ( E, \mathbf x _ { \odot } + s \mathbf u ),$$

where the 4 𝜋 normalization factor accounts for the isotropy of the decay processes. Finally, the expected number of events above a threshold energy follows from

$$n _ { \gamma } ( E ) = \int d \mathfrak { u } \int _ { > E } d E ^ { \prime } \, \mathcal { E } _ { \gamma } ( E ^ { \prime }, \mathfrak { n } ) J _ { \gamma } ( E ^ { \prime }, \mathfrak { u } ).$$

Assuming that the relic abundance of dark matter is saturated by super-heavy particles, constraints can be inferred in the plane ( 𝜏 𝑋 , 𝑀 𝑋 ) by requiring the number of gamma-rays expected from Eq. 4 to be less than the upper limits on the number of observed events accounting for the expected background. For generic two-body decay channels, the constraints generally lead to lifetime much longer than the age of the universe for 𝑀 &gt; 𝑋 10 9 GeV.

§ 2 Lifetime constraints in the perturbative domain. To comply with the lifetime constraints, some models postulate the existence of super-weak couplings between the dark and standard-model sectors. The lifetime 𝜏 𝑋 of the particles is then governed by the strength of the couplings 𝑔 𝑋 Θ (or reduced couplings 𝛼 𝑋 Θ = 𝑔 2 𝑋 Θ /( 4 𝜋 ) ) and by the mass dimension 𝑛 of the operator Θ standing for the standard model fields in the effective interaction [8]. Even without knowing the theory behind the decay of the dark-matter particle, generic constraints on 𝛼 𝑋 Θ and 𝑛 can be derived [9]. The effective interaction term that couples the field 𝑋 associated with the heavy particle to the standard-model fields is taken as

$$\mathcal { L } _ { \text{int} } = \frac { g _ { X \Theta } } { \Lambda ^ { n - 4 } } X \Theta,$$

a

X

Figure 2: Exclusion regions in the plane ( 𝛼 𝑋 Θ , 𝑀 𝑋 ) for several values of mass dimension 𝑛 of operators responsible for the perturbative decay of the super-heavy particle, and for an energy scale of new physics Λ= 10 16 GeV. From [9].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00795v1-figure-f30de431-e10f-433b-b564-73e319ec3729.jpg)

where Λ is an energy parameter typical of the scale of the new interaction. In the absence of further details about the operator Θ , the matrix element describing the decay transition is considered flat in all kinematic variables so that it behaves as |M| 2 ∼ 4 𝜋𝛼 𝑋 Θ / Λ 2 𝑛 -4 . On the basis of dimensional arguments, the lifetime of the particle 𝑋 is then given as

$$\tau _ { X \Theta } = \frac { V _ { n } } { 4 \pi M _ { X } \alpha _ { X \Theta } } \left ( \frac { \Lambda } { M _ { X } } \right ) ^ { 2 n - 8 },$$

where 𝑉 𝑛 is a phase space factor. As a proxy for this factor, we use the expression derived for 𝑛 -1 particles in the final state [10],

$$V _ { n } = \left ( \frac { 2 } { \pi } \right ) ^ { n - 1 } \Gamma ( n - 1 ) \Gamma ( n - 2 ),$$

with Γ ( 𝑥 ) the Euler gamma function. It is apparent that the coupling constant 𝛼 𝑋 Θ and the dimension 𝑛 have to take specific values for super-heavy particles to be stable enough. In practice, for a specific upper limit on 𝑛 𝛾 ( 𝐸 ) at one energy threshold 𝐸 , a lower limit of the 𝜏 𝑋 Θ parameter is derived for each value of mass 𝑀 𝑋 . The lower limit on 𝜏 𝑋 Θ is subsequently transformed into an upper limit on 𝛼 𝑋 Θ by means of Eq. (6). By repeating the procedure for each upper limit on 𝑛 𝛾 ( 𝐸 ) , a set of curves is obtained, reflecting the sensitivity of a specific energy threshold to some range of mass. The union of the excluded regions finally provides the constraints in the plane ( 𝛼 𝑋 Θ , 𝑀 𝑋 ) . In this manner we obtain the contour lines shown in Fig. 2 for several values of 𝑛 and for an emblematic choice of GUT Λ value. The scale chosen for 𝛼 𝑋 Θ ranges from 1 down to 10 -5 . It is observed that for the limits on photon fluxes to be satisfied, the mass of the super-heavy particle cannot exceed ≳ 10 9 GeV ( ≳ 10 11 GeV)for operators of dimension equal to or larger than 𝑛 = 8 ( 𝑛 = 10), while larger masses require an increase in 𝑛 . To approach the large masses while keeping operators of dimension relatively low, 'astronomically-small' coupling constants should be at work. The same conclusions hold for other choices of Λ . All in all, for perturbative processes to be responsible for the decay of super-heavy dark matter particles requires quite 'unnatural' fine-tuning.

§ 3 Constraints on instanton-induced decays. Stability of super-heavy dark matter particles is consequently calling for a new quantum number to be conserved in the dark sector so as to protect the particles from decaying. Nevertheless, even stable particles in the perturbative domain will in general eventually decay due to non-perturbative effects (instantons) in non-abelian gauge theories. Instanton-induced decay can thus make observable a dark sector that would otherwise be totally hidden by the conservation of a quantum number [11]. Assuming quarks and leptons carry this quantum number and so contribute to anomaly relationships with contributions from the dark sector, they will be secondary products in the decays of dark matter together with the lightest hidden fermion. The lifetime of the decaying particle is mainly driven by the instanton suppression factor that leads to [12]

$$\tau _ { X } \simeq M _ { X } ^ { - 1 } \exp { ( 4 \pi / \alpha _ { X } ) },$$

with 𝛼 𝑋 the reduced coupling constant of the hidden gauge interaction.

a

X

Figure 3: Upper limits at 95% C.L. on the coupling constant 𝛼 𝑋 of a hidden gauge interaction as a function of the mass 𝑀 𝑋 of a dark matter particle 𝑋 (assumed to compose 100% of the observed dark matter abundance) decaying into a dozen of 𝑞𝑞 ¯ pairs. From [13].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00795v1-figure-dc61b69d-0e01-499b-8152-d4b5293e239b.jpg)

Quite independently of the hidden gauge interaction, the exact content in instanton-induced decays of quarks and leptons, which will eventually produce hadrons decaying into photons and neutrinos, obeys selection rules that involve very large multiplicities. Assuming ten pairs of quarks and leptons in the final state, the energy spectrum 𝑑𝑁 𝛾 ( 𝐸, 𝑀 𝑋 )/ 𝑑𝐸 then follows from fragmentation effects [14, 15, 16, 17].

Eq. (8) provides us with a relationship connecting the lifetime 𝜏 𝑋 to the coupling constant 𝛼 𝑋 . In the same way as in the perturbative case above, upper limits on 𝛼 𝑋 can be obtained. They are shown as the shaded red area in Fig. 3. The coupling should be less than ≃ 0 09 for a wide range of masses. . Numerical factors in front of the exponential could however arise in Eq. (8) depending on the underlying model for the hidden gauge sector. Explicit constructions of the dark sector are, however, well beyond the scope of this contribution. Although the limits presented in Fig. 3 are hardly destabilized due to the exponential dependence in 𝛼 -1 𝑋 , we note that model-dependent numerical factors in front of the exponential such as 10 ± 𝑘 lead to a shift of ± 0 0013 . 𝑘 in 𝛼 𝑋 limits. We show in dotted and dashed lines the bounds for 𝑘 = 2 and 𝑘 = 4, respectively. These factors are by far the dominant systematic uncertainties.

§ 4 Constraints on dark matter coupled to sterile neutrinos. A superheavy particle that is metastable can also result from the coupling between a pseudo-scalar particle with sterile neutrinos embedded in an extended-seesaw framework [18]. In this beyond-standard-model extension, the dark-matter particle 𝑋 interacts only with sub-eV and superheavy (10 12 -14 GeV) sterile neutrinos, of masses 𝑚 𝑁 and 𝑀 𝑁 respectively, via Yukawa couplings 𝑦 𝑚 and 𝑦 𝑀 . In the mass-eigenstate basis, neutrinos 𝜈 1 and 𝜈 2 are then quasi-sterile or quasi-active respectively, depending on the mixture of active and sterile neutrinos governed by a small mixing angle 𝜃 𝑚 ≃ √ 2 𝑦 𝑚 𝑣 / 𝑚 𝜈 , with 𝑣 the electroweak scale (mass dimension) and 𝑚 𝜈 the mass of the known neutrinos. To leading order in 𝑦 𝑚 , quasi-active neutrinos are produced from quasi-sterile ones subsequent to the decay of 𝑋 . Consequently, the coupling 𝑦 𝑚 controls the dominant decay channels and allows for trading a factor ( 𝑀 𝑋 / 𝑀 P ) 2 (with 𝑀 P the Planck mass) for a ( 𝑚 𝜃 𝜈 𝑚 / ) 𝑣 2 one in the decay width of 𝑋 . This trading enables the reduction of the width by a factor ∼ 10 -25 𝜃 2 𝑚 for a benchmark value 𝑀 𝑋 = 10 9 GeV, leading to the required lifetimes well beyond the age of the universe.

To leading order in 𝜃 𝑚 , the total width Γ 𝑋 is dominated by three-body channels stemming from the interaction between active/sterile neutrinos and the Higgs isodoublet with Yukawa coupling 𝑦 𝑚 ≃ √ 2 𝜃 𝑚 𝑚 2 / 𝑣 . The channel

Figure 4: Top : Constraints on 𝜏 𝑋 as a function of 𝑀 𝑋 for a value of the couplings of sterile neutrinos with the inflationary sector 𝜆 𝑁 𝑚 = 10 -5 . The dotted line illustrates the constraints stemming from neutrino secondaries alone. Bottom : Constraints on 𝜃 𝑚 as a function of 𝑀 𝑋 for three different values of the coupling constant 𝛼 𝑋 , still for 𝜆 𝑁 𝑚 = 10 -5 . The hatched-red region 𝜃 𝑚 ≥ 9 × 10 -4 is excluded from the constraint on Δ 𝑁 eff (see text). From [6].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00795v1-figure-d17798e1-c813-4eb7-82e9-8f95d1552eda.jpg)

𝑋 → ℎ𝜈 𝜈 1 2, diagrammatically represented as

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00795v1-figure-56e2c473-78ac-49ea-a6ec-279364832317.jpg)

gives the most important contribution to the width [18]:

$$\Gamma _ { h \nu _ { 1 } \nu _ { 2 } } ^ { X } = \frac { \alpha _ { X } ^ { 2 } \theta _ { m } ^ { 2 } } { 1 9 2 \pi ^ { 3 } } \left ( \frac { M _ { X } } { M _ { P } } \right ) ^ { 2 } \left ( \frac { m _ { 2 } } { v } \right ) ^ { 2 } M _ { X }.$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00795v1-figure-d4e52050-76f3-40f9-9652-6fb1e66092e9.jpg)

X

Figure 5: Left: Constraints in the ( 𝐻 inf , 𝑀 𝑋 ) plane, where viable values are delineated by the blue lines for different values of reheating efficiency 𝜖 . Additional constraints from the non-observation of instanton-induced decay of super-heavy particles allow for excluding the mass ranges in the red-shaded regions, for the specified value of the dark-sector gauge coupling. From [13]. Right: Same, adding the possibility of a radiative production of dark matter in the inflaton decay.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00795v1-figure-dd5c9cc9-cda3-4964-8626-ac28c3946e07.jpg)

As announced, a factor ( 𝑀 𝑋 / 𝑀 P ) 2 expected from dimensional arguments is indeed traded for a factor ( 𝑚 2 / ) 𝑣 2 . Constraints in the planes ( 𝜏 𝑋 , 𝑀 𝑋 ) and ( 𝜃 𝑚 , 𝑀 𝑋 ) can be derived from the non-observation of ultra-high energy gamma-rays (and neutrinos) at the Observatory. First, 90% C.L. lower limits on the lifetime 𝜏 𝑋 are obtained as previously by setting, for a specific value of 𝑀 𝑋 , 𝑛 𝛾 ( 𝐸 ) or 𝑛 𝜈 ( 𝐸 ) to the 90% C.L. upper-limit numbers corresponding to the number of background-event candidates in the absence of signal. Subsequently, a scan in 𝑀 𝑋 is carried out. It leads to a curve in the plane ( 𝜏 𝑋 , 𝑀 𝑋 ) that pertains to the energy threshold 𝐸 considered. By repeating the procedure for several thresholds, a set of curves is obtained, reflecting the sensitivity of a specific energy threshold to some range of mass 𝑀 𝑋 . The union of the excluded regions finally provides the constraints in the ( 𝜏 𝑋 , 𝑀 𝑋 ) plane. Results are shown in Fig. 4 (top panel); lifetimes within the cross-hatched region are excluded. The region in full red pertains to a particular value of a Yukawa coupling 𝜆 𝑁 𝑚 = 10 -5 , the meaning of which will be explained below. To illustrate the contribution from each secondary at our disposal, we show as the dotted the contribution to the constraints stemming from neutrinos alone; an analysis of the IceCube exposure dedicated to the benchmark-scenario decay channels would likely provide better sensitivity for exploring masses 𝑀 𝑋 ≲ 10 8 5 . GeV. The lower limit on 𝜏 𝑋 is then transformed into an upper limit on 𝜃 𝑚 using the expressions of the total width of the particle 𝑋 . Results are shown in the bottom panel for separate values of 𝛼 𝑋 : color-coded regions pertain to 𝜆 𝑁 𝑚 = 10 -5 while their extension (in cross-hatched) would require smaller values of 𝜆 𝑁 𝑚 . Systematic uncertainties on 𝜃 𝑚 constraints amount overall to ≃ ± 15%; they are dominated by those on the neutrino exposure [19]. The restricted ranges of 𝑀 𝑋 for different 𝛼 𝑋 values come from the requirements not to overclose the universe with dark matter, while the exclusion hatched band comes from not altering the expansion history of the universe with the presence of ultra-light species such as sterile neutrinos 𝑁 𝑚 . We now briefly explain how these constraints are obtained.

§ 5 Cosmological constraints. Gravitational interaction alone may have been sufficient to produce the right amount of super-heavy dark matter particles at the end of the inflation era for a wide range of high masses, up to 𝑀 GUT, accounting for the production by annihilation of standard model particles (SM) [20] or of inflaton particles ( 𝜙 hereafter) [21] through the exchange of a graviton. In this scenario, the relic abundance of super-heavy dark matter particles can be estimated from the quite involved reheating dynamics [22, 23]. The time evolution of the 𝑋 -particle density 𝑛 𝑋 reads as

$$\mathbb { m } \, \infty \, \infty \, \frac { d n _ { X } ( t ) } { d t } + 3 H ( t ) n _ { X } ( t ) & \simeq \sum _ { i } \bar { n } _ { i } ^ { 2 } ( t ) \gamma _ { i },$$

where the sum in the right hand side stands for the contributions from the standard model [20] and inflationary [21] sectors to produce fermions. Introducing the dimensionless abundance 𝑌 𝑋 = 𝑛 𝑋 𝑎 3 / 𝑇 3 rh to absorb the expansion of the universe, with 𝑇 rh the reheating temperature, and using 𝑎𝐻 𝑎 𝑑𝑡 ( ) = 𝑑𝑎 from the definition of the Hubble parameter (with 𝑎 the scale factor), Eq. 11 becomes

$$\frac { d Y _ { X } ( a ) } { d a } \simeq \frac { a ^ { 2 } } { T _ { \text{th} } ^ { 3 } H ( a ) } \sum _ { i } \bar { n } _ { i } ^ { 2 } ( a ) \gamma _ { i },$$

which, using the dynamics of the expansion rate during reheating, yields the present-day dimensionless abundance 𝑌 𝑋, 0 assuming 𝑌 𝑋, inf = 0. The present-day relic abundance, Ω CDM, can then be related to 𝑀 𝑋 , 𝐻 inf , and 𝜖 = 𝑇 rh /( 0 25 . √ 𝑀 𝐻 P inf ) through [20] Ω CDM ℎ 2 = 9 2 10 24 . 𝜖 4 𝑀 𝑌 𝑋 𝑋, 0 / 𝑀 P. As a result, one interesting viable possibility in the ( 𝐻 inf , 𝑀 𝑋 ) parameter space is that 𝑋 particles with masses as large as the GUT energy scale could be sufficiently abundant to match the dark matter relic density, provided that the inflationary energy scale is high ( 𝐻 inf ∼ 10 13 GeV) and the reheating efficiency is high (so that reheating is quasi-instantaneous). This rules out values of the dark-sector gauge coupling greater than ≃ 0 085, as observed in the left panel of Fig. 5. . The mass values could however be smaller if the reheating temperature is not that high. In general, for high efficiencies 𝜖 (corresponding to short duration of the reheating era), the SM + SM → 𝑋 + 𝑋 reaction allows for a wide range of 𝑀 𝑋 values. For efficiencies below ≃ 0 01, the . 𝜙 + 𝜙 → 𝑋 + 𝑋 reaction allows for solutions in a narrower range of the ( 𝐻 inf , 𝑀 𝑋 ) plane close to 𝑀 𝑋 = 10 13 GeV, with in particular 𝑀 𝑋 ≤ 𝑀 𝜙 as a result of the kinematic suppression in the corresponding rate 𝛾 𝑖 [21].

Similar constraints can be drawn for the model invoking couplings of the 𝑋 particle with sterile neutrinos. In addition to its couplings to the dark-matter sector and to the standard-model one through the Higgs isodoublet, the sterile neutrino 𝑁 𝑚 is also coupled to an inflationary sector in the benchmark [18]. This coupling, governed by a unique Yukawa parameter 𝜆 𝑁 𝑚 for every 𝜈 1 neutrinos, yields to a 'radiative' production of 𝑋 via a diagram similar to that depicted in Eq. 9 (substituting 𝑋 by the inflaton Φ in the initial state, and ℎ and 𝜈 2 by 𝑋 and 𝜈 1 in the final states). Such a mechanism leads to a direct production of dark matter during the reheating period that can be sufficient, in general, to match the right amount of dark matter observed today [24]. In the benchmark [18], values for 𝜆 𝑁 𝑚 are then required to range preferentially around 10 -5 . To infer the dark matter density 𝑛 𝑋 produced mainly during the reheating epoch, we also consider the minimal setup of gravitational production of 𝑋 particles through the annihilation of standard-model (inflaton) particles as in [20] (as in [21]). In these conditions, 𝑋 particles can be produced as long as the collision rate of particles is larger than the expansion rate 𝐻 and/or as long as the inflaton field oscillates. By contrast, 𝑛 𝑋 is prohibitively low to allow any thermal equilibrium for dark matter. The collision term in the Boltzmann equation is then approximated as a source term only. Overall, the Boltzmann equation reads as

$$\frac { d n _ { X } ( t ) } { d t } + 3 H ( t ) n _ { X } ( t ) = \sum _ { i } \bar { n } _ { i } ^ { 2 } ( t ) \gamma _ { i } + \bar { n } _ { \phi } ( t ) \Gamma _ { X \nu _ { 1 } \nu _ { 1 } }.$$

As a result, viable couples of values for ( 𝐻 inf , 𝑀 𝑋 ) scale as 𝐻 inf ∝ 𝑀 2 𝑋 up to a maximum value for 𝑀 𝑋 , which depend on 𝜖 and 𝛼 𝑋 - see right panel of Fig. 5. This scaling is a consequence of the domination of the radiativeproduction process over the gravitational one as long as the allowed values of 𝐻 inf are too small for a given 𝑀 𝑋 value to generate significant particle production by gravitational interactions. For larger masses, the contribution from gravitational interactions added to the radiative production of 𝑋 leads to an overproduction of dark matter that overcloses the universe, and there is thus no longer solution. This explains why the color-coded regions extend up to some maximum values of 𝑀 𝑋 in Fig. 4, for a benchmark value of 𝜆 𝑁 𝑚 = 10 -5 . To the right of the regions shown in cross-hatched, 𝜆 𝑁 𝑚 would need to be smaller.

§ 6 Conclusions. We have shown that the data of the Pierre Auger Observatory provide stringent constraints on decaying super-heavy particles that would constitute dark matter. These constraints allow us to put a limit on instanton strength or on the angle 𝜃 𝑚 mixing sub-eV sterile and active neutrinos in the context of an extension to the SM that couples the sterile neutrinos to a superheavy DM candidate. Other physics providing mechanisms to stabilize super-heavy particles can be explored, such as supersymmetry broken at high scale with a tiny 𝑅 -parity violation. This will be the subject of a future study.

## References

- [1] Aab A et al. (Pierre Auger) 2015 Nucl. Instrum. Meth. A 798 172-213 ( Preprint 1502.01323 )
- [2] Abdul Halim A et al. (Pierre Auger) 2023 PoS ICRC2023 238
- [3] Abreu P et al. (Pierre Auger) 2022 Astrophys. J. 933 125 ( Preprint 2205.14864 )
- [4] Abdul Halim A et al. (Pierre Auger) 2024 ( Preprint 2406.07439 )
- [5] Abreu P et al. (Pierre Auger) 2023 JCAP 05 021 ( Preprint 2209.05926 )
- [6] Abdul Halim A et al. (Pierre Auger) 2024 Phys. Rev. D 109 L081101 ( Preprint 2311.14541 )
- [7] Jiao Y, Hammer F, Wang H, Wang J, Amram P, Chemin L and Yang Y 2023 Astron. Astrophys. 678 A208 ( Preprint 2309.00048 )
- [8] de Vega H J and Sanchez N G 2003 Phys. Rev. D 67 125019
- [9] Abreu P et al. (Pierre Auger) 2023 Phys. Rev. D 107 042002 ( Preprint 2208.02353 )
- [10] Kleiss R, Stirling W J and Ellis S D 1986 Comput. Phys. Commun. 40 359
- [11] Kuzmin V A and Rubakov V A 1998 Phys. Atom. Nucl. 61 1028 ( Preprint astro-ph/9709187 )
- [12] 't Hooft G 1976 Phys. Rev. Lett. 37 8-11
- [13] Abreu P et al. (Pierre Auger) 2023 Phys. Rev. Lett. 130 061001 ( Preprint 2203.08854 )
- [14] Aloisio R, Berezinsky V and Kachelriess M 2004 Phys. Rev. D 69 094023 ( Preprint hep-ph/0307279 )
- [15] Sarkar S and Toldra R 2002 Nucl. Phys. B 621 495-520 ( Preprint hep-ph/0108098 )
- [16] Barbot C and Drees M 2003 Astropart. Phys. 20 5-44 ( Preprint hep-ph/0211406 )
- [17] Bauer C W, Rodd N L and Webber B R 2021 JHEP 06 121 ( Preprint 2007.15001 )
- [18] Dudas E, Heurtier L, Mambrini Y, Olive K A and Pierre M 2020 Phys. Rev. D 101 115029 ( Preprint 2003.02846 )
- [19] Aab A et al. (Pierre Auger) 2019 JCAP 10 022 ( Preprint 1906.07422 )
- [20] Garny M, Sandora M and Sloth M S 2016 Phys. Rev. Lett. 116 101302 ( Preprint 1511.03278 )
- [21] Mambrini Y and Olive K A 2021 Phys. Rev. D 103 115009 ( Preprint 2102.06214 )
- [22] Chung D J H, Kolb E W and Riotto A 1999 Phys. Rev. D 60 063504 ( Preprint hep-ph/9809453 )
- [23] Giudice G F, Kolb E W and Riotto A 2001 Phys. Rev. D 64 023508 ( Preprint hep-ph/0005123 )
- [24] Kaneta K, Mambrini Y and Olive K A 2019 Phys. Rev. D 99 063508 ( Preprint 1901.04449 )
Olivier Deligny
Searches for ultra-high energy gamma-ray at the Pierre Auger Observatory and implications on super-heavy dark matter
The first interactions of photon-induced showers are of electromagnetic
nature, and the transfer of energy to the hadron/muon channel is reduced with
respect to the bulk of hadron-induced showers. This results in a lower number
of secondary muons. Additionally, as the development of photon showers is
delayed by the typically small multiplicity of electromagnetic interactions,
their maximum of shower development is deeper in the atmosphere than for
showers initiated by hadrons. These salient features have enabled searches for
photon showers at the Pierre Auger Observatory. They have led to stringent
upper limits on ultra-high-energy gamma-ray fluxes over four orders in
magnitude in energy. These limits are not only of considerable astrophysical
interest, but they also allow us to constrain beyond-standard-physics
scenarios. For instance, dark matter particles could be superheavy, provided
their lifetime is much longer than the age of the universe. Constraints on
specific extensions of the Standard Model of particle physics that meet the
lifetime requirement for a superheavy particle will be presented. They include
limits on instanton strength as well as on mixing angle between active and
sterile neutrinos.