Standard Model fails to explain neutrino oscillations, dark matter, and baryon asymmetry of the Universe. All these problems can be solved with three sterile neutrinos added to SM. Quite remarkably, if sterile neutrino masses are well below the electroweak scale, this modification—Neutrino Minimal Standard Model (MSM)—can be tested experimentally. We discuss a new experiment on search for decays of GeV-scale sterile neutrinos, which are responsible for the matter-antimatter asymmetry generation and for the active neutrino masses. If lighter than 2 GeV, these particles can be produced in decays of charm mesons generated by high energy protons in a target, and subsequently decay into SM particles. To fully explore this sector of MSM, the new experiment requires data obtained with at least incident protons on target (achievable at CERN SPS in future) and a big volume detector constructed from a large amount of identical single modules, with a total sterile neutrino decay length of few kilometers. The preliminary feasibility study for the proposed experiment shows that it has sensitivity which may either lead to the discovery of new particles below the Fermi scale—right-handed partners of neutrinos—or rule out seesaw sterile neutrinos with masses below 2 GeV.

1. Introduction

The discovery of neutrino oscillations provides an undisputed signal that the Standard Model (SM) of elementary particles is not complete. However, what kind of new physics it brings to us remains still unclear: we do not know yet the properties of new particles which are believed to be behind this phenomenon. An attractive possibility is the extension of the SM by three right-handed neutrinos, making the leptonic sector similar to the quark one, see Figure 1.

The masses and Yukawa couplings of new leptons remain largely unknown. Basically, can have any value between zero (corresponding to Dirac neutrinos) and  GeV, whereas can vary from (Dirac neutrino case) to   ~ (the onset of the strong coupling). The admitted region is sketched in Figure 2.

Independently on their mass, the new Majorana leptons can explain oscillations of active neutrinos. So, an extra input is needed to fix their mass range. It can be provided by the LHC experiments.

Suppose that the resonance found at the LHC by Atlas and CMS in the region 125-126 GeV is indeed the Higgs boson of the Standard Model. This number is remarkably close to the lower limit on the Higgs mass coming from the requirement of the absolute stability of the electroweak vacuum and from Higgs inflation, and to prediction of the Higgs mass from asymptotic safety of the Standard Model (see detailed discussion in [1] and in a proposal submitted to European High Energy Strategy Group by Bezrukov et al. [2]). The existence of the Higgs boson with this particular mass tells that the Standard Model vacuum is stable or metastable with the life-time exceeding that of the Universe. The SM in this case is a valid effective field theory up to the Planck scale, and no new physics is required above the Fermi scale from this point of view. Suppose also that the LHC finds no new particle and no deviations from the Standard Model. In this case the “naturalness paradigm”, leading the theoretical research over the last few decades will be much less attractive, as the proposals for new physics stabilizing the electroweak scale by existence of new particles in the TeV region and based on low energy supersymmetry, technicolor or large extra dimensions would require severe fine-tunings.

The solution of the hierarchy problem does not require in fact the presence of new particles or new physics above the Fermi scale. Moreover, the absence of new particles between the electroweak and Planck scales, supplemented by extra symmetries (such as the scale invariance) may itself be used as an instrument towards a solution of the problem of stability of the Higgs mass against radiative corrections (for detailed arguments see [35]).

Even regardless the hierarchy problem, it is clear that the Standard Model of elementary particles is not complete. It is in conflict with several observations. These are nonzero neutrino masses and oscillations discussed above, the excess of matter over antimatter in the Universe, and the presence of nonbaryonic dark matter. Any model of physics beyond the SM (BSM) should be able to deal with the experimental troubles of the SM.

Guided by the arguments steaming from alternative solutions to the hierarchy problem it is natural to ask whether the observational problems of the SM can be solved by new physics below the Fermi scale. And the answer is affirmative: an economic way to handle in a unified way the problems of neutrino masses, dark matter, and baryon asymmetry of the Universe is to add to the SM three Majorana singlet fermions with masses roughly of the order of masses of known quarks and leptons. This theory is called the MSM, for “Neutrino Minimal Standard Model” (for a review see [6]). The lightest of the three new leptons is expected to have a mass from keV to keV and plays the role of the dark matter particle (see detailed discussion in [6] and in a proposal submitted to European High Energy Strategy Group by Boyarsky et al. [7]). Two other neutral fermions are responsible for giving masses to ordinary neutrinos via the seesaw mechanism at the electroweak scale and to creation of the baryon asymmetry of the Universe. The masses of these particles and their couplings to ordinary leptons are constrained by particle physics experiments and cosmology. Two leptons should be almost degenerate, forming thus nearly Dirac fermion (this is coming from the requirement of successful baryogenesis). For comparison, we show in Figure 3 the summary of different possibilities for the masses of Majorana leptons.

