Next Article in Journal
Method of Wavelet-Decomposition to Research Cosmic Ray Variations: Application in Space Weather
Next Article in Special Issue
Searches for Exotic Interactions Using Neutrons
Previous Article in Journal
Developing an Optical Measuring System for Hole Saw Caps
Previous Article in Special Issue
Neutron Oscillations and the Parity Doubling Theorem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Possible Neutron-Antineutron Oscillation Experiment at PF1B at the Institut Laue Langevin

1
Department of Physics and Astronomy, University of South Carolina, Columbia, SC 29208, USA
2
Department of Physics, Technical University of Denmark, Anker Engelunds Vej 1, 2800 Lyngby, Denmark
3
Department of Physics, Stockholm University, Fysikum, 106 91 Stockholm, Sweden
4
Fysiska Institutionen, Lunds Universitet, 223 62 Lund, Sweden
5
Institut Max von Laue—Paul Langevin, 71 Avenue des Martyrs, 38042 Grenoble, France
6
Laboratoire de Physique Subatomique et de Cosmologie, UGA-CNRS/IN2P3, 38026 Grenoble, France
7
European Spallation Source, European Research Infrastructure Consortium, 224 84 Lund, Sweden
8
Department of Physics, Indiana University, 727 E. Third Str., Bloomington, IN 47405, USA
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(12), 2314; https://doi.org/10.3390/sym13122314
Submission received: 29 October 2021 / Revised: 23 November 2021 / Accepted: 26 November 2021 / Published: 3 December 2021

Abstract

:
We consider a possible neutron–antineutron ( n n ¯ ) oscillation experiment at the PF1B instrument at Institut Laue Langevin. It can improve the best existing constraint on the transition rate and also allow the testing of the methods and instrumentation which would be needed for a later larger-scale experiment at ESS. The main gain factors over the most competitive experiment, performed earlier at PF1 instrument at ILL, are: a more intense neutron beam and a new operating mode based on coherent n and n ¯ mirror reflections. The installation of such an experiment would need a temporary replacement of the existing ballistic neutron guide by a specially designed n / n ¯ guide with a gradually increasing cross section and a specially selected coating as well as the development and construction of an advanced n ¯ annihilation detector with a high efficiency and low background. The overall gain factor could reach up to an order of magnitude and depends on the chosen experiment configuration.

1. Introduction

An observation of neutron–antineutron oscillations ( n n ¯ ) 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 Δ B = 1 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 Δ B = 2 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 n n ¯ 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 n n ¯ 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 n ¯ mirror reflections. In Section 2, we present a general concept of this experiment. In Section 3, we analyze the interaction of n ¯ with the n / n ¯ guide walls and provide parameters for the design of the guide. In Section 4, we present a design of the n / n ¯ guide. In Section 5, we describe an advanced n ¯ 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 n n ¯ 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 n n ¯ experiment from PF1 to PF1B neutron facility. In Section 2.2, we describe its new operating mode based on coherent n and n ¯ mirror reflections, estimate systematic uncertainties associated with the interaction of n ¯ with the guide walls, and introduce a concept of the n / n ¯ guide.

2.1. Comparison of Statistical Sensitivity on Neutron Beams PF1 and PF1B at ILL

The best constraint on the n n ¯ oscillation time in experiments with free neutrons is equal τ n n ¯ > 0.86× 10 8 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
F PF 1 1.7× 10 11 n / s ,
the mean flight time was
t PF 1 0.11 s ,
and the total experiment duration was
T PF 1 2.4× 10 7 s .
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
F PF 1 B 1.0× 10 12 n / s .
It is higher than that at positions upstream of the beam, where a considered experiment can start. Not all n / n ¯ could reach the n ¯ annihilation detector. This is estimated below.
The mean flight time can be roughly calculated using the total experiment length and the mean neutron wavelength:
t ( 1 ) PF 1 B 0.06 s ,
t ( 2 ) PF 1 B 0.085 s .
For the length estimation, we consider here a possibility to temporarily replace the downstream sections of the ballistic PF1B neutron guide by the n n ¯ oscillations experiment. There are two options for doing this: (1) replacing all sections starting from Section 6 (see Figure 1) ( t ( 1 ) PF 1 B ), 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) ( t ( 2 ) P F 1 B ), 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 n n ¯ 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 :
T PF 1 B = 1.7× 10 7 s .
Using the formula below, we obtain a gain factor over the former PF1 position :
G a i n ( 1 ) = F PF 1 B F PF 1 · t P F 1 B t P F 1 2 · T PF 1 B T PF 1 1.2,
G a i n ( 2 ) = F PF 1 B F PF 1 · t PF 1 B t PF 1 2 · T PF 1 B T PF 1 2.5.
These estimations imply that a repetition of the n n ¯ 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 n / n ¯ guide, a longer experiment duration, a higher efficiency of n ¯ detection. We estimate these factors using direct simulations presented in Section 4.

