Telescope array bursts, radio pulses and axion quark nuggets

Telescope array (TA) experiment has recorded Abbasi et al. (Phys Lett A 381(32):2565–2572, 10.1016/j.physleta.2017.06.022, 2017), Okuda (J Phys Conf Ser 1181:012067, 10.1088/1742-6596/1181/1/012067, 2019) several short time bursts of air shower like events. These bursts are very distinct from conventional single showers, and are found to be strongly correlated with lightnings. In our previous work Zhitnitsky (J Phys G 48(6):065201, 10.1088/1361-6471/abd457, arXiv:2008.04325 [hep-ph], 2021) we proposed that these bursts represent the direct manifestation of the dark matter (DM) annihilation events within the so-called axion quark nugget (AQN) model. In the present work we suggest to test this proposal to search for the radio signals in frequency band \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu \in $$\end{document}ν∈ (0.5–200) MHz which must be synchronized with the TA bursts. We argued that the conventional lightning-induced radio emission can be easily discriminated from the AQN-induced radio pulses discussed in this work.


Introduction
This work is tightly linked to our recent proposal [3] interpreting the mysterious bursts observed by telescope array (TA) experiment [1,2] in terms of the AQN annihilation events under the thunderstorm. These events are very unusual and cannot be interpreted in terms of conventional single showers as reviewed below. In the present work we shall argue that the proposed mechanism [3] in terms of the AQN annihilation events inevitably predict the radio wave pulses which must be synchronized with TA bursts. Based on this prediction we suggest to test our proposal by searching for the radio signals in frequency band ν ∈ (0.5-200) MHz which must be synchronized with the TA bursts, which represents the main goal of the present work. a e-mail: xunyul@phas.ubc.ca b e-mail: arz@phas.ubc.ca (corresponding author) Such test would unambiguously support or refute the proposal (interpreting TA bursts as the AQN events) because the main arguments of the present work based on synchronization which represents pure geometrical property of the system. This is because the radio pulse and the TA burst are originated from the same location, emitted at the same instant, and propagate with the speed of light. This feature of synchronization is not sensitive to many assumptions and uncertainties which were inevitably present in estimates [3] interpreting the TA bursts as the AQN annihilation events. Furthermore, we also argue that the conventional lightninginduced radio emission is qualitatively different from the AQN -induced radio pulses discussed in this work. Therefore, these two different radio signals can be easily distinguished from each other, and we advocate to conduct such studies.
Finally, before we proceed one should emphasize from the start that AQN events are generically random events as they are related to DM particles. However, as mentioned in [3] the AQN hitting the thunderstorm area may serve as a trigger which sparks the lightning. Therefore, the AQNs play the dual role in our proposal. First, they play the same role of the CRs which are normally assumed to initiate the lightning processes. Secondly, they emit highly energetic particles being recorded as TASD bursts. This represents a key element why mysterious bursts [1,2] which are associated with lightning events can be indeed originated from the AQNs which happen to appear in the area under the thunderclouds. The event rate estimated in [3] is based on this assumption.
Our presentation is organized as follows. In next Sect. 2 we overview the basic ideas of the AQN model, in Sect. 3 we overview the TA bursts observations [1,2] and the basic ideas of the proposal [3] with the resolution of these puzzling bursts in terms of the AQN annihilation events. In Sects. 4 and 5 we argue that the emission of the radio pulse is inevitable consequence of the proposed mechanism, and we estimate its numerical characteristics. In Sects. 6 and 7 we estimate the radio pulse due to the Earth's magnetic field. In Sect. 8 we explain how the conventional radio signals induced by the lightnings could be easily discriminated from the AQNinduced radio pulses. Finally, Sect. 9 is our conclusion where we suggest to search for a synchronized radio pulse with TA burst.

