Multi-Modal Clustering Events observed by Horizon-10T and Axion Quark Nuggets

The Horizon-10T collaboration \cite{Beisembaev:2016cyg,2017EPJWC.14514001B,Beznosko:2019cI,2019EPJWC.20806002B,Beisembaev:2019nzd} have reported observation of Multi-Modal Events (MME) containing multiple peaks suggesting their clustering origin. These events are proven to be hard to explain in terms of conventional cosmic rays (CR). We propose that these MMEs might be result of the dark matter annihilation events within the so-called axion quark nugget (AQN) dark matter model, which was originally invented for completely different purpose to explain the observed similarity between the dark and the visible components in the Universe, i.e. $\Omega_{\rm DM}\sim \Omega_{\rm visible}$ without any fitting parameters. We support this proposal by demonstrating that the observations \cite{Beisembaev:2016cyg,2017EPJWC.14514001B,Beznosko:2019cI,2019EPJWC.20806002B,Beisembaev:2019nzd}, including the frequency of appearance, intensity, the spatial distribution, the time duration, the clustering features, and many other properties nicely match the emission characteristics of the AQN annihilation events in atmosphere. We list a number of features of the AQN events which are very distinct from conventional CR air showers. The observation (non-observation) of these features may substantiate (refute) our proposal.


I. INTRODUCTION
In this work we discuss two naively unrelated stories. The first one is the study of a specific dark matter (DM) model, the so-called axion quark nugget (AQN) model [6], see a brief overview of this model below. The second one deals with the recent puzzling observations [1][2][3][4][5] by the Horizon-10T (H10T) collaboration of the Multi-Modal Events (MME). We overview the corresponding events below in details. We also highlight some difficulties in interpretation of these events in terms of the standard CR air showers. The unusual features of MME include: 1. "clustering puzzle": Two or more peaks separated by ∼ 10 2 ns are present in several detection points, while entire event may last ∼ 10 3 ns. It can be viewed as many fronts separated by ∼ (10 2 − 10 3 ) ns, instead of a single front; 2. "particle density puzzle": The number density of particles recorded at different detection points apparently weakly dependent on distance from Extensive Air Showers (EAS) axis; 3. "pulse width puzzle": The width of each individual pulse is around (20 − 35) ns and apparently does not depend on distance from EAS axis; 4. "intensity puzzle": The observed intensity of the events (measured in units of a number of particles per unit area) is of order ρ ∼ (100 − 300)m −2 when measured at distances (300 − 800) m from EAS axis. Such intensity would correspond to the CR energy of the primary particle on the level E p > ∼ 10 19 eV which would have dramatically different event rate in comparison with observed rate recorded by H10T.
Before we proceed with our explanation of the proposal to view the MMEs as the AQN events one should mention that similar unusual features of the CR air showers have been noticed long ago for the first time by Jelley and Whitehouse [7] in 1953. Later, EAS exhibiting the unusual time structures were studied by several independent experiments, see e.g. [8]. We refer to ref. [2] for overview and references of the older observations and studies of such unusual events.
The considerable recent advances in rising the resolution (on the level of a few ns) has allowed the H10T collaboration dramatically improve the collection and analysis of the MMEs. In particular, during ∼ 3500 hours of operation the H10T collaboration collected more than 10 3 MMEs.
In what follows we overview the observation [1][2][3][4][5] by emphasizing the puzzling nature of these events if interpreted as conventional EAS events. At the same time the same observations can be explained in very natural way if interpreted in terms of the AQN annihilation events as we shall argue in this work. One should mention here that a similar conclusion has been also reached for a different type of unusual CR-like events. First, it has been argued in [9,10] that the Telescope Array (TA) "mysterious bursts" [11,12] can be naturally interpreted as the AQN events. Secondly, it has been also argued in [13] that the Antarctic Impulse Transient Antenna (ANITA) observation [14][15][16] of two anomalous events with noninverted polarity might be also a consequence of the AQN annihilation events. Important comment here that in all those cases the basic parameters of the AQN model remain the same as they have been fixed long ago from dramatically different observations in a very different context.
Our presentation is organized as follows. In next Section II we explain why the observation [1][2][3][4][5] listed as 1-4 puzzles are very mysterious events if interpreted as conventional EAS events. In Section III we give a brief overview of the AQN framework with emphasize on the key elements relevant for the present studies. In Section IV we formulate our proposal on identification of the unusual Multi-Modal Events with the upward moving AQN events. In Section V we estimate the event rate. Our main Section is VI where we estimate a variety of relevant time scales and the intensity of the events. In the same section we also confront our proposal with observations and argue that puzzles 1-4 can be naturally understood within the AQN framework. Section VII is our conclusion where we suggest several tests which may substantiate or refute our proposal on identification of the Multi-Modal Cluster Events with the upward moving AQNs.
II. CONVENTIONAL CR PICTURE CONFRONTS THE MME OBSERVATIONS [1][2][3][4] FIG. 1. top: Aerial view and geometry of H10T instrument with location of 10 detectors, adopted from [3]; bottom: A typical MME event recorded on April 6.2018 by H10T instrument at point #9, adopted from [3]. All pulses are recorded at a single detection point. Delay times, the width of each peak τ0.8 in ns, and the particle density ρ0.8 per m −2 within τ0.8 are also shown in the table. .
This section is devoted to the first part of our story where we describe the MME observations and argue why the observed events are inconsistent with conventional CR interpretation, while the next section III is devoted to the second part of our story, the AQN dark matter model.
We start by reviewing the conventional picture of the EAS. It is normally assumed that EAS can be thought as a disk -pancake with well defined EAS axis. It is also assumed that EAS represents an uniform, without any breaks structure. The standard picture also assumes that the particle density drops smoothly with distance when moving away from the core, while the thickness of the EAS pancake increases with the distance from the core. Now we want to see why this conventional picture is in dramatic contradiction with observations of the MME events.
1. Indeed, a typical MME is shown on Fig. 1 where a complicated temporal features are explicitly seen. Several peaks separated by ∼ 10 2 ns in a single detector represent the "clustering puzzle", listed above. In conventional EAS picture one should see a single pulse in FIG. 2. Solid lines: the particle density distribution ρ(R) in simulated EAS disk versus distance from axis for different energies shown by different colours, depending on energy of the CR. Blue dots: cumulative particle density for each bimodal pulse vs. distance to EAS axis, adopted from [3]. .

FIG. 3.
Solid lines: simulated EAS disk width versus distance from axis for different energies shown by different colours. Pulse width τ0.8 in ns for the first pulse (in blue) and second pulse (in red) for the bimodal pulses versus distance to the EAS axis, adopted from [3]. . each given detector with the amplitude which depends on the distance from the EAS axis. It is not what actually observed by H10T.
2. The manifestation of the "particle density puzzle" is as follows. On Fig. 2 we show the particle density distribution ρ(R) in simulated EAS disk versus distance from axis for different energies shown by solid lines with differ-ent colours, depending on energy of the CR. In particular, for energy of the primary particle on the level of 10 17 eV one should expect a strong suppression ∼ 10 3 when distance R to the EAS axis changes from R 100 m to R 700 m. It is not what actually observed by H10T: the density of the particles ρ(R) is not very sensitive to the distance to the EAS axis, and remains essentially flat for the entire region of observations. Furthermore, the magnitude of the density ρ(R) is much higher than it is normally expected for CR energies ∼ (10 17 − 10 18 ) eV.
3. The manifestation of the "pulse width puzzle" is as follows. On Fig. 3 we show the simulated EAS disk width versus distance from axis for different energies shown by solid lines in different colours. As explained above the thickness of the EAS pancake increases with the distance from the core. Therefore, the pulse width must also increase correspondingly as shown by sold lines on Fig. 3 for different energies. It is not what actually observed by H10T: all MME events show similar duration of the pulse width on the level of (20−35) ns irrespective to the distance to the AES axis. This observation is in dramatic conflict with conventional picture as outlined above.
4. The manifestation of the "intensity puzzle" is as follows. The Fig. 2 suggests that the charged particle density ρ(R) varies between (30 − 300) particles per m 2 at the distances (500 − 700) m from the EAS axis. This is at least factor 10 2 above the expected ρ(R) for the primary particle with energy in the interval (10 17 − 10 18 ) eV. Only the particles with energies well above 10 19 eV could generate such enormous particle density as shown on Fig. 2. However the frequency of appearance for such highly energetic particles is only once every few years. This represents the "intensity puzzle" when the intensity of the events estimated from ρ(R) for MME is several orders of magnitude higher than the energy estimated by the event rate. This shows a dramatic inconsistency between the measured intensity and observed event rate carried out by one and the same detector.
We conclude this section with the following comment. From the puzzles listed above it must be obvious that the events detected by H10T are not the conventional CR air showers as they demonstrate enormous inconsistency with standard CR interpretation. What it could be? Before we put forward our proposal we would like to briefly overview in next section III the second part of our story-the AQN dark matter model.

III. THE AQN DARK MATTER MODEL
We start with few historical remarks and motivation of the AQN model in subsection III A, while in subsection III B we overview recent observations of CR-like events (such as puzzling bursts observed by the TA experiment and ANITA's two anomalous events with non-inverted polarity) of some mysterious events which could be explained by the AQN events hitting the Earth. Finally, in section III C we briefly overview some specific features of the AQNs traversing the Earth (such as internal temperature, level of ionization, etc). These characteristics will be important for the present study interpreting the MMEs as the AQN events.

A. The basics
The AQN dark matter model [6] was invented long ago with a single motivation to explain in natural way the observed similarity between the dark matter and the visible densities in the Universe, i.e. Ω DM ∼ Ω visible without any fitting parameters. We refer to recent brief review [17] on the AQN model. Here we want to mention few key elements which are important for this work.
The AQN construction in many respects is similar to the Witten's quark nuggets, see [18][19][20], and review [21]. This type of DM is "cosmologically dark" not because of the weakness of the AQN interactions, but due to their small cross-section-to-mass ratio, which scales down many observable consequences of an otherwise stronglyinteracting DM candidate.
There are two additional elements in the AQN model compared to the original models [18][19][20][21]. First new element is the presence of the axion domain walls which are copiously produced during the QCD transition 1 . The axion field θ(x) plays a dual role in this framework: first it serves as an additional stabilization factor for the nuggets, which helps to alleviate a number of problems with the original nugget construction [18][19][20][21]. Secondly, the same axion field θ(x) generates the strong and coherent CP violation in the entire visible Universe. This is because the θ(x) axion field before the QCD epoch could be thought as classical CP violating field correlated on the scale of the entire Universe. The axion field starts to oscillate at the QCD transition by emitting the propagating axions. However, these oscillations remain coherent on the scale of the entire Universe. Therefore, the CP violating phase remains coherent on the same enormous scale.
Another feature of the AQN model which plays absolutely crucial role for the present work is that nuggets can be made of matter as well as antimatter during the QCD transition. Precisely the coherence of the CP violating field on large scale mentioned above provides a preferential production of one species of nuggets made of antimatter over another species made of matter. The 1 The axion field had been introduced into the theory to resolve the so-called the strong CP problem which is related to the fundamental initial parameter θ 0 = 0. This source of CP violation is no longer available at the present time as a result of the axion's dynamics in early Universe. One should mention that the axion remains the most compelling resolution of the strong CP problem, see original papers on the axion [22][23][24][25][26][27][28], and recent reviews [29][30][31][32][33][34][35][36][37].
preference is determined by the initial sign of the θ field when the formation of the AQN nuggets starts. The direct consequence of this feature along with coherent CP violation in entire Universe is that the DM density, Ω DM , and the visible density, Ω visible , will automatically assume the same order of magnitude densities Ω DM ∼ Ω visible without any fine tuning as they both proportional to one and the same dimensional parameter Λ QCD . We refer to the original papers [38][39][40][41] devoted to the specific questions related to the nugget's formation, generation of the baryon asymmetry, and survival pattern of the nuggets during the evolution in early Universe with its unfriendly environment, see also brief review [17] covering these topics.
One should emphasize that AQNs are absolutely stable configurations on cosmological scales. Furthermore, the antimatter which is hidden in form of the very dense nuggets is unavailable for annihilation unless the AQNs hit the stars or the planets.
However, when the AQNs hit the stars or the planets it may lead to observable phenomena. The strongest direct detection limit 2 is set by the IceCube Observatory's, see Appendix A in [43]: where we formulated the constraints in term of the AQN's baryon charge B rather than in terms of its mass M ≈ m p B. Similar limits are also obtainable from the ANITA and from geothermal constraints which are also consistent with (1) as estimated in [44].
While ground based direct searches offer the most unambiguous channel for the detection of quark nuggets the flux of nuggets is inversely proportional to the nugget's mass and consequently even the largest available conventional dark matter detectors are incapable to exclude the entire potential mass range of the nuggets. Instead, the large area detectors which are normally designed for analysing the high energy CR are much better suited for our studies of the AQNs as we discuss in next section III B.
B. When the AQNs hit the Earth... For our present work, however, the most relevant studies are related to the effects which may occur when the AQNs made of antimatter hit the Earth and continue to propagate in deep underground in very dense environment. In this case the most of the energy deposition will 2 Non-detection of etching tracks in ancient mica gives another indirect constraint on the flux of dark matter nuggets with mass M < 55g [42]. This constraint is based on assumption that all nuggets have the same mass, which is not the case as we discuss below. The nuggets with small masses represent a tiny portion of all nuggets in this model, such that this constraint is easily satisfied with any reasonable nugget's size distribution.
occur in the Earth's interior. The corresponding signals are very hard to detect as the photons, electrons and positrons will be quickly absorbed by surrounding dense material deep underground, while the emissions of the very weakly interacting neutrinos and axions are hard to recover. In this short subsection we want to mention several observed phenomena which could be related to the AQN annihilation events when the nuggets propagate in the Earth's atmosphere. We start with the Telescope Array (TA) experiment which has recorded [11,12] several short time bursts of air shower like events. These events are very unusual and cannot be interpreted in terms of conventional CR single showers. In particular, if one tries to fit the observed bursts (cluster events) with conventional code, the energy for CR events should be in 10 13 eV energy range to match the frequency of appearance, while the observed bursts correspond to (10 18 − 10 19 ) eV energy range as estimated by signal amplitude and distribution. Therefore, the estimated energy from individual events within the bursts is five to six orders of magnitude higher than the energy estimated by event rate [11,12]. It has been argued in [9,10] that these bursts represent the direct manifestation of the AQN annihilation events.
This feature, in many aspects, is very similar to the intensity puzzle listed in previous Section II, where analogous dramatic inconsistency (between the observed intensity and measured event rate) is recorded by H10T instrument.
Our next example is the observed seasonal variations of the X ray background in the near-Earth environment.
To be more precise, the XMM-Newton at 11σ confidence level [45] had recorded the seasonal variations in the 2-6 keV energy range at distances r > ∼ 8R ⊕ where the measurements are effectively had been performed. The authors [45] argue that conventional astrophysical sources have been ruled out. Furthermore, the XMM-Newton's operations exclude pointing at the Sun and at the Earth directly, diminishing a possible direct X ray background from the Earth and the Sun. It has been argued in [46] that this seasonal variation may be naturally explained by the AQNs exiting the Earth. The AQNs continue the emission of the X rays at r > ∼ 8R ⊕ , long after the nuggets exit the Earth's atmosphere. It has been also shown that the spectrum and the intensity computed in [46] almost identically match the observed spectrum [45].
The final example we want to mention here is the observation by ANITA [14][15][16] of two anomalous events with noninverted polarity. These two events had been identified as the Earth emergent upward going CR-like events with exit angles of −27 • and −35 • relative to the horizon. These anomalous events are in dramatic tension with the standard model because neutrinos are exceedingly unlikely to traverse through Earth at a distance of > ∼ 5×10 3 km with such ultrahigh energy, even accounting for the ν τ regeneration [14]. The analysis [47] reviewed the high-energy neutrino events from the IceCube Neutrino Observatory and inferred that the ν τ interpretation is excluded by at least 5σ confidence level. We advocated in [13] that these two anomalous events with noninverted polarity observed by ANITA could be interpreted as the upward going AQNs.
The applications [46] to the X ray emission in the near-Earth environment and the ANITA anomalous events [13] are especially relevant for the present work because both these phenomena are originated from the AQN upward going (Earth emergent) events when the AQNs traversed through the Earth interior and exit the Earth surface. As we shall argue in this work the MME events could be also the consequence of the same upward moving AQNs.
However, there is an additional technical challenging problem to be addressed for the application to MME analysis. The point is that the ANITA anomalous events are generated by AQNs which just crossed the surface and entered the Earth's atmosphere (upward-going events). Exactly at this instant large number of electrons instantaneously released into atmosphere and generated the radio pulse measured by ANITA. In contrast, the X rays are observed by XMM-Newton at very large distances r > ∼ 8R ⊕ from Earth. In this case the impact of the atmospheric material at r < ∼ 60 km can be ignored, and the cooling of the AQN is dominated by nugget's propagation in an empty space. The application to MME, which is the topic of the present work, deals with an intermediate stage between these two cases: it is not the first instant when the AQN enters the atmosphere emerging from the interior, and it is not the asymptotically far away when the AQN propagates in empty space. In other words, we need to know the emission pattern and the energy deposition rate for the AQNs when they propagate in the atmosphere, and the impact of the atmospheric surrounding material cannot be ignored. We formulate our proposal in Section IV where we argue that emission of the electrons by the AQN in form of the well-separated "bunches" may explain the unusual features of the MME observations.
However, before we present our arguments supporting this interpretation of the MME unusual events (as highlighted in section II) we should overview the basic characteristics of the AQNs traversing the Earth, which represents the topic of the next subsection III C.

C. Internal structure of the AQN
The goal here is to explain the basic features of the AQNs when they enter the dense regions of the surrounding material and annihilation processes start. The related computations originally have been carried out in [48] in application to the galactic environment with a typical density of surrounding visible baryons of order n galaxy ∼ 300 cm −3 in the galactic center, in dramatic contrast with dense region in the Earth's interior when n rock ∼ 10 24 cm −3 and atmosphere with n air ∼ 10 21 cm −3 . We review these computations with few additional elements which must be implemented in case of propagation in the Earth's atmosphere and interior when the density of the environment is much greater than in the galactic environment.
The total surface emissivity from electrosphere has been computed in [48] and it is given by where α ≈ 1/137 is the fine structure constant, m = 511 keV is the mass of electron, and T is the internal temperature of the AQN. One should emphasize that the emission from the electrosphere is not thermal, and the spectrum is dramatically different from black body radiation, see [48], see also Appendix A with more details.
A typical internal temperature of the AQNs for very dilute galactic environment can be estimated from the condition that the radiative output of Eq. (2) must balance the flux of energy onto the nugget where n represents the number density of the environment. The left hand side accounts for the total energy radiation from the AQN's surface per unit time as given by (2) while the right hand side accounts for the rate of annihilation events when each successful annihilation event of a single baryon charge produces ∼ 2m p c 2 ≈ 2 GeV energy. In Eq. (3) we assume that the nugget is characterized by the geometrical cross section πR 2 when it propagates in environment with local density n with velocity v AQN ∼ 10 −3 c. The factor κ is introduced to account for the fact that not all matter striking the AQN will annihilate and not all of the energy released by an annihilation will be thermalized in the AQNs by changing the internal temperature T . In particular, some portion of the energy will be released in form of the axions, neutrinos 3 . The high probability of reflection at the sharp quark matter surface lowers the value of κ. The propagation of an ionized (negatively charged) nugget in a highly ionized plasma (such as solar corona) will increase the effective cross section. As a consequence, the value of κ could be very large as discussed in [49] in application to the solar corona heating problem.
The internal AQN temperature had been estimated previously for a number of cases. It may assume dramatically different values, mostly due to the huge difference in number density n entering (3). In particular, for the galactic environment T galaxy ≈ 1 eV, while in deep Earth's interior it could be as high as T rock ≈ (100 − 200) keV. Precisely this value of T had been used as initial temperature of the nuggets in the proposal [46] explaining the seasonal variations of the X rays observed the XMM-Newton at 11σ confidence level [45] at distances r ∼ (6 − 10)R ⊕ from the Earth surface.
The crucial element for the study [46] was the emission rate by AQNs at very high temperatures T ≈ (100 − 200) keV. In this case the emission is not determined by simple Bremsstrahlung radiation given by (2) which is only valid for relatively low temperatures. In the high temperature regime a number of many-body effects in the electrosphere, that were previously ignored, become important. In particular, it includes: the generation of the plasma frequency in electrosphere, which suppresses the low frequency emission. Also, the the AQNs get ionized which dramatically decreases the density of positrons in electrosphere, and consequently suppresses the emission, among many others effects. All these complications have been carefully considered in [46] in applications to the seasonal variations of the X rays as observed the XMM-Newton, and we refer to that paper for the details. Here we quote the result of these studies. The rate of emission can be effectively represented as follows: where factor η(T, R) is a result of strong ionization of the electrosphere, which leads to the corresponding suppression of the emission. The details of this suppression are explained in Appendix A and given by (A12). This suppression is a direct consequence of high internal temperature T when a large number of weakly bound positrons are expanded over much larger distances order of R rather than distributed over much shorter distances of order m −1 around the nugget's core. This basically determines the suppression factor η ∼ (mR) −1 ∼ 10 −6 .
In both previously considered applications (to ANITA [13] and to XMM-Newton [46] experiments) when the upward moving AQNs with high T ≈ (100 − 200) keV play the key role in the explanations of the observations, we ignored an additional 4 annihilation events with atmospheric material as we already mentioned at the end of the previous subsection III B. Now we want to include the corresponding physics in our analysis.
First of all we want to estimate the number of direct head-on collisions of the atmospheric molecules with AQN per unit time. It can be estimated as follows: where N m 2.7 · 10 19 cm −3 is the molecular density in atmosphere when each molecule contains approximately 30 baryons such that the typical density of surrounding baryons in air is n air 30 · N m 10 21 cm −3 . The dominant portion of these collisions are the elastic scattering processes rather than successful annihilation events suppressed by parameters κ as discussed above. The energy being deposited to the AQN per unit time as a result of annihilating processes can be estimated as follows One should emphasize that this energy is being deposited into the AQN which is already at full capacity with very high temperature T (100 − 200) keV accumulated during a long journey in much denser Earth's interior.
What is the mechanism to release this extra energy? In other words: In what form the energy (6) could be emitted into the surrounding atmosphere as it cannot be easily transferred to the AQN's quark core (as it is already saturated, see footnote 4)?
There are two qualitatively different regimes which can be normally realized in such circumstances. If a typical time scale τ deposit to deposit a specific amount of energy (let us say, 1 GeV) determined by (6) is longer than the time scale τ cool to release the same amount of energy, i.e. τ deposit > ∼ τ cool than a continuous cooling process takes place and the temperature slowly decreases. This is a thermodynamically equilibrium behaviour which can be analyzed using conventional technical tools 5 .
The second option is realized when the time scale to deposit the energy τ deposit is much shorter than the time scale of cooling processes, i.e. τ deposit τ cool , in which case the burst-like (explosion like, eruption-like blasts) non-equilibrium processes must occur. These eruptionlike events will release the extra accumulated energy in form of very short pulses, which cannot be described as conventional cooling equilibrium processes. These short pulses must alternate with much longer periods of accumulation energy which inevitably end with subsequent eruption-like events. In other words, in this case, the nonequilibrium fast equilibration process can be thought as the eruption-like event 6 . 5 For example, the AQNs moving in empty space and slowly cooling by emitting the X rays as computed in [46] belongs to this class. 6 An analogy for such eruption-like event is the lightning flash under thunderclouds. The clouds accumulate the electric charge in form of the ionized molecule very efficiently. The corresponding time scales plays the role of τ deposit in our system. The neutralization of these ions is less efficient process, which is analogous to our τ cool . If some conditions are met (the so-called runaway breakdown conditions are satisfied, see [50,51] for review) the discharge occurs in form of the eruption which is the lightning strike in our analogy. The system is getting neutralized in form of the non-equilibrium lightning event (eruption) on the time scales which are much shorter than any time scales of the prob-Now we give a numerical estimation for the cooling rate determined by (4) to find out what regime is realized in given circumstances. The order of magnitude estimate can be expressed as follows: dE emiss dt ∼ 10 16 T 100 keV 17 4 GeV s .
Comparing this estimate for cooling rate with estimate for the deposition rate (6) one concludes that the rate of energy deposition is much shorter than the rate at which the system is capable to release the same amount of energy in form of the photon's emission. As a result, we arrive to the conclusion that second option when τ deposit τ cool is realized. This conclusion inevitably implies the eruption like events must occur, and these very short explosions must be alternated with much longer periods of time when the energy is being accumulated by the AQNs. These longer periods of the energy accumulations also must end with consequent eruption-like blasts. This alternating pattern continues as long as following condition holds: τ deposit τ cool ⇒ (short eruptions occur). (8) With these preliminary comments with estimates from the previous studies on the AQN features we are now in position to formulate the proposal interpreting the MME events observed by X10T collaboration in terms of the upward moving AQNs, which is the topic of the next section IV.

IV. MME AS THE AQN EVENT
In this Section we formulate the basic idea of our proposal on identification of the unusual Multi-Modal Events with the upward moving AQN events, while supporting estimates on the event rate, the intensity, and the variety of time scales will be presented in following sections V and VI.
The AQNs which propagate in the Earth's interior are very hot. Their temperature could be as hot as T ≈ (100−200) keV at the moment of exit on the Earth's surface as we discussed in [46] in the application to studies of the seasonal variations of the X rays observed by XMM-Newton at very large distances r ∼ (6 − 10)R ⊕ from the Earth. At such large distances the AQN's cooling can be described by conventional thermodynamically equilibrium processes with well defined spectrum, see footnote 4. It should be contrasted with our application to ANITA anomalous events [13] when the same upward lem. This analogy is in fact, quite deep, and will be used in section VI when we discuss different time scales of MMEs as observed by H10T detector.
moving AQNs just crossed the surface and entered the Earth's atmosphere. In this case the dominant cooling mechanism is realized in form of short burst-like events according to (8). The outcome of these short burst-like events is emission of numerous bunches (clumps) of relativistic electrons. These eruption like events are analogous to lightning flashes from footnote 6.
In our studies [13] we focused on an estimation of the number of emitted electrons N and their energy E ∼ 10 MeV which, according to the proposal [13], essentially determine the intensity and the spectral features of the ANITA anomalous events. These electrons will be always accompanied by the emission of much more numerous number of photons with typical energy of order ω ∼ T ≈ (100−200) keV, which in fact represent the dominant fast non-equilibrium mechanism of cooling of the system in these circumstances (8).
Precisely these multiple explosion-like emissions of photons and electrons will be identified with Multi-Modal Events with their unusual features reviewed in Section II. The mean free path of the photons with ω ∼ T is relatively short, measured in meters. Therefore, they cannot be directly observed. In contrast, the mean free path of the energetic electrons at sea level with E ∼ 10 MeV is around several kilometres in atmosphere (and even longer at the elevation of 3346 m where H10T is located). The propagation of these ultra-relativistic electrons takes place in the background of the geomagnetic field B ∼ 0.5 G which determines the instantaneous curvature ρ ∼ 3 km of the electron's trajectories with such energies, see section V with relevant estimates. Due to the large component of magnetic field parallel to the Earth's surface the electrons may move in downward direction even if they are initially emitted in upward direction. Our proposal is that precisely these downward moving energetic electrons could mimic the CR and could be responsible for the Multi-Modal pulses as detected by H10T instrument.
Any precise computation during a short period of time of emission is very hard problem of non-equilibrium dynamics, which is beyond the scope of the present work. Fortunately, the observable intensity, the event rate, the time delays between the pulses are determined in the AQN framework by the basic characteristics of the model and are not very sensitive to the details of this nonequilibrium mechanism of production and its corresponding time scales. This is because all these observables depend on several parameters such as typical energy E ∼ 10 MeV of the emitted electrons, typical number of emitted electrons N , the deposition rate (6) and the cooling rate (7) which had been previously estimated for very different purposes in different circumstances. We shall use the same parameters in the present work.
In particular, the number of electrons N ∼ (10 8 − 10 9 ) emitted by AQNs is very hard to compute from the first principles. However it can be fixed by the observed intensity and spectral features of the anomalous ANITA radio pulses assuming, of course, that these observed ra-dio pulses with inverted polarities are due to the upward moving AQNs [13]. The typical energy E ∼ 10 MeV of these electrons is also consistent with duration of the observed pulses [13]. The next section V is devoted to the estimation of the event rate of MME, while section VI is devoted to explanations of the puzzles which demonstrate a dramatic deviation from conventional CR picture as formulated in section II. We shall argue that the observed MME features are consistent with our AQN-based interpretation.

V. EVENT RATE OF MME
This section is devoted to estimate of the event rate of MMEs within the AQN framework. We anticipate that this evaluation is expected to be very qualitative estimation due to large uncertainties in parameters and rare occurrence of the observed MMEs. In particular, it is known that key parameter, the DM density locally may dramatically deviate from the well established average global value ρ DM ≈ 0.3 GeV cm −3 . Nevertheless we would like to present such estimate to demonstrate that our interpretation of MMEs as an outcome of upwardgoing AQNs is at least a self-consistent proposal. In what follows we use the same formulae we used for the estimates of the frequency of appearance of the ANITA Anomalous Events [13] such that the numerous uncertainties related to the local DM density or/and the AQN size distribution would not dramatically affect our estimates below if they are normalized to the ANITA event rate.
We start with the same formula from [13] for the expected number N of the MMEs assuming that they are induced by the AQNs: where ∆Ω ≈ 2π for the isotropic AQN flux, and the expression for the local rate of upward-going AQNs per unit area is given by: where ρ DM is the local density of DM and R ⊕ = 6371 km is the radius of the Earth and Φ is the total hit rate of AQNs on Earth [43]. In formula (9) the T is the time of operation while A eff is the effective area to be estimated below.
The estimation of the effective area A eff is a complicated task as it is not simply determined by the area of the detectors similar to standard CR analysis. Instead, it is determined by the area along the AQN's path where it can emit the bunches of electrons in form of short pulses which can mimic the CR as outlined in previous section IV. The area A eff ∼ L AQN · D AQN can be thought as a strip of length L AQN and width D AQN .
To estimate these parameters let us consider an AQN moving in upward direction with angle θ AQN with respect to the Earth's surface such that upward velocity component is (v AQN · sin θ AQN ). This vertical velocity component determines the maximal height h of the atmosphere where the atmospheric density is still sufficiently high such that condition (8) holds. We estimate h ≈ (30 − 50) km where the atmospheric density falls by two order of magnitude such that conventional X-ray emission (which had been used in computations [46] to explain the seasonal variations observed by XMM Newton) becomes the dominant cooling process at higher altitudes. This condition determines the maximal length L AQN and time scale τ AQN when the AQNs can emit the bunches of electrons in form of pulses which we identify with MMEs and which can mimic the CR air showers as outlined in previous section IV. Numerically these parameters are estimated as follows: where we used θ AQN 30 0 for an order of magnitude numerical estimates.
To estimate the width of the strip D AQN we recall that the geomagnetic field B in location of the H10T instrument can be characterized by two numbers: it has strong component in up to down direction B down 0.49 G, and strong component from south to north direction B north 0.24 G, see e.g. [52].
Assuming that the energetic electrons with E ∼ 10 MeV will be emitted along the v AQN -direction one can estimate instantaneous radius of curvature ρ, see e.g. Jackson [53]: where B ≈ B north ≈ 0.24 G is the local magnetic field strength, which changes the direction of the electrons from upward to downward moving such that they can mimic the CR air showers. The angle θ B is the angle between the particle electron velocity v and magnetic direction. We choose θ B ≈ 30 • in (12) for the numerical estimates. Now we estimate the relevant scale λ which determines the survival pattern of the electron's bunches. It is mostly determined by the Coulomb elastic scattering with cross section σ Coul The electrons with θ > ∼ 1/2 may strongly deviate from their main paths and cease to stay with majority of particles forming the bunch (which eventually becomes the pulse being interpreted as MME). The corresponding length scale λ(h) at altitude h (accounting for the air-density n air (h) variation) is estimated as follows The physical meaning of λ is the length distance particles propagate at which the majority of the particles in the bunch remain within the bunch before the dispersing to much larger distances when they cease to be a part of the pulse (burst). Important comment here is that λ > ∼ ρ such that the majority of electrons within the bunch do not strongly re-scatter, and therefore survive the change in orientation due to the geomagnetic field B (from upward to downward direction).
The width of the strip D AQN now can be roughly estimated as follows: Combining all the estimates above one arrives to the following order of magnitude estimation for the expected number of MME events: This order of magnitude estimate should be compared with observed frequency of appearance of the MMEs. According to [2] the observations from start until August 2016 during T ≈ 3, 500h ≈ 0.4yr of operation the H10T instrument had detected N obs 10 3 MMEs. According to [3] the observations from Feb 15, 2018 to May 12, 2018 (which corresponds to T ≈ 0.25 yr) the H10T instrument had detected N obs = 217 MMEs.
The AQN based estimations (16) suggest that number of events N is almost one order of magnitude lower than N obs events observed by H10T. Nevertheless, we consider this order of magnitude estimation (16) being consistent with our proposal due to many uncertainties which enter this estimate. There are several reasons for this optimistic view.
First, as we already mentioned the parameters entering (10) are in fact not precisely known. The fundamental parameters such as ρ DM and v AQN only have accuracy up to order one as the local flux distribution of DM and size distribution of AQN remain unknown to date as the DM density locally may dramatically deviate from the well established average global value ρ DM ≈ 0.3 GeV cm −3 , see introduction in [54] for the references and details. Furthermore, the size distribution factor entering (10) in form B −1 had been fixed from dramatically different physics (including solar corona heating puzzle) and can easily deviate by large factor. The estimation of the effective area A eff also suffers enormous uncertainty.
Furthermore, the numerical suppression factor ∼ 0.1 (between computed N and observed N obs values) which appears in our order of magnitude estimates (16) is very similar to suppression factor ∼ 0.1 which occurred in analogous estimates for the mysterious TA bursts [9] and the ANITA anomalous events [13]. This similarity hints on a common origin for all these phenomena. Therefore, if one normalize N to the observed ANITA anomalous events (assuming that all three phenomena originated from the same AQN based physics) one arrives to correct order of magnitude estimate as the suppression factor ∼ 0.1 is approximately the same for all three phenomena estimated for the mysterious TA bursts in [9], the ANITA anomalous events in [13], and Multi-Modal Events estimated in (16).
Our main arguments, however, leading to the identification of the MMEs with the AQN events are not based on an order of magnitude estimate of the frequency of appearance (16) which is presented in this work exclusively for illustrative purposes. Rather, our main arguments are based on specific qualitative features (formulated as puzzles) which have been listed in Section I, and which cannot be explained in terms of the conventional CR air showers as highlighted in Section II. We consider a qualitative agreement between the observations and our theoretical estimates (to be discussed in next section VI) as the strong arguments supporting our identification.

VI. AQN PROPOSAL CONFRONTS THE MME
OBSERVATIONS [1][2][3][4][5] We start this section by estimating the particle number density ρ(R) which appears in formulation of the puzzles in Section II. The particle number density in the AQN framework is determined by the number of particles which could be observed at the detector site. Assuming that the total number of electrons being emitted in form of a pulse (as a result of erupted burst) is similar to the number of particles which was used in our analysis [13] of the ANITA anomalous events N (10 8 − 10 9 ) we arrive to the following estimate for ρ(R): where the area which is hit by the bunch of particles from a single burst is estimated as (λ∆θ) 2 with λ(h) given by (14) at the detector site with h 4 km. It is assumed that the particles are propagating within the cone with angle ∆θ < ∼ 1/2. Our next task is the estimation of the time delay between the pulses within the same cluster. The corresponding time delay can be estimated from the condition that the energy deposited into the AQN during time τ delay must be released in form of the eruption when N (10 8 − 10 9 ) electrons being emitted in form of a single pulse. As we mentioned previously, the dominant portion of the energy in such eruption is released in form of the photons with energy ω ∼ T which always accompany the electron's emission as a result of these bursts.
The relative ratio r between these two components is estimated in Appendix B, see eq. (B1). Numerically, the estimate can be presented as follows: where for numerical estimates we use T 200 keV. We must emphasize that (18) is really an order of magnitude estimate as mechanism of emission is determined by a very complex physics as highlighted in Appendix B. Furthermore, the corresponding eruption event is not a thermodynamically equilibrium process as we already mentioned previously. The dominant contribution in form of the photon's emission during a single pulse can be estimated in terms of these parameters as follows E γ ∼ N · (2m) · r −1 ∼ 1.4 · (10 10 − 10 11 ) GeV, (19) where we use the numerical value for N ∈ (10 8 − 10 9 ) extracted from analysis of the ANITA anomalous events [13] and parameter r is given by (18). In this estimate we assumed that the energy of the pair E e + e − ∼ 2mN computed at the moment of the pair creation. The released energy during a single pulse (19) determines the time delay τ delay between the pulses within the same cluster: where we used numerical value for (dE deposit /dt) from (6) with density of atmosphere at h 4 km where the detector H10T is located. Our final, the order of magnitude estimate for the time delay τ delay between the pulses can be represented as follows where we use κ 0.1 in (21) for numerical estimates which was previously extracted from studies of the emission in the dilute galactic environment, see footnote 3.
While this is very crude, an order of magnitude estimate, the important point here is that the time scale τ delay may dramatically vary from event to event as it strongly depends on many parameters entering the problem. In particular, it obviously depends on intensity (number of emitted particles N ) of a previous pulse. Essentially the parameter (21) should be interpreted as a preparation time the system requires for the next eruption.
The strong variation of the τ delay should be contrasted with another parameter which is the time duration of a single pulse τ pulse . In the AQN framework this parameter must not vary much from one event to another as it is entirely determined by internal dynamics of the AQN during the eruption irrespective to the intensity of the previous events, and irrespective to the prehistory of the AQN propagation as long as condition (8) is met 7 . In different words, our proposal suggests that where the constant cannot be computed from first principles as it is entirely determined by the non-equilibrium dynamics of the AQN at the instant of eruption 8 . Now we are prepared to confront the basic consequences of our proposal (identifying the MMEs with the AQN annihilation events) with observations. 1. "clustering puzzle": The multiple number of events is a very generic feature of the system as explained at the very end of subsection III C as long as condition (8) is met. The time delays between the pulses τ delay dramatically fluctuate between different events. The window for these variations is huge according to (21). In fact, it could be even well outside of this window. The time delays τ delay may considerably vary even for different events within the same cluster at the same detection point, which is consistent with observed cluster shown on Fig. 1. Furthermore, the AQN itself remains almost at the same location as the displacement ∆l AQN during entire cluster of events is very tiny which implies that all individual bunches making the cluster are likely to be emitted by the AQN along the same direction, and can be recorded and classified as MME by H10T. Each event can be viewed as an approximately uniform front as mere notion of the "EAS axis" does not exist in this framework, see also next item. However, each individual event may appear to arrive from slightly different direction due to the inherent spread of the emitted electrons at the moment of eruption.
2. "particle density puzzle": Particle density distribution ρ(R) in the AQN framework is estimated by (17) and it shows strong fluctuation from one event to another event. These variations are mostly related to the intensity of the individual bursts being expressed by the number of electrons in the bunch N . However, the distinct feature of the distribution in the AQN framework is that it does not depend on R, as the notion of the "EAS axis" does not exist in this framework as we already mentioned. All these generic features of the AQN framework are perfectly consistent with H10T observations as presented on Fig. 2. However, these observed features are in dramatic contradiction with conventional EAS prediction shown by solid lines on Fig. 2 with different colours, depending on energy of the CR.
Furthermore, the magnitude of the density ρ(R) in the AQN framework (which is mostly determined by parameter N ) is much higher than one normally expects for CR energies ∼ (10 17 − 10 18 ) eV. The corresponding parameter N representing the number of electrons in the bunch was not fitted for the present studies to match the observations. Instead, it was extracted from different experiment in dramatically different circumstances (in proposal [13] to explain the ANITA anomalous events as the AQN events).
3. "pulse width puzzle": In the AQN framework the width of the pulse (22) cannot vary much from one event to another. It is a fundamental feature of the framework because the duration of the pulse is entirely determined by internal dynamics of the AQN during the blast as long as condition (8) is met, see footnote 7 with a comment. This feature is in perfect agreement with observations [3] τ pulse ≈ (20 − 35)ns ⇐ [observations] (24) for all recorded MMEs. At the same time, this feature is in dramatic conflict with conventional picture when the duration of the pulse must depend on the distance to the EAS axis as shown by sold lines on Fig. 3. This basic prediction of the conventional CR analysis is due to increase of the thickness of the EAS pancake with the distance from the EAS axis. As we already mentioned the mere notions such as the "EAS axis" and the "thickness of the EAS pancake" do not exist in the AQN framework. 4. "intensity puzzle": Particle density distribution ρ(R) in the AQN framework is estimated by (17). The corresponding event to event fluctuations do not depend on the distance to the EAS axis as we already mentioned. Such intensity of the events as given by (17) in the AQN framework is consistent with observations shown on Fig.  2. However, the observations are in dramatic conflict with conventional CR analysis when such huge intensity could be generated by a primary particle with energy well above 10 19 eV with dramatically lower event rate on the level of once every few years. The frequency of appearance in the AQN framework is estimated in section V, and it is consistent with observed event rate.
We conclude this section with the following comment. All our formulae presented in this section are the order of magnitude estimations at the very best, as they include many inherent uncertainties which are inevitable features of any composite system (such as AQN) propagating in a complex environment (such as Earth atmosphere) with very large Mach number M ≡ v AQN /c s 1 where c s is speed of sound.
Nevertheless, the emergent picture suggests that all the puzzles formulated in Sections I and II can be naturally understood within the AQN framework as explained above in items 1-4. Needless to say that all phenomenological parameters used in the estimates had been fixed long ago for dramatically different observations in different circumstances for different environments as overviewed in Section III.

VII. CONCLUSION AND FUTURE DEVELOPMENT
Our basic results can be summarized as follows. We have argued that all the puzzles formulated in Sections I and II can be naturally understood within the AQN framework as explained in items 1-4 in previous Section, and we do not need to repeat these arguments again in this Conclusion. Instead, we want to discuss the drastic differences between the events induced by conventional CR showers and the AQNs. These dramatic distinct features can be tested by the future experiments, such that our proposal can be discriminated from any other proposals and suggestions. We list below the following typical features of the AQN events and contrast them with any other possible mechanisms which could be responsible for MMEs.
1. The events which are generated by the bunches of electrons as a result of eruption of the propagating AQN in the Earth atmosphere suggests an enormous number of possibilities to generate different clusters when each event within a given cluster may have very different intensity from a previous and consequent events with very different time delays between the events. In other words, the AQN proposal suggests that there should be large variety of shapes and delays between the events with very different patterns due to the complexity of the AQN system. It should be contrasted, for example, with hypothesis of "delayed particles" (which was originally suggested to explain the MMEs) in which case all clusters must be the same as they should be determined by a specific pattern of decaying fundamental particle of unknown nature.
2. A "rule of thumb" suggests that a typical number of charged particles (mostly electrons and positrons) in CR air shower is E CR /GeV, which implies that N ≈ (10 8 − 10 9 ) for the energy of the primary particle E CR ≈ (10 17 − 10 18 ) eV. This estimate suggests that any detector which is designed to study the EAS with energies E CR ≈ (10 17 − 10 18 ) eV are, in principle, capable to study MMEs if the resolution of the detectors is in ∼ 10-ns level, similar to H10T, see also item 4 below as an alternative option to properly select and discriminate the MMEs.
3. In particular, we expect that the extension of the H10T detector would produce more multiple pulses (at each given detector) instead of simple bimodal pulses. We also expect that more detectors in the area will be recording MMEs because the area covered by each individual pulse is relatively large (few kilometres) according to (17), which is well above the present size of H10T instrument. 4. A large number of charged particles N ≈ (10 8 −10 9 ) in the background of the geomagnetic field B ∼ 0.5 gauss will produce the radio pulse in both cases: the CRinduced radio pulse [55,56] as well as AQN-induced radio pulse [13]. However, these pulses can be easily discriminated from each other as argued in [13].
The main reason for the dramatic differences between these two radio pulses is that the AQN event could be viewed as an (approximately) uniform front of size λ∆θ ∼km with a constant width, while EAS is characterized by central axis. In different words, the number of particles per unit area ρ(R) in the AQN case does not depend on the distance from the central axis, in huge contrast with conventional CR air showers when ρ(R) strongly depends the distance from the central axis. The width of the "pancake" in CR air shower also strongly depends on R. As a result, the effective number of coherent particles contributing to the radio pulse is highly sensitive to the width of the "pancake" when it becomes close to the wavelength of the radio pulse. These distinct features lead to very different spectral properties of the radio pulses in these two cases, which can be viewed as an independent characteristic of MMEs. In fact, this unique feature can be used in future studies for purpose of discrimination and proper selection of the Multi-Modal Clustering Events.
If future studies and tests (including the detecting of the synchronized radio pulses with MMEs as suggested above) indeed substantiate our proposal it would be a strong argument supporting the AQN nature of the MMEs.
We conclude this work with the following comment. We estimated the event rates for three dramatically different puzzling CR-like events: for mysterious TA bursts [9], for the ANITA anomalous events [13], and finally for the Multi-Modal Events in this work as estimated in (16). All three puzzling phenomena are proportional to one and the same AQN flux (10). The self-consistency between all three estimates hints on a common nature for these puzzling CR-like events. We interpret this selfconsistency in the event rates as an additional indirect argument supporting the AQN nature for all three mysterious phenomena, while our direct arguments are presented in section VI and listed as item 1-4. We finish on this optimistic note. with where n(z) is the positron number density of the electrosphere, z = 0 corresponds to the surface of the nuggets, and n(z = 0) reproduces an approximate formula for the plasma density in the Boltzmann regime at the temperature T as computed in [48]. The 1d approximation formulated in terms of distance z from the surface was more than sufficient sufficient to study the low temperature behaviour. However, after an AQN crosses the Earth, it acquires a very high internal temperature of the order T (100 − 200) keV, and approximation (A1) is not sufficient. This is because temperature's rise will cause the electrosphere to expand well beyond the thin layer ∼z surrounding the nugget surface. Some positrons will leave the system but the majority will stay in close vicinity of the moving, negatively charged nugget core. Consequently, the nugget will acquire a negative charge of approximately −|e|Q with the number of positrons Q estimated as: The distance ρ at which the positrons remain attached to the nugget is given by the capture radius R cap (T ), determined by the Coulomb attraction: In order to estimate η entering (4) we start our analysis with an estimation of the positron density n(ρ, T ) when the AQNs enter the Earth atmosphere moving in upward direction. Formula (A3) shows that the capture radius R cap (T ) could be much larger than R, i.e. we have R < ∼ ρ < ∼ R cap . The expression (A1), which is valid when z = |ρ − R| R, in close vicinity to the core, breaks down for ρ > ∼ R. In that case the curvature of the nugget surface cannot be neglected and one should use a truly 3-dimensional formula for n(ρ, T ) instead of the 1-dimensional approximation (A1). We will assume that the density n(ρ, T ) behaves as a power law for ρ > ∼ R with exponent p: where n 0 (T ) ≡ n(ρ = R, T ) is a normalization factor and p a free parameter. Eq. (A4) is consistent with our previous numerical studies [57] of the electrosphere with p 6. It is also consistent with the Big Bang Nucleosyntheses (BBN) study discussed in [58] and with the conventional Thomas-Fermi model 9 at T = 0 [59]. While keeping p as a free parameter, we will show below that our main claim is not very sensitive to the precise value of p. The normalization n 0 (T ) can be estimated from the condition that the finite portion of the positrons satisfying Eq. (A3) is such that the total number of positrons surrounding the nugget is approximately equal to the ionization charge Q determined by (A2). Therefore, we arrive to the following expression for the normalization factor n 0 (T ): Note that the integrant ρ 2 n(ρ, T ) is mostly dominated by the inner shells ρ ∼ R such that the intgration can be extended to infinity with very high accuracy, instead of cutting off at R cap . The resulting estimate for n 0 (T ) assumes the form: which, as expected, is smaller than n(z = 0) (mT ) 3/2 because the positrons now are distributed over a distance of order R from the nugget core surface, rather than over a distance of orderz. This suppression is convenient to represent in terms of the dimensionless ratio: The next step is the calculation of the Bremsstrahlung emissivity by the nuggets due to the strong suppression of the plasma density (A7) as a result of the electrosphere's expansion. The spectral surface emissivity is denoted as dF/dω = dE/dtdAdω, representing the energy emitted by a single nugget per unit time, per unit area of nugget 9 In [59], the dimensionless function χ(x) behaves as χ ∼ x −3 at large x. The potential φ = χ(x)/x behaves as φ ∼ x −4 . The density of electrons in Thomas-Fermi model scales as n ∼ φ 3/2 ∼ x −6 at large x.
surface and per unit frequency. For low temperature, when approximation (A1) for the flat geometry is justified, the corresponding expression assumes the form [48]: where n(z) is the local density of positrons at distance z from quark nugget surface and function K(ω) does not depend on z and describes the spectral dependence of the system: (A9) The dimensionless function h ω T in (A9) is a slowly varying logarithmic function of ω/T which was explicitly computed in [48]. In order to calculate the emissivity with the positron density (A4), when the nugget core is ionized, we have to replace the integral dz in (A8) by the spherical integration over ρ: where the normalization n 0 (T ) is determined by (A6).
Using the positrons distributed according to Eq. (A4) the nugget surface emissivity can be calculated as follows: The key result here is the strong suppression factor R −1 which was not present in the no-ionization case (A8) when the positrons are localized in a thin layer. It is instructive to compare the emissivity given by eq. (A11) accounting for ionization of the nugget core with the original formula (A8). The dimensionless ratio is given by The total emissivity integrated over all frequencies in case of low temperature without ionization has been computed in [48] and it is given by (2). In case of strong ionization, Eq. (A12) implies: The critical element leading to the suppression factor (A12) is the very small quantity (mR) −1 ∼ 10 −6 . It is a direct result of the expanding electrosphere and high level of ionization at very high temperature T . The expression (A13) with numerical suppression (A12) is precisely the formula (4) which we used in main body of the text to proceed with our key arguments.
The main goal of this Appendix is to overview the photon emission and the e + e − production at high temperature T > ∼ 10 2 keV from high density QCD phases. The corresponding studies [60][61][62][63][64][65][66] have been carried out in the past in context of the quark stars. In context of the present work all the key ingredients relevant for e + e − production are also present in the system. Indeed, the AQN is characterized by very high temperature T > ∼ 10 2 keV, the quark core is assumed to be in CS dense phase, and a strong internal electric field is also present in the system. However, we cannot literally use the results from the previous studies obtained in context of the quark stars because the size of the AQN is much smaller than the relevant mean free paths for all elementary processes as discussed in details in [13]. As a result the thermal equilibration cannot be achieved in the AQN system, and entire physics is determined by non-equilibrium dynamics in high temperature regime. It should be contrasted with large size quark stars where the thermal equilibrium is maintained.
Nevertheless, it is very instructive to review the relevant results from the previous studies [60][61][62][63][64][65][66] on quark stars due to the following reasons. First, it explicitly shows the role of the main ingredients of the system, such as temperature and the large electric field. Secondly, it demonstrates the complexity of the problem when even a much simpler case of the bare quark star remains to be a matter of debates.
The idea on possibility of the e + e − emission at high temperature from quark stars was originally suggested in [60,61]. The temperature range considered in [60,61] includes the typical temperatures T > ∼ 10 2 keV which is expected to occur in our case when the AQN exits the Earth's surface as mentioned in Section III C. In refs. [62,63] the authors argue that bremsstrahlung radiation from the electrosphere could be much more important than e + e − emission. It has been also argued that a number of effects such as the boundary effects, inhomogeneity of the electric field, and the Landau-Pomeranchuk-Migdal (LPM) suppression may dramatically modify the emission rate. In [64] it has been argued that the Pauli blocking will strongly suppress the bremsstrahlung emission. Finally, in refs. [65,66] it has been argued that the so-called mean field bremsstrahlung could be the dominant mechanism.
It is not the goal of the present work to critically analyze all these suggested mechanisms of the emission. Rather, our goal is to demonstrate that even a relatively simple system of the bare quark star remains to be a matter of debates. Our case of the emission from AQN at high temperature is even more complicated as it is determined by non-equilibrium dynamics.
In the present work we do not even attempt to solve this very ambitious problem of estimating the absolute intensity rate of the e + e − and γ emissions. Rather, the absolute number N (10 8 − 10 9 ) of the e + e − pairs in a single pulse was extracted from ANITA anomalous events as explained in Section VI. In particular, this number N enters the estimate (17) for the density of the observed particles by H10T experiment [1][2][3][4]. Our goal here is much less ambitious as we try to estimate the relative ratio r between the e + e − pair production and the γ radiation assuming that the emission of both components is mostly originated from the region in electrosphere where the Boltzmann regime is justified, and the Pauli blocking suppression is not very dramatic and can be ignored.
In this case, the relative ratio r between these two components contributing to the emission can be estimated as follows. It is assumed that the dominant contribution to the γ emission comes from the Bremsstrahlung radiation (e + e + → e + e + + γ) resembling the low temperature case considered in [48]. The e + e − pair is produced by a similar mechanism through the virtual photon (e + e + → e + e + + γ * ) which consequently decays to the pair: γ * → e + e − . If this is the dominant mechanism of radiation the pair production is expected to be suppressed by a factor α/2π as a result of conversion of the virtual photon γ * → e + e − to pair 10 . It must be also suppressed by a factor exp(−2m/T ) describing the suppression in the density distribution because the virtual photons must have sufficient energy > ∼ 2m to produce e + e − pair. This oversimplified estimate leads to the following expression where for numerical estimates we use T 200 keV, and ignored a possible complicated function which could depend on parameter T /m ∼ 0.4, which is order of unity for our case. We must emphasize that (B1) is really an order of magnitude estimate as mechanism of emission is not a thermodynamically equilibrium process as we already previously mentioned. Formula (B1) is used in the main body of the text in (18) for our estimate (21) for τ delay which is measured by H10T experiment [1][2][3][4] as it describes a typical time delay between subsequent pulses representing the Multi-Modal Clustering Events.