2.2. The New Experimental Approach Based on a n / n ¯ Guide

Using the new idea of a guide for both n and n ¯ [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 n ¯ with the walls of a short guide are small. However, they should be taken into account.
The interaction of n ¯ with the wall is described in terms of the n ¯ -nucleus ( n ¯ A ) optical potential in full analogy with the well-known n-nucleus ( n A ) interaction [19]. It assumes almost complete reflection of the n ¯ 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 n ¯ 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 n ¯ A scattering length.
One uncertainty in this description is associated with the imaginary part of the n ¯ A scattering length. It accounts for the annihilation of subcritical n ¯ in the guide walls. To minimise this effect, materials are preferred with a large critical velocity for n ¯ and a guide shape that decreases normal velocities of n / n ¯ as explained in more detail in Section 4. Virtually all subcritical n ¯ produced in the n / n ¯ guide would reach the annihilation detector. The terms “subcritical” and “above-critical” n ¯ 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 n ¯ is approximately equal to the ratio
t PF 1 B 2 τ n ¯ 1.5% ,
where τ n ¯ is the storage time of n ¯ in the guide, typically ∼2 s. A better knowledge of n ¯ annihilation rates reduces further this uncertainty. A conservative theoretical analysis of the existing models of n ¯ A interaction is presented in ref. [20]. Additional experimental information can be obtained from measurements of antiproton–nucleus ( p ¯ A ) interaction at low energies. Such experiments can be performed at CERN on the n ¯ facility AD, as was recently discussed in ref. [21]. Conversion of the p ¯ A interaction to the n ¯ A 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 n ¯ are subcritical at all collisions with the n / n ¯ 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 n ¯ . For n ¯ , 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 n / n ¯ velocities in the initial super-mirror section of the PF1B neutron guide is much larger than the critical velocity of the n / n ¯ guide material for n ¯ , 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 n n ¯ experiment, the neutron guide has a rectangular cross section with the vertical size of h in = 20 cm and the horizontal size of w in = 9 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 n / n ¯ guide, it is larger and can be simulated.
The gravitational energy of a n ¯ falling from the top of the n / n ¯ guide to its bottom is about equal or even larger than the Fermi potential of a n / n ¯ guide material. If the guide is long enough, it should be subdivided into a few superimposed guides in order to decrease gravitational effects on n / n ¯ trajectories. If the length of the diverging part of the n / n ¯ guide appears to be sufficiently long to shape the perpendicular components of n / n ¯ velocities, then further upstream sections would be straight.
The actual shape of the n / n ¯ guide is optimized in this work by direct simulations described in Section 4. Before doing this, we consider in Section 3 the interaction of n ¯ with the guide walls, provide parameters of this interaction and estimate lifetimes of n ¯ in this experiment.

3. n ¯ Interaction with a Straight Guide in n n ¯ Oscillation Experiments

The interaction of n ¯ with a wall is described using an optical n ¯ A 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 ( N ¯ ) with nuclei. At very low energies, it can be presented in terms of scattering length:
Re   a ( N ¯ A ) = ( 1.54± 0.03) A 0.311± 0.005 fm , Im   a ( N ¯ A ) = 1.00± 0.04 fm .
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 n ¯ 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 N ¯ A interaction.
Using Equation (11), one can calculate the n ¯ Cu Fermi potential:
V F ( n ¯ C u ) V 0 i W = ( 104 ± 2 ) i ( 22 ± 1 ) neV .
The real part corresponds to the critical velocity of ∼4.5 m/s. The imaginary part is responsible for n ¯ annihilation in Cu.
In this simple approach, n ¯ 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 n ¯ . A real configuration of the experiment will have to be analyzed using direct simulations.
For a short n / n ¯ guide, the most important parameter is the real part of Fermi potential. As subcritical n ¯ 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 5 σ and also increase its imaginary part by 5 σ . Thus, the most conservative estimation of the Cu Fermi potential is:
V F ( n ¯ C u ) V 0 i W = 94 i 27 neV ,
and the form of the well potential is:
V ( x ) = V F , x < L 2 , 0 , L 2 x L 2 , V F , x > L 2 .
The vertical potential is modified by gravity:
V ( x ) = V F m g H 2 , x < H 2 , m g x , H 2 x H 2 , V F + m g H 2 , x > H 2 .
Note that m g H = 32 neV and m V F H 2 and m V F L 2 are very large numbers of the order 10 7 10 8 .
Strictly speaking, potential inside the bulk is a sum of Fermi and optical potentials. On the left side x = H / 2 , n ¯ 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:
R 2 X 2 = X tan X   and   R 2 X 2 = X cot X ,
with X 2 m L 2 2 2 E and R 2 m L 2 2 2 V F . 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
Γ = 4 2 m L E V 0 E 3 / 2 W .
For simplicity we assumed W V 0 E ; we also noted E E n ( 0 ) the real part of the quantum level energy. Let us note that Γ 0 when E 0 .
For larger values of W, the last expression is more complicated:
Γ = 4 2 m L E V 0 E 3 / 2 W 2 1 + η 2 1 η 1 + η 2 ,
with η = W / ( V 0 E ) .

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:
V ( x ) = V 0 , x < 0 , m g x , x 0 .
The energy spectrum of this problem can be found within the WKB-approximation [24], which is very suitable for the linear potential as usual:
E ( V E ) 1 / 2 ε 3 / 2 = 2 16 E 3 / 2 ε 3 / 2 cot E 3 / 2 ε 3 / 2 + π 4 ,
where ε = 2 m d 2 1.2 peV and d = 2 2 m g 2 1 / 3 are the characteristic energy and distance of the problem.
Both the real and imaginary parts of potential (12) are larger than ε 1.2 peV. One can thus solve the last equation using the small parameter ε V 0 , W . In the zero order approximation, one obtains the well known pure real expression
E n ( 0 ) = ε 3 2 8 π n 1 4 2 / 3 , n = 1 , 2 ,
The imaginary part (width) appears in the next order. One obtains:
Γ = 1 2 2 ε V 0 E 3 / 2 W ,
where, for simplicity, we set W V 0 E , and we note E E n ( 0 ) . This is the same behavior as a function of V as in (8) in [16].
For larger values of W, the last expression is:
Γ = 1 2 2 ε V 0 E 3 / 2 W 2 1 + η 2 1 η 1 + η 2 ,
with η = W / ( V 0 E ) .
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
Γ = 4.2× 10 16 eV ,
which corresponds to the n ¯ survival time in the lower quantum states
τ n ¯ = Γ 1.6   s .
For the highest energy accessible in this approximation E = m g H , one obtains Γ = 7.3× 10 16 eV and τ n ¯ = 0.9 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:
2 m 0 E / m g ( E m g x ) d x = ( n + C L + C R ) π ,
where C L and C R are the constants related to the potential form at the classical turning points.
For a “linear” function, their values are equal to 1 4 , for an “infinite” wall to 1 2 . For a finite abrupt wall, in the leading term on E / V , one can rewrite
C L , R = 1 2 1 π arctan k in κ out ,
where k in and κ out 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 κ out = 2 m ( V E ) and k in = 2 m E .
Let us note that expression (25) covers the limits of a “linear” function ( k in = κ out , C L , R = 1 4 ) and of an “infinite” wall ( k in / κ out = 0 , C L , R = 1 2 ).
For particles of higher energies, the turning points are those of box boundaries x = 0 and x = h . Thus, the quantification condition can be rewritten in the form:
2 m E / 2 m g E / 2 m g ( E m g x ) d x = ( n + C L + C R ) π ,
and the integral can be easily calculated
2 m 2 3 E + m g H 2 3 / 2 E m g H 2 3 / 2 = ( n + C L + C R ) π .
For a box with the infinite potential well (“a quantum bouncer in a closed court’’), C L = C R = 1 2 ) , the problem was studied in detail (see [25] and references therein). In particular, it is shown that for high energies E > m g H , the spectrum finds its usual n 2 -behavior corresponding to the spectrum in a pure box.
For a box with a finite height and complex potential (which makes k i n and κ o u t complex), the energies also become complex with the imaginary part, which is directly calculated from (27):
Γ = 2 g H m g H 2 E + m g H 2 E m g H 2 × Im arctan E + m g H 2 V 0 i W E + m g H 2 + arctan E m g H 2 V 0 i W E m g H 2 ,
In the limit of low energies, E V 0 , W V 0 , and m g H E , one obtains again Equation (16).
After some algebra calculations, one can write Equation (28) in a more explicit form:
Γ = 2 g H m g H 2 E + m g H 2 E m g H 2 φ ( ν + , ω + ) + φ ( ν , ω ) ,
φ ( ν , ω ) = 1 4 ln 1 + ν 2 + ω 2 + 2 ν 2 + ω 2 ν 1 + ν 2 + ω 2 2 ν 2 + ω 2 ν ,
ν ± = V 0 E ± m g H 2 1 , ω ± = W E ± m g H 2 .
Let us note that the width tends to a constant value when the energy is approaching the Fermi potential corrected by gravity ( ν + 0 or E V 0 m g H 2 ). 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 n ¯ , which separates the ranges of “subcritical” and “above-critical” n ¯ . With the parameters from (13) and H = 30 cm, the annihilation time for the highest energy E = V 0 m g H 2 , or perpendicular velocity ∼3.9 m/s, appears to be τ n ¯ 0.077 s, which is slightly larger than the time of flight of n / n ¯ 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 n / n ¯ guide as described in Section 4. “Fine tuning” of the effective critical velocity/energy of the guide wall for n ¯ 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 n / n ¯ guides, where the second straight section of the guide provides the dominant contribution to the experiment sensitivity. Note that a much longer n / n ¯ would slightly decrease the value of the effective critical velocity as defined here.

4. Design of the n / n ¯ 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 n / n ¯ oscillation experiment, starting from 11.6 m from the source with a special diverging guide with length of 75 m or starting from 31.6 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 n / n ¯ with perpendicular-to-surface velocities at each bounce below V c u t o f f ∼ 3.9 m/s to be reflected and reach the annihilation detector and we consider lost all n / n ¯ with perpendicular-to-surface velocities above this cut-off value. This is a conservative estimate since some n / n ¯ 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 × 1010 n/cm2s) is higher then the results of actual measurements at nominal reactor power (2.2 × 1010 n/cm2s). Therefore, a correction factor of 2.2/ 2.6= 0.81 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 Q c value of the guide to 0.56 · Q c N i = 0.0122 Å−1 ( Q c N i corresponds to v N i = 6.9 m/s) with a hard cut-off of the reflectivity for scattering vectors greater then Q c . The reflectivity is hence computed by McStas using the empirical formula derived from experimental data in [28]:
I ( Q ) = R = R 0 Q Q c R = 1 2 R 0 1 tanh ( Q m Q c ) / W 1 α Q Q c Q > Q c
where R 0 = 0.99 , m = 1 , slope α = 3.2 Å and width of cut-off W = 0.0015 Å−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 L d i v , the dimensions of the guide at the entry w i n , h i n (respectively, the width and the height), at the exit w o u t , h o u t , the divergence angle for the vertical a v and the horizontal a h plane are given.
    Symmetry 13 02314 i001
  • Assuming n sections with length L d i v / n and halving divergence at each step, the height (as well width) of each section is defined as:
    h 1 i = h 2 i 1 h 2 i = h 1 i + 2 L d i v n tan a v 2 i .
  • The following relation (shown for the vertical plane, but true for the horizontal too) links the starting and the ending section.
    h 2 o u t = h 1 i n + 2 L d i v n tan a v + 2 L d i v n tan a v 2 + + 2 L d i v n tan a v 2 n 1 h 1 i n + 2 L d i v n a v i = 0 n 1 1 2 i = h 1 i n + L d i v a v 4 2 2 n n
  • Let us impose h 1 i n = h i n and then the condition:
    h 2 o u t h o u t .
    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 n the guide is simply a straight guide. The interesting solution is the minimum integer n, which would give the biggest h 2 out allowed;
  • The same applies to w 2 out . It is important to notice that n m in h and n m in v 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 n m in h and n m in v . For example, if the 75 m-guide requires n m i n h = 1 and n m i n v = 2 , then it will end up generating a 2-sections guide, where a v is halved once after 37.5 m, while the vertical mirrors keep the same a h for the whole length;
  • If L d i v 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 1.1 m [14] diameter surrounded by the detector. Therefore, the first interesting guide exit to study is a square of 1 × 1 m 2 . In addition, two smaller (0.4 × 0.4 m 2 and 0.8 × 0.8 m 2 ) and one bigger (1.2 × 1.2 m 2 ) exit windows were also considered. The cross section at the end of the n / n ¯ guide defines the size of the annihilation detector to be developed (see Section 5). The parameters left to optimize are L d i v and the divergences a v and a h . 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:
F O M = i N n i · t n i 2 ,
where for velocity spectrum bins i, N n i is the number of neutrons per unit time reaching the annihilation detector after t n i 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 2 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 2 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 t 2 factor in the calculation of the FOM (see Equation (36)).

5. Design of the n ¯ Annihilation Detector

Using the n / n ¯ 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 n ¯ 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 n ¯ 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 d E d x . 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 π 0 γ γ . 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 Δ t = t 2 t 1 , defined as the difference between the timing of the first ( t 1 ) and last ( t 2 ) 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, Δ t 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 n n ¯ 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 n ¯ 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 n ¯ 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 n / n ¯ neutron guide.
All the following estimations are preliminary, and a more precise future analysis can change them slightly. The estimated statistical sensitivity for the n n ¯ 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 n ¯ 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 n ¯ annihilation detector, a longer n / n ¯ guide (a factor of ∼1.6, see Table 1), a more accurate consideration of n and n ¯ 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 n / n ¯ guide. This would lead to a nearly quadratic gain with regard to the increased length in the experimental sensitivity, as far as n ¯ 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 n ¯ from the walls of the n / n ¯ 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 n ¯ 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.

Author Contributions

Conceptualization, B.M., D.M., V.V.N., K.V.P., V.S. and W.M.S.; methodology, E.K., B.M., D.M., V.V.N., K.V.P., N.R. and V.S.; software, B.M. and N.R.; validation, V.G., E.K., R.W. and W.M.S.; formal analysis, V.V.N., K.V.P. and N.R.; investigation, B.M., V.V.N., K.V.P., N.R. and S.-C.Y.; writing, B.M., D.M., V.V.N., K.V.P., N.R., V.S. and W.M.S.; visualization, B.M. and N.R.; supervision, E.K., B.M., D.M. and V.S. All authors have read and agreed to the published version of the manuscript.