The AQN DM model
The AQN model was invented two decades ago [4] to explain the observed similarity between the dark matter and the visible matter densities in the Universe. 1 In this model, DM is made out of macroscopic lumps of quarks (or antiquarks) in colour superconducting (CS) phase with characteristic mass and size of order grams and 0.1 mm respectively. In the AQN framework the baryogenesis is actually a charge segregation (rather than charge generation) process when the global baryon number of the universe remains zero at all times. Formation of AQNs relies on the coherent C P-violating axion field in early Universe, and consequently it leads to an asymmetry between matters and antimatter nuggets, which results in asymmetry between visible matter and antimatter in the observable Universe.
Unlike many conventional DM candidates such as weakly interacting massive particles (WIMPs) and axion, the AQN model naturally explains the observed similarity between the dark and visible densities in the Universe, i.e. Ω DM ∼ Ω vis , with no fitting parameters, as both DMs and visible matters share the same QCD origin. As its name suggests, the AQN is similar to Witten's quark nugget, see [5][6][7], and review [8], in many respects. This type of DM is in fact stronglyinteracting by virtue of its macroscopic size, but it serves as the DM due to its diminutive number density.
In contrast to Witten's original proposal, formation of AQNs does not rely on a first order QCD phase transition. Rather, accumulation and compression of baryon (or antibaryon) charge come from oscillations of axion domain wall (DW) bubbles during the QCD transition in early Universe which can squeeze the quark matter inside the bubbles. Due to the surface tension of the DW, ultimately the baryon (or antibaryon) charge trapped in a bubble is squeezed into a quark (or antiquark) nugget in CS phase. Because the CS phase is far more stable than the conventional hadronic phase, AQNs do not suffer from the evaporation problem, in contrast with the original Witten's proposal [5][6][7][8]. We refer to the original papers [9][10][11][12] and a brief review [13] devoted to the specific questions related to the AQN's formation, generation of the baryon asymmetry, and survival pattern in early Universe. We also refer to independent analysis [14] supporting the basic elements on the formation and survival pattern of the AQNs during the early stages of the evolution, including the BBN and CMB epochs.
For the present studies, however, we take the agnostic viewpoint, and assume that such nuggets made of antimatter are present in our Universe today irrespective to their formation mechanism. This assumption is consistent with all presently available cosmological, astrophysical and terrestrial constraints as long as the average baryon charge of the nuggets is sufficiently large as we review below.
The AQN DM model is consistent with all presently available cosmological, astrophysical and terrestrial constraints. The strongest direct detection limit is set by the IcuCube observatory, see Appendix A in [15]: (1) Similar limits are from the ANITA experiment and geothermal constrains which are also consistent with (1) as estimated in [16]. One soft constraint comes from the indirect nondetection of etching tracks in ancient mica [17], as it gives a more stringent limit in the available parameter space of DM nuggets with mass M > 55 g, or correspondingly B 10 25 in the AQN model. However this constraint is based on an oversimplified assumption that all nuggets have equal mass, which is invalid in the AQN model since the AQN's formation mechanism implies a reasonably broad distribution of nugget's size. It is also claimed in [18] that the AQNs cannot account for more than 20% of the DM density according to the neutrino flux limits in the 20-50 MeV range observed from Super-Kamiokande (SuperK). However, as pointed out in [19], the claim [18] is based on an incorrect assumption that AQNs produce a neutrino spectrum similar to the conventional baryon-antibaryon annihilation events. 2 Rather, the annihilation processes in CS phase are dramatically different from the conventional scenarios [19] because the energy scale of the lightest pseudo Goldstone mesons (pions and kaons) is in 20 MeV range, rather than 140 MeV in hadronic phase. The resulting neutrino spectrum features in CS phase computed in [19,20] are consistent with the observations. The presence of the antimatter nuggets in the system implies that there will be annihilation events leading to large number of observable effects on different scales: from galactic scales to the terrestrial rare events. In fact, there are many hints suggesting that such annihilation events may indeed took place in early Universe as well as they are happening now in present epoch. In particular, the AQNs might be responsible for a resolution of the "Primordial Lithium Puzzle" [21] during BBN epoch. The AQNs may also alleviate the tension between standard model cosmology and the recent EDGES observation of a stronger than anticipated 21 cm absorption feature as argued in [22]. The AQNs may also explain a recently observed "exotic" diffuse UV radiation in our galaxy [23] with very puzzling features, which are very hard to interpret within conventional astrophysical models [24]. The AQNs might be also responsible for famed long standing problem of the "Solar Corona Mystery" [25,26] when the so-called "nanoflares" conjectured by Parker long ago [27] are identified with the annihilation events in the AQN framework. The AQNs could be also responsible for other mysterious and anomalous CR like events (along with [3] we already mentioned). It includes a mysterious anomalous events with noninverted polarity observed by the Antarctic Impulse Transient Antenna (ANITA) collaboration [28], and Multi Modal Clustering anomalous events [29] observed by the HORIZON 10T collaboration.
We conclude this short overview of the AQN model by emphasizing that the model is consistent with all present constraints as long as average baryon charge of the nuggets satisfies the relation (1). Furthermore, with the same set of parameters within the same framework it may explain a large number of mysterious phenomena listed above, which apparently suggest that the DM might be indeed in form of the mater and antimatter quark nuggets. In what follows we assume that such antimatter nuggets exist and together with nuggets saturate the DM density today.

