• ### The 17 MeV Anomaly in Beryllium Decays and $U(1)$ Portal to Dark Matter(1609.07198)

April 13, 2018 hep-ph
The experiment of Krasznahorkay \textit{et al} observed the transition of a $\rm{^{8}Be}$ excited state to its ground state and accompanied by an emission of $e^{+}e^{-}$ pair with 17 MeV invariant mass. This 6.8$\sigma$ anomaly can be fitted by a new light gauge boson. We consider the new particle as a $U(1)$ gauge boson, $Z'$, which plays as a portal linking dark sector and visible sector. In particular, we study the new $U(1)$ gauge symmetry as a hidden or non-hidden group separately. The generic hidden $U(1)$ model, referred to as dark $Z$ model, is excluded by imposing various experimental constraints. On the other hand, a non-hidden $Z'$ is allowed due to additional interactions between $Z'$ and Standard Model fermions. We also study the implication of the dark matter direct search on such a scenario. We found the search for the DM-nucleon scattering excludes the range of DM mass above 500 MeV. However, the DM-electron scattering for MeV-scale DM is still allowed by current constraints for non-hidden $U(1)$ models. It is possible to test the underlying $U(1)$ portal model by the future Si and Ge detectors with $5e^{-}$ threshold charges.
• ### Constraining the Mass Scale of a Lorentz-Violating Hamiltonian with the Measurement of Astrophysical Neutrino-Flavor Composition(1704.04027)

Jan. 12, 2018 hep-ph, hep-ex, astro-ph.HE
We study Lorentz violation effects on flavor transitions of high energy astrophysical neutrinos. It is shown that the appearance of Lorentz violating Hamiltonian can drastically change the flavor transition probabilities of astrophysical neutrinos. Predictions of Lorentz violation effects on flavor compositions of astrophysical neutrinos arriving on Earth are compared with IceCube flavor composition measurement which analyzes astrophysical neutrino events in the energy range between $25~{\rm TeV}$ and $2.8~{\rm PeV}$. Such a comparison indicates that the future IceCube-Gen2 will be able to place stringent constraints on Lorentz violating Hamiltonian in the neutrino sector. We work out the expected sensitivities by IceCube-Gen2 on dimension-$3$ CPT-odd and dimension-$4$ CPT-even operators in Lorentz violating Hamiltonian. The expected sensitivities can improve on the current constraints obtained from other types of experiments by more than two orders of magnitudes for certain range of the parameter space.
• ### Letter of Intent: Jinping Neutrino Experiment(1602.01733)

June 21, 2016 hep-ex, physics.ins-det
Jinping Neutrino Experiment (Jinping) is proposed to significantly improve measurements on solar neutrinos and geoneutrinos in China Jinping Laboratory - a lab with a number of unparalleled features, thickest overburden, lowest reactor neutrino background, etc., which identify it as the world-best low-energy neutrino laboratory. The proposed experiment will have target mass of 4 kilotons of liquid scintillator or water-based liquid scintillator, with a fiducial mass of 2 kilotons for neutrino-electron scattering events and 3 kilotons for inverse-beta interaction events. A number of initial sensitivities studies have been carried out, including on the transition phase for the solar neutrinos oscillation from the vacuum to the matter effect, the discovery of solar neutrinos from the carbon-nitrogen-oxygen (CNO) cycle, the resolution of the high and low metallicity hypotheses, and the unambiguous separation on U and Th cascade decays from the dominant crustal anti-electron neutrinos in China.
• ### Probing Neutrino Mass Hierarchy by Comparing the Charged-Current and Neutral-Current Interaction Rates of Supernova Neutrinos(1603.00692)

June 13, 2016 hep-ph, astro-ph.HE
The neutrino mass hierarchy is one of the neutrino fundamental properties yet to be determined. We introduce a method to determine neutrino mass hierarchy by comparing the interaction rate of neutral current (NC) interactions, $\nu(\hat{\nu}) + p\rightarrow\nu(\hat{\nu}) + p$, and inverse beta decays (IBD), $\bar{\nu}_e + p\rightarrow n + e^+$, of supernova neutrinos in scintillation detectors. Neutrino flavor conversions inside the supernova are sensitive to neutrino mass hierarchy. Due to Mikheyev-Smirnov-Wolfenstein effects, the full swapping of $\bar{\nu}_e$ flux with the $\bar{\nu}_x$ ($x=\mu,~\tau$) one occurs in the inverted hierarchy, while such a swapping does not occur in the normal hierarchy. As a result, more high energy IBD events occur in the detector for the inverted hierarchy than the high energy IBD events in the normal hierarchy. By comparing IBD interaction rate with the mass hierarchy independent NC interaction rate, one can determine the neutrino mass hierarchy.
• ### Determination of the total absorption peak in an electromagnetic calorimeter(1603.04433)

A physically-motivated function was developed to accurately determine the total absorption peak in an electromagnetic calorimeter and to overcome biases present in many commonly used methods. The function is the convolution of a detector resolution function with the sum of a delta function, which represents the complete absorption of energy, and a tail function, which describes the partial absorption of energy and depends on the detector materials and structures. Its performance was tested with the simulation of three typical cases. The accuracy of the extracted peak value, resolution, and peak area was improved by an order of magnitude on average, relative to the Crystal Ball function.
• ### Heavy-Flavor-Conserving Hadronic Weak Decays of Heavy Baryons(1512.01276)

Feb. 21, 2016 hep-ph, hep-ex
More than two decades ago, we studied heavy-flavor-conserving weak decays of heavy baryons within the framework that incorporates both heavy-quark and chiral symmetries. In view of the first observation of $\Xi_b^-\to\Lambda_b^0\pi^-$ by LHCb recently, we have reexamined these decays and presented updated predictions. The predicted rates for $\Xi_b^-\to\Lambda_b^0\pi^-$ in the MIT bag and diquark models are consistent with experiment. The major theoretical uncertainty stems from the evaluation of baryon matrix elements. The branching fraction of $\Xi_c\to\Lambda_c\pi$ is predicted to be of order $10^{-4}$. It is suppressed relative to $Br(\Xi_b\to\Lambda_b\pi)$ owing to the shorter lifetime of $\Xi_c$ relative to $\Xi_b$ and the destructive nonspectator $W$-exchange contribution. The kinematically accessible weak decays of the sextet heavy baryon $\Omega_Q$ are $\Omega_Q\to\Xi_Q\pi$. Due to the absence of the $B_6-B_{\bar 3}$ transition in the heavy quark limit and the $B_6-B_6$ transition in the model calculations, $\Omega_Q\to\Xi_Q\pi$ vanish in the heavy quark limit.
• The Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton multi-purpose underground liquid scintillator detector, was proposed with the determination of the neutrino mass hierarchy as a primary physics goal. It is also capable of observing neutrinos from terrestrial and extra-terrestrial sources, including supernova burst neutrinos, diffuse supernova neutrino background, geoneutrinos, atmospheric neutrinos, solar neutrinos, as well as exotic searches such as nucleon decays, dark matter, sterile neutrinos, etc. We present the physics motivations and the anticipated performance of the JUNO detector for various proposed measurements. By detecting reactor antineutrinos from two power plants at 53-km distance, JUNO will determine the neutrino mass hierarchy at a 3-4 sigma significance with six years of running. The measurement of antineutrino spectrum will also lead to the precise determination of three out of the six oscillation parameters to an accuracy of better than 1\%. Neutrino burst from a typical core-collapse supernova at 10 kpc would lead to ~5000 inverse-beta-decay events and ~2000 all-flavor neutrino-proton elastic scattering events in JUNO. Detection of DSNB would provide valuable information on the cosmic star-formation rate and the average core-collapsed neutrino energy spectrum. Geo-neutrinos can be detected in JUNO with a rate of ~400 events per year, significantly improving the statistics of existing geoneutrino samples. The JUNO detector is sensitive to several exotic searches, e.g. proton decay via the $p\to K^++\bar\nu$ decay channel. The JUNO detector will provide a unique facility to address many outstanding crucial questions in particle and astrophysics. It holds the great potential for further advancing our quest to understanding the fundamental properties of neutrinos, one of the building blocks of our Universe.
• ### Complementary Test of the Dark Matter Self-Interaction by Direct and Indirect Detections(1505.03781)

Sept. 28, 2015 hep-ph
The halo dark matter (DM) can be gravitationally captured by the Sun. For self-interacting DM (SIDM), we show that the number of DM trapped inside the Sun remains unsuppressed even if the DM-nucleon cross section is negligible. We consider a SIDM model where $U(1)$ gauge symmetry is introduced to account for the DM self-interaction. Such a model naturally leads to isospin violation for DM-nucleon interaction, although isospin symmetry is still allowed as a special case. We show that the detection of neutrino signature from DM annihilation in the Sun can probe those SIDM parameter ranges not reachable by direct detections. Those parameter ranges are either the region with a very small $m_{\chi}$ or the region opened up due to isospin violations.
• ### Thermal transport of the solar captured dark matter and its impact on the indirect dark matter search(1508.05263)

Sept. 1, 2015 hep-ph
We study the thermal transport occurring in the system of solar captured dark matter (DM) and explore its impact on the DM indirect search signal. We particularly focus on the scenario of self-interacting DM (SIDM). The flows of energies in and out of the system are caused by solar captures via DM-nucleon and DM-DM scatterings, the energy dissipation via DM annihilation, and the heat exchange between DM and solar nuclei. We examine the DM temperature evolution and demonstrate that the DM temperature can be higher than the core temperature of the Sun if the DM-nucleon cross section is sufficiently small such that the energy flow due to DM self-interaction becomes relatively important. We argue that the correct DM temperature should be used for accurately predicting the DM annihilation rate, which is relevant to the DM indirect detection.
• ### The dark matter self-interaction and its impact on the critical mass for dark matter evaporations inside the sun(1412.6739)

Dec. 21, 2014 hep-ph
We study the capture, annihilation and evaporation of dark matter (DM) inside the Sun. It has been shown that the DM self-interaction can increase the DM number inside the Sun. We demonstrate that this enhancement becomes more significant in the regime of small DM mass, given a fixed DM self-interaction cross section. This leads to the enhancement of neutrino flux from DM annihilation. On the other hand, for DM mass as low as as a few GeVs, not only the DM-nuclei scatterings can cause the DM evaporation, DM self-interaction also provides non-negligible contributions to this effect. Consequently, the critical mass for DM evaporation (typically 3 ~ 4 GeV without the DM self-interaction) can be slightly increased. We discuss the prospect of detecting DM self-interaction in IceCube- PINGU using the annihilation channels $\chi\chi\rightarrow\nu\bar{\nu},\:\tau^{-}\tau^{+}$ as examples. The PINGU sensitivities to DM self-interaction cross section $\sigma_{\chi\chi}$ are estimated for track and cascade events.
• ### Probing Dark Matter Self-Interaction in the Sun with IceCube-PINGU(1408.5471)

Nov. 20, 2014 hep-ph, astro-ph.HE
We study the capture, annihilation and evaporation of dark matter (DM) inside the Sun. It has been shown that the DM self-interaction can increase the DM number inside the Sun. We demonstrate that this enhancement becomes more significant in the regime of small DM mass, given a fixed DM self-interaction cross section. This leads to the enhancement of neutrino flux from DM annihilation. On the other hand, for DM mass as low as a few GeVs, not only the DM-nuclei scatterings can cause the DM evaporation, DM self-interaction also provides non-negligible contributions to this effect. Consequently, the critical mass for DM evaporation (typically 3 $\sim$ 4 GeV without the DM self-interaction) can be slightly increased. We discuss the prospect of detecting DM self-interaction in IceCube-PINGU using the annihilation channels $\chi\chi \rightarrow \tau^{+}\tau^{-}, \nu\bar{\nu}$ as examples. The PINGU sensitivities to DM self-interaction cross section $\sigma_{\chi\chi}$ are estimated for track and cascade events.
• ### Probing the coupling of heavy dark matter to nucleons by detecting neutrino signature from the Earth's core(1409.3094)

Sept. 10, 2014 hep-ph, astro-ph.HE
We argue that the detection of neutrino signature from the Earth's core can effectively probe the coupling of heavy dark matter ($m_{\chi}>10^{4}$ GeV) to nucleons. We first note that direct searches for dark matter (DM) in such a mass range provide much less stringent constraint than the constraint provided by such searches for $m_{\chi}\sim 100$ GeV. Furthermore the energies of neutrinos arising from DM annihilation inside the Sun cannot exceed a few TeVs at the Sun surface due to the attenuation effect. Therefore the sensitivity to the heavy DM coupling is lost. Finally, the detection of neutrino signature from galactic halo can only probe DM annihilation cross sections. We present neutrino event rates in IceCube and KM3NeT arising from the neutrino flux produced by annihilation of Earth-captured DM heavier than $10^{4}$ GeV. The IceCube and KM3NeT sensitivities to spin independent DM-proton scattering cross section $\sigma_{\chi p}$ in this mass range are presented for both isospin symmetric and isospin violating cases.
• ### Probing the coupling of heavy dark matter to nucleons by detecting neutrino signature from the Earth core(1404.0446)

April 2, 2014 hep-ph
We argue that the detection of neutrino signature from the Earth core is an ideal approach for probing the coupling of heavy dark matter ($m_{\chi}>10^{4}$ GeV) to nucleons. We first note that direct searches for dark matter (DM) in such a mass range do not provide stringent constraints. Furthermore the energies of neutrinos arising from DM annihilations inside the Sun cannot exceed a few TeV at the Sun surface due to the attenuation effect. Therefore the sensitivity to the heavy DM coupling is lost. Finally, the detection of neutrino signature from galactic halo can only probe DM annihilation cross sections. After presenting the rationale of our studies, we discuss the event rates in IceCube and KM3NeT arising from the neutrino flux produced by annihilations of Earth-captured DM heavier than $10^{4}$ GeV. The IceCube and KM3NeT sensitivities to spin independent DM-proton scattering cross section $\sigma_{\chi p}$ and isospin violation effect in this mass range are presented. The implications of our results are also discussed.
• ### Probing Neutrino Flavor Transition Mechanism with Ultra High Energy Astrophysical Neutrinos(1308.1828)

Aug. 30, 2013 hep-ph, astro-ph.HE
Observation of ultra-high energy astrophysical neutrinos and identification of their flavors have been proposed for future neutrino telescopes. The flavor ratio of astrophysical neutrinos observed on the Earth depends on both the initial flavor ratio at the source and flavor transitions taking place during propagations of these neutrinos. The flavor transition mechanisms are well-classified with our model-independent parametrization. We find a new parameter R={\phi}_e/({\phi}_{\mu} + {\phi}_{\tau}) can probe directly the flavor transition in the framework of our model-independent parametrization, without the assumption of the {\nu}_{\mu}-{\nu}_{\tau} symmetry. A few flavor transition models are employed to test our parametrization with this new observable. The observational constraints on flavor transition mechanisms by the new observable is discussed through our model-independent parametrization.
• ### Constraining dark matter capture and annihilation cross sections by searching for neutrino signature from the Earth core(1308.1782)

Aug. 8, 2013 hep-ph
We study the sensitivity of IceCube/DeepCore detector to dark matter annihilations in the Earth core. We focus on annihilation modes $\chi\chi\to \nu \bar{\nu}, \, \tau^+ \tau^-, \, b \bar{b}$, and $W^+W^-$. Both track and cascade events are considered in our analysis. By fixing the dark matter annihilation cross section $\langle \sigma\upsilon\rangle$ at some nominal values, we study the sensitivity of IceCube/DeepCore detector to dark matter spin-independent cross section $\sigma_p^{\rm SI}$ for $m_{\chi}$ ranging from few tens of GeV to 10 TeV. This sensitivity is compared with the existing IceCube 79-string constraint on the same cross section, which was obtained by searching for dark matter annihilations in the Sun. We compare this sensitivity to dark matter direct detection results as well, in particular the XENON100 (2012) limit and the parameter regions preferred by DAMA and CRESST-II experiments. We also present IceCube/DeepCore sensitivity to $\langle \sigma\upsilon \rangle$ as a function of $m_{\chi}$ by fixing $\sigma_p^{\rm SI}$ at XENON100 (2012) and XENON1T limits, respectively. This sensitivity is compared with the preferred dark matter parameter range derived from the combined fitting to PAMELA and AMS02 positron fraction data. We conclude that the search for dark matter annihilations in the Earth core provides competitive constraints on $\sigma_p^{\rm SI}$ and $\langle \sigma\upsilon \rangle$ in the case of low-mass dark matter. Particularly, the expected constraint on $\sigsip$ for 5 years of data taking in IceCube/DeepCore is more stringent than the current IceCube 79-string limit mentioned above.
• ### Sensitivities of IceCube DeepCore Detector to Signatures of Low-Mass Dark Matter in the Galactic Halo(1209.6226)

March 4, 2013 hep-ph, astro-ph.CO
We discuss the event rate in DeepCore array due to neutrino flux produced by annihilations and decays of galactic dark matter. This event rate is calculated with a 10 GeV threshold energy, which is smaller than the threshold energy taken in previous works. Taking into account the background event rate due to the atmospheric neutrino flux, we evaluate the sensitivity of DeepCore array for probing dark matter annihilation cross section and decay time. The sensitivity studies include the annihilation modes $\chi\chi\to b\bar{b}, \ \tau^+ \tau^-$, \ $\mu^+\mu^-$, and $\nu\bar{\nu}$, and decay modes $\chi\to b\bar{b}, \ \tau^+ \tau^-$, \ $\mu^+\mu^-$, and $\nu\bar{\nu}$. We compare our results with corresponding constraints derived from observations of WMAP, ACT and Fermi-LAT.
• ### Probing annihilations and decays of low-mass galactic dark matter in IceCube DeepCore array: Track events(1105.5719)

Feb. 11, 2012 hep-ph, astro-ph.HE
The deployment of DeepCore array significantly lowers IceCube's energy threshold to about 10 GeV and enhances the sensitivity of detecting neutrinos from annihilations and decays of light dark matter. To match this experimental development, we calculate the track event rate in DeepCore array due to neutrino flux produced by annihilations and decays of galactic dark matter. We also calculate the background event rate due to atmospheric neutrino flux for evaluating the sensitivity of DeepCore array to galactic dark matter signatures. Unlike previous approaches, which set the energy threshold for track events at around 50 GeV (this choice avoids the necessity of including oscillation effect in the estimation of atmospheric background event rate), we have set the energy threshold at 10 GeV to take the full advantage of DeepCore array. We compare our calculated sensitivity with those obtained by setting the threshold energy at 50 GeV. We conclude that our proposed threshold energy significantly improves the sensitivity of DeepCore array to the dark matter signature for $m_{\chi}< 100$ GeV in the annihilation scenario and $m_{\chi}<300$ GeV in the decay scenario.
• ### Determination of the Source Flavor Ratio of Ultrahigh Energy Neutrinos(1004.5154)

Oct. 11, 2011 hep-ph, hep-ex
We discuss the reconstruction of neutrino flavor neutrino at a distant source in the very high en- ergy regime. This reconstruction procedure is relevant to the confirmation of detecting cosmogenic neutrinos, for example. To facilitate such a reconstruction, it is imperative to achieve effective flavor discriminations in terrestrial neutrino telescopes. We note that, for energies beyond few tens of PeV, a tau-lepton behaves like a track similar to a muon. Hence, while it is rather challenging to separate {\nu}{\mu} from {\nu}{\tau} in this case, one can expect to isolate {\nu}e from the rest by a distinctive shower signature. We present the result of flavor ratio reconstruction given the anticipated accuracies of flavor measurement in neutrino telescopes and current uncertainties of neutrino mixing parame- ters. It is shown that the further separation between {\nu}{\mu} and {\nu}{\tau} events does not improve the flavor reconstruction due to the approximate {\nu}{\mu} - {\nu}{\tau} symmetry.
• ### Flavor Transition Mechanisms of Propagating Astrophysical Neutrinos -A Model Independent Parametrization(1004.1583)

Nov. 27, 2010 hep-ph, astro-ph.HE
One of the important goals for future neutrino telescopes is to identify the flavors of astrophysical neutrinos and therefore determine the flavor ratio. The flavor ratio of astrophysical neutrinos observed on the Earth depends on both the initial flavor ratio at the source and flavor transitions taking place during propagations of these neutrinos. We propose a model independent parametrization for describing the above flavor transitions. A few flavor transition models are employed to test our parametrization. The observational test for flavor transition mechanisms through our parametrization is discussed.
• ### Neutrino Flavor Ratio on Earth and at Astrophysical Sources(1004.1030)

April 7, 2010 hep-ph
We present the reconstruction of neutrino flavor ratios at astrophysical sources. For distinguishing the pion source and the muon-damped source to the 3$\sigma$ level, the neutrino flux ratios, $R\equiv\phi(\nu_\mu)/(\phi(\nu_e)+\phi(\nu_\tau))$ and $S\equiv\phi(\nu_e)/\phi(\nu_\tau)$, need to be measured in accuracies better than 10%.
• ### Determination of the Neutrino Flavor Ratio at the Astrophysical Source(0905.4003)

Nov. 3, 2009 hep-ph, astro-ph.HE
We discuss the reconstruction of neutrino flavor ratios at astrophysical sources through the future neutrino-telescope measurements. Taking the ranges of neutrino mixing parameters $\theta_{ij}$ as those given by the current global fit, we demonstrate by a statistical method that the accuracies in the measurements of energy-independent ratios $R\equiv\phi (\nu_{\mu})/(\phi (\nu_{e})+\phi (\nu_{\tau}))$ and $S\equiv\phi (\nu_e)/\phi (\nu_{\tau})$ among integrated neutrino flux should both be better than 10% in order to distinguish between the pion source and the muon-damped source at the $3 \sigma$ level. The 10% accuracy needed for measuring $R$ and $S$ requires an improved understanding on the background atmospheric neutrino flux to a better than 10% level in the future. We discuss the applicability of our analysis to practical situations that the diffuse astrophysical neutrino flux arises from different types of sources and each point source has a neutrino flavor ratio varying with energies. We also discuss the effect of leptonic CP phase on the flavor-ratio reconstruction.
• ### GeoSynchrotron Radiation from Earth Skimming Tau Neutrino Shower(0909.1504)

Sept. 8, 2009 astro-ph.CO, astro-ph.HE
Cosmogenic neutrinos are expected from ultrahigh energy cosmic rays undergoing the GZK process and anticipated to be observed by detecting air showers from the decays of tau leptons. We use CORSIKA simulated shower structure to calculate the coherent geosynchrotron radio emissions of the tau decay showers above $10^{17}$eV. We present the pattern and spectrum of radio waves and discuss their detections by radio antennae.
• ### A New Type of Plasma Wakefield Accelerator Driven by Magnetowaves(0808.3036)

We present a new concept for a plasma wakefield accelerator driven by magnetowaves (MPWA). This concept was originally proposed as a viable mechanism for the "cosmic accelerator" that would accelerate cosmic particles to ultra high energies in the astrophysical setting. Unlike the more familiar Plasma Wakefield Accelerator (PWFA) and the Laser Wakefield Accelerator (LWFA) where the drivers, the charged-particle beam and the laser, are independently existing entities, MPWA invokes the high-frequency and high-speed whistler mode as the driver, which is a medium wave that cannot exist outside of the plasma. Aside from the difference in drivers, the underlying mechanism that excites the plasma wakefield via the ponderomotive potential is common. Our computer simulations show that under appropriate conditions, the plasma wakefield maintains very high coherence and can sustain high-gradient acceleration over many plasma wavelengths. We suggest that in addition to its celestial application, the MPWA concept can also be of terrestrial utility. A proof-of-principle experiment on MPWA would benefit both terrestrial and celestial accelerator concepts.
• ### Ultra High Energy Cosmic Ray Puzzle and the Plasma Wakefield Acceleration(0808.0223)

Aug. 4, 2008 astro-ph
Magnetowave induced plasma wakefield acceleration (MPWA) in a relativistic astrophysical outflow has been proposed as a viable mechanism for the acceleration of cosmic particles to ultra high energies. Here we present simulation results that demonstrate the viability of this mechanism. We invoke the high frequency and high speed whistler mode for the driving pulse. The plasma wakefield so induced validates precisely the theoretical prediction. This mechanism is shown capable of accelerating charged particles to ZeV energies in Active Galactic Nuclei (AGN).
• ### GZK Horizons and the Anisotropy of Highest-energy Cosmic Ray Sources(0807.3391)

July 22, 2008 astro-ph
Motivated by recent Pierre Auger result on the correlation of the highest-energy cosmic rays with the nearby active galactic nuclei, we explore possible ultrahigh energy cosmic ray (UHECR) source distributions and their effects on GZK horizons. Effects on GZK horizons by local over-density of UHECR sources are examined carefully with constraints on the degree of local over-density inferred from the measured UHECR spectrum. We include the energy calibration effect on the Pierre Auger data in our studies. We propose possible local over-densities of UHECR sources which are testable in the future cosmic ray astronomy.