Funding

W.M.S. acknowledges support from US National Science Foundation grant PHY-1913789 and the Indiana University Center for Spacetime Symmetries. This work is also funded by the HighNESS project at the European Spallation Source. HighNESS is funded by the European Framework for Research and Innovation Horizon 2020, under grant agreement 951782. The authors also gratefully acknowledge project grant support from Vetenskapsrådet. V.G. acknowledge the support of the U.S. Department of Energy Office of Science, Office of Nuclear Physics program under Award No. DE-SC0020687.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to Torsten Soldner and Michael Kreuz for their very important help in this analysis.

Conflicts of Interest

The authors declare no conflict of interest..

References

  1. Sakharov, A.D. Violation of CP Invariance, C asymmetry, and baryon asymmetry of the universe. Pisma Zh. Eksp. Teor. Fiz. 1967, 5, 32, Erratum in Usp. Fiz. Nauk 1991, 161, 61. [Google Scholar]
  2. Smith, G.L.; Hoyle, C.D.; Gundlach, J.H.; Adelberger, E.G.; Heckel, B.R.; Swanson, H.E. Short range tests of the equivalence principle. Phys. Rev. D 2000, 61, 022001. [Google Scholar] [CrossRef]
  3. Schlamminger, S.; Choi, K.Y.; Wagner, T.A.; Gundlach, J.H.; Adelberger, E.G. Test of the equivalence principle using a rotating torsion balance. Phys. Rev. Lett. 2008, 100, 041101. [Google Scholar] [CrossRef] [Green Version]
  4. Dolgov, A. Non-GUT baryogenesis. Phys. Rep. 1992, 222, 309. [Google Scholar] [CrossRef]
  5. Hooft, G. Computation of the Quantum Effects Due to a Four-Dimensional Pseudoparticle. Phys. Rev. D 1976, 14, 3432, Erratum in Phys. Rev. D 1978, 18, 2199. [Google Scholar] [CrossRef]
  6. Kuzmin, V.A.; Rubakov, V.A.; Shaposhnikov, M.E. Anomalous Electroweak Baryon Number Nonconservation and GUT Mechanism for Baryogenesis. Phys. Lett. B 1987, 191, 171. [Google Scholar] [CrossRef]
  7. Babu, K.S.; Mohapatra, R.N.; Nasri, S. Post-Sphaleron Baryogenesis. Phys. Rev. Lett. 2006, 97, 131301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Babu, K.; Dev, P.S.B.; Mohapatra, R. Neutrino mass hierarchy, neutron-antineutron oscillation from baryogenesis. Phys. Rev. D 2009, 79, 015017. [Google Scholar] [CrossRef] [Green Version]
  9. Babu, K.S.; Dev, P.S.B.; Fortes, E.C.F.S.; Mohapatra, R.N. Post-Sphaleron Baryogenesis and an Upper Limit on the Neutron-Antineutron Oscillation Time. Phys. Rev. D 2013, 87, 115019. [Google Scholar] [CrossRef] [Green Version]
  10. Dev, P.S.B.; Mohapatra, R.N. TeV scale model for baryon and lepton number violation and resonant baryogenesis. Phys. Rev. D 2015, 92, 016007. [Google Scholar] [CrossRef] [Green Version]
  11. Allahverdi, R.; Dev, P.S.B.; Dutta, B. A simple testable model of baryon number violation: Baryogenesis, dark matter, neutron–antineutron oscillation and collider signals. Phys. Lett. B 2018, 779, 262. [Google Scholar] [CrossRef]
  12. Phillips, D.G., II; Snow, W.M.; Babu, K.; Banerjee, S.; Baxter, D.V.; Berezhiani, Z.; Bergevin, M.; Bhattacharya, S.; Brooijmans, G.; Castellanos, L.; et al. Neutron-Antineutron Oscillations: Theoretical Status and Experimental Prospects. Phys. Rep. 2016, 612, 1–45. [Google Scholar] [CrossRef] [Green Version]
  13. Addazi, A.; Anderson, K.; Ansell, S.; Babu, K.S.; Barrow, J.L.; Baxter, D.V.; Bentley, P.M.; Berezhiani, Z.; Bevilacqua, R.; Biondi, R.; et al. New high-sensitivity searches for neutrons converting into antineutrons and/or sterile neutrons at the HIBEAM/NNBAR experiment at the European Spallation Source. J. Phys. G 2021, 35, 070501. [Google Scholar] [CrossRef]
  14. Baldo-Ceolin, M.; Benetti, P.; Bitter, T.; Bobisut, F.; Calligarich, E.; Dolfini, R.; Dubbers, D.; El-Muzeini, P.; Genoni, M.; Gibin, D.; et al. A new experimental limit on neutron-antineutron oscillations. Z. Phys. C 1994, 63, 409. [Google Scholar] [CrossRef] [Green Version]
  15. Abele, H.; Dubbers, D.; Häse, H.; Klein, M.; Knöpfler, A.; Kreuz, M.; Lauera, T.; Märkischa, B.; Munda, D.; Nesvizhevsky, V.; et al. Characterization of a ballistic supermirror neutron guide. Nucl. Instr. Meth. A 2006, 562, 407. [Google Scholar] [CrossRef] [Green Version]
  16. Nesvizhevsky, V.V.; Gudkov, V.; Protasov, K.V.; Snow, W.M.; Voronin, A.Y. Experimental approach to search for free neutron-antineutron oscillations based on coherent neutron and antineutron mirror reflection. Phys. Rev. Lett. 2019, 122, 221802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Nesvizhevsky, V.V.; Gudkov, V.; Protasov, K.V.; Snow, W.M.; Voronin, A.Y. A new operating mode in experiments searching for free neutron-antineutron oscillations based on coherent neutron and antineutron mirror reflections. Eur. Phys. J. Web Conf. 2018, 191, 01005. [Google Scholar] [CrossRef] [Green Version]
  18. Nesvizhevsky, V.V.; Gudkov, V.; Protasov, K.V.; Snow, W.M.; Voronin, A.Y. Comment on B.O. Kerbikov, “The effect of collisions with the wall on neutron-antineutron transitions”. Phys. Lett. B 2019, 795, 362. [Google Scholar]
  19. Fermi, E. Sul moto dei neutroni nelle sostanze idrogenate. Rice Sci. 1936, 7, 13. [Google Scholar]
  20. Protasov, K.V.; Gudkov, V.; Kupriyanova, E.A.; Nesvizhevsky, V.V.; Snow, W.M.; Voronin, A.Y. Theoretical analysis of antineutron-nucleus data needed for antineutron mirrors in neutron-antineutron oscillation experiments. Phys. Rev. D 2020, 102, 075025. [Google Scholar] [CrossRef]
  21. Aghai-Khozani, H.; Barna, D.; Corradini, M.; De Salvador, D.; Hayano, R.S.; Hori, M.; Lealide, M.; Lodi-Rizzinide, E.; Mascagnadei, V.; Prest, M.; et al. Limits on antiproton-nuclei annihilation cross sections at ∼125 keV. Nucl. Phys. A 2021, 1009, 122170. [Google Scholar] [CrossRef]
  22. Batty, C.J.; Friedman, E.; Gal, A. Unified optical model approach to low energy antiproton annihilation on nuclei and to antiprotonic atoms. Nucl. Phys. A 2001, 669, 721. [Google Scholar] [CrossRef] [Green Version]
  23. Nesvizhevsky, V.V.; Börner, H.G.; Petukhov, A.K.; Abele, H.; Baeßler, S.; Rueß, F.J.; Stöferle,T. ; Westphal, A.; Gagarski, A.M.; Petrov, C.A.; et al. Quantum states of neutrons in the Earth’s gravitational field. Nature 2002, 415, 297. [Google Scholar] [CrossRef] [PubMed]
  24. Amthong, A. WKB approximation for abruptly varying potential wells. Eur. J. Phys. 2014, 35, 065009. [Google Scholar] [CrossRef]
  25. Robinett, R.W. WKB energy quantization and first order perturbation theory. Am. J. Phys. 1997, 65, 320. [Google Scholar] [CrossRef]
  26. Lefmann, K.; Nielsen, K. McStas, a general software package for neutron ray-tracing simulations. Neutron News 1999, 10, 1044–8632. [Google Scholar] [CrossRef]
  27. PaNOSC-ViNYL/McStasScript. Available online: https://github.com/PaNOSC-ViNYL/McStasScript (accessed on 29 October 2021).
  28. Clausen, K.N.; McMorrow, D.F.; Lefmann, K.; Aeppli, G.; Mason, T.E.; Schröder, A.; Issikiie, M.; Noharae, M.; Takagi, H. The RITA spectrometer at Risø—Design considerations and recent results. Phys. B Condens. Matter 1997, 241–243, 50–55. [Google Scholar] [CrossRef]
  29. Yiu, Sze-Chun.; Meirose, Bernhard; Barrow, Joshua; Bohm, Christian; Brooijmans, Gustaaf; Dunne, Katherine; Milstead, David; Nepomuceno, André; Oskarsson, Anders; Silverstein, Samuel; et al. Status of the design of an annihilation detector to observe neutron-antineutron conversions. Symmetry 2021. in preparation. [Google Scholar]
  30. Barrow, J.; Brooijmans, G.; Damian, J.I.M.; DiJulio, D.; Dunne, K.; Frost, M.; Golubeva, E.; Kamyshkov, Y.; Kittelmann, T.; Klinkby, E.; et al. Computing and Detector Simulation Framework for the HIBEAM/NNBAR Experimental Program at the ESS. EPJ Web Conf. 2021, 251, 02062. [Google Scholar] [CrossRef]
  31. Agostinelli, S.; Allison, J.; Amako, K.A.; Apostolakis, J.; Araujo, H.; Arce, P.; Asaigai, M.; Axenit, D.; Banerjeebil, S.; Barrand, G.; et al. GEANT4—A simulation toolkit. Nucl. Instrum. Meth. A 2003, 506, 250. [Google Scholar] [CrossRef] [Green Version]
  32. Allison, J.; Amako, K.; Apostolakis, J.; Araujo, H.; Dubois, P.A.; Asai, M.; Barrand, G.; Capra, R.; Chauvie, S.; Chytracek, R.; et al. Geant4 developments and applications. IEEE Trans. Nucl. Sci. 2006, 53, 270. [Google Scholar] [CrossRef] [Green Version]
  33. Allison, J.; Apostolakis, J.; Lee, S.B.; Amako, K.; Chauvie, S.; Mantero, A.; Shin, J.I.; Toshito, T.; Truscott, P.R.; Yamashita, T.; et al. Recent developments in Geant4. Nucl. Instrum. Meth. A 2016, 835, 186. [Google Scholar] [CrossRef]
  34. Golubeva, E.S.; Barrow, J.L.; Ladd, C.G. Model of n annihilation in experimental searches for n¯ transformations. Phys. Rev. D 2019, 99, 035002. [Google Scholar] [CrossRef] [Green Version]
  35. Barrow, J.L.; Golubeva, E.S.; Paryev, E.; Richard, J.M. Progress and simulations for intranuclear neutron-antineutron transformations in Ar 18 40 . Phys. Rev. D 2020, 101, 036008. [Google Scholar] [CrossRef] [Green Version]
  36. Hagmann, C.; Lange, D.; Wright, D. Cosmic-ray shower generator (CRY) for Monte Carlo transport codes. In Proceedings of the 2007 IEEE Nuclear Science Symposium Conference Record, Honolulu, HI, USA, 26 October–3 November 2017; Volume 2, pp. 1143–1146. [Google Scholar]
Figure 1. Graphical representations lateral view (not-to-scale) of the 55 m (top) and 75 m (bottom) guide with a 1 × 1 m 2 exit. The short guide starts at 31.6 m from the neutron source and has w i n = 0.09 m, h i n = 0.2 m, horizontal and vertical initial divergences of 0.01 rad and 0.009 rad, respectively. There are two sections for a total divergence length of 55 m, w o u t = 0.915 m and h o u t = 0.942 m. Similarly, the long guide starts at 11.6 m from the neutron source and has w i n = 0.077 m, h i n = 0.2 m, and horizontal and vertical initial divergences of 0.011 rad and 0.009 rad, respectively. There are three sections for a total divergence length of 75 m, w o u t = 0.852 m and h o u t = 0.987 m.
Figure 1. Graphical representations lateral view (not-to-scale) of the 55 m (top) and 75 m (bottom) guide with a 1 × 1 m 2 exit. The short guide starts at 31.6 m from the neutron source and has w i n = 0.09 m, h i n = 0.2 m, horizontal and vertical initial divergences of 0.01 rad and 0.009 rad, respectively. There are two sections for a total divergence length of 55 m, w o u t = 0.915 m and h o u t = 0.942 m. Similarly, the long guide starts at 11.6 m from the neutron source and has w i n = 0.077 m, h i n = 0.2 m, and horizontal and vertical initial divergences of 0.011 rad and 0.009 rad, respectively. There are three sections for a total divergence length of 75 m, w o u t = 0.852 m and h o u t = 0.987 m.
Symmetry 13 02314 g001
Figure 2. From top to bottom: neutron wavelength distribution at the guide entrance, wavelength and velocity distribution at the guide exit, weighted by the FOM, for the configuration with 1 × 1 m 2 exit window and 75 m-long guide shown in Figure 1. Notably, the absence in the distribution of a wavelength cut-off suggests that the contribution of low energy neutrons stays relevant even when the source absolute intensity drastically drops.
Figure 2. From top to bottom: neutron wavelength distribution at the guide entrance, wavelength and velocity distribution at the guide exit, weighted by the FOM, for the configuration with 1 × 1 m 2 exit window and 75 m-long guide shown in Figure 1. Notably, the absence in the distribution of a wavelength cut-off suggests that the contribution of low energy neutrons stays relevant even when the source absolute intensity drastically drops.
Symmetry 13 02314 g002
Figure 3. Event displays with the ILL detector showing (top) a signal event with five pions and (bottom) a cosmic muon.
Figure 3. Event displays with the ILL detector showing (top) a signal event with five pions and (bottom) a cosmic muon.
Symmetry 13 02314 g003aSymmetry 13 02314 g003b
Figure 4. Timing quantity, Δ t , for signal and cosmic muon background for the ILL detector.
Figure 4. Timing quantity, Δ t , for signal and cosmic muon background for the ILL detector.
Symmetry 13 02314 g004
Figure 5. True and detector level multiplicities for (top left) charged, (top right) neutral and (bottom) all pions.
Figure 5. True and detector level multiplicities for (top left) charged, (top right) neutral and (bottom) all pions.
Symmetry 13 02314 g005aSymmetry 13 02314 g005b
Figure 6. (Left) the invariant mass distributions from the signal and cosmic background events. Both truth level and detector level results are shown. (Right) the expected distribution of sphericity for signal events at the generator level and detector level and for cosmic ray muons at the detector level.
Figure 6. (Left) the invariant mass distributions from the signal and cosmic background events. Both truth level and detector level results are shown. (Right) the expected distribution of sphericity for signal events at the generator level and detector level and for cosmic ray muons at the detector level.
Symmetry 13 02314 g006
Table 1. Optimal FOM value for each guide cross section at the exit and guide total length ( 55 m and 75 m indicated as short and long, respectively), along with the quasi-free TOF expectation value, intensity at the exit and gain factor for one-year-long experiment defined in Section 2.1. Gravitational fall was considered at all stages of the study.
Table 1. Optimal FOM value for each guide cross section at the exit and guide total length ( 55 m and 75 m indicated as short and long, respectively), along with the quasi-free TOF expectation value, intensity at the exit and gain factor for one-year-long experiment defined in Section 2.1. Gravitational fall was considered at all stages of the study.
Exit Square Side, m
0.40.811.2
ShortLongShortLongShortLongShortLong
TOF (s)7.98 × 10 2 1.20 × 10 1 6.94 × 10 2 1.02 × 10 1 6.94 × 10 2 9.49 × 10 2 6.95 × 10 2 9.27× 10 2
I (n/s)4.66 × 10 11 3.92 × 10 11 8.89 × 10 11 6.85 × 10 11 1.02 × 10 12 8.73 × 10 11 1.11 × 10 12 9.79 × 10 11
FOM4.68 × 10 9 8.45 × 10 9 6.69 × 10 9 1.14 × 10 10 7.47 × 10 9 1.26 × 10 10 8.13 × 10 9 1.37 × 10 10
Gain1.683.202.434.042.784.463.044.77
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gudkov, V.; Klinby, E.; Meirose, B.; Milstead, D.; Nesvizhevsky, V.V.; Protasov, K.V.; Rizzi, N.; Santoro, V.; Snow, W.M.; Wagner, R.; et al. A Possible Neutron-Antineutron Oscillation Experiment at PF1B at the Institut Laue Langevin. Symmetry 2021, 13, 2314. https://doi.org/10.3390/sym13122314