TA bursts as the AQN annihilation events under thunderclouds
The TA is designed for detection of extensive air showers induced by ultrahigh energy cosmic rays (CRs) consisting of 507 ground surface particle detectors (SDs) and 3 atmospheric fluorescence telescope stations. The TA collaboration have reported 10 excessively unlikely bursts of CR-like events are observed by the SDs [1,2]. Comparing to conventional CRs, the TA burst events: 1. have a much smaller dispersion in reconstructed air shower fronts [see Fig. 3 and Fig. 4 in [1]), and do not have sharp edges in waveforms, the so-called "curvature" puzzle; 2. are temporally clustered within 1 ms that is highly unlikely (chance coincidence less than 10 −4 for five-year observation) for ultrahigh energy CRs in the fitted energy range (10 18 -10 19 ) eV, the so-called "clustering" puzzle; 3. are all recorded under thunderstorm, and most of them are "synchronized" (less than 1 ms) or "related" (less than 200 ms) with the lightning events, the so-called "synchronization" puzzle, see original articles [1,2] and short reviews in [3,13] for technical details. Proposal [3] suggested the TA bursts can be a natural consequence of an antimatter AQN traversing through the thundercloud, where about 10 9 weakly bound positrons are emitted from the AQN instantaneously in presence of the strong intracloud electric field (∼ kV/cm) at altitude above 10 km.
The key mechanism in proposal [3] may be summarized as follows. First, the electric field E in the thunderclouds is characterized by the following parameters [30,31] where l a is the so-called avalanche length. The electric field is sufficiently strong to ionize an antimatter AQN and induce liberation of weakly bound positrons: where E bound ∼ keV is the binding energy of the positrons, and R cap ∼ 2 cm is the typical distance (from the nugget's core where positrons reside [3]. These liberated positrons will be accelerated to MeV energies in the background of electric field characterized by typical length scale l a 100 m according to (2): Thereafter, the positrons exit the region of strongly fluctuating electric field which is known to be present under thunderclouds.
Comparison to conventional CR air shower, the positron flux induced by AQN has a smaller dispersion angle Δα and spatial spread Δs: [3] Δs r Δα cos α 1 km cos α r 10 km where we choose the transverse component (with respect to the electric field) of the velocity v ⊥ √ 2ΔE/m 0.1c, see Fig. 1 for precise definitions of the parameters. This behaviour 3 is consistent with feature 1 coined as the "curvature" puzzle.
A TA burst in the AQN framework represents the cluster of events related to one and the same AQN which traverses through the distance L burst 250 m during Δt 1 ms at relativistic speed, rather than a "coincidental" combination of  Fig. 1 The positrons move along the cone with angle Δα and inclination angle α with respect to the vertical direction. The angular spread Δα α is assumed to be small. The spatial spread on the surface is determined by Δs, while the additional travelling path is determined by Δr , see estimates in the text. The altitude is assumed to be within conventional range h (4-12) km. Instant direction of the electric field E at the moment of exit of the positrons is also shown. Figure adapted from [3] independent CR events. It also resolves the dramatic inconsistency (when interpreted in terms of the CR air showers) for the burst events when the event rate favours an energy range in 10 13 eV, but the intensity of the events suggests (10 18 -10 19 ) eV, which explains the item 2 coined as the "clustering" puzzle.
Lastly, an AQN can serve as a trigger of runaway breakdown avalanche as discussed in [3]. 4 This explains the strong correlation between bursts and lightning events within AQN interpretation. The AQN-induced burst may or may not initiate lightnings depending on presence of other required ingredients for runaway breakdown avalanche to start.
To rephrase it, an AQN plays the dual role when it propagates in the thunderclouds: it emits the very energetic positrons, and it may also trigger the lightning, similar to conventional CR. The energetic positrons is the source for TASD bursts, while the feature of triggering explains the observed correlation between the bursts and lightnings in the AQN framework. It explains feature 3 coined as the "synchronization" puzzle, as well as it answers the question of why most but not all bursts are correlated with lightnings.
To summarize this section: the mysterious bursts (with highly unusual features 1-3 as listed at the beginning of the section) CR-like events observed by TASD [1,2] are naturally interpreted as the cluster events generated by the AQNs propagating in thunderstorm environment. Some features of these events (such as intensity and the basic normalization factors) suffer inevitable uncertainties due to very complex dynam-ical properties of the system (the AQN propagation under thunderstorm with large Mach number). However, there are some features of the system such as given by Eq. (5) which are not sensitive to these uncertainties and represent almost model-independent consequences of the proposal [3].

Radio pulse and TA burst as synchronized events
This section is devoted to studies of the radio signals which always accompany TA bursts when interpreted in terms of the AQN annihilation events under the thunderstorm as presented in previous Sect. 3. We shall argue below that the emission of the radio pulse is inevitable consequence of the proposed mechanism. Furthermore, the radio pulse must be synchronized with TA burst irrespective of whether the bursts are related or unrelated to the lightning events. This synchronization must be within 10 µs from every TA event as the positrons and radio waves originated from the same location at the same instant and both propagate with the speed of light. As the typical time and spatial spreads of the TA events do not exceed 10 µs, the synchronization must be within the same temporal and spatial range.
This section is organized as follows. In next Sect. 4.1 we study the trajectory of the liberated positrons in a typical electric field characterized by parameters (2). We analyze the total intensity and the angular distribution of the radio emission in Sect/ 4.2, while the spectral properties of the pulse are studied in Sect. 4.3. Finally, in Sect. 4.4 we study the properties of the electric field's pulse which can be detected at the surface detector area.

Positron's acceleration in electric field under thunderstorm
Our goal here is to describe the dynamics of the liberated positrons which are emitted with initial energies in keV range according to Eq. (3). To simplify things we assume that the electric field E directed along z-axis is uniform and characterized by parameters (2). In this case the solution is well known [37] and can be described as follows. At the very initial moment of acceleration at t t 0 the particle moves with non-relativistic velocities according to conventional formulae: where p 0 z ∼ p 0 ⊥ is the initial positron's momentum at the moment of liberation from AQN. Numerically, this nonrelativistic portion of the positron's journey is very short as t 0 ∼ 10 −8 s. It will be ignored in all our discussions which follow.
The portion of the positron's journey which plays the key role in our discussions is the motion with constant acceleration eE /m and the velocity close to c. The corresponding motion is determined by the formulae: where E 0 (mc 2 + ΔE) is the total initial energy of positrons at the moment of liberation.
The most important observation here is that the positrons move with ultra-relativistic velocities for most of the time, and the displacement in transverse direction (with respect to the orientation of the electric field) is very modest such that the majority of the particles remain in the system and continue their acceleration for the entire time interval τ E ∼ 0.3 µs determined by parameter l a according to (2). At this point most of the particles assume very high energy close to 10 MeV according to Eq. (4) while transverse momentum distribution is characterized by (v ⊥ /c) 0.1. The elastic scattering of the positrons off the atmospheric molecules during this short journey plays very minor role as the kinetic energy of positrons being in keV range at the initial instant assumes MeV energies very quickly such that the corresponding cross section never becomes sufficiently large to modify the positron's dynamics as computed above in vacuum.

Intensity and the angular distribution
Our next task is to present conventional formulae for the intensity and the angular distribution of the radio emission due to the relativistic positrons (7) accelerating in electric field E with parameters (2). As we discussed above the velocity v of the positrons is mostly oriented along the electric field while v ⊥ is small. Therefore, we assume that E v to simplify our formulae which follow. In this simplified geometry the angular distribution for the intensity assumes the following form, see e.g. [37,38]: where N is the number of coherent positrons participating in the radio wave emission, to be estimated below. The angle θ in this expression is defined as usual as the angle between the velocity v of the charged particle at the moment of emission and the direction n of the observer, i.e. v · n = v cos θ .
Formula (8) explicitly shows that the emission mostly occurs along E direction as E v. This is precisely the same direction where burst events are recorded as shown on Fig. 1.
Integration over all angles leads to well known expression for the total intensity: Formula (9) represents the well known feature of the emission that the intensity is strongly enhanced for ultra-relativistic particles, which is obviously the case for positrons accelerated up to 10 MeV energies according to Eq. (4). This implies that the relativistic enhancement factor entering Eq. (9) is enormous: (E exit /m) 6 ∼ 10 8 .

Spectral characteristics of the radio emission
The main goal of this subsection is to understand the spectral characteristics of the radio emission. It will also allow to estimate the coherence factor N entering Eqs. (8) and (9) as this factor obviously depends on the frequency of radiation and corresponding wave length. The starting point for these studies is spectral property of the electric field E ω emitted by the accelerating positrons, see e.g. [37,38]: where all quantities at the right hand side of Eq. (10) must be computed as the retarded times t : where we assumed that the positrons move with approximately constant time-independent velocity |v| ≈ c. This assumption allows us to represent the Fourier component E ω in the following form where ω and a ω are defined as follows which precisely the combination entering the Fourier transform (10). Using Fourier component for E ω as given by Eq. (12) and magnetic component B ω = i(k ω × E ω ) one can compute the energy d E nω emitted into solid angle dΩ with frequency interval dω [37,38]: One can simplify this expression by assuming that the a v as we previously discussed. Furthermore, one can carry out the integration over dΩ by integrating over dω for a given ω as these two variables are related according to (13). Indeed, the emission is mostly concentrated along the cone with where we expanded cos θ ≈ (1−θ 2 /2) and represented small factor (1−v/c) ≈ 1/(2γ 2 ) in terms of the conventional combination γ which is valid approximation for ultra-relativistic positrons with |v| c. Now, integration over angles dΩ can be replaced by integration over dω as these two variables are related by (15). Therefore, such that the energy d E ω emitted within the frequency interval dω assumes the form where the low limit of integration is determined by θ = 0 in relation (15). The same expression can be thought as the angular distribution of the emission where |a ω | should be expressed in terms of θ according to relation (15). Our next task is to model the acceleration a(t) for a typical thunderstorm electric field with parameters (2). The corresponding Fourier transform (13) determines a ω which enters expression for d E ω /dω as given by Eq. (17). Our simplest possible choice is as follows: while a(t) ≈ 0 for t being outside of this interval. The corresponding expression for a ω defined by (13) assumes the form Now we can perform the integration dω in Eq. (17). For our simple model for acceleration (19) the integral can be approximately computed for small (ω τ E ) 1. This leads us to the order of magnitude estimate for d E ω /dω: The same expression can be derived from Eq. (18) by integrating over the angles dΩ. This spectral density holds as long as ω is sufficiently small. To be more precise: which represents the dominant contribution to the integral (17) for small ω. The emission with higher frequencies will be power-suppressed as ω −2 . Few comments are in order. The spectral density (21) integrated over all angles dΩ approximately constant and does not depend on ω as long as condition (22) is satisfied. This is of course a well-known feature of the constant acceleration (19) which we used as a simplified model for the electric field. In reality, the electric field under thunderstorm obviously fluctuates in time and space. This will obviously modify the spectral features (21). However, we expect that a typical frequency of the emission as given by (22) will hold because it is basically determined by the typical time scale of the problem τ E and expected relativistic factor γ which is not very sensitive to the details of the fluctuating electric field E .
The crucial observation here is that the typical frequency of emission is not simply given by the scale τ −1 E as one could naively expect based on dimensional arguments. Rather, the emission extends to much broader region due to the large relativistic factor γ 2 as inequality (22) states.
We postpone for the detail numerical estimates to Sect. 5. Now we want to make few simple numerical estimates supporting the main claim of this work that the frequency of emission is in the radio band ν ∈ (0.5-200) MHz. Furthermore, the emission is mostly oriented along the same direction where TA bursts are recorded. Therefore, one should anticipate a strong synchronization between the TA bursts and the radio pulses as both emission occurs at the same location at the same instant, and propagate to the corresponding detectors with the same speed of light. These features are not very sensitive to details of the dynamics of the AQNs, nor specific features of electric field under thunderstorm. Rather, all these features are inevitable consequences of the basic AQN framework along with pure geometrical properties of these observables.
With these comments in mind we can represent the typical band as follows which we expect to hold irrespective of specific features of the AQN model as it is entirely determined by well established typical time scales of the electric field under the thunderstorm (2). One could naively think that the typical frequency of emission could be extrapolated to very low ν as the spectral density (21) apparently does not depend on frequency. This is not quite correct conclusion though as we shall discuss in next Sect. 5.

Radio pulse of the electric field
The goal of this section is to estimate the intensity of the electric field (10) at very large distances R where it can be potentially detected. The orientation of E field is determined by cross product (10) where one can assume (for the simplicity of the numerical estimates) that v a, similar to our previous estimates for the spectral intensity of the emission (17). The absolute value for |E| at large distances can be estimated as follows where R is the distance from the emission site to the area where the electric field |E| could be recorded. It should not be confused with parameter r which enters all formulae from Sect. 3 and describes the distance from the emission site to the SD site where the energetic positrons could be detected. Numerically these parameters are similar, of course. However, the radio signal can be observed in a different location from the local area where the positrons hit the TASD. This formula explicitly shows that the time duration of the radio pulse and its properties are unambiguously determined by the features of the electric field E under thunderstorm. For simple model (19) this implies that the radio pulse lasts τ E while the absolute value of the pulse is estimated as It is important to emphasize that the duration of the pulse τ E is not very sensitive to the details of the AQN model as it is entirely determined by well established typical time scales of the electric field under the thunderstorm (2). This feature is very similar to our previous arguments regarding a typical frequency of the radio emission (23). At the same time the intensity of the pulse is highly sensitive to the details of the AQN model as it depends on the number of positrons N participating in the emission. In this respect this feature of insensitivity to any specific details of the AQN model for the pulse duration τ E and the frequency ν as given by (23) is similar to our studies of the TA burst events when such features as the "curvature", the timing and the spatial spread of the bursts are not sensitive to the details of the AQN model, but entire determined by the geometry as reviewed in Sect. 3. This should be contrasted with estimations [3] of the intensity of the TA bursts which are highly sensitive to specific features of the model, which are hard to carry out in a quantitative manner.

Numerical estimates
Our goal here is to present some numerical results using parameters which had been used in our previous estimations [3] related to puzzling TASD bursts as reviewed in Sect. 3. One should emphasize from the very beginning that the corresponding estimates suffer from huge uncertainties, and should be considered as the order of magnitude estimates, at the very best. Our main arguments of the present work are not based on these highly model dependent estimates. Rather, our main arguments are based on essentially model-independent features of the radio emission such as typical frequency band (23) and strong synchronization between TASD burst events and the radio emission. Nevertheless, we think that an order of magnitude estimates for the amplitude of the electric field as presented below could be useful as they demonstrate the consistency of the proposal. The same estimates also demonstrate the principle feasibility to detect such radio signals.
We start with numerical estimation of the strength of the electric field of the pulse (25) in conventional units (mV/m): where factor N is the number of the coherent positrons participating in the emission, to be estimated below. The orientation of the electric field is determined by the cross product as given by (10), and it is approximately (up to small angle θ ) points along a which represents the direction of the electric field under thundercloud at the instant of emission. The temporal shape of the pulse (bipolar, unimodal or even more complicated form) is also determined by the same fluctuating electric field E at the moment of emission. The intensity of the field is strongly peaks along n with typical angle θγ 1, see Fig. 2 for a precise angular distribution of |E(θ )|. Our next task is the estimation of parameter N entering Eq. (26). Our estimates will be entirely based on the observations, rather than on theoretical computations which inevitably suffer from huge uncertainties as reviewed in Sect. 3. Therefore, the logic of our estimation of parameter N is as follows. Let us assume that the TASD bursts which have been recorded are indeed due to the positrons emitted by AQNs as proposed in [3]. In this case the number of detected particles must be proportional to the density of particle detectors. This number is known for presently operating 507 SD detectors with area 3 m 2 each. These detectors are covering 680 km 2 total area. The number of detected particles must be also proportional to the area ∼ π(r Δα) 2 where particles had been recorded by TASD during a single event within the burst. Therefore, the number of positrons N at the emission site can be estimated from the following relation: In our estimate (27) the number of detected particles being recorded over the surface area π(Δs) 2 [left hand side of (27)] varies for different events within a single burst, and it normally varies in the range (2-6) ·10 2 . In our estimate (27) we inserted the suppression factor exp(−r/λ) ∼ 0.1 introduced in [3] to account for some attenuation of the positrons travelling a distance r ∼ 10 km or so when the mean free path λ for positrons with few MeV energy is order of kilometre at the sea level and several kilometres at higher altitudes. One should also note that π(r Δα) 2 on the right hand side of Eq. (27) is not identically the surface area ∼ π(Δs) 2 where particles are being recorded. These areas are directly related to each other through angle α according to Eq. (5), see Fig. 1 for notations. The relation (27) suggests that number N entering Eq. (26) for (r Δα) ∼ 1 km (which corresponds to a typically observed spatial spread for the burst events) can be estimated as N ≈ (0.3-1) · 10 9 for (r Δα) ∼ 1 km. (28) This estimate for N corresponds to the amplitude of the electric field (26) on the level 20 (mV/m) measured at distance R 10 km from the source of the emission within the angle range (θ γ ) 1 where the most of the radio emission occurs. Our last task in this section is establishing the lower frequency bound for the radio emission. The spectral density  (21) which is approximately a constant, naively suggests that the spectrum extends to arbitrary low frequencies. In fact, it is a premature conclusion as formula (21) was derived assuming a constant acceleration. In reality, a typical acceleration time is determined by (2) such that one should expect that ν (2πτ E ) −1 ≈ 0.5 MHz. Therefore, the AQN-induced ∼ 0.3µs pulse represents a radio emission in the bandwidth ν ∈ (0.  MHz with the amplitude of the corresponding electric field of the order |E(θ, t)| ∼ 20 (mV/m) at distance R ∼ 10 km as estimated above. Important feature of this spectrum is that it is very broad. Indeed, the intensity of the radiation changes by factor two or so when the frequency varies by two orders of magnitude. Another comment which follows from Fig. 3 is that the portion of the positron's energy being converted into the radio waves is very tiny as initial positrons energy can be estimated as 10N MeV when N is estimated in Eq. (28). This implies that the initial energy is many orders of magnitude greater than the radio wave energy shown on Fig. 3.
The crucial point of our analysis can be formulated as follows. This radio pulse must be synchronized with an event within TA burst within 10 µs because the TA burst and the radio emission are originated from the same location emitted at the same instant and both signals propagate with the speed of light along the same direction to be detected within ∼ km 2 area.
It is important to emphasize that such strong synchronization must hold even if the TA bursts are not related to any lightning events and can be easily discriminated from much stronger radio emission which always accompany lightnings. Therefore, such synchronization between the radio pulses and the TA bursts, if observed, can be easily distinguished from any other possible correlation between radio emission and conventional lightning events, see Sect. 8 for the details.

Geosynchrotron radiation
Another significant radio emission is the geosynchrotron radiation which is due to presence of the Earth's magnetic field B ∼ 0.5 gauss. This effect becomes important near the end of radio emission under thunderstorm. At this instant all the positrons are ultrarelativistic (γ ∼ 20) but they are no longer accelerating as they already left the region of the external electric field E . Trajectory of positrons are helical in presence of the magnetic field of Earth B ∼ 0.5 gauss. The acceleration in this case is perpendicular to the direction of velocity v ⊥ a B , so it is convenient to set up the following coordinate with respect to the direction of observer n: v · n = cos θ,â B · n = sin θ cos φ (29) Radius of the helical trajectory can be estimated as where θ B is defined to be the angle between direction of B and the velocity of positrons v. Radiation from an ultrarelativistic particle has a narrow emission angle θ γ −1 , similar to our previous discussions. However, the pulse of radiation for geosynchrotron radiation is much shorter [37,38]: This interval is much shorter than the time scale τ E ∼ 0.3 µs discussed previously due to extra suppression factor ∼ γ −3 . The energy spectrum integrated over angle dΩ is well known [37][38][39]: where K 5/3 (x) is the modified Bessel function, and ω c is the critical angular frequency beyond which radiation becomes negligible. 32 The electric field of the geosynchrotron radiation follows from Eq. (10): The estimation of electric field strength holds as long as the positron current remains coherent. The coherent length l coh can be estimated in what follows. The dispersion in velocity is δv ∼ v ⊥ 0.1c from Eq. (5). The condition of coherence requires the dispersion length δl to be much less than the radius ρ of trajectory: It implies l coh ∼ ρ that is insensitive to details of local parameters such as γ and B, and the geosynchrotron radiation is only significant within the first cycle of rotation.

Numerical estimates on geosynchrontron radiation
Similar to Sect. 5, we make numerical estimation for the electric field strength and the frequency bandwidth. Numerical estimation for the strength of the electric field of pulse Eq. (33) in conventional units (mV/m) is given as: see Fig. 4 for a precise angular distribution with specific parameters. The upper limit of frequency band is determined by the critical frequency ω c : and the lower frequency cutoff is given by the fundamental frequency: Therefore, the geosynchrotron induces a radio pulse with time duration ∼ 0.14 ns in the bandwidth ν ∈ (0.05-600) MHz (see Fig. 5) and the amplitude of the corresponding electric field is of order ∼ 50 mV/m at distance R ∼ 10 km.
While numerical values of the frequency bandwidth ν ∼ (0.1-10 2 ) MHz and the maximal strength of electric field of the radio wave [|E| ∼ (20-50) mV/m] are similar in both cases of radio pulses (which we term as E -type and B-type correspondingly) there is an important distinct feature, which is manifested in dramatically different time durations in these two cases. To be more precise, a B-type pulse induced by geosynchrotron radiation (τ B ∼ 0.14 ns) is much shorter than E -type pulse with (τ E ∼ 0.3 µs).
This distinct feature is in fact an absolutely crucial element for possible potential observations of the emission. This difference between τ B and τ E is the main reason for dramatically different energy injection which is measured in units MeV · MHz −1 in case of E -type emission to be contrasted with much lower keV · MHz −1 scale in case of B-type emission, see Figs. 3 and 5 correspondingly.
Therefore, the total energy of the radio emission integrated over the pulse time is much greater for the E -type radiation in comparison with B-type radiation. Furthermore, the required resolution to detect the very short B-type pulse with τ B ∼ ns requires much more demanding instruments than recording of the conventional E -type pulses with duration time of order τ E ∼ µs. These arguments strongly suggest that E -type radiation plays the dominant role in emission processes. Therefore, any future studies (including possible observations) should be focused on E -type radiation analyzed in Sect. 5.

Can conventional physics be responsible for proposed effects?
The question we address in this section can be formulated as follows. It has been known for quite sometime that the thunderstorm lightnings are always accompanied by the radio emissions, see e.g. reviews [30,31]. Can these lightninginduced radio pulses be responsible for the effects studied in this work? Before we answer this question we first recall [3] the basic features of the TA bursts interpreted as the AQN events, see item 1. Then, in item 2 we overview the main features of the radio emission which always accompanies the thunderstorm lightning events. After that in item 3 we explain the dramatic qualitative differences between the radio pulses associated with thunderstorm lightnings and the radio pulses studied in this work.

1.
As we previously mentioned, lightnings are correlated with AQN events when an AQN hits the region of the thunderclouds and potentially may initiate lightnings serving as a trigger [3]. The mechanism is similar to the role of conventional CRs that are considered to be one of the dominant triggers of lightnings, see footnote 4 with related comments and references. This is the basic reason why the bursts being interpreted as random AQN events are nevertheless associated with the lightnings.
Further to this point, it is important to emphasize that in the AQN framework the direction of the electric field E characterized by (2) (also shown on Fig. 1) and the direction of the lightning current are not correlated. This is because the AQN event may serve as a trigger for the lightning strike. The trigger should not be confused with the initial stepped lightning leader which determines the direction of the current of the lightning being developed. In different words, the direction of the observed current of a lightning does not determine the sign of the fluctuating electric field E which itself defines the direction of the positrons to be recorded as TASD unusual events.
To rephrase it, an AQN plays the dual role when it propagates in the thunderclouds: it emits the very energetic positrons, and it also triggers the lightning, similar to conventional CR, see footnote 4 with relevant references. The former phenomenon is recorded as a TASD burst, while the latter one explains the observed correlation between the bursts and lightnings.
As a next remark: as we mentioned, all mysterious bursts occur at the very initial moment of the lightning flashes or even earlier than lightning. Some bursts are not related to the lightning flashes at all (but occur under the thunderclouds). In original interpretation [1] this feature is explained in terms of the initial stepped lightning leaders propagating to the ground on a time scale of ten or tens of milliseconds. In the AQN framework this feature is also automatically satisfied as the AQNs can serve as the triggers of lightning flashes, which may or may not initiate the lightnings. The particles which are recorded as TASD events are produced directly by AQN itself at this moment, not at the initial stage of the strike determined by the lightning leader.
2. Having discussed the basic features of the TASD events within AQN framework we now turn to our next item on radio emission features as a result of lightnings. The correlation between lightning and radio emission during thunderstorms is well known and well documented generic feature of the the lightning discharges, see [40] with large number of references on observations. In particular, the lightning discharges are characterized by a very large number of radio pulses which last in total for about 1 s. Each pulse is characterized by full width (0.2-0.3) µs with electric field strength which could be as large as |E| ∼ 10 3 mV/m. Most of the pulses, though, show the strength of the electric field in the |E| ∼ (100-200) mV/m range. Another important feature of the radio emission: the gaps between pulses are in the range (10 − −10 2 ) µs. Therefore, total number of pulses could be very large 10 3 during a single lightning event. Finally, the typical frequency of the radiation is strongly peaked in few MHz bands, while it completely diminishes for ν 10 MHz.
3. These features must be contrasted with the AQNinduced radio pulses studied in this work. While a typical duration of an individual pulse τ E and the strength of the corresponding electric field |E| assume similar orders of magnitude in both cases, some other key characteristics are dramatically different. It allows an easy discrimination between these two cases.
Indeed, the number of clustered radio pulses associated with the mysterious bursts must be very few (it ought to be 3+ corresponding to the number of clustered events in a single burst). It should be contrasted with ∼ 10 3 in case of the lightning-induced radio pulses. The total durations of the radio emissions are also very different for these two cases.
The most important distinct feature which discriminates the different sources of the emission is that the frequency bands of the radiation are dramatically different in these two cases. The lightning-induced radio emission is strongly peaked in few MHz bands, while AQN-induced radio pulse is characterized by the flat spectrum with ν 200 MHz according to (23). The basic reason for this dramatic difference is that the electric current responsible for the lightning is represented by the particles with γ ∼ 1 while for the AQNinduced case the positrons are characterized by γ ∼ 20. This difference is translated into dramatic modification of the frequency bands according to (23)  Next, we want to mention that there are some dramatic differences between the radio pulses induced by conventional CR showers [41,42] and the AQN-induced radio pulses. These differences have been discussed in details in [28]. Here we want to mention that these dramatic differences are related to very different structures of the showers. In case of conventional CR the radio emission is based on the picture when a CR air shower is characterized by the "pancake" of particles and its central axis. The density of the particles strongly depend on the distance to the central axis such that the spectral density of the radio pulse is highly sensitive to the width of the "pancake", which becomes much thicker further away from the axis. It is very different from the AQN-induced radio signal as the notions of the shower axis and the "pancake" do not exist in our case such that the density of the particles is approximately the same irrespective of the distance to the central axis.
Additionally, conventional CRs are well studied even being heavily distorted by electric field in thundercloud, and in fact it is nowadays an efficient probe to study the electric field in thunderstorms by analysing the radio emission [43]. 5 The dramatic difference with the radio pulse studied in the present work is that the radio emission due to the CR is strongly correlated with the Earth's magnetic field, in contrast with dominant E -type emission studied in the present work. As a result the width of the corresponding signal is in nanosecond range, in contrast with the AQN induced pulse when the width is ∼ 0.3 µs. Furthermore, the polarization pattern for the pulses related to CRs are dramatically different from the AQN-induced emission. Therefore, there should be no technical difficulty to distinguish an AQN-induced radio pulse from a conventional CR under the thunderstorm. The corresponding studies could support or refute our proposal.
There are many other distinct features between the AQN and conventional CR air showers. In particular, the AQNs mostly emit photons in X ray bands, such that a signal cannot be observed by a fluorescence detector which is designed to detect the visible and UV light. 6 Next, we note that gamma rays can be produced by an avalanche of relativistic runaway electrons during initiation of lightning [46]. Bursts of gamma rays initiated in the Earth's atmosphere, commonly referred to as terrestrial gamma ray flashes, can be generated in a lightning leader system [47,48]. The terrestrial gamma ray flashes can trigger atmospheric photonuclear reactions that produce neutrons and positrons [49]. Evidently these gamma rays and their byproducts do not interfere with our proposed radio detection, as they are excessively more energetic and only appear at frequency band above EHz. Another important point here is that this powerful emission occurs in a latter stage of a lightning, while an AQN-induced event appears before or at the initial stage as it serves as a potential trigger of a lightning strike.
Based on these dramatic differences in frequencies and timings of the radio emissions we conclude that the signals due to the thunderstorm lightning events and conventional CR air showers can be easily discriminated from the AQNinduced radio pulses studied in this work. Therefore, we suggest to study the corresponding radio signal to support or refute our proposal to interpret the TA bursts as the AQNinduced events.
Lastly, one can assume that the source of the TASD bursts is entirely due to the complicated and not well understood physics of the lightning strikes (which may or may not be recorded by the system on the surface). In addition to many problems with explanation of the observed intensity, timing, clustering features, and geometry, there is an additional problem of the very low observed event rate (which is 10) in comparison with recorded ∼ 10 4 number of lightnings in the same area. At the same time the estimated event rate [3] based on assumption that the observed TASD bursts is a result of the AQN-induced events is consistent with observations.
To conclude: we are not aware of any studies which could explain the observed intensity, timing, clustering features, and geometry of the signal in form of the mysterious bursts being produced during initial stage of the lightning as recorded, while proposal [3] naturally explains all these features within AQN framework.

Conclusion
As we stated in the Introduction the main goal of the present work is to test the proposal [3] by searching for the radio signals in frequency band ν ∈ (0.5-200) MHz which must be synchronized with the TA bursts. Such test would unambiguously support or refute the proposal (interpreting TA bursts as the AQN events) as the radio signals due to the AQN annihilation events can be easily discriminated from conventional radio pulses which always accompany the thunderstorm. This is precisely the goal of this work: we want to eliminate all common objections which essentially state that a thunderstorm is very complicated system 7 such that everything is possible, including TA bursts as a result of flashes. The proposal of the present work is to study the radio signals in frequency band ν ∈ (0.5-200) MHz. The corresponding results would unambiguously answer these rhetoric questions, see Sect. 8 with more comments on this. Now we summarize the results of our studies. As we argued at the very end of previous Sect. 7 the E -type radiation plays the dominant role in radio emission. The corresponding results have been summarized at the very end of Sect. 5 where it has been argued that the radio pulse can be recorded if proper instruments are designed and built at the TASD site. The pulse must be synchronized with TASD bursts. Such a signal cannot be confused with any other spurious and noise signals as a result of this synchronization. Therefore, observing (not observing) such synchronized signals can confirm, substantiate or refute our proposal. We further argued in Sect. 8 that the AQN-induced radio signals can be easily discriminated from the pulses generated by the thunderstorm lightning events.
One should comment here that the strength (26) of the amplitude |E| on the level 20 mV/m with the duration time of the radio pulse on the scale ∼ 0.3 µs could be recorded with existing technology as discussed in [40] in application to thunderstorm lightning events where wide-band radio interferometry (0.1-30 MHz) has been used to detect such short signals with amplitudes of electric field on the level 10 2 mV/m, which is very close to what is required for the purposes of the present work. Therefore, the predicted radio signals (26) can be in principle measured if proper instruments are designed and built at the TASD site.
Furthermore, the short radio pulses can be, in principle, recorded by any (sufficiently sensitive) radio telescope outside of the TASD site as the corresponding radio pulses will be emitted in form of the clustering events, similar to TAbursts. Such clusters of individual short radio pulses can be discriminated from any spurious signals representing the radio noise. It can be also discriminated from radio emission occurring as a result of the thunderstorm lightning events as argued in Sect. 8.
Our present proposal suggests that the TA bursts with very unusual features will be synchronized with radio pulses. If this synchronization is observed and the interpretation of TA bursts as the AQN annihilation events is confirmed by future studies it would be the direct (non-gravitational) evidence which reveals the nature of the DM, in contrast with large number of indirect hints mentioned in Sect. 2. distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .