Detection of ultra high energy cosmic rays and neutrinos with lunar orbital radio telescope

Particle cascades induced by ultra-high-energy (UHE) cosmic rays and neutrinos impacting on the lunar regolith usually radiate Cherenkov radio emissions due to the presence of excess negative charge, which is known as Askaryan effect. Several experiments have been carried out to detect the Cherenkov radio emissions in the lunar regolith. To prepare for future lunar Ultra-Long Wavelength (ULW, frequencies below 30 MHz) radio astronomy missions, we study the detection of the Cherenkov radio emissions with the ULW radio telescope that are operating at the lunar orbit. We have carried out instrument modelling and analytic calculations for the analysis of aperture, flux and event rate, and the analyses show the detectability of the Cherenkov radiation. Based on the properties of the Cherenkov radiation, we have demonstrated that the cosmic ray and neutrino events could be reconstructed with the three ULW vector antennas onboard the lunar satellites via measurements of the Askaryan radio pulse intensity, polarizations, etc. The results obtained by this study would be useful for future lunar radio explorer mission, where the detections of UHE cosmic rays and neutrinos could be successfully attempted.


Introduction
Studies on the behaviour of ultra-high energy (UHE) particles and their source of origins have been a major topic of interest in modern high energy astrophysics and particle physics.However, the origin of ultra-high-energy (> 10 18 eV) cosmic rays (UHECR) still remains unknown to date.It is mainly because of their deflected trajectories caused by cosmic magnetic fields.Above a threshold energy of ∼ 10 19.6 eV, cosmic rays usually interact with cosmic microwave background to produce pions and lose a substantial amount of energy within a distance of tens of mega parsec (Mpc).This is well known as Greisen-Zatsepin-Kuzmin (GZK) effect ( [28], [54]).The few cosmic rays that observed above this threshold energy must originate from the distant sources.In order to determine the source origin of such UHE cosmic rays, the detection of UHE neutrinos, those produced together with the cosmic rays or during the GZK interaction, has usually been employed as one of the method.Since neutrinos are chargeless and weakly interacting particles, the propagating UHE neutrinos (UHECν) usually remain unaffected over large cosmic distances, and therefore, their arrival directions carry direct information of their source of origins.
The chief problem of both UHECR and UHECν event detection come from their extreme rarity.The flux at GZK energy is around one particle per square kilometer per century, which calls for huge detectors (thousands km 2 ) to collect a significant amount of events at ultra-high energies.For instance, even the Pierre Auger Observatory (PAO) with a collecting area of 3000 km 2 can only detect of order 30 particles above the GZK energy per year [37].
However, for the detection of particles above 10 21 eV, the detectors are needed to have a thousand-fold increase in their collecting area.Thus, space-based detection systems are, therefore, appeared to be a feasible approach that is currently being attempted and implemented.
The Moon acts as a huge natural detector with large collecting area and is first proposed to be a target to detect cosmic ray and neutrino showers in [20].The detection is based on the Askaryan effect, which is pointed out by G.A. Askaryan from the Lebedev Physical Institute in 1961 [5].When high-energy particle interaction occurs in a denser medium like ice, rock salt, lunar regolith, or the atmosphere, it attracts electrons from the surrounding medium and transfer them into the shower disk.With the annihilation of showering positrons in flight, a net excess of electrons are produced.The excess electrons rapidly appear and disappear to produce a radio pulse that could be observed by a sensitive telescope.Using the above technique several experiments have been proposed and performed to detect high-energy cosmic particles, such as Parkes ( [32], [33]), GLUE ( [26], [25]), LUNASKA [38], [39], NuMoon [13], RESUN [36], LaLuna [52], LOFAR [50], FAST [37], all of which employ the terrestrial telescopes to carry out the observations.In the literatures ( [31], [53], [4]), the scheme of a lunar orbital satellite system was presented and analyzed to detect the UHE cosmic rays and neutrinos.
The ultra-long wavelength (ULW) (≥ 10 m) range, the last unexplored region of the electromagnetic (EM) spectrum, is presently one of the new promising areas of radio astronomy.Strengthening the science cases for extending the radio astronomy into the unexplored ULW window has been significantly increasing.Due to the reflection of the ULW signals by the Earth's ionosphere and the presence of the strong terrestrial anthropogenic radio frequency interference (RFI), the far side of Moon can be an ideal site for carying out the ULW radio observations.To date, various Moon-based telescopes have been made and launched for exploring the ULW region, such as OLFAR [22], FARSIDE [47], DARE [14], LRX [44], DSL ( [12], [18]), NCLE [11], LFRS [41], etc.The lunar ULW radio telescopes provide ample opportunities to detect the Cherenkov radiation usually caused by the UHNCR or UHECν on the Moon.
In this paper we present the radio detection of the UHNCR or UHECν with the lunar orbital ULW radio antennas.The paper is organized as follows.In Section 2 we describe the model for the detection of the lunar Cherenkov emission.The aperture and flux limit for the lunar UHECR and UHEC events are estimated and analyzed in Section 3. In Section 4 we discuss the characteristic of the UHNCR or UHECν with the multi-satellite detections.Finally, summary and conclusions are provided in Section 5.

Model for radio emission detection
When energetic particles interact in the denser medium like lunar regolith, intense coherent radio pulses are produced in the frequency range from MHz to GHz.The radio emissions that escaped from the lunar surface could be detected by a Moon-based or terrestrial radio telescope as shown in Figure 1, which strongly depends on the sensitivity of the telescope.Though the terrestrial radio telescope have a huge collecting area and the best sensitivity, the intensity of the lunar Askaryan radio emission is weak due to the long-distance attenuation.A lunar radio telescope may provide a reasonable detection capability for the UHECR and UHECν events since it is thousand times closer to the lunar surface than to the ground-based telescopes.For the detection of the UHE particles, the effective aperture is another vital factor which is directly dependent on the physical lunar surface illuminated by the antenna.The physical surface area is determined by the antenna beamwidth and antenna altitude from the lunar surface, and the altitude affects the detected electric field of the Askaryan radio emission.In order to study the lunar UHECR and UHECν detections, the electromagnetic properties of Cherenkov radiation have been modelled including the system parameters of the radio detector.
where, R is the distance from the source point to the detector, E s is the shower energy in eV.E 0 = 10 20 eV, ν 0 = 2.32GHz ( [39], [9]).The angular distribution of the radio emissions around the Cherenkov angle θ c is crucial for the detection of lunar UHNCR and UHECν events.Due to the loss of coherence, the radio emission decreases in the direction away from θ c .[48] discussed the angular spread of the Cherenkov radiation around the Cherenkov angle.Through comparing several different formulae of intensity distribution, Scholten concluded that the gaussian parametrization formula is more accurate since it agrees with both the Monte Carlo simulations at high frequencies and the analytic results at low frequencies.In the present study, we use the formula given in [48] for angular spread of the Cherenkov radiation as follows, The spreading of the radiation intensity is determined by ∆ c = 0.0754[L(E 0 )/(νL(E s )] (in radians), which is inversely proportional to the shower length and the radiation frequency.Here L(E is the shower length determined by the energy, L(E) = 12.7 + 2/3 log(E/E 0 ) [48].
In the detection analysis of lunar UHE particle events, a point-source is treated for the coherent Cherenkov radiation.When the Cherenkov radio emissions escape from the lunar surface, it will be reflected and transmitted at the interface.According the Fresnel law, the transmission coefficient for parallel polarization can be expressed as where θ s is the polar angle which is uniquely related to the angle r and i. r is the angle of refraction relative to the normal (outside the Moon) as the rays pass through the lunar surface into free space, i is the angle of incidence as shown in Figure 1.r and i have the relation as sin(i) = sin(r)/n, n is the refraction index of lunar regolith, and chosen to be 1.73.The detected radiation from a lunar regolith shower, with energy E, at a frequency ν and an angle θ, can be expressed as Here the absorption of the radio waves before exiting the lunar surface has not been discussed, and it will be taken into account in the aperture calculations.

Detection with lunar ULW radio antenna
The radio emission from the hadronic cascades induced by the UHE cosmic rays or neutrinos in the lunar regolith covers a broad frequency band from MHz to GHz, and it peaks in the GHz regime [39].Therefore, most of the past lunar radio experiments to UHECR or UHECν operate at GHz frequency band ( [32], [26], [38], [36]).With the decreasing frequency the peak intensity of the radiation decreases, however, the increasing angular spread allows the radio emission detected by the lunar radio telescope to escape from the lunar surface within a broader solid angle, which makes it increasingly efficient to detect the Cherenkov radiation at low frequencies.Furthermore, the surface roughness does not play an important role in the detection analysis at lower frequencies since a considerable fraction of the radiation will penetrate the surface.According to these factors, some experiments running at low frequencies(∼ 150 MHz) for the lunar particle detection have been carried out or planned with the Westerbork Synthesis Radio Telescope (WSRT), NuMoon [13], the LOw Frequency ARray (LOFAR), and the Square Kilometer Array (SKA-low).Although the ULW band ≤ 30MHz is beyond the optimum frequency window where the wavelength is comparable with the shower length [48], and the bandwidth is also limited at ULW band, the increased particle energy could well compensate the loss of electric field strength at lower frequency.The net effect is that the detection probability could increase reasonably at adequately high shower energies.
In the present study, we discuss the UHECR and UHECν detection with lunar orbital ULW radio telescope.
For the space ULW radio observations, a tripole antenna has been selected as the receiving element in most of the radio missions as DSL ( [12], [18]), NCLE [11], LFRS [41], and OLFAR [15], since it has specific advantages of measuring the three-dimension electric field, estimating the directions of arrival (DOA) of incident signals with single antenna ( [17], [16]), and protecting the desired signals from interference signals [19].In the lunar ULW mission, the Dark Ages is undoubtedly one of the greatly interesting sciences, it requires the observation of the entire 1 − 30 MHz band, which could be extend up to 50 MHz for cross-referencing with the terrestrial radio facilities such as LOFAR.Therefore, the tripole antenna of the lunar ULW radio telescope is supposed to operate at the frequency range up to 40 or 50 MHz ( [41], [12]).
For a radio telescope, its sensitivity can be determined by the system equivalent flux density (SEFD) F = kT sys /A ef f in Jansky (10 −26 Wm −2 Hz −1 ), .Here k is the Boltzmann's constant, T sys the system temperature and A eff the effective collecting area.However, a electric field strength (V/m/Hz) of a coherent radio pulse is usually utilized to characterize the Cherenkov radiation.In order to analyze the sensitivity of a radio telescope for the Askaryan pulse detection, it is required to express the telescope sensitivity in terms of the electric field strength.
If a radio telescope has a flat system noise spectrum over the observed bandwidth, the equivalent root mean square (RMS) spectral electric field [9] can be written as follows, where E rms is the RMS electric field strength over the bandwidth ∆ν in one polarization, and Z 0 is the impedance of free space.Then the sensitivity of an antenna to detect a coherent radio pulse (the minimum detectable electric field) can be defined as ε min = N σ ε rms , where N σ is the minimum number of standard deviations needed to reject statistical noise pulses.
At the ULW band, the galactic synchrotron emission dominates the sky radio foreground, its spectrum varies greatly with the frequencies, the brightness temperature rises from about 10 4 K at 30 MHz to about 2.6 × 10 7 K at around 1 MHz; Furthermore, the antenna effective collecting area changes significantly with the frequencies too.
Therefore, the Equation ( 4) cannot be applied directly to the electric field calculation for a ULW radio antenna.
In order to truly characterize the equivalent spectral electric field of the telescope, the RMS electric field strength E rms should be re-computed by integrating the SEFD over the frequency.According to the studies in Ref. [16], the sky noise temperature is absolutely dominant in the total system temperature for the ULW radio antennas of different lengths, then E rms can be re-written as, Where ν L and ν H are the lower and upper limit of the frequency respectively.Using a sky temperature model [40], the sensitivity of a ULW radio antenna longer than 5 meters [16] to detect the Askaryan radio pulse has been calculated with a nominal value N σ = 3 for different center frequencies and bandwidths, as shown in Table  that the antenna length has less influence on the sensitivity since the sky noise temperature is absolutely dominant.
In order to improve the detection probability of lunar UHECR and UHECν events, we use the central frequency of 30 MHz and the bandwidth of 40 MHz to obtain the best sensitivity in our analysis.

Detection analysis for lunar UHECR and UHECν events
We used the method presented in Ref. [31] to analyze the detection of lunar UHECR and UHECν events.In order to decide the cascade event rates for the UHECR and UHECν, one way is in term of the aperture size, which is multiplied by the isotropic cosmic ray or neutrino flux in the given energy bin to obtain the detection rate.We calculate the analytic aperture by integrating the angular aperture for specified energy E over the angles (θ s , ϕ s ) at the lunar surface S, and over the angles (θ n , ϕ n ) of the cascade axis relative to the normal to this surface, as shown in Figure 1.
For a volume element near the lunar surface (z ≃ 0) with the polar angle θ s , (see Figure 1) radio emissions from it received by the detector has a given refraction angle r, sin(r , is the distance between the surface element and the detector.The maximum value of the polar angle At low frequencies, the lunar surface roughness will not have a major influence on the detection sensitivity, and thus can be neglected [48], which could be another advantage of UHE particle detection with a lunar ULW radio antenna.In our analysis, the Moon is assumed to be an approximately smooth sphere.

UHE cosmic rays
For the UHE cosmic rays, they cannot penetrate through the lunar regolith, and cosmic ray events occur in the lunar regolith near the surface.Therefore, only the upper hemisphere (cos θ n > 0) contributes to the aperture calculation for the cosmic ray cascades.The total aperture A CR (E) is defined by Equation ( 6) and (7) as the integral of the angular aperture ∆Ω CR (E, θ s ) over the surface S [31].
where, A 0 = 2πR 2 M is the physical area of lunar hemisphere, the unit step function Θ(ε − ε rms ) only selects the cascade events that produce radio signals above the threshold value ε min .The radio attenuation length of lunar soil depends on the fraction of its composition.For the lunar regolith we use the value ℓ = 60λ (λ, the wavelength in the vacuum) given in the literatures ( [39], [9]), it is a quite large values at low frequencies.Thus the attenuation of radio wave prior to their escape from the lunar surface can be ignored since the cosmic ray cascades exist only few meters deep.The lunar subsurface soil consists of several different strata according to the detection of Lunar Penetrating Radar.Therefore, the reflected radio waves by the strata interfaces that can be detected by the lunar orbital radio antenna should also be taken into account besides the direct radio waves emitted by the hadronic shower, as shown in Figure 2. Due to the very small contribution, the secondary reflections at the interfaces are neglected.In order to simplify the analysis, the lunar soil structure is modelled as only two layers, lunar regolith and lunar sub-regolith.
The regolith is assumed to be 12 meter deep [45], and the refractive index of lunar regolith (n = 1.73) and subregolith (n = 2.5) ( [39], [9]) have been used in the aperture calculation.Using analytic integrations of the angular aperture over the Moon surface, the total aperture are computed for the lunar orbital ULW radio antenna with different altitudes.As shown in Figure 3, with the higher altitude of lunar ULW radio antenna, the physical aperture illuminated by the antenna increases, however, the detected Cherenkov radiation on the antenna becomes weaker due to the distance-decay effect.As a result, the minimum (threshold) cascade energy at which the radio emission intensity on the antenna is equal to the threshold intensity ε min increases when the altitude rises.For all the altitudes, the aperture reaches around the peak above the energy of ∼ 10 22 eV.
Based on the model-independent estimate [46], for the zero detected events, the upper limit for differential (in energy) flux J(E) of primary particles in the observation time T, can be obtained by the following relation, Here the Poisson factor S up = 2.3 for a limit with 90% confidence level.For one-year observation with the lunar orbital ULW radio antenna, the flux limits are calculated for different antenna altitudes.As shown in Figure 4, only the cosmic ray with the energy above a threshold can be detected by the lunar orbital antenna, it is 4 × 10 19 eV for the antenna at the 50 km high orbit, and 4 × 10 20 eV for the 500 km high orbit.The results in Figure 4 also show that the flux limit of single lunar orbital antenna will be competitive with the estimated limits of large terrestrial radio telescopes such as LOFAR [51] and SKA-Low [10].Note that the flux limit of SKA-Low is calculated according to the SKA array for Phase I.
For an isotropic cosmic ray or neutrino flux J(E), during the observation time T the number of detected events in the given energy interval between E 1 and E 2 can be obtained by N = E2 E1 T J(E)A(E)dE, where A(E) is the total aperture as a function of energy.
Using the flux parametrization of the Auger Collaboration in Ref. [1], the expected event rates of UHE cosmic rays are calculated in the energy range between 10 18 and 10 21 , as shown in Figure 5, for the lunar ULW radio antenna with different altitudes.It can be seen that the total detected cosmic ray events decreases for the antenna Figure 5: The number of detected (per year) cosmic ray events for the energy E ≥ 10 18 eV with the flux parametrization of the Auger Collaboration [1].Top, the total number detected events as a function of lunar antenna orbital altitude; Bottom, the detected events per energy interval for the lunar antenna of 50 km-high orbit (left) and 200 km-high orbit (right).with higher altitudes.The detected event number reaches to about 400 for the 10 km-high orbit, and it drops rapidly from tens at the altitude of 200 km to zeros at the attitude above 500 km.The simulation results of detected cosmic ray events per energy interval show that the events mostly lie in the higher energy range with the increasing of orbital altitude.

UHE neutrinos
The cascades initiated by the UHE neutrinos are different from that produced by the cosmic rays.Due to the long attenuation length for neutrinos in lunar regolith, most of the UHE neutrinos induce cascades very deep inside the Moon.Therefore, the attenuation of radio wave propagation prior to their escape from the lunar surface can not be neglected, and the contribution of both the upper and lower lunar hemisphere have to be taken into account while calculating the total aperture.Considering the neutrino absorption on the path L(z, θ n ), we integrate the angular aperture over z and obtain the contributions of the upper and lower hemisphere respectively [31] as follows, here, it is assumed that only 20% of the energy of the original neutrino goes into the hadronic particle cascade.ℓ is the radio wave attenuation length.Since the showers produced by the UHE neutrinos occur in the lunar soil over a large range of depths from the surface, we use the value ℓ = 29λ of lunar sub-regolith ( [39], [9]) for the aperture calculation.The neutrino-nucleon interaction length was L νN (km) = 122km(E 0 /E ν ) 1/3 , where E 0 = 10 20 eV ([23], [9]).For the upper hemisphere cos θ n > 0, we assume l/L νN ≪ 1 [31], which is valid for most of cases.
Figure 6: Direct and reflected radio emissions initiated by the UHE neutrinos propagate between the lunar soil and lunar orbital radio antenna.The lunar soil is modelled with two layers, the regolith and the sub-regolith.
For the UHE neutrinos, the multi-strata structure of the lunar subsurface soil has also been taken into account in the aperture calculation.Since the cascade showers generated by the UHE neutrinos mostly occur deep inside the lunar soil, that is in the lunar sub-regolith, we consider the refractions at both the lunar surface and the subsurface strata interface, as shown in Figure 6.Here, the reflected signals of the Askaryan radio emissions at the lunar surface is ignored to reduce the complexity of aperture analysis.Similar to the cosmic rays, the total aperture of the UHE neutrinos are calculated using analytic methods for the lunar orbital ULW radio antenna at different altitudes.As shown in Figure 7, the aperture increases with the rising of antenna attitude.However, the increase becomes slower above the attitude of 300 km.Being different from the cosmic rays, the aperture of the UHE neutrinos increase over the entire energy range from 10 20 eV to 10 24 eV without reaching peaks for all the four altitudes.For 50 km distance, the energy of detectable events begin at 2 × 10 20 eV.With the increasing of orbital attitudes, only higher energetic neutrino events become detectable, for 500 km attitude the threshold energy of detectable events increases to ∼ 2 × 10 21 eV.show the predicted fluxes from the GZK [21], and topological defects (TD) [49].The green dashed lines show the flux limits set by the past experiments of ANITA-II [24], IceCube [3] and Auger [2].
Figure 8 shows the flux limits for one-year observation with the lunar orbital antennas of different altitudes.
The flux limits of the SKA-Low for a 1000-hour observation, the LOFAR for a 90-day observation, the past experiments of ANITA-II [24], IceCube [3] and Auger [2] are also plotted for comparison, as well as the predicted fluxes from the models of the GZK [21] and topological defects (TD) [49].The results indicate that the lunar orbital observations in one year will set the comparable flux limits with the estimated fluxes of the SKA-Low and the LOFAR core (even lower than the energy of ∼ 3 × 10 22 eV), which signifies the higher detection efficiency of the lunar ULW radio antenna for the UHE neutrinos.
Figure 9: The total number of detected (per year) neutrino events with the flux limits set by the ANITA-II experiment, and the predicted fluxes for the GZK neutrinos [21] and topological defects (TD) [21], as a function of the lunar antenna altitude in the energy E ≥ 10 18 eV.
Based on the predicted fluxes from the models of the GZK neutrinos, topological defects (TD), the neutrino flux limit set by the experiment of ANITA-II, the expected event rates per year are estimated and presented in Figure 9 as a function of attitudes for the lunar orbital ULW radio antenna.It is evident from Figure 9 that the GZK neutrinos almost cannot be detected for all the altitudes.In order to increase the detection sensitivity of the GZK neutrinos, the aperture must be improved by increasing the effective receiving area of the radio telescope, for instance, employing a radio array.For the TD neutrino detection, the event rate depends on the altitude, and reaches the peak of about 70 per year at the attitude of about 100 km, which could be the optimum attitude for a lunar orbital ULW antenna.With the neutrino flux limit of ANITA-II, the expected neutrino events increase exponentially at the altitude range from 1 km to 1000 km, and thousands of events per year could be detected at most by the lunar ULW antenna.

Characteristic analysis of UHE cosmic rays and neutrinos
For the detection of the UHE cosmic rays and neutrinos, the ultimate goal is to reconstruct the events to obtain all the parameters characterizing the cosmic rays or neutrinos, including the source, the primary particle energy, etc.
In the literatures ([29], [30]), a random search method was employed to determine the UHECR event parameters for the LORD experiment under the assumption that the cosmic ray flux is known.Due to the lack of direct information, the reconstructed parameters have large uncertainties even for the detection with two satellites.Compared with the LORD experiment, the future lunar ULW radio missions like DSL ( [12], [18]), OLFAR [22] have more advantages in the UHE particle detection.They consist of many antennas onboard micro-satellites to form a large radio array, which make it feasible to detect the UHECR and UHECν with multi satellites simultaneously.Furthermore, in lunar ULW radio observations, the tripole antenna can measure the three-dimension electric fields, which make it possible to measure the three-dimension polarization of the Askaryan radio emission.In the detection of the UHECR and UHECν, if the Askaryan radio pulses can be observed by multi antennas on-board lunar satellites, the time delays can be measured between the radio pulses sensed by different antennas, and the particle cascade location (P) on the lunar surface can be solely determined with the time delays among at least three non-linearly distributed antennas, as shown in Figure 10.We know the Askaryan radio emission is highly linearly polarized, which is always in the same plane of the Poynting vector and the shower axis ( [37], [34], [25]).
Therefore, using the three-dimension polarizations of the electric fields (E 1 , E 2 , E 3 ) measured by the tripole antennas onboard three orbital satellites, the show axis can be determined by the intersection line of the three planes 10, which means the original direction of the UHECR or UHECν is known now.For a linear ULW radio array like DSL, it is necessary to note that the cascade location deduced by the time delays could have two candidates symmetric with respect to the satellite orbit plane.However, since only the right location can account with the fact that the shower axis, the Poynting vectors and the polarizations of the Askaryan radio emissions are in the same plane, the measured polarizations of the coherent radio pulses on three tripole antennas will definitely exclude one of the candidate locations except for the case that the shower axis is parallel to the satellite orbit plane.
Once the shower axis and the location of the UHE particle cascade are confirmed, the distance R s between the particle cascade and the radio detector, the three angles θ n , r, i shown in Figure 1 will also be known according to the simple geometrical relationship.Using the formulae (1),( 2) and (3), the particle energy E can be calculated with the intensities of electric fields induced on the lunar orbital antennas.The results of the three satellites will be cross-checked with each other to improve the accuracy and reliability.

Summary and conclusions
In the present study, using the lunar Askaryan technique, the detection of the UHE cosmic rays and neutrinos has been analytically analyzed for the future lunar ULW radio missions.Our results indicate that the single lunar ULW radio antenna could detect the cosmic ray and neutrino flux as low as that observed by the present or future experiments in the energy range E > 10 20 eV.With the known flux of the UHE cosmic rays, the simulated detectable events per year has found to increase when the antenna altitude is lowered, and it reaches the peak at the altitude of ∼ 10 km.For the UHE neutrinos, during the observations in one year, it is shown that the single lunar ULW radio antenna could detect ∼ 70 topological defect neutrinos at an optimal altitude of 100 km, while the GZK neutrinos could not be detected.
Investigating the properties of the lunar Cherenkov radiation, it has been demonstrated that the UHECR or UHECν events could be reconstructed directly using the radio observations with at least three tripole antennas onboard the lunar satellites that makes it feasible to detect the UHECR and UHECν with the lunar ULW radio mission.The method could also be used to build a dedicated lunar radio telescope for the detections of the UHE cosmic ray and the UHE neutrinos.where, z max = 2ℓ cos i ln(ε/ε min ).For the ultra-high energy neutrino, when ℓ/L νN ≪ 1, formula A.6 will be simplified to the formula (10).The angular aperture is then calculated numerically.

Figure 1 :
Figure 1: Schematic view of lunar antenna detections of Askaryan radio emissions produced by the UHECR or UHECν in the lunar soil.

Figure 2 :
Figure 2: Askaryan radio emissions initiated by UHE cosmic rays propagate between the lunar soil and lunar orbital radio antenna.The lunar soil is modelled with two layers, the regolith and the sub-regolith.

Figure 3 :
Figure 3: Total aperture of UHE cosmic rays for a lunar orbital ULW radio antenna at different altitudes of 50 km, 100 km, 300 km, and 500 km with a sensitivity of 1.34 uV/m/MHz.

Figure 4 :
Figure 4: UHE cosmic ray flux of a lunar orbital ULW radio antenna at different altitudes of 50 km, 100 km, 300 km, and 500 km for one-year observation.The results are compared with the flux limits of both SKA Low for a 1000-hour observation between 100 MHz and 350 MHz (reproduced from Bray et al., 2014) and LOFAR for a 90-day observation between 125 MHz and 175 MHz (reproduced from Singh et al., 2009).

Figure 7 :
Figure 7: Total aperture of the UHE neutrinos for a lunar orbital ULW radio antenna at different altitudes of 50 km, 100 km, 300 km, and 500 km with sensitivity of 1.34 uV/m/MHz.

Figure 8 :
Figure 8: The UHE neutrino flux of a lunar orbital ULW radio antenna at different altitudes of 50 km, 100 km, 300 km, and 500 km for one-year observation.The results are compared with the flux limits of the SKA Low for a 1000-hour observation between 100 and 350 MHz (reproduced from Bray et al., 2014), LOFAR for a 90-day observation between 125 MHz and 175 MHz (reproduced from Singh et al., 2009).The blue chain-dotted lines show the predicted fluxes from the GZK[21], and topological defects (TD)[49].The green dashed lines show the flux limits set by the past experiments of ANITA-II[24], IceCube[3] and Auger[2].

Figure 10 :
Figure 10: Scheme of Askaryan pulse detections in the lunar soil with three antennas onboard micro-satellites.

Figure 11 :
Figure 11: Angular aperture ∆Ω CR for the cosmic ray detection with a lunar orbital ULW radio antenna at different altitudes of 50km, 100km and 200km, and for the initial cosmic ray energy W = 10 21 eV .