1. Introduction
An observation of neutron–antineutron oscillations (
) would be a major scientific discovery with fundamental implications for particle physics and cosmology. This process would violate baryon number (
B) by two units. Although baryon number violation has not yet been seen in any laboratory experiment, it is the most obvious necessary ingredient in any attempt to explain the matter–antimatter asymmetry of the universe in terms of the Big Bang theory, as Sakharov explained long ago [
1]. Unlike electric charge, the conservation of which is intimately associated with its role within the Abelian gauge theory of electromagnetism, there is no experimental evidence for any similar gauge interaction associated with baryon number, which would automatically lead to baryon number conservation. If present, the long range that any such weakly-interacting “baryphoton” should possess would ruin the stringent experimental tests of the equivalence principle [
2,
3]. The spectacular observation of gravitational waves with properties as predicted by general relativity, the theory founded in part on this principle, arguably leaves even less room to expect baryphotons than before.
Given the above, we therefore expect that the baryon number is violated based on our present understanding of both particle physics and cosmology [
4]. If this is the case, then it is “only” a matter of finding it. One should of course look in any system that supports a sensitive search, which is consistent with all other known constraints. The present Standard Model of particles and interactions can exhibit baryon number and lepton number violation through nonperturbative electroweak gauge field configurations with a nontrivial topological winding number [
5]. The finite action of these field configurations leads to the exponential suppression of these amplitudes in our vacuum today. However, thermally-activated tunneling in the hot early universe can activate them [
6], and serious discussions about the possibility of finding these “sphalerons” in next-generation high energy colliders have begun. The other alternatives are baryon number violation by one or two units.
From the cosmological point of view, there is an important distinction to be made among all cosmological baryogenesis models, namely, does the baryogenesis occur at an energy scale above, at, or below the electroweak crossover/phase transition, which is one of the most natural candidates for the out-of-equilibrium dynamics required by another one of Sakharov’s baryogenesis conditions. If the sphaleron dynamics expected from the Standard Model are a feature of nature, they can erase any baryon number violation from high energy scale processes of the type expected by dimensional analysis from the
operators responsible for proton decay. It can also convert high-scale lepton number violation into baryon number violation, as suggested in leptogenesis models. By contrast, the scales associated with
operators can be slightly below the electroweak scale and still be consistent with present experiments, leading to the idea of “post-sphaeleron” baryogenesis [
7,
8,
9,
10,
11] (PSB), so it is still possible to assert that sphaleron dynamics need not influence the cosmological baryon to photon ratio, which a successful theory of baryogenesis must explain. The long-term reach of the experiment to search for free neutron–antineutron oscillations described below has a potential for either a discovery of fundamental significance or the ability to narrow the phase space for PSB models so much that this option becomes very unlikely if the result is null, in which case it would be possible to conclude from experiment that sphaleron dynamics must (not merely could) be relevant for baryogenesis within the Sakharov paradigm. Opportunities for one laboratory experiment to close a loophole and make a qualitatively nontrivial statement of this nature through a null result about a process as difficult to access and as fascinating as baryogenesis are few and far between and should be seized. Neither searches for proton decay nor for CP violation in neutrino oscillations with foreseeable sensitivity can make a comparably crisp claim about the nature of baryogenesis even if one or both of these processes are discovered. A detailed overview of the theoretical status and experimental prospects of searches for
oscillations could be found in ref. [
12].
In this note, we explore the feasibility of an experiment at the PF1B instrument at the ILL to search for
oscillations, which has a potential to explore beyond the existing constraint. The main gain factors over the best experiment performed earlier at PF1 are: a more intense neutron beam and a new operating mode based on coherent
n and
mirror reflections. In
Section 2, we present a general concept of this experiment. In
Section 3, we analyze the interaction of
with the
guide walls and provide parameters for the design of the guide. In
Section 4, we present a design of the
guide. In
Section 5, we describe an advanced
annihilation detector. Such an experiment could take place before and would be complementary to the proposed HIBEAM/NNBAR program of neutron conversion searches at the European Spallation Source [
13], at which an ultimate improvement of sensitivity to neutron–antineutron oscillations of three orders of magnitude compared to the last search with free neutrons [
14] is expected. In all such experiments, vacuum along the flight path should be good enough and magnetic shielding around the experiment efficient enough to satisfy the “quasi-free” condition, which means that the probability of
oscillations is not suppressed.
2. A General Concept of the Proposed Experiment
In
Section 2.1, we estimate a possible statistical sensitivity gain due to the move of the
experiment from PF1 to PF1B neutron facility. In
Section 2.2, we describe its new operating mode based on coherent
n and
mirror reflections, estimate systematic uncertainties associated with the interaction of
with the guide walls, and introduce a concept of the
guide.
2.1. Comparison of Statistical Sensitivity on Neutron Beams PF1 and PF1B at ILL
The best constraint on the
oscillation time in experiments with free neutrons is equal
s [
14]. It was obtained using the PF1 neutron facility at ILL.
We start from comparing characteristic parameters of PF1 and PF1B neutron beams and their effect on the statistical sensitivity of this experiment; here, we still assume the experimental configurations to be the same.
The total neutron flux at PF1 at the entrance to the experiment was
the mean flight time was
and the total experiment duration was
Since the last experiment, a more intense neutron beam has been built, and the PF1 instrument was moved to a new position called PF1B. Its parameters at the exit of the guide have been measured [
15]; the guide modernization mentioned in this paper has been done. The beam parameters at any intermediate point of the guide have not been directly measured, but could be calculated to reasonable accuracy. The total neutron particle flux at the exit of the PF1B neutron guide is
It is higher than that at positions upstream of the beam, where a considered experiment can start. Not all could reach the annihilation detector. This is estimated below.
The mean flight time can be roughly calculated using the total experiment length and the mean neutron wavelength:
For the length estimation, we consider here a possibility to temporarily replace the downstream sections of the ballistic PF1B neutron guide by the
oscillations experiment. There are two options for doing this: (1) replacing all sections starting from
Section 6 (see
Figure 1) (
), with the guide length of 55 m guide plus the size of the annihilation detector, (2) replacing all sections starting from
Section 4 (see
Figure 1) (
), with the guide length of 75 m plus the size of the annihilation detector. Although technically possible, the second option is more difficult and expensive to realize. We do not consider further continuation of the experiment downstream the PF1B experimental zone because this modification would affect or even exclude the operation of some other ILL instruments. The mean neutron wavelength at the PF1B neutron guide exit is ∼4.5 A. However, it depends on the procedure of measurement and is different at the location upstream of the beam where a
experiment can take place. It is simulated below.
A typical duration of major particle physics experiments at PF1B is 4 reactor cycles; there are 50 days per reactor cycle. This results in the experiment duration :
Using the formula below, we obtain a gain factor over the former PF1 position :
These estimations imply that a repetition of the
experiment at the PF1B neutron beam would give a factor of 1.2 or 2.5 improvement in the sensitivity. However, there are also other gain factors to be considered and taken into account. These are an eventual higher neutron flux and a different spectrum at the entrance to the experiment, better transport of slow neutrons due to the
guide, a longer experiment duration, a higher efficiency of
detection. We estimate these factors using direct simulations presented in
Section 4.
2.2. The New Experimental Approach Based on a Guide
Using the new idea of a guide for both
n and
[
16,
17,
18], the experiment can become more compact in the transverse directions and less expensive and thus more feasible. The theoretical uncertainties associated with the interaction of
with the walls of a short guide are small. However, they should be taken into account.
The interaction of
with the wall is described in terms of the
-nucleus (
) optical potential in full analogy with the well-known
n-nucleus (
) interaction [
19]. It assumes almost complete reflection of the
from the wall if the component of the neutron velocity normal to the wall is smaller than a certain critical value. In the opposite case the
annihilates in the wall with a high probability. This model of critical reflection is based on basic quantum mechanics and some knowledge of the complex
scattering length.
One uncertainty in this description is associated with the imaginary part of the
scattering length. It accounts for the annihilation of subcritical
in the guide walls. To minimise this effect, materials are preferred with a large critical velocity for
and a guide shape that decreases normal velocities of
as explained in more detail in
Section 4. Virtually all subcritical
produced in the
guide would reach the annihilation detector. The terms “subcritical” and “above-critical”
in this context are intuitively clear by analogy with normal
n; they will be defined more rigorously at the end of
Section 3. A fraction of annihilated subcritical
is approximately equal to the ratio
where
is the storage time of
in the guide, typically ∼2 s. A better knowledge of
annihilation rates reduces further this uncertainty. A conservative theoretical analysis of the existing models of
interaction is presented in ref. [
20]. Additional experimental information can be obtained from measurements of antiproton–nucleus (
) interaction at low energies. Such experiments can be performed at CERN on the
facility AD, as was recently discussed in ref. [
21]. Conversion of the
interaction to the
one contains limited theoretical uncertainties
Another uncertainty is associated with the real part of the scattering length. It accounts for the validity of condition that are subcritical at all collisions with the guide walls. To decrease it, one essentially needs to meet the same conditions as above.
First, the guide wall material should provide a large enough critical velocity for both
n and
. For
, this condition favors materials with a large atomic number [
16], for instance, copper, which is a material routinely used for early neutron guides (before the super-mirror era). Several other materials could be also considered with roughly the same performance. We limit the present analysis to only copper for simplicity.
Second, the guide shape has to be specially designed. As the spread of perpendicular components of velocities in the initial super-mirror section of the PF1B neutron guide is much larger than the critical velocity of the guide material for , the guide cross-section has to be gradually increased along the guide. It follows from the Liouville theorem that the increase in size, in the adiabatic case, would be equal to the decrease in the perpendicular velocity spread. In the realistic non-adiabatic case, this factor is a bit larger.
More details about the PF1B guide and neutron beam are given in ref. [
15]. At the entrance to the proposed
experiment, the neutron guide has a rectangular cross section with the vertical size of
cm and the horizontal size of
cm, or smaller, depending on the exact point where the experiment starts. The spread of perpendicular neutron velocities at the exit of the PF1B guide is ∼15 m/s to any direction. At the entrance to the
guide, it is larger and can be simulated.
The gravitational energy of a falling from the top of the guide to its bottom is about equal or even larger than the Fermi potential of a guide material. If the guide is long enough, it should be subdivided into a few superimposed guides in order to decrease gravitational effects on trajectories. If the length of the diverging part of the guide appears to be sufficiently long to shape the perpendicular components of velocities, then further upstream sections would be straight.
The actual shape of the
guide is optimized in this work by direct simulations described in
Section 4. Before doing this, we consider in
Section 3 the interaction of
with the guide walls, provide parameters of this interaction and estimate lifetimes of
in this experiment.
3. Interaction with a Straight Guide in Oscillation Experiments
The interaction of
with a wall is described using an optical
Fermi potential which can be evaluated using the optical model approach described in ref. [
22] and developed to describe the low energy interaction of antinucleons (
) with nuclei. At very low energies, it can be presented in terms of scattering length:
Annihilation in such systems is very strong, so this A-dependence allows a simple geometrical interpretation: the real part of the scattering length corresponds to simply the scattering of on a black disc with about the nucleus radius, whereas the imaginary part is about the same for all nuclei and proportional to the diffuseness of the interaction.
Using Equation (
11), one can calculate the
Cu Fermi potential:
The real part corresponds to the critical velocity of ∼4.5 m/s. The imaginary part is responsible for annihilation in Cu.
In this simple approach, moves in a one dimension square well with a complex Fermi potential. In the following, we perform calculations for a guide with a constant cross section as this approach allows us to get an analytical solution. A diverging guide effectively decreases perpendicular velocities, thus leading to even smaller losses of . A real configuration of the experiment will have to be analyzed using direct simulations.
For a short guide, the most important parameter is the real part of Fermi potential. As subcritical have no time to annihilate, main losses are associated with above-critical reflections.
To be conservative, we reduce the estimated real part of Cu Fermi potential (
12) by
and also increase its imaginary part by
. Thus, the most conservative estimation of the Cu Fermi potential is:
and the form of the well potential is:
The vertical potential is modified by gravity:
Note that neV and and are very large numbers of the order –.
Strictly speaking, potential inside the bulk is a sum of Fermi and optical potentials. On the left side
,
could tunnel through if their energy is very close to Fermi potential and escape from the guide. However, this tunneling effect is very small and is important only for the energies of a few peV close to the top of the Fermi potential (of the same order as energies of bound states of neutrons in a gravity field [
23]). Our conservative estimation of Fermi potential renders this a negligible effect.
We are not interested in the energy spectrum very close to Fermi potential, so we neglect all these modification of the potential outside the box.
3.1. Square Well Problem
A solution for the square well potential is well known and can be found in textbooks. It can be easily generalized to the complex potential. As usually, there two families (symmetric and anti-symmetric) of wave functions obeying standard transcendental equations:
with
and
. For the states far from Fermi potential, these equations can be solved analytically. The width of the levels due to imaginary part
W of Fermi potential is equal to
For simplicity we assumed ; we also noted the real part of the quantum level energy. Let us note that when .
For larger values of
W, the last expression is more complicated:
with
.
3.2. Step-Linear Problem
For the lowest states in the gravitation plus box potential (
15), one can use an approximation neglecting the right wall of the box:
The energy spectrum of this problem can be found within the WKB-approximation [
24], which is very suitable for the linear potential as usual:
where
peV and
are the characteristic energy and distance of the problem.
Both the real and imaginary parts of potential (
12) are larger than
peV. One can thus solve the last equation using the small parameter
. In the zero order approximation, one obtains the well known pure real expression
The imaginary part (width) appears in the next order. One obtains:
where, for simplicity, we set
, and we note
. This is the same behavior as a function of
V as in (8) in [
16].
For larger values of
W, the last expression is:
with
.
Let us note that contrary to the case of a pure box, this width tends to a constant value when energy tends to zero. In our case (
13), this constant value is equal to
which corresponds to the
survival time in the lower quantum states
For the highest energy accessible in this approximation
, one obtains
eV and
s. Let us be reminded that the time of flight of neutron in the guide is as small as 0.060–0.085 s (
5).
3.3. Linear Potential in a Box
To obtain a more general expression for a linear potential in a finite size box (
15) (An exercise for a box with infinite potential is usually called “quantum bouncer in a closed court’’). Let us note that quantification condition (
19) can be rewritten in a more general way which relates the Bohr–Sommerfelfeld integral to the obtained one at the turning points:
where
and
are the constants related to the potential form at the classical turning points.
For a “linear” function, their values are equal to
, for an “infinite” wall to
. For a finite abrupt wall, in the leading term on
, one can rewrite
where
and
are the momenta (real and complex) inside and outside the well on its left and right side.
For instance, in the previous example, on the left abrupt turning point and .
Let us note that expression (
25) covers the limits of a “linear” function (
,
) and of an “infinite” wall (
,
).
For particles of higher energies, the turning points are those of box boundaries
and
. Thus, the quantification condition can be rewritten in the form:
and the integral can be easily calculated
For a box with the infinite potential well (“a quantum bouncer in a closed court’’),
, the problem was studied in detail (see [
25] and references therein). In particular, it is shown that for high energies
, the spectrum finds its usual
-behavior corresponding to the spectrum in a pure box.
For a box with a finite height and complex potential (which makes
and
complex), the energies also become complex with the imaginary part, which is directly calculated from (
27):
In the limit of low energies,
,
, and
, one obtains again Equation (
16).
After some algebra calculations, one can write Equation (
28) in a more explicit form:
Let us note that the width tends to a constant value when the energy is approaching the Fermi potential corrected by gravity ( or ). This is due to the contribution of the imaginary part of potential to the reflection.
We will use this formalism to conservatively estimate the value of the effective critical velocity of the guide walls for
, which separates the ranges of “subcritical” and “above-critical”
. With the parameters from (
13) and
cm, the annihilation time for the highest energy
, or perpendicular velocity ∼3.9 m/s, appears to be
s, which is slightly larger than the time of flight of
through the “short” guide (0.06 s) and slightly smaller than the time of flight through the “long” guide (0.086 s) (
5). Thus, the natural choice for the maximum allowed perpendicular velocity is ∼3.9 m/s corresponding to ∼75 neV. We use it for the design of the
guide as described in
Section 4. “Fine tuning” of the effective critical velocity/energy of the guide wall for
is not important at the stage of this feasibility study. Moreover, this value would not depend significantly on the parameters of the guide. It can be done later, as soon as the guide configuration is fixed. This “fine tuning” is expected to slightly increase the experiment sensitivity. Additionally, we will use this formalism in further developments of this analysis, in particular for the cases of long
guides, where the second straight section of the guide provides the dominant contribution to the experiment sensitivity. Note that a much longer
would slightly decrease the value of the effective critical velocity as defined here.
4. Design of the Guide
The neutron guide H113 and beam parameters that feed the PF1B facility are described in ref. [
15]. The guide consists of a long ballistic middle section that is enclosed by straight entry and exit sections. There are two different possibilities, given the constraints of the facility, to install the ideal guide for performing the
oscillation experiment, starting from
m from the source with a special diverging guide with length of 75 m or starting from
m from the source with a guide of 55 m.
Based on these constraints, the ideal layout of the guide for the proposed experiment has been optimized. In the simulations, we consider all with perpendicular-to-surface velocities at each bounce below ∼ 3.9 m/s to be reflected and reach the annihilation detector and we consider lost all with perpendicular-to-surface velocities above this cut-off value. This is a conservative estimate since some with velocities above this value could still be able to reach the annihilation detector and contribute to the sensitivity. However, systematic uncertainties associated with this contribution would start increasing rapidly with increasing the perpendicular component of velocity. Therefore, we ignore this small sensitivity gain at this stage of the feasibility study.
The design of the guide was performed using McStas 2.7 [
26], a popular neutron ray tracing code ideal for this kind of study. As a starting point, we use the previously developed McStas instrument files describing the ILL cold source, the H113 guide and the PF1B experimental area (instrument file provided in a private communication by Torsten Soldner). The flux produced by the source component of the simulation (2.6 × 10
10 n/cm
2s) is higher then the results of actual measurements at nominal reactor power (2.2 × 10
10 n/cm
2s). Therefore, a correction factor of
is applied to the output intensity of the simulations. Notably, the over-prediction is not necessarily related to the source brightness since several imperfections of the H113 guide (e.g., waviness, alignment imperfections) are not included in the simulation, and they are likely to be the dominant source of the disagreement. Hence, even though we are replacing the H113 guide with a new guide, this factor was still applied to keep a conservative approach. The code used to generate the different geometries was written in Python, exploiting the interface provided by McStasScript, which is the McStas API for creating and running instruments from python scripting [
27]. The gravitational fall is taken into account throughout the whole optimization process.
The aforementioned perpendicular-to-surface velocity condition on the copper guide walls corresponds to setting the
value of the guide to
Å
−1 (
corresponds to
m/s) with a hard cut-off of the reflectivity for scattering vectors greater then
. The reflectivity is hence computed by McStas using the empirical formula derived from experimental data in [
28]:
where
,
, slope
Å and width of cut-off
Å
−1. The need for a gradual increase in the guide cross-section required by the Liouville theorem, along with the necessity of keeping the height small, lead the optimization problem in the direction of a “fixed entry and exit size” approach, in which, given the parameters that define the geometry of the guide, the number of sections is calculated such that the cross-section constraints at the entry and at the exit of the guide are respected. The approach used can be described as follows.
The total divergent length
, the dimensions of the guide at the entry
(respectively, the width and the height), at the exit
, the divergence angle for the vertical
and the horizontal
plane are given.
Assuming
n sections with length
and halving divergence at each step, the height (as well width) of each section is defined as:
The following relation (shown for the vertical plane, but true for the horizontal too) links the starting and the ending section.
Let us impose
and then the condition:
From (
34), it is clear that one can always find a value of
n large enough for (
35) to be true, where in the limit
the guide is simply a straight guide. The interesting solution is the minimum integer
n, which would give the biggest
allowed;
The same applies to . It is important to notice that and can be considerably different, hence, the algorithm determines the final number of sections of the guide in such a way to halve independently the angles based on and . For example, if the 75 m-guide requires and , then it will end up generating a 2-sections guide, where is halved once after m, while the vertical mirrors keep the same for the whole length;
If is less then the total length of the guide, a straight guide is inserted for the remaining distance.
The 99.7% of the beam of the previous experiment at PF1 was contained within the target of a
m [
14] diameter surrounded by the detector. Therefore, the first interesting guide exit to study is a square of 1 × 1 m
. In addition, two smaller (0.4 × 0.4 m
and 0.8 × 0.8 m
) and one bigger (1.2 × 1.2 m
) exit windows were also considered. The cross section at the end of the
guide defines the size of the annihilation detector to be developed (see
Section 5). The parameters left to optimize are
and the divergences
and
. The figure of merit (FOM) for the optimization used in this analysis is the same used for the optimization of the NNBAR experiment at ESS [
13] and is given by the following quantity:
where for velocity spectrum bins
i,
is the number of neutrons per unit time reaching the annihilation detector after
seconds of flight.
The parameter space to be explored by the simulations was chosen to be wide enough to also include its surroundings. In
Table 1, the optimal FOM value for each guide exit and guide total length ( 55 m and 75 m), along with the quasi-free time of flight (TOF) expectation value, the intensity at the exit and gain factor for a one-year-long experiment (
Section 2.1) are summarized. For convenience, the guide parameters that produced the optimal FOM values for all the different designs of
Table 1 are not reported in this work, but in
Figure 1 two graphical representations of the 55 m and 75 m guide with 1 × 1 m
exit are shown. Overall, we observe that the requirement of a smaller guide cross section at the exit of the beamline, produces a higher number of sections as well as a lower divergence in both the vertical and the horizontal plane. The optimal divergent length, instead, is in general always close to the maximum allowed, but hits earlier a plateau for small exit sizes where the high number of sections makes the guide almost straight.
In
Figure 2, we show the neutron wavelength and velocity distribution at the guide exit, both weighted using the FOM defined in Equation (
36), for the configuration with the exit window of 1 m
and a 75 m-long guide. The important feature to notice is that no wavelength cut-off is present in the distribution and the contribution of low energy neutrons stays relevant even when the source absolute intensity drastically drops, as expected by the
factor in the calculation of the FOM (see Equation (
36)).
5. Design of the Annihilation Detector
Using the guide simulations and the geometrical constraints in the experimental zone of PF1B, we estimate the external size of the annihilation detector in both directions to be below ∼2.8 m.
The maximum sensitivity of the experiment is achieved when the expected background is well below one event for the duration of the complete experiment and the efficiency of detection of
is maximized. A conservative estimation is that a new annihilation detector should at least achieve the 52% detection efficiency of the previous ILL [
14] experiment, but most likely, due to the use of new technologies, it is expected to surpass it.
The detector must be sensitive to the characteristic antineutron–nucleon annihilation signal. The final state consists mainly of charged pions and photons from neutral pion decays. The detector consists of a thin (∼100 m) carbon foil in which the would annihilate a tracking chamber, which will allow particle identification as well as determination of the primary vertex, and a calorimeter.
The tracking chamber will consist mainly of a Time Projection Chamber (TPC), which will provide three-dimensional tracking and a measurement of the mean energy loss
. The calorimeter comprises arrays of plastic scintillators and lead-glass modules. The primary function of the scintillators will be to identify hits originating outside the inner detector volume, which will be important for particle identification via range determinations and for the rejection of cosmic ray background events. Electromagnetic calorimetry is provided by lead-glass modules and uses the Cerenkov effect. A high precision electromagnetic calorimeter is needed to identify neutral pion production via the decay
. The calorimeters would be position sensitive, with a segmentation to be determined by simulation. This design is therefore in essence the same as being planned for the HIBEAM experiment at the European Spallation source [
29]. A simulation and analysis software framework [
30], which is based on detector simulation using
Geant-4 [
31,
32,
33], can also be used for a search at the ILL. A complete analysis with simulated datasets is beyond the scope of this work. Here, distributions of sensitive observables in signal and one of the major sources of background (cosmic ray muons) are shown to demonstrate that a feasible detector design concept exists.
Figure 3 (top) shows a signal event with five final-state pions in the ILL detector with a nuclear fragment from the carbon target. The pions’ kinetic energies range from around 220–320 MeV. The antineutron–carbon annihilation signal was calculated with the model in Refs. [
34,
35].
Figure 3 (bottom) shows a cosmic muon of kinetic energy 495 MeV impinging the ILL detector. The cosmic muon enters from the top detector module and leaves the detector from the bottom. This was made using the CRY [
36] cosmic ray program.
Figure 4 shows the quantity
, defined as the difference between the timing of the first (
) and last (
) signal in the scintillators. The spectra are shown for an annihilation event and a cosmic muon background event. Each cosmic event contains one charged cosmic muon passing through the ILL detector. Since the cosmic muons cross the top and bottom of the detector,
for cosmic background is expected to be larger than the signal. As expected, clear separation between the two distributions is observed.
Figure 5 shows the multiplicities of neutral and charged pions for 1500 annihilation events. Measurements of pion multiplicity represent important evidence that an annihilation event has occurred.
Figure 6 (Left) shows the invariant mass distributions annihilation events at truth level and at detector level, i.e., using information available from the detector such as energy loss and particle range. Detector-level background predictions for cosmic muons misidentified as pions are also shown for single muon events. The truth level invariant mass has a peak 1.88 GeV while the distribution of the reconstructed invariant is spread broadly and has a peak around 1.45 GeV in the ILL detector. The cosmic muon event distribution has a far lower invariant mass (typically around 500 MeV).
Figure 6 (Right) shows the sphericity distribution calculated for the same sets of events as shown in
Figure 6 (Left). As expected, the pure cosmic ray events have small values of sphericity, closer to zero, while signal (and signal with cosmic) events have larger sphericity.
6. Conclusions
We explore the feasibility of an experiment to search for oscillations at the PF1B instrument at ILL. The main gain factors over the best experiment performed earlier at PF1 instrument at ILL are: a stronger neutron beam and a new operating mode based on coherent n and mirror reflections. We show that the overall configuration is feasible. Due to the relatively short length available at PF1B, systematic uncertainties of the method are negligible. Virtually all subcritical would be transported without losses to the annihilation detector. A major fraction of all the initial n could be converted to subcritical ones in a special diverging neutron guide.
All the following estimations are preliminary, and a more precise future analysis can change them slightly. The estimated statistical sensitivity for the
transition rate is up to an order of magnitude higher than that in the best performed experiment [
14]. As a conservative estimate, the gain factor is ∼2.8 for the “short” neutron guide and the middle-size
detector of 1 m, see
Table 1. Potentially, it can be further improved by a longer measuring time (say, a factor of ∼2), a more efficient
annihilation detector, a longer
guide (a factor of ∼1.6, see
Table 1), a more accurate consideration of
n and
transport and the account of their interaction with the guide walls (say, a factor of ∼1.3), giving all together an additional gain factor of ∼4.1), or the total gain factor of up to ∼10.
This is large enough to provide a discovery potential. If, however, the actual transition rate is beyond the sensitivity of this experiment, it would be a significant step towards a future more sensitive larger-scale experiment at the ESS [
13]. It would allow one to test the main experimental approaches and components of a future experimental setup. In such an experiment, the length can be increased using a guide with a constant cross section that is put after the diverging part of the
guide. This would lead to a nearly quadratic gain with regard to the increased length in the experimental sensitivity, as far as
annihilation losses are not too large. Such a saturation of the sensitivity corresponds to experiment lengths much larger, or neutron spectra much softer than those considered in this paper.
The new proposed experimental scheme based on the coherent reflection of both n and from the walls of the guide, allows us to reduce the transverse size of the experiment, and therefore to reduce its cost and simplify the limitations associated with the need to provide very low magnetic fields and high vacuum, as well as to develop the annihilation detector.
The gain factor estimated in this article compares the sensitivity of the experiment [
14] at PF1 and the possible sensitivity of a future experiment at PF1B. However, the same idea of coherent reflection of
n and
can be also applied to the analysis of the already performed experiment [
14]. The fact is that, as noted by one of the referees of this article, it contained an initial focusing section of a neutron guide 33.6 m long. It effectively increases the length of the experiment and therefore increases the sensitivity of the experiment. This result will be presented separately when the corresponding calculations are performed.