A lot of experimental efforts were devoted to the direct search of Majorana neutral leptons in the past [816]. No new particles were found, but several constraints on their mixing angles with ordinary leptons were derived. The interest to these searches declined considerably at nineties, most probably due to the theoretical prejudice that the masses of Majorana leptons should be associated with the Grand Unified scale GeV [1720], making their direct search impossible. As a result, no dedicated searches of relatively light neutral leptons were performed in the last decade, with an exception of several experiments [1315] related to the so-called Karmen anomaly [21]. Today, the reason for the negative result became clear: in the domain of parameters accessible to most of these experiments the neutral leptons interact too strongly to produce the baryon asymmetry of the Universe, not satisfying one of the Sakharov conditions for baryogenesis.

The recent theory developments provide a guideline for the required experimental sensitivity [22], which was absent at the time when the early experiments were made. The aim of the present work is to discuss the experimental signatures of the existence of these particles and estimate the parameters of the experiment of beam-target type required to fully explore the model parameter space if sterile neutrinos are lighter than -mesons (see a proposal submitted to European High Energy Strategy Group by Gorbunov and Shaposhnikov [23]).

The paper is organized as follows. In Section 2 we give the essential features of the model and discuss the constraints on new particles. Section 3 is devoted to general description of experiments on searches for sterile neutrinos of masses in GeV range. In Section 4 we present an estimate of sensitivity required to test the model and suggest the experimental setup to accomplish the ultimate task. Section 5 contains conclusions.

2. The Model and Constraints on Properties of New Particles

The MSM is described by the most general renormalizable Lagrangian of all the SM particles and 3 singlet fermions. For the purpose of the present discussion we take away from it the lightest singlet fermion (the “dark matter sterile neutrino”), which is coupled extremely weakly to the ordinary leptons. In addition, we take and degenerate in mass, . The approximate equality of the masses of comes from requirement of resonance production of the baryon asymmetry of the Universe. Then a convenient parameterization of the mass term and the interaction of 's with the leptons of SM is given in terms of three parameters (, , ): where and are the combinations of the left-handed doublets , and , with . In (2.1) GeV is the vacuum expectation value of the Higgs field , , and are the active neutrino masses. Note that one of them is negligibly small in comparison with others in the MSM [6], so that with  eV being the atmospheric neutrino mass difference, for normal (inverted) hierarchy of neutrino masses. Mixing parameters and can be expressed through the parameters of the active neutrino mixing matrix (explicit relations can be found in [24]). The parameter (by definition, ) and the CP-breaking phase cannot be determined by neutrino oscillation experiments.

If the value of mass parameter is fixed, smaller yields stronger interactions of singlet fermions to the SM leptons. This leads to equilibration of these particles in the early Universe above the electroweak temperatures, and, therefore, to erasing of the baryon asymmetry. In other words, the mixing angle between neutral leptons and active neutrinos must be small, explaining why these new particles have not been seen previously. For small , The region, where baryogenesis is possible in plane is shown in Figure 4. We also plot there the exclusion regions coming from different experiments such as CHARM [10], NuTeV [13], and CERN PS191 experiment [11, 12] (see also discussion of different experiments in [25, 26]). Only CERN PS191 significantly entered into the cosmologically interesting part of the parameter space of the MSM for sterile neutrinos lighter than kaon. In Figure 5 we present the expected lifetime of the singlet fermions in an experimentally interesting region GeV.

As we see, the mass of these particles is limited from below by current experiments and cosmological considerations by ~ MeV, while no known solid upper bound, better than the electroweak scale, can be applied. At the same time, various considerations indicate that their mass may be in GeV region [24].

3. General Considerations for the Experiment

The most efficient mechanism of sterile neutrino production is associated with weak decays of heavy mesons and baryons, see left panel of Figure 6 for examples of relevant two- and three-body decays. Therefore, from experimental point of view, the distinct mass ranges are associated with the masses of mesons in which sterile neutrinos can be created: below  MeV (-meson), between  MeV and GeV (-mesons), between and GeV (-mesons), and above GeV. We will concentrate here on the masses below the charm threshold, the higher values will be discussed shortly in Section 5.

The mechanism works as follows. Heavy mesons (and baryons) can be produced by energetic protons scattering off the target material. As a reference machine we take CERN SPS with beam energy of 400 GeV. Here the relevant heavy hadrons contain charm as the heaviest quark: , , ,, With typical weak lifetime (in the rest frame) of about  s these mesons mostly decay before further interaction in the target, no matter how thick it is. Sterile-active neutrino mixing gives rise to sterile neutrino production in weak decays of the heavy mesons (note that searches for those heavy hadron decays is one of the possible strategy of exploring the MSM parameter space in full, which at the present level of experimental technique seems unrealistic, see discussion in Section 5). For charmed mesons typical branching ratios expected for the MSM parameter space are at the level of [22] (referring to upper and lower limits in Figure 4, resp.).

These sterile neutrinos further weakly decay to the SM particles due to mixing with active neutrinos. In the MSM mixing angles are small, see Figure 4, hence the sterile neutrinos live much longer (see Figure 5) as compared to that of the weakly decaying SM particles of similar masses. Before the decay (and with account of -factor) relativistic sterile neutrinos would cover quite a large distance significantly exceeding ten kilometers. Hence, sterile neutrino decays into SM particles due to mixing with active neutrino can be searched for in the near detector, see Figure 7. Two examples of the interesting decay modes are presented in right panel of Figure 6. More examples and typical patterns of sterile neutrino decay branching ratios expected within the MSM can be found in [22].

Let us first make rough estimates of length-scale of the detector needed to fully explore the MSM parameter space for sterile neutrino masses in the interesting range –2GeV, where sterile neutrinos are dominantly produced in charmed hadron decays. Let and be average transverse and longitudinal momenta of produced sterile neutrinos. We start with detector of a special geometry designed to cover entirely the solid angle ~ around the beam axis. Then all the sterile neutrinos travel through the fiducial volume of the detector which has a conical form. For 400 GeV proton beam operating at SPS the charm production cross section was measured by LEBC/EHS collaboration revealing b [28]. Thus for protons incident on the target per year one expects charmed hadrons (mb [29] is the total proton cross section at SPS energy). With charm branching ratio to sterile neutrino given above (3.1) one expects about sterile neutrinos crossing the detector each year of operating.

Then, for reference sterile neutrino mass of 1 GeV one estimates the decay length of  km (see Figure 5) and hence decay events per year for 1 km-length detector. Note, that with ~ [28] and heavy neutrino following the parent meson travel direction, the detector at a distance of 1 km has to cover the area of about m to accomplish the task.

Considering a more realistic setup, where detector is placed on the Earth surface and its height is of about 5 m (cf. the size of CHARM detector [10]) one concludes that the detector has to be about 5 times longer to achieve the same sensitivity to the sterile neutrinos as in the previous conical configuration.

To register all the expected sterile neutrino decays it seems appropriate to install many detectors of a reasonable m size rather than one large detector. These small detectors may operate separately by independent groups of experimentalists but together cover the same volume as the large detector, see Figure 8.

4. Preliminary Study for the Feasibility of the Proposed Experiment

The searches for Majorana leptons were undertaken in the past at PS191, BEBC, CHARM, and NuTeV and led to a negative result, not surprising in view of the cosmological constraints on the properties of sterile neutrinos, that were derived later on. The decays of the MSM sterile neutrinos are very rare events, and their observation presents a challenge for the detector design and performance. The sensitivity of the search has dependence on the mixing stength, and in order to increase it substantially as compared to the previous searches, one has to use, (i)higher intensity of the proton beam and larger amounts of pot’s, (ii)shorter distance to the proton target,(iii)larger decay fiducial volume, (iv)longer time of data accumulation, (v)better efficiency of events detection and reconstruction, background rejection, and so forth. To compromise the cost of the detector and its performance, we focus mainly on discussions of the experimental setup to search for the decays in a detector volume filled with a (diluted) gas, not in vacuum. The setup with vacuum requirements would provide possibly better sensitivity, but is more complicated and costly.

The main components of the experimental setup to search for the decays of sterile neutrino are schematically illustrated in Figure 8. The full detector is composed of single identical modules. Each module consists of a decay volume tank (DV) with diameter of about 5–7 meters and length of 20–30 meters. To minimize cost and amount of passive material, the decay tank could be, for example, a rubber balloon filled with helium. Each DV is followed by a high precision tracker, which could be a layers of straw tube (ST), see, for example, [30]. The entire detector is assembled from 50–100 such modules and is followed by a calorimeter (E/HCAL) which has the electromagnetic and hadronic parts. The detector is protected against charged secondaries from neutrino interactions in the dump by a Veto (V).

The design of a single module is similar to that one used in the PS191, or CHARM experiments, that searched for decays of heavy neutrinos in the mass range from 10 MeV to 1.8 GeV originated from decays , and charmed mesons decays [1012]. The later experiment, specifically designed to search for neutrino decays in a high-energy neutrino beam, was performed by using 400 GeV protons from the CERN Super Proton Synchrotron (SPS) with the total number of protons on (Cu) target (pot). The CHARM decay detector, located at the distance of 480 m from the target, consists of decay volume of , three chambers modules located inside the volume to detect charged tracks and followed by a calorimeter. The decay volume was essentially an empty region to substantially reduced the number of ordinary neutrino interactions. The signature of the heavy neutrino decay were events originating in the decay region at a small angle with respect to the neutrino beam axis with one or two separate electromagnetic showers in the calorimeter [10]. No such events were observed and limits were established on the mixing strength as a function of the mass.

4.1. Production and Decay of Sterile Neutrinos

If the decays exist, one expects a flux of high energy 's from the SPS target, since neutral mesons and are abundantly produced in the forward direction by high energy protons in the target. If is a relatively long-lived particle, this flux would penetrate the downstream shielding without significant attenuation and would be observed in the proposed detector via the decay into a high energy or pair, as schematically illustrated in Figure 7. The experimental signature of the sterile neutrino decays is an observation of either or pairs originated from a common vertex located in empty space. This signature is clean and the signal events are expected to be selected in the detector with a small background. The occurrence of decays would appear as an excess of or pairs in the detector above those expected from standard neutrino interactions. Up to a mass of 2 GeV, the mass difference between the and the , 's can originate both from the produced directly in the decay and from the produced indirectly in the subsequent decay of produced in decays. For a mass below than 1.8 GeV ’s produced in decay can contribute. However these indirect 's have a lower acceptance resulting in a smaller probability to observe the decays in the detector. Therefore we focus mainly on directly produced ’s.

As discussed previously, the mixing between the sterile neutrino and the muon neutrino results in the decay , as illustrated in Figure 8. The meson, which is normally decays into a and , might instead decay to a and a heavy neutrino . For the mass interval GeV the chirality-flip is mostly due to sterile neutrino mass which results in In the SM, meson decays leptonically via annihilation of the and quarks through a virtual . The decay rate of this process is given by where the is the meson mass, is the mass of the charged lepton, is the decay constant, is the Fermi constant, and is a Cabibbo-Kobayashi-Maskawa matrix element which value equals 0.97334 [31]. The decay rate (4.2) is suppressed by the lepton mass squared, since the very leptonic decay is due to chirality-flip. As follows from (4.1), in neutrino scattering experiments the decay signal rate is proportional to (the mixing appears twice, through the heavy neutrino production, and through the its decays in the detector) and, thus is more suppressed. The branching ratio of a particular decay mode, for example, (or ) is given by where the function is calculated in [22], and the total rate is dominated by the decay channel.

To make quantitative estimates, we performed simplified simulations of the production in the inclusive reaction with the emission of a GeV-scale sterile neutrino subsequently decaying semileptonically or leptonically in the detector, as illustrated in Figures 7 and 8. In the simulations it is assumed, that the is a long-lived particle and the flux of 's penetrates the downstream shielding without significant attenuation. The decay cannot be distinguished from the antineutrino decay or and the obtained result therefore refers to the sum of all these decay modes.

The flux of sterile neutrinos from decays of ’s and ’s produced in the target by primary protons can be expressed as follows: where is the number of pot, is the meson production cross-section, is the decay mode branching fraction [31], and is the decay phase space factor, respectively.

To perform quantitative we used simulations of the process shown in Figure 7 by taking into account the relative normalization of the yield of charmed meson species and from the original publications. The invariant cross-section of a charm meson production can be expressed as [31] where is the transverse momentum of the particle, is its rapidity, and in the last equality integration over the full 2 azimuthal angle is performed. For the production cross-sections of and charm mesons we used the Bourquin-Gaillard (B-G) formula from [32]: where with , being the rapidity, is the Mandelstam variable in , and is the transverse energy (momentum) in GeV. The coefficient in (4.7) is a normalization factor that was tuned to obtain the cross sections of the charm production in pp collisions at 400 GeV. The spectra of ’s produced in the Be target by primary protons could be calculated by using the approach reported in [33] (see also [34]). However, the Bourquin-Gaillard (B-G) approach gives the parametric form of (4.6) for the production of many different hadrons, including charm, which is in good agreement with data over the full phase-space in high-energy hadronic collisions and with a precision which is satisfactory at the level of our accuracy. The total and production cross sections in -Be collisions were calculated from its linear extrapolation to the target atomic number. The contribution of protons not interacting in the Be target and interacting in the SPS beam dump was not taken into account.

The calculated fluxes and energy distributions of produced in the Be target were used to predict the flux of ’s, as a function of its mass. An example of 2D-flux distribution calculated for 1 GeV sterile neutrino as a function of , the angle between the proton beam axis and momentum of , and the energy for 400 GeV proton incident on target is shown in Figure 9. In these calculation the ’s are produced through the decay assuming mixing strength . One can see, that for detector located at the distance 100 m from the target its fiducial area for efficient detection should be of the order .

4.2. Background

There are several expected sources of background in the experiment, which come from (i)neutrino interactions in helium, (ii)neutrino interactions in surrounding materials, (iii)cosmics, (iv)accidentals,(v)punch-through energetic neutral particles, such as or . Because of the high number of pot’s and many possible sources of background the detail study for the feasibility of the proposed search would require the simulation of a very large number of events resulting in a prohibitively large amount of computer time. Consequently, only small fraction of the required statistics for inelastic reactions and other background components, for example, such as coherent and quasielastic reactions, and so forth…, were considered for preliminary estimate and combined with numerical calculations.

The estimate shows that the dominant background is expected from the primary interactions in passive material, which correspond roughly to / pot events per 1 ton of matter. The largest contribution to background is expected from neutrino interactions yielding a single pair with little hadronic activity in the final state. Neutrino interactions that occur in the surrounding material, for example, upstream of the DV fiducial volume or in material of the decay tanks, also may produce an isolated pairs if the accompanying particles are not detected. These background events from surrounding material could be rejected with an effective method developed in [15], by introducing the collinearity variable , where is the angle between, for example, the pair momentum direction and the neutrino beam axis. The main idea is that if the pair originates from the decay chain , at high energies the pair momentum extrapolated back should point to the production target. An example of the distribution of the variable obtained in the NOMAD experiment, searching for the decay , is shown in Figure 5 of [15]. One can see that potentially one could measure the pair momentum direction with precision . By using this variable, two independent techniques could be used for the background estimation in the signal region. The first is obviously based on the MC. The second method relies on the data themselves. By extrapolating the observed events to the signal region with the shape of the fully simulated MC events a background estimate can be obtained from the data. Using this technique would result in an additional rejection of background events.

4.3. Sensitivity Estimate

For a given flux of 's the expected number of decays occurring within the fiducial length of the detector located at a distance from the neutrino target is given by where , and are the energy, momentum, and the lifetime at rest, respectively, and (%) is the () pairs reconstruction efficiency. The acceptance was calculated by tracing 's produced in the Be-target to the detector taking the relevant momentum and angular distributions into account. In this estimate the average momentum is GeV and %.

In the case of no signal observation, the background free experiment would result in stringent constraints on existence of the sub-GeV sector of the MSM. An example of the exclusion region calculated in the () plane for  pot and the detector assembled from modules, assuming at 90% CL, is shown in Figure 10. For the mass range GeV these limits correspond to the bounds on tne mixing strength which are in the range and exclude seesaw sterile neutrinos with masses below 2 GeV.

Note, that the uncertainty for the production at low does not significantly affect the sensitivity estimate. For example, the variation of the yield in (4.9) by 30% results in the corresponding variation of the limits on mixing strength of order 5%. This is because the sensitivity of the search is proportional to the . Indeed, in (4.9) the branching fraction of (4.2) and the decay rate of (4.1) both are proportional to .

5. Conclusions

The energy scale of new physics is not known. If it exists at energies above the Fermi scale (examples include supersymmetry, large or warped extra dimensions, and models with dynamical electroweak symmetry breaking), the search for new particles can be carried out in direct experiments, such as ATLAS and CMS at LHC. In addition, new hypothetical heavy particles inevitably appear as virtual states, leading to different rare processes, absent in the SM. These effects can be found at experiments such as LHCb and are competitive with the direct searches. If no new physics is found at the LHC, this may indicate the absence of new physics all the way between the Fermi and Planck scales.

Quite paradoxically, the largely unexplored up to now domain of energies where the new physics can be hidden is related to physics below the Fermi scale. If the new particles are light and weakly interacting (as in the MSM, discussed here) then the search for rare processes is superior to high energy experiments. It provides a unique possibility for discovery of new physics, not accessible by any of the LHC experiments.

The present proposal deals with the searches of new relatively light Majorana leptons which “give” the masses to ordinary neutrinos and produce baryon asymmetry of the Universe. We discussed here the parameters of the beam-target experiment which can explore completely the mass range below GeV. Due to the weak coupling of these particles to other fermions these experiments are very challenging and require construction of a kilometer scale detector. However, an experiment at much smaller scales, with 100 meters detectors, would be able to explore the dominant fraction of the parameter space.

Below we describe shortly the prospect for other experiments and higher range of sterile neutrino masses.

We start with remark that signal statistics in beam-target experiment proposed in this paper is determined by the fourth power of sterile-active neutrino mixing, , since both sterile neutrino production in meson decays and rate of subsequent neutrino decay within a detector are proportional to . Experiments on searches of meson decays to sterile neutrinos requires (potentially) less statistics of mesons, since the signal statistics are proportional only to the second power of sterile-active neutrino mixing .

With expected for MSM branching ratios (3.1) the relevant experiment has to collect not less than hadrons in order to fully explore the entire part of parameter space at GeV. This statistics is far above what has been obtained at previous - and -factories (CLEO(-), BaBar, Belle) or at proposed future - (Novosibirsk) and -factories (Tsukuba, Frascati), where total number of -meson is planned at the level of , see, for example, [35]. Only at LHC the -meson production rate is enough to perform this search, however, limits from data taking rate, QCD-background, vertex resolution, particle misidentification etc, prevent measurement of -meson branching ratios at the required level (3.1).

As we have already mentioned, there are no solid-motivated upper limits on sterile neutrino mass of subelectroweak scale. Sterile neutrinos can be heavier than -meson, and then charmed hadrons are out of use. If sterile neutrinos are lighter than 5 GeV, they can be produced in beauty hadron decays. This mass range is favorable [27] for a particular mechanism of dark matter production (resonant oscillations in the primordial plasma with lepton asymmetry) in the early Universe. There are other options for dark matter production, which are insensitive to the 1 GeV scale sterile neutrinos, see, for example, [36, 37]. Typical branching ratios of -mesons to sterile neutrino expected within MSM are [22] by one-two orders of magnitude smaller, than those of -mesons, as both upper and lower bounds on mixing decrease with increase of sterile neutrino masses, see Figure 4. Proposed upgrade of LHCb [38] will be able of entering the interesting region of MSM parameter space by searching for -meson decays into sterile neutrinos, but cannot cover the entire region. Once can further discuss an idea of using LHC-beam (4 TeV, 7 TeV, higher) for a beam-target experiment of the kind proposed in this paper, in order to fully explore the MSM parameter space, if sterile neutrino mass belongs to the the interval 2–5 GeV.

The preliminary analysis shows that the quoted sensitivity could be obtained with a proton beam and a detector optimized for several its properties. Namely, (i) the primary beam intensity and the total number of pot’s, (ii) the fiducial decay volume of the detector, (iii) the composition of the material of decay tank and surrounding material and (iv) the efficiency and precision of the tracking system, (v) the E/HCAL granularity, energy resolution, and particle identification, and (vi) different sources of background a are of importance.

As far as the cost concerned, each module may be relatively cheap, empty-space with simple tracker system either inside or behind the decay volume, and the calorimeter at the far end. Its design may repeat the design of the experiment on precision neutrino measurements and searches for sterile neutrino decays at FNAL [30]. The detector of the 100 length has a sensitivity to worth by a factor only . Thus even one (the first) section will be able to deepen considerably into the parameter range interesting for cosmology, see Figure 4.


The work was supported in part by the Grant of the President of the Russian Federation NS-5590.2012.2, by the SCOPES program (D.G. and M.S.), and by the RFBR Grant 11-02-01528a (DG).