AMA Style

Gudkov V, Klinby E, Meirose B, Milstead D, Nesvizhevsky VV, Protasov KV, Rizzi N, Santoro V, Snow WM, Wagner R, et al. A Possible Neutron-Antineutron Oscillation Experiment at PF1B at the Institut Laue Langevin. Symmetry. 2021; 13(12):2314. https://doi.org/10.3390/sym13122314

Chicago/Turabian Style

Gudkov, Vladimir, Esben Klinby, Bernhard Meirose, David Milstead, Valery V. Nesvizhevsky, Konstantin V. Protasov, Nicola Rizzi, Valentina Santoro, William Michael Snow, Richard Wagner, and et al. 2021. "A Possible Neutron-Antineutron Oscillation Experiment at PF1B at the Institut Laue Langevin" Symmetry 13, no. 12: 2314. https://doi.org/10.3390/sym13122314

APA Style

Gudkov, V., Klinby, E., Meirose, B., Milstead, D., Nesvizhevsky, V. V., Protasov, K. V., Rizzi, N., Santoro, V., Snow, W. M., Wagner, R., & Yiu, S. -C. (2021). A Possible Neutron-Antineutron Oscillation Experiment at PF1B at the Institut Laue Langevin. Symmetry, 13(12), 2314. https://doi.org/10.3390/sym13122314

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop